Now sometimes we over sample by a lot, for example validation/preikszas.py spends way too much time calculating the required derivatives. Maybe a good first guess would be to take the sampling equal to the mesh size.
It looks like we first need a more generic Solution (or maybe Field) class. The sample function only has access to the current solution tuple which does not contain the mesh size...
Now sometimes we over sample by a lot, for example
validation/preikszas.py
spends way too much time calculating the required derivatives. Maybe a good first guess would be to take the sampling equal to the mesh size.