Added the optimized first version of FLUSI.
It is operational and can:
 - time steppers: AB2, EI1, RK2
 - volume penalization for fixed obstacles
 - the routines are optimized for efficiency