0024002: Overall code and build procedure refactoring -- automatic
[occt.git] / src / BRepBlend / BRepBlend_AppSurface.cxx
1 // Created on: 1996-11-26
2 // Created by: Philippe MANGIN
3 // Copyright (c) 1996-1999 Matra Datavision
4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
5 //
6 // This file is part of Open CASCADE Technology software library.
7 //
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.
13 //
14 // Alternatively, this file may be used under the terms of Open CASCADE
15 // commercial license or contractual agreement.
16
17
18 #include <Approx_SweepFunction.hxx>
19 #include <BRepBlend_AppSurface.hxx>
20 #include <Standard_DomainError.hxx>
21 #include <Standard_OutOfRange.hxx>
22 #include <StdFail_NotDone.hxx>
23 #include <TColgp_Array1OfPnt.hxx>
24 #include <TColgp_Array1OfPnt2d.hxx>
25 #include <TColgp_Array1OfVec.hxx>
26 #include <TColgp_Array1OfVec2d.hxx>
27 #include <TColStd_Array1OfReal.hxx>
28
29 BRepBlend_AppSurface::BRepBlend_AppSurface(
30                                Handle(Approx_SweepFunction)& Func,
31                                const Standard_Real First, 
32                                const Standard_Real Last,
33                                const Standard_Real Tol3d, 
34                                const Standard_Real Tol2d,
35                                const Standard_Real TolAngular,
36                                const GeomAbs_Shape Continuity,
37                                const Standard_Integer Degmax, 
38                                const Standard_Integer Segmax) :
39                                approx(Func)
40 {
41    Standard_Integer Nb2d = Func->Nb2dCurves();
42    Standard_Integer NbPolSect, NbKnotSect, udeg;
43    GeomAbs_Shape   continuity = Continuity;
44
45 // (1) Verification de la possibilite de derivation
46   if (continuity != GeomAbs_C0) {
47     if (Nb2d == 0) Nb2d =1;
48     Func->SectionShape(NbPolSect, NbKnotSect, udeg);
49     TColStd_Array1OfReal W  (1, NbPolSect);
50     TColgp_Array1OfPnt   P  (1, NbPolSect);
51     TColgp_Array1OfPnt2d P2d(1, Nb2d);
52     TColgp_Array1OfVec   V  (1, NbPolSect);
53     TColgp_Array1OfVec2d V2d(1, Nb2d);
54     Standard_Boolean Ok;
55     if (continuity == GeomAbs_C2) {
56         Ok = Func->D2( First, First, Last, P, V, V, P2d, V2d, V2d, W, W, W);
57         if (!Ok) { continuity = GeomAbs_C1; }
58       }
59     if (continuity == GeomAbs_C1) {
60       Ok = (Func->D1(First, First, Last, P, V, P2d, V2d, W, W));
61       if (!Ok) { continuity = GeomAbs_C0; }
62     }
63   }
64   
65
66 // (2) Approximation
67    approx.Perform(First, Last, 
68                   Tol3d,  Tol3d,  Tol2d,  TolAngular, 
69                   continuity, Degmax, Segmax);    
70  }
71
72 void BRepBlend_AppSurface::SurfShape (Standard_Integer& UDegree,
73                                   Standard_Integer& VDegree,
74                                   Standard_Integer& NbUPoles,
75                                   Standard_Integer& NbVPoles,
76                                   Standard_Integer& NbUKnots,
77                                   Standard_Integer& NbVKnots) const
78 {
79   approx.SurfShape(UDegree, VDegree, 
80                    NbUPoles, NbVPoles, 
81                    NbUKnots,NbVKnots); 
82 }
83
84
85 void BRepBlend_AppSurface::Surface(TColgp_Array2OfPnt& TPoles,
86                                TColStd_Array2OfReal& TWeights,
87                                TColStd_Array1OfReal& TUKnots,
88                                TColStd_Array1OfReal& TVKnots,
89                                TColStd_Array1OfInteger& TUMults,
90                                TColStd_Array1OfInteger& TVMults) const
91
92 {
93   approx.Surface(TPoles, TWeights, TUKnots,TVKnots, TUMults,TVMults);
94 }
95
96
97 Standard_Real BRepBlend_AppSurface::MaxErrorOnSurf() const 
98 {
99   return approx.MaxErrorOnSurf();
100 }
101
102
103 void BRepBlend_AppSurface::Curves2dShape(Standard_Integer& Degree,
104                                    Standard_Integer& NbPoles,
105                                    Standard_Integer& NbKnots) const
106 {
107   approx.Curves2dShape( Degree, NbPoles, NbKnots);
108 }
109
110 void BRepBlend_AppSurface::Curve2d(const Standard_Integer Index,
111                                TColgp_Array1OfPnt2d& TPoles,
112                                TColStd_Array1OfReal& TKnots,
113                                TColStd_Array1OfInteger& TMults) const
114 {
115   approx.Curve2d(Index, TPoles, TKnots, TMults);
116 }
117
118 Standard_Real BRepBlend_AppSurface::Max2dError(const Standard_Integer Index) const 
119 {
120   return approx.Max2dError(Index);
121 }
122
123
124 Standard_Real BRepBlend_AppSurface::TolCurveOnSurf(const Standard_Integer Index) const 
125 {
126   return approx.TolCurveOnSurf(Index);
127 }
128
129 inline void  BRepBlend_AppSurface::TolReached (Standard_Real& Tol3d,
130                                          Standard_Real& Tol2d) const
131 {
132   Tol3d = approx.MaxErrorOnSurf();
133   Tol2d = 0;
134   for (Standard_Integer ii=1; ii<=approx.NbCurves2d(); ii++) {
135     Tol2d = Max(Tol2d, approx.Max2dError(ii));
136   }
137 }
138
139 void BRepBlend_AppSurface::Dump(Standard_OStream& o) const 
140 {
141   approx.Dump(o);
142 }