Some major refactoring of refode and refpde, possibly buggy
Some further changes to resode.
Larger class of initial data, still doesn't work for d>5
Testing the initialization values for derivatives when calling ddassl
Added linpack to the repository
Fixed compatibility between initial data and boundary conditions in harmonic.f90
Changed all real(8) and double precision to real
New makefiles, less garbage on output
Trivial change of the residuum function and it works like a charm
Removed the debug information and (ut, uxt) initialization from movcol
Makefile now compiles linpack with less output
Added an equivariance parameter
Going back to double precision
Movcol now supports multiple output files, harmonic now outputs data to
subfolders depending on d, k and npts
Added a simple python script to visualize results
Fixed the formula for k-equivariant equation
Nicer harmonic.py
Removed a print instruction from ddassl.f90
Some minor experiments with harmonic
I reworked the directory structure, updated makefiles