almost finished the linear solver

we now have a solver called colddae which has solve_causal_ddae and
solve_noncausal_ddae as subsolvers