def METHOD_NAME(solver, state):
    """
    Compare each derrivatives to the finite differences equivalent.
    """
    t = 0
    N = solver.num_collapse
    a = solver.drift
    solver.set_state(t, state)
    assert _data.norm.l2(solver.drift(t, state) - solver.a()) < 1e-6
    for i in range(N):
        b = lambda *args: solver.diffusion(*args)[i]
        assert b(t, state) == solver.bi(i)
        for j in range(N):
            _check_equivalence(
                L(solver, j, b), solver.Libj(j, i), (t, state)
            )
    _check_equivalence(L0(solver, a), solver.L0a(), (t, state))
    for i in range(N):
        b = lambda *args: solver.diffusion(*args)[i]
        _check_equivalence(L0(solver, b), solver.L0bi(i), (t, state))
        _check_equivalence(L(solver, i, a), solver.Lia(i), (t, state))
        for j in range(i, N):
            for k in range(j, N):
                _check_equivalence(
                    LL(solver, k, j, b), solver.LiLjbk(k, j, i), (t, state)
                )