Skip to content
Snippets Groups Projects
Commit 9ed42317 authored by Emilie Marchandise's avatar Emilie Marchandise
Browse files

No commit message

No commit message
parent 522567f5
No related branches found
No related tags found
No related merge requests found
...@@ -873,7 +873,7 @@ void Centerline::cutMesh(){ ...@@ -873,7 +873,7 @@ void Centerline::cutMesh(){
MVertex *v2 = lines[j]->getVertex(1); MVertex *v2 = lines[j]->getVertex(1);
SVector3 pt(v1->x(), v1->y(), v1->z()); SVector3 pt(v1->x(), v1->y(), v1->z());
SVector3 dir(v2->x()-v1->x(),v2->y()-v1->y(),v2->z()-v1->z()); SVector3 dir(v2->x()-v1->x(),v2->y()-v1->y(),v2->z()-v1->z());
printf("-> cut length (AR=%g) split %d\n", AR, nbSplit); printf("-> cut length (AR=%g) nbSplit=%d\n", AR, nbSplit);
cutByDisk(pt, dir, edges[i].maxRad); cutByDisk(pt, dir, edges[i].maxRad);
nbSplit--; nbSplit--;
lc = 0.0; lc = 0.0;
...@@ -916,7 +916,7 @@ void Centerline::cutByDisk(SVector3 &PT, SVector3 &NORM, double &maxRad){ ...@@ -916,7 +916,7 @@ void Centerline::cutByDisk(SVector3 &PT, SVector3 &NORM, double &maxRad){
double b = NORM.y(); double b = NORM.y();
double c = NORM.z(); double c = NORM.z();
double d = -a * PT.x() - b * PT.y() - c * PT.z(); double d = -a * PT.x() - b * PT.y() - c * PT.z();
printf("cut disk (R=%g)= %g %g %g %g \n", maxRad, a, b, c, d); //printf("cut disk (R=%g)= %g %g %g %g \n", maxRad, a, b, c, d);
const double EPS = 0.007; const double EPS = 0.007;
std::set<MEdge,Less_Edge> allEdges; std::set<MEdge,Less_Edge> allEdges;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment