double ang2 = fabs(VN2.Angle(VI2));
double coef1 = 0.5 - (( ang1 < M_PI/3. ) ? cos(ang1)*0.25 : 0 );
double coef2 = 0.5 - (( ang2 < M_PI/3. ) ? cos(ang2)*0.25 : 0 ); // cos(ang2) ?
-// double coef2 = 0.5;
-// if(ang2<PI/3)
-// coef2 -= cos(ang1)*0.25;
+ // double coef2 = 0.5;
+ // if(ang2<PI/3)
+ // coef2 -= cos(ang1)*0.25;
VN1.Scale(coef1);
VN2.Scale(coef2);
MergeAdjacent( PrmI, nodesToMove );
MergeAdjacent( PrmJ, nodesToMove );
+ apexI = PrmI->GetNode( PYRAM_APEX ); // apexI can be removed by merge
+
} // end if(hasInt)
- } // loop on suspectPyrams
- } // loop on 4 base nodes of PrmI
+ } // loop on 4 base nodes of PrmI
+ } // loop on suspectPyrams
} // loop on all pyramids