Deep Equilibrium Models

params:

# initialize
y_obs: Array["Dx Dy"] = ...
y0: Array["Dx Dy"] = zeros_like(y_obs)
solver: Solver = GaussNewton()

# solve fixed point method
y_pred: Array["Dx Dy"] = fixed_point_solver(
    y=y_obs, y0=y0,
    model=model, params=params,
    solver=solver
)