Added support of non-cubic domains
added support of different boundary conditions. FREE_SLIP not implemented yet.
Treating the boundary is slow.