diff --git a/Numeric/gmshCrossConf.cpp b/Numeric/gmshCrossConf.cpp index 6fd73a4340027ef82e87423fdf7d10d29b9aaf97..32ef8a02ddee944c6b0a76057c1512509e8e9c5b 100644 --- a/Numeric/gmshCrossConf.cpp +++ b/Numeric/gmshCrossConf.cpp @@ -40,17 +40,10 @@ void gmshCrossConfTerm::elementMatrix(MElement *e, gmshMatrix<double> &m) const } SVector3 N (jac[2][0],jac[2][1],jac[2][2]); - //SVector3 N (jac[0][2],jac[1][2],jac[2][2]); - - double norm = N[0]*N[0]+ N[1]*N[1]+ N[2]*N[2]; - if (norm != 1.0){ - printf("Normal N=%g %g %g NORM = %g\n", N[0], N[1], N[2], norm); - //exit(1); - } for (int j = 0; j < nbNodes; j++){ for (int k = 0; k <= j; k++){ - m(j, k) += dot(crossprod(Grads[j],Grads[k]),N) * weight * detJ * _diff; + m(j, k) += 1.0*dot(crossprod(Grads[j],Grads[k]),N) * weight * detJ * _diff; } } } diff --git a/benchmarks/3d/Torus_GEO.geo b/benchmarks/3d/Torus_GEO.geo index b208bc85cc33d8186bdee675f2ac0b851102fa5d..c0dc7f1354e026c6701d2a155880563e5feb4fd9 100644 --- a/benchmarks/3d/Torus_GEO.geo +++ b/benchmarks/3d/Torus_GEO.geo @@ -1,5 +1,4 @@ -Mesh.CharacteristicLengthFactor=0.2; - +Mesh.CharacteristicLengthFactor=1.0; Merge "Torus2_CLASS.msh"; CreateTopology;