Initial import into google code, version 218.
Added space-dependend external field, requested by M. Curcic.
Fixed  bug in space-dependend external field, field was not added to h.
Backup of version used for simulations of symmetry-breaking due to roughness.
Added ToplayerHelix Homeomorphism, an efficient Helix-shaped sample geometry.
Added possibility to load saved configurations.
3D rendering now resolves the Faces' sub-triangles.
Bug-fixes in Translated/Transformed Homeomorphisms; added Outerlayer, Toplayer, which wrap other Homeomorphisms and limit them to the top/outer layers of cells to improve speed.
Added Gradient and Gyrofield data models.
Backup before implementing amusolver.
Implemented a series of AmuSolvers (0-5), which grow towards different steps for the different field contributions. AmuSolver5 employs three levels of time steps: a root time step for HSmooth, smaller steps for HFast=Hex+Hdemag, and even smaller substeps for the precssion.

Dynamically chaning the large time steps, based on the accuracy of the predictor step still needs to be implemented.
Usable christmas version with failsafe solver.
Some optimalizations in the overhead of solver AmuSolver5. The overhead seems to be the current bottleneck. Therefore, at least AmuSolver5.substeps() should be parallellized. The performance should be massively improved when only substeps are taken for the cells that really require them.
Added Unit Testing methods.

Also reverted to uncached Extrapolator2 as workaround for a bug in the cached version (moved to Extrapolator2Cached), which will be superseeded by SplineExtrapolator anyway. Such bugs are now checked for by the unit tests.
Fixed SplineExtrapolator and Extrapolator2Cached.
Adaptive mesh step 1: charge-free cells (inside magnetic domains, e.g.) can be ignored in FMM calculation. Use Problem.setChargeCutoff(). 0.01 seems a reasonable value.
Added Runge-Kutta solvers
Added Runge Kutta Fehlberg 45 and 54 solvers.
Still solves standard problem 4 perfectly.
Backup before implementing adaptive mesh stage 2: periodic re-wiring with opening angle weighed by the cell multipole. Should allow to keep the same accuracy by lowering the max angle, while still gaining performance.
Implemented dynamic rewiring. Overhead of wiring+recalculating shift factors and kernel could be substantially reduced, but would probably require more coding than a full adaptive mesh.