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