Added MIC preconditioner, added boundaries

- Added modified incomplete cholesky preconditioner to CG solver
- Removed unstable timestep control
- Fixed out-of-bound access in matrix multiply
- Removed unnecessary comments
- Added solid boundary framework
- Added variational pressure solve