[Getfem-commits] r4778 - /trunk/getfem/src/getfem_error_estimate.cc

[Yves . Renard]

Author: renard
Date: Tue Oct 14 14:46:00 2014
New Revision: 4778

URL: http://svn.gna.org/viewcvs/getfem?rev=4778&view=rev
Log:
minor fix

Modified:
    trunk/getfem/src/getfem_error_estimate.cc

Modified: trunk/getfem/src/getfem_error_estimate.cc
URL: 
http://svn.gna.org/viewcvs/getfem/trunk/getfem/src/getfem_error_estimate.cc?rev=4778&r1=4777&r2=4778&view=diff
==============================================================================
--- trunk/getfem/src/getfem_error_estimate.cc   (original)
+++ trunk/getfem/src/getfem_error_estimate.cc   Tue Oct 14 14:46:00 2014
@@ -252,15 +252,15 @@
          gmm::add(gmm::scaled(E, 2*mu), S1);  
          gmm::mult(S1, up, sig); // sig = sigma(u)n
          sign = gmm::vect_sp(sig,up);// sign = sigma_n(u)
-        Un = gmm::vect_sp(U1,up);// un = u_n
+          Un = gmm::vect_sp(U1,up);// un = u_n
          scal = Un-gamma*sign;
          if (scal<0)
            Pr = sign;
          else
            Pr = (scal/gamma + sign);
          
-         ERR[v.cv()] += coefficient*sqrt(gamma) *Pr*Pr; 
-         eta4 +=  coefficient*sqrt(gamma)* Pr*Pr;
+         ERR[v.cv()] += coefficient*gamma *Pr*Pr; 
+         eta4 +=  coefficient*gamma*Pr*Pr;
          
          gmm::copy(up,sigt);
          gmm::scale(sigt, - sign);