Provide Functional Building Blocks

[danielkelshaw]

At the moment the library is very much centred around solving the Kolmogorov flow. There is no reason why the same tool can not be used to solve problems such as non/linear convection-diffusion. By providing building blocks such as:

• phi_dot_nabla_zeta()
• laplacian_phi

The user will very easily be able to define their own dynamics() function and construct their own solver. In order to do this the structure of the project should be re-thought a little bit to ensure that this is a natural extension.

The Kolmogorov solver can then be written as a sum of these parts - ideally making the dynamics() much easier to read as well.