new solver colddae_noncausal2

the new solver uses a slightly different selection scheme of the
regularized system

also the option IsConst is now available for colddae_noncausal(_2) in
order to handle DDAEs with constant coeeficients and delay faster