else
{
Curve2d.Nullify(); // Pas de solution continue
-#if DEB
+#if PLATE_DEB
cout << "BuildPlateSurace :: Pas de projection continue" << endl;
#endif
}
if (Projector.NbCurves() != 1) {
HProjector.Nullify(); // Pas de solution continue
-#if DEB
+#if PLATE_DEB
cout << "BuildPlateSurace :: Pas de projection continue" << endl;
#endif
}
else
{
HProjector.Nullify(); // Pas de solution continue
-#if DEB
+#if PLATE_DEB
cout << "BuildPlateSurace :: Pas de projection complete" << endl;
#endif
}
myInitOrder->SetValue(l,l);
if (!CourbeJointive(myTol3d))
{// Standard_Failure::Raise("Curves are not joined");
-#ifdef DEB
+#ifdef PLATE_DEB
cout<<"WARNING : Courbes non jointives a "<<myTol3d<<" pres"<<endl;
#endif
}
do
{
-#if DEB
+#if PLATE_DEB
if (Affich && NbBoucle) {
cout<<"Resultats boucle"<< NbBoucle << endl;
cout<<"DistMax="<<myG0Error<<endl;
Fini = VerifSurface(NbBoucle);
if ((NbBoucle >= myNbIter)&&(!Fini))
{
-#ifdef DEB
+#ifdef PLATE_DEB
cout<<"Warning objectif non atteint"<<endl;
#endif
Fini = Standard_True;
VerifPoints(di,an,cu);
}
} while (!Fini); // Fin boucle pour meilleur surface
-#ifdef DEB
+#ifdef PLATE_DEB
if (NTLinCont != 0)
{ cout<<"======== Resultats globaux ==========="<<endl;
cout<<"DistMax="<<myG0Error<<endl;
Disc2dContour ( const Standard_Integer nbp,
TColgp_SequenceOfXY& Seq2d)
{
- if (nbp!=4) cout<<"nbp doit etre egal a 4 pour Disc2dContour"<<endl;
+#ifdef PLATE_DEB
+ if (nbp!=4)
+ cout<<"nbp doit etre egal a 4 pour Disc2dContour"<<endl;
+#endif
// initialisation
Seq2d.Clear();
const Standard_Integer iordre,
TColgp_SequenceOfXYZ& Seq3d)
{
- if (nbp!=4) cout<<"nbp doit etre egal a 4 pour Disc3dContour"<<endl;
- if (iordre!=0&&iordre!=1) cout<<"iordre incorrect pour Disc3dContour"<<endl;
+#ifdef PLATE_DEB
+ if (nbp!=4)
+ cout<<"nbp doit etre egal a 4 pour Disc3dContour"<<endl;
+ if (iordre!=0&&iordre!=1)
+ cout<<"iordre incorrect pour Disc3dContour"<<endl;
+#endif
// initialisation
Seq3d.Clear();
// echantillonnage en "cosinus" + 3 points sur chaque intervalle
} //if (isHalfSpace)
if (!isHalfSpace)
{
+#ifdef PLATE_DEB
cout<<endl<<"Normals are not in half space"<<endl<<endl;
+#endif
myIsLinear = Standard_False;
nopt = 2;
}
nopt = 1; //Calcul par la methode du plan d'inertie
else if (!CourbeJoint || NTLinCont != myNbBounds)
{// Standard_Failure::Raise("Curves are not joined");
-#ifdef DEB
+#ifdef PLATE_DEB
cout<<"WARNING : Courbes non jointives a "<<myTol3d<<" pres"<<endl;
#endif
nopt = 1;
}
}
}
-#if DEB
+#if PLATE_DEB
if (! myIsLinear)
cout <<"Metrics are too different :"<< Ratio<<endl;
#endif
{ int2d = Intersection.Point(k);
myLinCont->Value(i)->D0(int2d.ParamOnFirst(),P1);
myLinCont->Value(j)->D0(int2d.ParamOnSecond(),P2);
-#if DEB
+#if PLATE_DEB
if (Affich> 1)
{
cout << " Intersection "<< k << " entre " << i
if (A1>(PI/2))
A1= PI - A1;
if (Abs(Abs(A1)-PI)<myTolAng) Tol = 100000 * myTol3d;
-#if DEB
+#if PLATE_DEB
if (Affich) cout <<"Angle entre Courbe "<<i<<","<<j
<<" "<<Abs(Abs(A1)-PI)<<endl;
#endif
if (A1 > PI/2)
A1= PI - A1;
if (Abs(Abs(A1) - PI) < myTolAng) Tol = 100000 * myTol3d;
-#if DEB
+#if PLATE_DEB
if (Affich) cout <<"Angle entre Courbe "<<i<<","<<j
<<" "<<Abs(Abs(A1)-PI)<<endl;
#endif
{
coin = Ci.Resolution(Tol);
coin *= Angle / myTolAng * 10.;
+#if PLATE_DEB
cout<<endl<<"coin = "<<coin<<endl;
+#endif
Standard_Real Par1 = int2d.ParamOnFirst() - coin;
Standard_Real Par2 = int2d.ParamOnFirst() + coin;
// Stockage de l'intervalle pour la courbe i
{
coin = Cj.Resolution(Tol);
coin *= Angle / myTolAng * 10.;
+#if PLATE_DEB
cout<<endl<<"coin = "<<coin<<endl;
+#endif
Standard_Real Par1 = int2d.ParamOnSecond() - coin;
Standard_Real Par2 = int2d.ParamOnSecond() + coin;
// Stockage de l'intervalle pour la courbe j
Append( int2d.ParamOnSecond() + tolint);
}
-#ifdef DEB
+#ifdef PLATE_DEB
cout<<"Attention: Deux points 3d ont la meme projection dist="
<<Dist<<endl;
#endif
NbPtInter= PntInter->Value(i).Length();
NbPtG1G1= PntG1G1->Value(i).Length();
-#if DEB
+#if PLATE_DEB
if (Affich > 1) {
cout << "Courbe : " << i << endl;
cout << " NbPnt, NbPtInter, NbPtG1G1 :" << NbPnt_i << ", "
mySurface = AppPlate.Surface(1);
myAppError = AppPlate.MaxError(3,1);
myCritError = AppPlate.CritError(3,1);
-#if DEB
+#if PLATE_DEB
cout<<"Approximation results"<<endl;
cout<<" Approximation error : "<<myAppError<<endl;
cout<<" Criterium error : "<<myCritError<<endl;
Standard_Real seuil = Tol3d;
if (CritOrder==0&&Tol3d<10*dmax) {
seuil=10*dmax;
-#if DEB
+#if PLATE_DEB
cout<<"Seuil G0 choisi trop faible par rapport au contour. On prend "<<seuil<<endl;
#endif
}
if (CritOrder==1&&Tol3d<10*dmax) {
seuil=10*dmax;
-#if DEB
+#if PLATE_DEB
cout<<"Seuil G1 choisi trop faible par rapport au contour. On prend "<<seuil<<endl;
#endif
}
mySurface = AppPlate.Surface(1);
myAppError = AppPlate.MaxError(3,1);
myCritError = 0.;
-#if DEB
+#if PLATE_DEB
cout<<"Approximation results"<<endl;
cout<<" Approximation error : "<<myAppError<<endl;
#endif
mySurface = AppPlate.Surface(1);
myAppError = AppPlate.MaxError(3,1);
myCritError = AppPlate.CritError(3,1);
-#if DEB
+#if PLATE_DEB
cout<<"Approximation results"<<endl;
cout<<" Approximation error : "<<myAppError<<endl;
cout<<" Criterium error : "<<myCritError<<endl;
mySurface = AppPlate.Surface(1);
myAppError = AppPlate.MaxError(3,1);
myCritError = AppPlate.CritError(3,1);
-#if DEB
+#if PLATE_DEB
cout<<"Approximation results"<<endl;
cout<<" Approximation error : "<<myAppError<<endl;
cout<<" Criterium error : "<<myCritError<<endl;