Hi and thank you for this package! Is there a change you can extend it to solve boundary value DAEs w/ COLDAE?