What are you trying to do and how would you want to do it differently?

I would like to implement a linear Verlet solver which uses matrices.

Is it something you currently cannot do? Is this related to an issue/problem?

Coupling of matrix and point based methods.

Implementation details

No response

Would you be willing to submit a PR?

Yes