diff --git a/Numeric/gmshLaplace.cpp b/Numeric/gmshLaplace.cpp index 3eea76f3a8c6060e2efde2ffc8373311181a12de..28b06b0c54a4ee43c029c86d8b29c7f62380b2f3 100644 --- a/Numeric/gmshLaplace.cpp +++ b/Numeric/gmshLaplace.cpp @@ -57,10 +57,10 @@ void gmshLaplaceTerm::elementMatrix(MElement *e, gmshMatrix<double> &m) const //check for positive scheme // //printf("ELEM\n"); - for (int j = 0; j < nbNodes; j++){ - double sum = m(j,0)+m(j,1)+m(j,2); - if (sum < 0.0){ - printf("Sum LINE gmshLaplace Term NEG <0 %g %g %g\n", m(j,0),m(j,1),m(j,2)); +// for (int j = 0; j < nbNodes; j++){ +// double sum = m(j,0)+m(j,1)+m(j,2); +// if (sum < 0.0){ +// printf("Sum LINE gmshLaplace Term NEG <0 %g %g %g\n", m(j,0),m(j,1),m(j,2)); // for (int k = 0; k < nbNodes; k++) m(j,k) = -1.; // m(j,j) = (nbNodes-1); } diff --git a/benchmarks/stl/PelvisHEMI_CLASS_GEO.geo b/benchmarks/stl/PelvisHEMI_CLASS_GEO.geo index 7fe21a3b92136d13312e0b976c17b2fcc0b34374..6fdd4b88975437d05a63d6547ab048b9c5b8219b 100644 --- a/benchmarks/stl/PelvisHEMI_CLASS_GEO.geo +++ b/benchmarks/stl/PelvisHEMI_CLASS_GEO.geo @@ -1,4 +1,4 @@ -Mesh.CharacteristicLengthFactor=0.01; +Mesh.CharacteristicLengthFactor=0.03; Merge "PelvisHEMI_CLASS.msh"; CreateTopology;