before branching
working on explicit solve of advection terms...
before erasing stuff
advection contribution complete.
before modigying viscouse matrix
viscosity matrix and load vector assemble complete.
viscosity solve boundary condition complete
jacobi iteration implemented for viscous solve
boost:multi_array issue...
boost::multi_array no used any more.
seems that its good up to viscosity part
before presentation.
grid_matrix probably not right
seems working properly now!
before merge