Repeating small tests many times was order of magnitude slower. The reason is the computation of
Lagrange polynomials was both inefficient and performed every time an RegularInterpLagrange instance was
created. This causes overhead which is not relevant for large problem sizes but makes running many small
tests slow. The fix for now is to cache the creation of the polynomials. Now they are created only onnce
for each interpolation order, but not once per instance.