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

[mathieu . fabre]

Author: fabremathieu
Date: Mon Nov 10 14:30:28 2014
New Revision: 4808

URL: http://svn.gna.org/viewcvs/getfem?rev=4808&view=rev
Log:
final estimator

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=4808&r1=4807&r2=4808&view=diff
==============================================================================
--- trunk/getfem/src/getfem_error_estimate.cc   (original)
+++ trunk/getfem/src/getfem_error_estimate.cc   Mon Nov 10 14:30:28 2014
@@ -39,7 +39,7 @@
                               scalar_type mu,
                               scalar_type gamma0,
                               scalar_type f_coeff,
-                              scalar_type vertical_force,
+                              scalar_type vertical_force,
                               base_vector &ERR) {
    
     
@@ -263,8 +263,8 @@
          gmm::copy(up,sigt);
          gmm::scale(sigt, - sign);
          gmm::add(sig,sigt);
-         ERR[v.cv()] += coefficient *gmm::vect_norm2_sqr(sigt); 
-         eta3 +=  coefficient *gmm::vect_norm2_sqr(sigt);
+         ERR[v.cv()] += coefficient *radius*gmm::vect_norm2_sqr(sigt); 
+         eta3 +=  coefficient *radius*gmm::vect_norm2_sqr(sigt);
        }                 
        //  cout << "Erreur en contact sur " << v.cv() << " ERR[v.cv()] = " << 
ERR[v.cv()]-ee << endl;