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