Add a few simple examples, that can also be used as test cases of "full" programs.
This should test quite a few cases for all methods.
Current important cases to consider:
ALL_Staggered.cpp
already.Also todo is to create a way to check the tests for validity. The proposal here is to generate a known good output and commit that. Then run the example and compare its output with the known good solution. For this we need a method to compare (with a certain degree of leeway to account for floating point rounding) the two sets of data. This shold not be difficult to code, once the output format is decided upon.