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);
_______________________________________________
Getfem-commits mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-commits
