optimisation of predictor corrector

NonLinearSolver/internalPoints/ipNonLocalDamageGurson.cpp

@@ -24,7 +24,7 @@ IPNonLocalDamageGurson::IPNonLocalDamageGurson() : IPVariableMechanics(),  _elas
   _Fp(1,1)=1.;
   _Fp(2,2)=1.;
   _DFvStarDStrain=0;
-  _logsqrtCe *= 0.; 
+  _logsqrtCe = 0.; 
 };
 
 IPNonLocalDamageGurson::IPNonLocalDamageGurson(double fVinitial, const J2IsotropicHardening *j2IH,
@@ -58,7 +58,7 @@ IPNonLocalDamageGurson::IPNonLocalDamageGurson(double fVinitial, const J2Isotrop
   _Fp(1,1)=1.;
   _Fp(2,2)=1.;
   
-  _logsqrtCe *=0.;
+  _logsqrtCe =0.;
   
   _DFvStarDStrain=0;
 };

NonLinearSolver/materialLaw/mlawNonLocalDamageGurson.cpp

@@ -679,9 +679,13 @@ void mlawNonLocalDamageGurson::predictorCorrector(double tildefVstar, const STen
         for(int l=0; l<3; l++)
           S(i,j)+= Le(i,j,k,l)*kcor(k,l);
   
-  P_=Fe1;
-  P_*=S;
-  P_*=Fp1inv.transpose();
+  //P_=Fe1;
+  //P_*=S;
+  //P_*=Fp1inv.transpose();
+  
+  static STensor3 FeS;
+  STensorOperation::multSTensor3(Fe1,S,FeS);
+  STensorOperation::multSTensor3SecondTranspose(FeS,Fp1inv,P_);
 }