From 2c4d57009d38376d888ddfd902f4db4c29d57a7e Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Jean-Fran=C3=A7ois=20Remacle=20=28students=29?=
<jean-francois.remacle@uclouvain.be>
Date: Thu, 10 Sep 2009 08:02:06 +0000
Subject: [PATCH] *** empty log message ***
---
Numeric/gmshCrossConf.cpp | 9 +--------
benchmarks/3d/Torus_GEO.geo | 3 +--
2 files changed, 2 insertions(+), 10 deletions(-)
diff --git a/Numeric/gmshCrossConf.cpp b/Numeric/gmshCrossConf.cpp
index 6fd73a4340..32ef8a02dd 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 b208bc85cc..c0dc7f1354 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;
--
GitLab