Commit f004189c authored by Alexandre Halbach's avatar Alexandre Halbach

Check geo

parent 79bf0e57
......@@ -125,13 +125,16 @@ void rawarc::mesh(void)
///// Rotate the arc plane to have it parallel to the xy plane:
std::vector<double> xyrot = geotools::getplaneangles(pccoords, p1coords, p2coords);
sons[0]->rotate(-xyrot[1],xyrot[0],0);
sons[1]->rotate(-xyrot[1],xyrot[0],0);
mycenterpoint->rotate(-xyrot[1],xyrot[0],0);
if (xyrot[0] != 0 || xyrot[1] != 0)
{
sons[0]->rotate(-xyrot[1],xyrot[0],0);
sons[1]->rotate(-xyrot[1],xyrot[0],0);
mycenterpoint->rotate(-xyrot[1],xyrot[0],0);
p1coords = *(sons[0]->getcoords());
p2coords = *(sons[1]->getcoords());
pccoords = *(mycenterpoint->getcoords());
p1coords = *(sons[0]->getcoords());
p2coords = *(sons[1]->getcoords());
pccoords = *(mycenterpoint->getcoords());
}
///// Rotated
......@@ -163,6 +166,7 @@ void rawarc::mesh(void)
///// Rotate everything back:
rotate(xyrot[1],-xyrot[0],0);
if (xyrot[0] != 0 || xyrot[1] != 0)
rotate(xyrot[1],-xyrot[0],0);
}
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment