Implement C-level kernel and solver
Want to be sure that the Python-C-Python machinery works, and the current "Sparse BLAS kernel" and "boost solver" are not sufficient for this purpose because they are too complicated.
Instead we should:
-
Implement a C implementation of a Kernel, with sparse Matvec explicitly coded, and create its python wrapper (CKernel) -
Grab the DOPRI5 code and implement it as a C-level solver.
Edited by Joseph Weston