def linearize(self, params, unknowns, resids):

        J = {}

        # Output equation
        J[('y', 'x')] = np.array([1.0])
        J[('y', 'z')] = np.array([2.0])

        # State equation
        J[('z', 'z')] = np.array([params['x'] + 1.0])
        J[('z', 'x')] = np.array([unknowns['z']])

        return J
Analytical derivatives.
dataset/ETHPy150Open OpenMDAO/OpenMDAO/examples/implicit_ext_solve.py/SimpleImplicitComp.linearize