Hi all,
For electrophysiology problems, I’m using the PointIntegralSolver within cbcbeat:
The FEniCS PointIntegralSolvers are special-purpose solvers dedicated to solving multistage
integration schemes defined for each vertex of a mesh. Such problems typically occur in connection with spatially varying systems of ODEs, for instance in electrophysiology.
And from the Dolfin documentation:
It only includes Point integrals with piecewise linear test functions. Such problems are disconnected at the vertices and can therefore be solved locally.
Unfortunately, it does not work with higher-order test functions. Does anybody know a work-around to solve sets of ODEs at all DOFs?
Thank you very much!