print also max gamma_f now

src/elasticity/assemble/IElasticity.cpp

@@ -83,13 +83,17 @@ void IElasticity::Finalize(Vectors &vecs) {
   GetInvariant(vecs[0], vecs[2]);
 }
 
-void IElasticity::CalculateAndPrintMinGammaF() const{
+void IElasticity::CalculateAndPrintMinMaxGammaF() const{
   if (gammaEvaluation.empty()) THROW("gammaEvaluation is empty, should not be used!")
+
+  double max_gamma = gammaEvaluation[0][0];
   double min_gamma = gammaEvaluation[0][0];
   for (VectorField v : gammaEvaluation){
     min_gamma = std::min(min_gamma, v[0]);
+    max_gamma = std::max(max_gamma, v[0]);
   }
   vout(0) << "minimal value of gamma_f is " << min_gamma << endl ;
+  vout(0) << "maximal value of gamma_f is " << max_gamma << endl ;
 }
 
 void IElasticity::PointEvaluationGamma(const std::string &varname, const vector<Point> &point)  {