Optimization

[mnottoli]

Work in progress:

• working on the parallel efficiency of the matrix vector products and of the forces.
• precomputing the grid points which bring nonzero contributions to avoid if statements in the loops.