PointIntegralSolver with higher-order test functions

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!