1 // Created on: 1998-03-23
2 // Created by: Jean Yves LEBEY
3 // Copyright (c) 1998-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 #include <gp_Vec2d.hxx>
18 #include <TopoDS.hxx>
19 #include <TopExp.hxx>
20 #include <TopOpeBRepTool_2d.hxx>
21 #include <BRep_Tool.hxx>
22 #include <Geom_Curve.hxx>
23 #include <ProjLib_ProjectedCurve.hxx>
24 #include <Geom_Surface.hxx>
25 #include <BRepAdaptor_Surface.hxx>
26 #include <BRepAdaptor_Curve.hxx>
27 #include <GeomAdaptor_Curve.hxx>
28 #include <Geom_TrimmedCurve.hxx>
29 #include <TopOpeBRepTool_CurveTool.hxx>
30 #include <TopOpeBRepTool_DataMapOfShapeListOfC2DF.hxx>
31 #include <TopOpeBRepTool_C2DF.hxx>
32 #include <TopOpeBRepTool_ListOfC2DF.hxx>
33 #include <TopOpeBRepTool_ListIteratorOfListOfC2DF.hxx>
34 #include <TopOpeBRepTool_tol.hxx>
35 #include <TopOpeBRepTool_EXPORT.hxx>
36 #include <TopOpeBRepTool_TOOL.hxx>
37
38 #ifdef DRAW
39 #include <TopOpeBRepTool_DRAW.hxx>
40 #include <Geom2d_Curve.hxx>
41 #endif
42
43 #ifdef OCCT_DEBUG
debc2dnull(void)44 void debc2dnull(void) {}
45 Standard_EXPORT Standard_Boolean TopOpeBRepTool_GettraceC2D();
46 #endif
47
48 // structure e -> C2D/F
49 static TopOpeBRepTool_DataMapOfShapeListOfC2DF *GLOBAL_pmosloc2df = NULL;
50 static Standard_Integer GLOBAL_C2D_i = 0; // DEB
51
52 // structure ancetre
53 static TopTools_IndexedDataMapOfShapeListOfShape *GLOBAL_pidmoslosc2df = NULL;
54 static TopoDS_Face *GLOBAL_pFc2df = NULL;
55 static TopoDS_Shape *GLOBAL_pS1c2df = NULL;
56 static TopoDS_Shape *GLOBAL_pS2c2df = NULL;
57
58 Standard_EXPORT Handle(Geom2d_Curve) MakePCurve(const ProjLib_ProjectedCurve& PC);
59
60 // ------------------------------------------------------------------------------------
FC2D_FancestorE(const TopoDS_Edge & E)61 static const TopoDS_Face& FC2D_FancestorE(const TopoDS_Edge& E)
62 {
63 if (GLOBAL_pmosloc2df == NULL) GLOBAL_pmosloc2df = new TopOpeBRepTool_DataMapOfShapeListOfC2DF();
64 Standard_Integer ancemp = (*GLOBAL_pidmoslosc2df).Extent();
65 if (ancemp == 0) {
66 TopExp::MapShapesAndAncestors(*GLOBAL_pS1c2df,TopAbs_EDGE,TopAbs_FACE,(*GLOBAL_pidmoslosc2df));
67 TopExp::MapShapesAndAncestors(*GLOBAL_pS2c2df,TopAbs_EDGE,TopAbs_FACE,(*GLOBAL_pidmoslosc2df));
68 }
69 Standard_Boolean Eb = (*GLOBAL_pidmoslosc2df).Contains(E);
70 if ( !Eb ) return *GLOBAL_pFc2df;
71 const TopTools_ListOfShape& lf = (*GLOBAL_pidmoslosc2df).FindFromKey(E);
72 if (lf.IsEmpty()) return *GLOBAL_pFc2df;
73 const TopoDS_Face& F = TopoDS::Face(lf.First());
74 return F;
75 }
76
77 // ------------------------------------------------------------------------------------
FC2D_Prepare(const TopoDS_Shape & S1,const TopoDS_Shape & S2)78 Standard_EXPORT int FC2D_Prepare(const TopoDS_Shape& S1,const TopoDS_Shape& S2)
79 {
80 if (GLOBAL_pmosloc2df == NULL) GLOBAL_pmosloc2df = new TopOpeBRepTool_DataMapOfShapeListOfC2DF();
81 GLOBAL_pmosloc2df->Clear();
82 GLOBAL_C2D_i = 0;
83
84 if (GLOBAL_pidmoslosc2df == NULL) GLOBAL_pidmoslosc2df = new TopTools_IndexedDataMapOfShapeListOfShape();
85 GLOBAL_pidmoslosc2df->Clear();
86
87 if (GLOBAL_pFc2df == NULL) GLOBAL_pFc2df = new TopoDS_Face();
88 GLOBAL_pFc2df->Nullify();
89
90 if (GLOBAL_pS1c2df == NULL) GLOBAL_pS1c2df = new TopoDS_Shape();
91 *GLOBAL_pS1c2df = S1;
92
93 if (GLOBAL_pS2c2df == NULL) GLOBAL_pS2c2df = new TopoDS_Shape();
94 *GLOBAL_pS2c2df = S2;
95
96 return 0;
97 }
98
99 // ------------------------------------------------------------------------------------
FC2D_HasC3D(const TopoDS_Edge & E)100 Standard_EXPORT Standard_Boolean FC2D_HasC3D(const TopoDS_Edge& E)
101 {
102 TopLoc_Location loc; Standard_Real f3d,l3d; Handle(Geom_Curve) C3D = BRep_Tool::Curve(E,loc,f3d,l3d);
103 Standard_Boolean b = (!C3D.IsNull());
104 return b;
105 }
106
107 // ------------------------------------------------------------------------------------
FC2D_HasCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F)108 Standard_EXPORT Standard_Boolean FC2D_HasCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F)
109 {
110 Handle(Geom2d_Curve) C2D;
111 Standard_Boolean hasold = FC2D_HasOldCurveOnSurface(E,F,C2D);
112 Standard_Boolean hasnew = FC2D_HasNewCurveOnSurface(E,F,C2D);
113 Standard_Boolean b = hasold || hasnew;
114 return b;
115 }
116
117 // ------------------------------------------------------------------------------------
FC2D_HasOldCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F,Handle (Geom2d_Curve)& C2D,Standard_Real & f2d,Standard_Real & l2d,Standard_Real & tol)118 Standard_EXPORT Standard_Boolean FC2D_HasOldCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Handle(Geom2d_Curve)& C2D,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol)
119 {
120 Standard_Boolean hasold = Standard_False;
121 tol = BRep_Tool::Tolerance(E);
122 C2D = BRep_Tool::CurveOnSurface(E,F,f2d,l2d);
123 hasold = (!C2D.IsNull());
124 return hasold;
125 }
FC2D_HasOldCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F,Handle (Geom2d_Curve)& C2D)126 Standard_EXPORT Standard_Boolean FC2D_HasOldCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Handle(Geom2d_Curve)& C2D)
127 { Standard_Real f2d,l2d,tol; Standard_Boolean hasold = FC2D_HasOldCurveOnSurface(E,F,C2D,f2d,l2d,tol); return hasold;
128 }
129
130 // ------------------------------------------------------------------------------------
FC2D_PNewCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F)131 static TopOpeBRepTool_C2DF* FC2D_PNewCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F)
132 {
133 TopOpeBRepTool_C2DF* pc2df = NULL;
134 if (GLOBAL_pmosloc2df == NULL) return NULL;
135 Standard_Boolean Eisb = GLOBAL_pmosloc2df->IsBound(E);
136 if (!Eisb) return NULL;
137 TopOpeBRepTool_ListIteratorOfListOfC2DF it(GLOBAL_pmosloc2df->Find(E));
138 for(;it.More();it.Next()) {
139 const TopOpeBRepTool_C2DF& c2df = it.Value();
140 Standard_Boolean isf = c2df.IsFace(F);
141 if (isf) {
142 pc2df = (TopOpeBRepTool_C2DF*)&c2df;
143 break;
144 }
145 }
146 return pc2df;
147 }
FC2D_HasNewCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F,Handle (Geom2d_Curve)& C2D,Standard_Real & f2d,Standard_Real & l2d,Standard_Real & tol)148 Standard_EXPORT Standard_Boolean FC2D_HasNewCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Handle(Geom2d_Curve)& C2D,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol)
149 {
150 const TopOpeBRepTool_C2DF* pc2df = FC2D_PNewCurveOnSurface(E,F);
151 Standard_Boolean hasnew = (pc2df != NULL);
152 if (hasnew) C2D = pc2df->PC(f2d,l2d,tol);
153 return hasnew;
154 }
FC2D_HasNewCurveOnSurface(const TopoDS_Edge & E,const TopoDS_Face & F,Handle (Geom2d_Curve)& C2D)155 Standard_EXPORT Standard_Boolean FC2D_HasNewCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Handle(Geom2d_Curve)& C2D)
156 { Standard_Real f2d,l2d,tol; Standard_Boolean b = FC2D_HasNewCurveOnSurface(E,F,C2D,f2d,l2d,tol); return b;
157 }
158
159 // ------------------------------------------------------------------------------------
FC2D_AddNewCurveOnSurface(Handle (Geom2d_Curve)C2D,const TopoDS_Edge & E,const TopoDS_Face & F,const Standard_Real & f2d,const Standard_Real & l2d,const Standard_Real & tol)160 Standard_Integer FC2D_AddNewCurveOnSurface(Handle(Geom2d_Curve) C2D,const TopoDS_Edge& E,const TopoDS_Face& F,const Standard_Real& f2d,const Standard_Real& l2d,const Standard_Real& tol)
161 {
162 if (C2D.IsNull()) return 1;
163 TopOpeBRepTool_C2DF c2df(C2D,f2d,l2d,tol,F);
164 if (GLOBAL_pmosloc2df == NULL) return 1;
165 TopOpeBRepTool_ListOfC2DF thelist;
166 GLOBAL_pmosloc2df->Bind(E, thelist);
167 TopOpeBRepTool_ListOfC2DF& lc2df = GLOBAL_pmosloc2df->ChangeFind(E);
168 lc2df.Append(c2df);
169 return 0;
170 }
171
172 // ------------------------------------------------------------------------------------
173 static Handle(Geom2d_Curve) FC2D_make2d(const TopoDS_Edge& E,const TopoDS_Face& F,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol,const Standard_Boolean trim3d = Standard_False);
FC2D_make2d(const TopoDS_Edge & E,const TopoDS_Face & F,Standard_Real & f2d,Standard_Real & l2d,Standard_Real & tol,const Standard_Boolean trim3d)174 static Handle(Geom2d_Curve) FC2D_make2d(const TopoDS_Edge& E,const TopoDS_Face& F,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol,const Standard_Boolean trim3d)
175 {
176 Handle(Geom2d_Curve) C2D = BRep_Tool::CurveOnSurface(E,F,f2d,l2d);
177 if (!C2D.IsNull()) return C2D;
178
179 // pas de 2D
180 Standard_Real f3d,l3d;TopLoc_Location eloc;
181 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E,eloc,f3d,l3d);
182 Standard_Boolean hasC3D = (!C1.IsNull());
183
184 if ( hasC3D ) {
185 Standard_Boolean elocid = eloc.IsIdentity();
186 Handle(Geom_Curve) C2;
187 if (elocid ) C2 = C1;
188 else C2 = Handle(Geom_Curve)::DownCast(C1->Transformed(eloc.Transformation()));
189 Standard_Real f=0.,l=0.;
190 if (trim3d) {f=f3d; l=l3d;}
191 C2D = TopOpeBRepTool_CurveTool::MakePCurveOnFace(F,C2,tol,f,l);
192 f2d = f3d; l2d = l3d;
193 return C2D;
194 }
195 else {
196 // E sans courbe 2d sur F, E sans courbe 3d
197 // une face accedant a E : FE
198 const TopoDS_Face& FE = FC2D_FancestorE(E);
199 if (FE.IsNull()) return C2D;
200 Standard_Boolean compminmaxUV = Standard_False;
201 BRepAdaptor_Surface BAS(F,compminmaxUV);
202 Handle(BRepAdaptor_Surface) BAHS = new BRepAdaptor_Surface(BAS);
203 BRepAdaptor_Curve AC(E,FE);
204 Handle(BRepAdaptor_Curve) AHC = new BRepAdaptor_Curve(AC);
205 Standard_Real tolin; FTOL_FaceTolerances3d(F,FE,tolin);
206 ProjLib_ProjectedCurve projcurv(BAHS,AHC,tolin);
207 C2D = MakePCurve(projcurv);
208 Standard_Real f,l; BRep_Tool::Range(E,f,l);
209 f2d = f; l2d = l;
210 }
211
212 #ifdef OCCT_DEBUG
213 if (TopOpeBRepTool_GettraceC2D() && C2D.IsNull()) {
214 std::cout<<"#FC2D_make2d1 --> PCurve IsNull"<<std::endl;
215 }
216 #endif
217
218 return C2D;
219 } // make2d1
220
221 // ------------------------------------------------------------------------------------
222 //modified by NIZHNY-MZV Mon Oct 4 10:37:36 1999
Handle(Geom2d_Curve)223 Standard_EXPORT Handle(Geom2d_Curve) FC2D_MakeCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Standard_Real& f,Standard_Real& l,Standard_Real& tol,const Standard_Boolean trim3d)
224 {
225 #ifdef DRAW
226 if (TopOpeBRepTool_GettraceC2D()) {
227 std::cout<<"\n#FC2D_MakeCurveOnSurface : "<<std::endl;
228 TCollection_AsciiString fnam("c2df");GLOBAL_C2D_i++;fnam=fnam+GLOBAL_C2D_i;FDRAW_DINS("",F,fnam,"");
229 TCollection_AsciiString enam("c2de");GLOBAL_C2D_i++;enam=enam+GLOBAL_C2D_i;FDRAW_DINE(" ",E,enam,"\n");
230 std::cout.flush(); debc2dnull();
231 }
232 #endif
233
234 Handle(Geom2d_Curve) C2D = FC2D_make2d(E,F,f,l,tol,trim3d);
235 FC2D_AddNewCurveOnSurface(C2D,E,F,f,l,tol);
236 return C2D;
237 }
238
239 // ------------------------------------------------------------------------------------
Handle(Geom2d_Curve)240 Standard_EXPORT Handle(Geom2d_Curve) FC2D_CurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Standard_Real& f,Standard_Real& l,Standard_Real& tol,const Standard_Boolean trim3d)
241 {
242 Handle(Geom2d_Curve) C2D;
243 Standard_Boolean hasold = FC2D_HasOldCurveOnSurface(E,F,C2D,f,l,tol);
244 if (hasold) {
245 return C2D;
246 }
247 Standard_Boolean hasnew = FC2D_HasNewCurveOnSurface(E,F,C2D,f,l,tol);
248 if (hasnew) {
249 return C2D;
250 }
251 C2D = FC2D_MakeCurveOnSurface(E,F,f,l,tol,trim3d);
252 return C2D;
253 }
254
255 // ------------------------------------------------------------------------------------
Handle(Geom2d_Curve)256 Standard_EXPORT Handle(Geom2d_Curve) FC2D_EditableCurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,Standard_Real& f,Standard_Real& l,Standard_Real& tol,const Standard_Boolean trim3d)
257 {
258 Standard_Boolean hasold = Standard_False;
259 {
260 Handle(Geom2d_Curve) C2D;
261 hasold = FC2D_HasOldCurveOnSurface(E,F,C2D,f,l,tol);
262 if (hasold) {
263 Handle(Geom2d_Curve) copC2D = Handle(Geom2d_Curve)::DownCast(C2D->Copy());
264 return copC2D;
265 }
266 }
267 Standard_Boolean hasnew = Standard_False;
268 {
269 Handle(Geom2d_Curve) newC2D;
270 hasnew = FC2D_HasNewCurveOnSurface(E,F,newC2D,f,l,tol);
271 if (hasnew) {
272 return newC2D;
273 }
274 }
275 Handle(Geom2d_Curve) makC2D = FC2D_MakeCurveOnSurface(E,F,f,l,tol,trim3d);
276 return makC2D;
277 }
278
279 // ------------------------------------------------------------------------------------
FC2D_translate(Handle (Geom2d_Curve)C2D,const TopoDS_Edge &,const TopoDS_Face & F,const TopoDS_Edge & EF)280 static void FC2D_translate(Handle(Geom2d_Curve) C2D,
281 // const TopoDS_Edge& E,
282 const TopoDS_Edge& ,
283 const TopoDS_Face& F,
284 const TopoDS_Edge& EF)
285 {
286 TopLoc_Location sloc; const Handle(Geom_Surface)& S1 = BRep_Tool::Surface(F,sloc);
287 Standard_Boolean isperio = S1->IsUPeriodic() || S1->IsVPeriodic();
288 gp_Dir2d d2d; gp_Pnt2d O2d; Standard_Boolean isuiso,isviso;
289 Standard_Boolean uviso = TopOpeBRepTool_TOOL::UVISO(C2D,isuiso,isviso,d2d,O2d);
290 Standard_Boolean EFnull = EF.IsNull();
291
292 if (isperio && uviso && !EFnull) {
293 // C2D prend comme origine dans F l'origine de la pcurve de EF dans F
294 TopoDS_Face FFOR = F;
295 FFOR.Orientation(TopAbs_FORWARD);
296 gp_Pnt2d p1,p2; BRep_Tool::UVPoints(EF,FFOR,p1,p2);
297 Standard_Real pEF = isuiso ? p1.X() : p1.Y();
298 Standard_Real pC2D = isuiso ? O2d.X() : O2d.Y();
299 Standard_Real factor = pEF - pC2D;
300 Standard_Boolean b = (Abs(factor) > 1.e-6);
301 if ( b ) {
302 gp_Vec2d transl(1.,0.); if (isviso) transl = gp_Vec2d(0.,1.);
303 transl.Multiply(factor);
304 C2D->Translate(transl);
305 }
306 }
307 }
308
309 // ------------------------------------------------------------------------------------
310 static Handle(Geom2d_Curve) FC2D_make2d(const TopoDS_Edge& E,const TopoDS_Face& F,const TopoDS_Edge& EF,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol,const Standard_Boolean trim3d = Standard_False);
FC2D_make2d(const TopoDS_Edge & E,const TopoDS_Face & F,const TopoDS_Edge & EF,Standard_Real & f2d,Standard_Real & l2d,Standard_Real & tol,const Standard_Boolean trim3d)311 static Handle(Geom2d_Curve) FC2D_make2d(const TopoDS_Edge& E,const TopoDS_Face& F,const TopoDS_Edge& EF,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol,const Standard_Boolean trim3d)
312 {
313 Handle(Geom2d_Curve) C2D = BRep_Tool::CurveOnSurface(E,F,f2d,l2d);
314 if (!C2D.IsNull()) return C2D;
315
316 // pas de 2D
317 Standard_Real f3d,l3d;TopLoc_Location eloc;
318 Handle(Geom_Curve) C1 = BRep_Tool::Curve(E,eloc,f3d,l3d);
319 Standard_Boolean hasC3D = (!C1.IsNull());
320
321 if ( hasC3D ) {
322 Standard_Boolean elocid = eloc.IsIdentity();
323 Handle(Geom_Curve) C2;
324 if (elocid ) C2 = C1;
325 else C2 = Handle(Geom_Curve)::DownCast(C1->Transformed(eloc.Transformation()));
326 Standard_Real f=0.,l=0.;
327 if (trim3d) {f=f3d; l=l3d;}
328 C2D = TopOpeBRepTool_CurveTool::MakePCurveOnFace(F,C2,tol,f,l);
329 f2d = f3d; l2d = l3d;
330 FC2D_translate(C2D,E,F,EF);
331 return C2D;
332 }
333 else {
334 // E sans courbe 2d sur F, E sans courbe 3d
335 // une face accedant a E : FE
336 const TopoDS_Face& FE = FC2D_FancestorE(E);
337 if (FE.IsNull()) return C2D;
338 Standard_Boolean compminmaxUV = Standard_False;
339 BRepAdaptor_Surface BAS(F,compminmaxUV);
340 Handle(BRepAdaptor_Surface) BAHS = new BRepAdaptor_Surface(BAS);
341 BRepAdaptor_Curve AC(E,FE);
342 Handle(BRepAdaptor_Curve) AHC = new BRepAdaptor_Curve(AC);
343 Standard_Real tolin; FTOL_FaceTolerances3d(F,FE,tolin);
344 ProjLib_ProjectedCurve projcurv(BAHS,AHC,tolin);
345 C2D = MakePCurve(projcurv);
346 Standard_Real f,l; BRep_Tool::Range(E,f,l);
347 f2d = f; l2d = l;
348 FC2D_translate(C2D,E,F,EF);
349 }
350
351 #ifdef OCCT_DEBUG
352 if (TopOpeBRepTool_GettraceC2D() && C2D.IsNull()) {
353 std::cout<<"#FC2D_make2d2 --> PCurve IsNull"<<std::endl;
354 }
355 #endif
356
357 return C2D;
358 } // make2d2
359
360 // ------------------------------------------------------------------------------------
Handle(Geom2d_Curve)361 Standard_EXPORT Handle(Geom2d_Curve) FC2D_CurveOnSurface(const TopoDS_Edge& E,const TopoDS_Face& F,const TopoDS_Edge& EF,Standard_Real& f2d,Standard_Real& l2d,Standard_Real& tol,const Standard_Boolean trim3d)
362 {
363 Handle(Geom2d_Curve) C2D;
364
365 Standard_Boolean hasold = FC2D_HasOldCurveOnSurface(E,F,C2D,f2d,l2d,tol);
366 if (hasold) return C2D;
367
368 TopOpeBRepTool_C2DF* pc2df = FC2D_PNewCurveOnSurface(E,F);
369 if (pc2df != NULL) {
370 C2D = pc2df->PC(f2d,l2d,tol);
371 FC2D_translate(C2D,E,F,EF);
372 pc2df->SetPC(C2D,f2d,l2d,tol);
373 return C2D;
374 }
375
376 #ifdef DRAW
377 if (TopOpeBRepTool_GettraceC2D()) {
378 std::cout<<"\n#FC2D_CurveOnSurface : "<<std::endl;
379 TCollection_AsciiString fnam("c2df");GLOBAL_C2D_i++;fnam=fnam+GLOBAL_C2D_i;FDRAW_DINS("",F,fnam,"");
380 TCollection_AsciiString enam("c2de");GLOBAL_C2D_i++;enam=enam+GLOBAL_C2D_i;FDRAW_DINE(" ",E,enam,"\n");
381 std::cout.flush();debc2dnull();
382 }
383 #endif
384
385 C2D = FC2D_make2d(E,F,EF,f2d,l2d,tol,trim3d);
386 FC2D_AddNewCurveOnSurface(C2D,E,F,f2d,l2d,tol);
387 return C2D;
388 }
389