- Build in CG tolerance
- Changed from Astar r * Astar r to r * Astar r in CG method
- Began to add stability terms parameter, but didn't get it added to the rest of the assembly routine
- Changed all fluid assembly for the trilinear term since several terms were missing
- Moved the count on total_solves
- Implemented rudimentary preconditioner