Added rhs_augmentation field to AugmentedLagrangian_basic class to avoid
repeated allocation/deallocation of temporary

I think in augmented_lagrangian_rhs() function the expression
div_ThUh*TminusaG creates a temporary field which should be deallocated
at every call to function.
Just looks that Augmented Lag with inertia is giving reasonable results
Before creating a new branch for refactoring
Major refactoring in classes