1 // Created on: 1996-09-09
2 // Created by: Herve LOUESSARD
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
6 // This file is part of Open CASCADE Technology software library.
8 // This library is free software; you can redistribute it and/or modify it under
9 // the terms of the GNU Lesser General Public License version 2.1 as published
10 // by the Free Software Foundation, with special exception defined in the file
11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
12 // distribution for complete text of the license and disclaimer of any warranty.
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
18 #include <Geom2d_Curve.hxx>
19 #include <Geom_Surface.hxx>
20 #include <GeomLProp_SLProps.hxx>
23 #include <gp_Pnt2d.hxx>
25 #include <LocalAnalysis_SurfaceContinuity.hxx>
26 #include <StdFail_NotDone.hxx>
28 /*********************************************************************************/
29 /*********************************************************************************/
30 void LocalAnalysis_SurfaceContinuity::SurfC0 (const GeomLProp_SLProps& Surf1,
31 const GeomLProp_SLProps& Surf2
33 { myContC0=(Surf1.Value()).Distance(Surf2.Value());
36 /*********************************************************************************/
38 void LocalAnalysis_SurfaceContinuity::SurfC1( GeomLProp_SLProps& Surf1,
39 GeomLProp_SLProps & Surf2)
40 { gp_Vec V1u, V2u, V1v, V2v;
41 Standard_Real norm1u, norm2u, norm1v, norm2v,angu,angv;
49 norm1u = V1u.Magnitude();
50 norm2u = V2u.Magnitude();
51 norm1v = V1v.Magnitude();
52 norm2v = V2v.Magnitude();
54 if ((norm1u>myepsnul )&&(norm2u>myepsnul )&&(norm1v>myepsnul)
56 { if (norm1u >= norm2u )
57 myLambda1U= norm2u / norm1u;
58 else myLambda1U = norm1u / norm2u;
59 if (norm1v >= norm2v )
60 myLambda1V= norm2v / norm1v;
61 else myLambda1V = norm1v / norm2v;
63 if (angu>M_PI/2) myContC1U=M_PI-angu;
66 if (angv>M_PI/2) myContC1V=M_PI-angv;
67 else myContC1V=angv; }
68 else {myIsDone = Standard_False;
69 myErrorStatus=LocalAnalysis_NullFirstDerivative;}
72 /*********************************************************************************/
74 void LocalAnalysis_SurfaceContinuity::SurfC2( GeomLProp_SLProps& Surf1,
75 GeomLProp_SLProps& Surf2)
77 { gp_Vec V11u, V12u, V21u, V22u, V11v, V12v, V21v, V22v;
78 Standard_Real norm11u, norm12u, norm21u, norm22u, norm11v, norm12v, norm21v, norm22v;
84 norm11u = V11u.Magnitude();
85 norm12u = V12u.Magnitude();
86 norm21u = V21u.Magnitude();
87 norm22u = V22u.Magnitude();
89 if ((norm11u>myepsnul)&&(norm12u>myepsnul))
90 { if( (norm21u>myepsnul)&&(norm22u>myepsnul))
91 { if (norm11u >= norm12u )
92 {myLambda1U= norm12u / norm11u;
93 myLambda2U = norm22u /norm21u;}
94 else {myLambda1U = norm11u / norm12u;
95 myLambda2U = norm21u / norm22u;}
97 if(ang>M_PI/2) myContC2U=M_PI-ang;
100 { myIsDone=Standard_False;
101 myErrorStatus=LocalAnalysis_NullSecondDerivative;}
104 else { myIsDone=Standard_False;
105 myErrorStatus=LocalAnalysis_NullFirstDerivative;}
111 norm11v = V11v.Magnitude();
112 norm12v = V12v.Magnitude();
113 norm21v = V21v.Magnitude();
114 norm22v = V22v.Magnitude();
116 if ((norm11v>myepsnul)&&(norm12v>myepsnul))
117 {if ((norm21v>myepsnul)&&(norm22v>myepsnul))
118 {if ( norm11v >= norm12v )
119 {myLambda1V= norm12v / norm11v;
120 myLambda2V= norm22v / norm21v;}
121 else{ myLambda1V = norm11v / norm12v;
122 myLambda2V = norm21v / norm22v;}
123 ang= V21v.Angle(V22v);
124 if (ang>M_PI/2) myContC2V=M_PI-ang;
127 else{ myIsDone= Standard_False;
128 myErrorStatus=LocalAnalysis_NullSecondDerivative; }
130 else{ myIsDone=Standard_False;
131 myErrorStatus=LocalAnalysis_NullFirstDerivative;}
134 /*********************************************************************************/
135 void LocalAnalysis_SurfaceContinuity::SurfG1( GeomLProp_SLProps& Surf1,
136 GeomLProp_SLProps& Surf2)
137 { if (Surf1.IsNormalDefined()&&Surf2.IsNormalDefined())
138 { gp_Dir D1 = Surf1.Normal();
139 gp_Dir D2 = Surf2.Normal();
140 Standard_Real ang=D1.Angle(D2);
141 if (ang>M_PI/2) myContG1= M_PI-ang;
144 else{ myIsDone=Standard_False;
145 myErrorStatus=LocalAnalysis_NormalNotDefined;}
148 /*********************************************************************************/
150 void LocalAnalysis_SurfaceContinuity::SurfG2 ( GeomLProp_SLProps& Surf1,
151 GeomLProp_SLProps & Surf2)
152 { gp_Dir DMIN1, DMIN2, DMAX1, DMAX2;
153 Standard_Real RMIN1,RMIN2,RMAX1,RMAX2;
154 Standard_Real x1, x2, y1, y2, z1, z2;
156 if ( Surf1.IsCurvatureDefined() && Surf2.IsCurvatureDefined())
158 Surf1.CurvatureDirections(DMIN1, DMAX1);
159 Surf2.CurvatureDirections(DMIN2, DMAX2);
160 DMIN1.Coord(x1, y1, z1);
161 DMAX1.Coord(x2, y2, z2);
162 gp_Dir MCD1( (Abs(x1)+Abs(x2))/2, (Abs(y1)+Abs(y2))/2, (Abs(z1)+Abs(z2))/2 );
163 DMIN2.Coord(x1, y1, z1);
164 DMAX2.Coord(x2, y2, z2);
165 gp_Dir MCD2( (Abs(x1)+Abs(x2))/2, (Abs(y1)+Abs(y2))/2, (Abs(z1)+Abs(z2))/2 );
167 myAlpha = MCD1.Angle( MCD2 );
168 RMIN1 = Surf1.MinCurvature();
169 RMAX1 = Surf1.MaxCurvature();
170 RMIN2 = Surf2.MinCurvature();
171 RMAX2 = Surf2.MaxCurvature();
172 myETA1 = (RMIN1+RMAX1)/2;
173 myETA2 = (RMIN2+RMAX2)/2;
174 myETA = (myETA1+myETA2)/2;
175 myZETA1 = (RMAX1-RMIN1)/2;
176 myZETA2 = (RMAX2-RMIN2)/2;
177 myZETA = (myZETA1+myZETA2)/2;
178 Standard_Real DETA,DZETA;
179 DETA = ( myETA1- myETA2)/2;
180 DZETA =( myZETA1- myZETA2)/2;
181 myGap= Abs(DETA) + sqrt( DZETA*DZETA*Cos(myAlpha)*Cos(myAlpha)
182 + myZETA*myZETA*Sin(myAlpha)*Sin(myAlpha));
184 else {myIsDone = Standard_False;
185 myErrorStatus=LocalAnalysis_CurvatureNotDefined;}
188 LocalAnalysis_SurfaceContinuity::LocalAnalysis_SurfaceContinuity(const Standard_Real EpsNul,
189 const Standard_Real EpsC0,
190 const Standard_Real EpsC1,
191 const Standard_Real EpsC2,
192 const Standard_Real EpsG1,
193 const Standard_Real Percent,
194 const Standard_Real Maxlen)
220 myIsDone = Standard_True;
222 void LocalAnalysis_SurfaceContinuity::ComputeAnalysis(GeomLProp_SLProps& Surf1,
223 GeomLProp_SLProps& Surf2,
224 const GeomAbs_Shape Order)
225 { myTypeCont = Order;
227 { case GeomAbs_C0 : {
228 SurfC0(Surf1, Surf2);
232 SurfC0(Surf1, Surf2);
233 SurfC1(Surf1, Surf2);}
236 SurfC0(Surf1, Surf2);
237 SurfC1(Surf1, Surf2);
238 SurfC2(Surf1, Surf2);}
241 SurfC0(Surf1, Surf2);
242 SurfG1(Surf1, Surf2);}
245 SurfC0(Surf1, Surf2);
246 SurfG1(Surf1, Surf2);
247 SurfG2(Surf1, Surf2);}
254 /*********************************************************************************/
256 LocalAnalysis_SurfaceContinuity::LocalAnalysis_SurfaceContinuity( const Handle(Geom_Surface)& Surf1,
257 const Standard_Real u1, const Standard_Real v1,
258 const Handle(Geom_Surface)& Surf2, const Standard_Real u2,
259 const Standard_Real v2, const GeomAbs_Shape Ordre,
260 const Standard_Real EpsNul,
261 const Standard_Real EpsC0 ,
262 const Standard_Real EpsC1 ,
263 const Standard_Real EpsC2,
264 const Standard_Real EpsG1,
265 const Standard_Real Percent,
266 const Standard_Real Maxlen )
285 { myTypeCont = Ordre;
293 myIsDone = Standard_True;
295 { case GeomAbs_C0 : {
296 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 0, myepsnul);
297 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 0, myepsnul);
298 SurfC0(Surfa1, Surfa2);
302 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 1, myepsnul);
303 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 1, myepsnul);
304 SurfC0(Surfa1, Surfa2);
305 SurfC1(Surfa1, Surfa2);}
308 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 2, myepsnul);
309 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 2, myepsnul);
310 SurfC0(Surfa1, Surfa2);
311 SurfC1(Surfa1, Surfa2);
312 SurfC2(Surfa1, Surfa2);}
315 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 1, myepsnul);
316 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 1, myepsnul);
317 SurfC0(Surfa1, Surfa2);
318 SurfG1(Surfa1, Surfa2);}
321 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 2, myepsnul);
322 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 2, myepsnul);
323 SurfC0(Surfa1, Surfa2);
324 SurfG1(Surfa1, Surfa2);
325 SurfG2(Surfa1, Surfa2);}
331 /*********************************************************************************/
333 LocalAnalysis_SurfaceContinuity::LocalAnalysis_SurfaceContinuity(const Handle(Geom2d_Curve)& curv1,
334 const Handle(Geom2d_Curve)& curv2, const Standard_Real U,
335 const Handle(Geom_Surface)& Surf1 ,
336 const Handle(Geom_Surface)& Surf2,
337 const GeomAbs_Shape Ordre,
338 const Standard_Real EpsNul,
339 const Standard_Real EpsC0 ,
340 const Standard_Real EpsC1 ,
341 const Standard_Real EpsC2,
342 const Standard_Real EpsG1,
343 const Standard_Real Percent,
344 const Standard_Real Maxlen )
363 { Standard_Real pard1, parf1, pard2, parf2, u1, v1, u2, v2;
373 myIsDone = Standard_True;
375 pard1 = curv1->FirstParameter();
376 pard2 = curv2->FirstParameter();
377 parf1 = curv1->LastParameter();
378 parf2 = curv2->LastParameter();
381 if (!(((U <= parf1) && (U >= pard1)) &&((U <= parf2) && (U >= pard2))) ) myIsDone = Standard_False;
383 { gp_Pnt2d pt1 = curv1->Value(U);
384 gp_Pnt2d pt2 = curv2->Value(U);
389 { case GeomAbs_C0 : {
390 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 0, myepsnul);
391 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 0, myepsnul);
392 SurfC0(Surfa1, Surfa2);
396 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 1, myepsnul);
397 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 1, myepsnul);
398 SurfC0(Surfa1, Surfa2);
399 SurfC1(Surfa1, Surfa2);}
402 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 2, myepsnul);
403 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 2, myepsnul);
404 SurfC0(Surfa1, Surfa2);
405 SurfC1(Surfa1, Surfa2);
406 SurfC2(Surfa1, Surfa2);}
409 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 1, myepsnul);
410 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 1, myepsnul);
411 SurfC0(Surfa1, Surfa2);
412 SurfG1(Surfa1, Surfa2);}
415 GeomLProp_SLProps Surfa1 ( Surf1, u1, v1, 2, myepsnul);
416 GeomLProp_SLProps Surfa2 ( Surf2, u2, v2, 2, myepsnul);
417 SurfC0(Surfa1, Surfa2);
418 SurfG1(Surfa1, Surfa2);
419 SurfG2(Surfa1, Surfa2);}
426 /*********************************************************************************/
428 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsC0() const
429 { if (!myIsDone) { throw StdFail_NotDone();}
430 if ( myContC0 <= myepsC0 )
431 return Standard_True;
432 else return Standard_False;
435 /*********************************************************************************/
437 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsC1() const
438 { if (!myIsDone) { throw StdFail_NotDone();}
439 if ( IsC0 () && (myContC1U <= myepsC1) && (myContC1V <= myepsC1))
440 return Standard_True;
441 else return Standard_False;
444 /*********************************************************************************/
446 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsC2() const
447 { Standard_Real eps1u, eps1v, eps2u, eps2v;
449 if (!myIsDone) { throw StdFail_NotDone();}
451 { eps1u = 0.5*myepsC1*myepsC1*myLambda1U;
452 eps1v = 0.5*myepsC1*myepsC1*myLambda1V;
453 eps2u = 0.5*myepsC2*myepsC2*myLambda2U;
454 eps2v = 0.5*myepsC2*myepsC2*myLambda2V;
455 if ((myContC2U < myepsC2) && (myContC2V < myepsC2))
456 { if (Abs(myLambda1U*myLambda1U-myLambda2U) <= (eps1u*eps1u+eps2u))
457 if (Abs(myLambda1V*myLambda1V-myLambda2V) <= (eps1v*eps1v+eps2v))
458 return Standard_True;
459 else return Standard_False;
460 else return Standard_False;
462 else return Standard_False;
464 else return Standard_False;
468 /*********************************************************************************/
470 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsG1() const
472 if (!myIsDone) { throw StdFail_NotDone();}
473 if ( IsC0 () &&( myContG1 <= myepsG1))
474 return Standard_True;
475 else return Standard_False;
478 /*********************************************************************************/
480 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsG2()const
481 { Standard_Real EPSNL;
482 Standard_Integer itype;
484 if (!myIsDone) { throw StdFail_NotDone();}
486 EPSNL= 8*myepsC0/(mymaxlen*mymaxlen);
488 { if ( ( Abs(myETA)< EPSNL) && ( Abs(myZETA)< EPSNL))
489 return Standard_True;
490 if ( ( Abs(myZETA1)< EPSNL) && ( Abs(myZETA2)< EPSNL))
492 else if ( ( Abs(myETA1)< EPSNL) && ( Abs(myETA2)< EPSNL))
494 else if ((Abs(Abs(myZETA)-Abs(myETA))) < EPSNL)
496 else if ((myETA1<myZETA1)&&(myETA2<myZETA2))
498 else if ((myETA1>myZETA1)&&(myETA2>myZETA2))
503 if (( myETA >= (2*myZETA))&&(myGap<=(myperce*(myETA-myZETA)))) return Standard_True;
504 if (( myZETA>= myETA) && ( myGap<= (myperce*myZETA))) return Standard_True;
505 if (( myZETA<= myETA) && ( myETA <= (2*myZETA)) && (myGap <= (myperce * myETA)))
506 return Standard_True;
507 else return Standard_False;
510 else return Standard_False;
512 else return Standard_False;
517 /*********************************************************************************/
519 GeomAbs_Shape LocalAnalysis_SurfaceContinuity::ContinuityStatus() const
521 if (!myIsDone) { throw StdFail_NotDone();}
525 /*********************************************************************************/
527 Standard_Real LocalAnalysis_SurfaceContinuity::C0Value() const
529 if (!myIsDone) { throw StdFail_NotDone();}
533 /*********************************************************************************/
535 Standard_Real LocalAnalysis_SurfaceContinuity::C1UAngle() const
537 if (!myIsDone) { throw StdFail_NotDone();}
538 return ( myContC1U );
541 /*********************************************************************************/
543 Standard_Real LocalAnalysis_SurfaceContinuity::C1VAngle() const
545 if (!myIsDone) { throw StdFail_NotDone();}
546 return ( myContC1V );
549 /*********************************************************************************/
551 Standard_Real LocalAnalysis_SurfaceContinuity::C2UAngle() const
553 if (!myIsDone) { throw StdFail_NotDone();}
554 return ( myContC2U );
557 /*********************************************************************************/
559 Standard_Real LocalAnalysis_SurfaceContinuity::C2VAngle() const
561 if (!myIsDone) { throw StdFail_NotDone();}
562 return ( myContC2V );
565 /*********************************************************************************/
567 Standard_Real LocalAnalysis_SurfaceContinuity::G1Angle() const
569 if (!myIsDone) { throw StdFail_NotDone();}
573 /*********************************************************************************/
575 Standard_Real LocalAnalysis_SurfaceContinuity::C1URatio() const
577 if (!myIsDone) { throw StdFail_NotDone();}
578 return ( myLambda1U );
581 /*********************************************************************************/
583 Standard_Real LocalAnalysis_SurfaceContinuity::C2URatio() const
585 if (!myIsDone) { throw StdFail_NotDone();}
586 return ( myLambda2U );
589 /*********************************************************************************/
591 Standard_Real LocalAnalysis_SurfaceContinuity::C1VRatio() const
593 if (!myIsDone) { throw StdFail_NotDone();}
594 return ( myLambda1V );
597 /*********************************************************************************/
599 Standard_Real LocalAnalysis_SurfaceContinuity::C2VRatio() const
601 if (!myIsDone) { throw StdFail_NotDone();}
602 return ( myLambda2V );
605 /*********************************************************************************/
607 Standard_Real LocalAnalysis_SurfaceContinuity::G2CurvatureGap() const
609 if (!myIsDone) { throw StdFail_NotDone();}
613 /*********************************************************************************/
615 Standard_Boolean LocalAnalysis_SurfaceContinuity::IsDone() const
616 { return ( myIsDone );
619 /*********************************************************************************/
620 LocalAnalysis_StatusErrorType LocalAnalysis_SurfaceContinuity::StatusError() const
622 return myErrorStatus;