Can treat advection implicitly

v0.7

Many small cosmetic improvements.

Computes nodal velocities

This release computes nodal velocities (computed in function "nodal_uvw") which will be used for Lagrangian particle tracking.

Periodicity seems to work well too, but I didn't test it for collocated variable arrangement.

Matrix-vector product uses boundary values

Added function "mat_vec_bnd", matrix-vector product with boundary values. It is quite fundamental, really. It is an important landmark towards the implementation of periodic boundary condition, it may prove to be a step towards parallelization (boundary values might be buffer cells), it simplified function "create_matrix" a lot, and the code doesn't depend any longer on addition source term with boundary values.

Removed the notorious "all.py" from the code

v0.4

Removed the notorious "all.py" module everywhere, and replaced it wit…

Re-indented all sources to four

v0.3

Framed the main function.

Tests are in a loop and in its own directory

v0.2

Added a little script to clean the "__pycache__" directories.

Has its own linear solvers

v0.1

SciPy's sparse module is not used any more.