Skip to content
Snippets Groups Projects
Commit 8c1478ff authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

No commit message

No commit message
parent 3168e463
No related branches found
No related tags found
No related merge requests found
...@@ -203,6 +203,8 @@ static void copyMesh(GFace *source, GFace *target) ...@@ -203,6 +203,8 @@ static void copyMesh(GFace *source, GFace *target)
SVector3(mean_target.a,mean_target.b,mean_target.c)); SVector3(mean_target.a,mean_target.b,mean_target.c));
LINE = myLine(PLANE_SOURCE, PLANE_TARGET); LINE = myLine(PLANE_SOURCE, PLANE_TARGET);
// FIXME: this fails when the 2 planes have a common edge (= rotation axis)
// LINE is the axis of rotation // LINE is the axis of rotation
// let us compute the angle of rotation // let us compute the angle of rotation
count = 0; count = 0;
...@@ -227,7 +229,9 @@ static void copyMesh(GFace *source, GFace *target) ...@@ -227,7 +229,9 @@ static void copyMesh(GFace *source, GFace *target)
ANGLE = myAngle(t1, t2, LINE.t); ANGLE = myAngle(t1, t2, LINE.t);
else { else {
double ANGLE2 = myAngle(t1, t2, LINE.t); double ANGLE2 = myAngle(t1, t2, LINE.t);
if (fabs (ANGLE2-ANGLE) > 1.e-8) rotation = false; if (fabs (ANGLE2 - ANGLE) > 1.e-8){
rotation = false;
}
} }
count++; count++;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment