2726     if (OutSide|| !Gh.subdomains || !Gh.NbSubDomains )  in FindSubDomain()
2824 	if (Gh.NbSubDomains != NbSubDomains && subdomains)  in FindSubDomain()
2827 	    subdomains = new SubDomain[ Gh.NbSubDomains];  in FindSubDomain()
2828 	NbSubDomains =Gh.NbSubDomains;  in FindSubDomain()
2841 	    GeometricalEdge &eg = *Gh.subdomains[i].edge;  in FindSubDomain()
2842 	    subdomains[i].ref = Gh.subdomains[i].ref;  in FindSubDomain()
2846 	    Edge &e = *GeometricalEdgetoEdge[Gh.Number(eg)];  in FindSubDomain()
2850 	    int sens = Gh.subdomains[i].sens;  in FindSubDomain()
2879 			     << ": Bad sens  " << Gh.Number(eg) <<" "<< sens <<  endl;  in FindSubDomain()
2911 			<< " edge=" << Gh.Number(eg) << " " << sens << endl;  in FindSubDomain()
3225 …identity = new char[ lg + strlen(buf) + strlen(countbuf)+ 2  + 10 + ( Gh.name ? strlen(Gh.name) + …  in PreInit()
3230       if (Gh.name)  in PreInit()
3231 	strcat(strcat(identity,"G="),Gh.name);  in PreInit()
3260   Gh.NbRef++;// add a ref to Gh  in GeomToTriangles1()
3281   assert(&BTh.Gh == &Gh);  in GeomToTriangles1()
3293   int * bcurve = new int[Gh.NbOfCurves]; //  in GeomToTriangles1()
3311   for ( i=0;i<Gh.nbv;i++)  in GeomToTriangles1()
3312     if (Gh[i].Required()) NbVerticesOnGeomVertex++;  in GeomToTriangles1()
3323   for (i=0;i<Gh.nbv;i++)  in GeomToTriangles1()
3324     if (Gh[i].Required()) {//Gh  vertices Required  in GeomToTriangles1()
3325         vertices[nbv] =  Gh[i];  in GeomToTriangles1()
3327         Gh[i].to = vertices + nbv;// save Geom -> Th  in GeomToTriangles1()
3328         VerticesOnGeomVertex[nbv]= VertexOnGeom(vertices[nbv],Gh[i]);  in GeomToTriangles1()
3331     else Gh[i].to=0;  in GeomToTriangles1()
3348   Gh.UnMarkEdges();  in GeomToTriangles1()
3357   for (int i=0;i<Gh.NbOfCurves;i++)  in GeomToTriangles1()
3374 	         ei.on==Gh.curves[nc].be    &&  in GeomToTriangles1()
3375 …      (GeometricalVertex *) *ei[je].on == &(*Gh.curves[nc].be)[Gh.curves[nc].kb] //  same extremity  in GeomToTriangles1()
3384    assert( bfind==Gh.NbOfCurves);  in GeomToTriangles1()
3396       Gh.UnMarkEdges();  in GeomToTriangles1()
3407        for (int icurve=0;icurve<Gh.NbOfCurves;icurve++)  in GeomToTriangles1()
3411               if( ! Gh.curves[icurve].master) continue; // we skip all equi curve  in GeomToTriangles1()
3425 		  for(Curve * curve= Gh.curves+icurve;curve;curve= curve->next)  in GeomToTriangles1()
3428 		    int icurveequi= Gh.Number(curve);  in GeomToTriangles1()
3502 			  ongequi = Gh.ProjectOnCurve(eeequi,se,*A1,*GA1);  in GeomToTriangles1()
3557 			      << " Gh.vertices " << Gh.vertices  in GeomToTriangles1()
3559 			cerr << ee[k1].on->gv-Gh.vertices << endl;  in GeomToTriangles1()
3604 		    for(Curve * curve= Gh.curves+icurve;curve;curve= curve->next)  in GeomToTriangles1()
3675   Gh.NbRef++;// add a ref to GH  in GeomToTriangles0()
3683   Gh.InitDraw() ;  in GeomToTriangles0()
3697   for (i=0;i<Gh.nbv;i++)  in GeomToTriangles0()
3698     if (Gh[i].Required() && Gh[i].IsThe() ) NbVerticesOnGeomVertex++;  in GeomToTriangles0()
3706   for (i=0;i<Gh.nbv;i++)  in GeomToTriangles0()
3707     if (Gh[i].Required()&& Gh[i].IsThe()  ) {//Gh  vertices Required  in GeomToTriangles0()
3709 	vertices[nbv] = Gh[i];  in GeomToTriangles0()
3710       Gh[i].to = vertices + nbv;// save Geom -> Th  in GeomToTriangles0()
3711       VerticesOnGeomVertex[nbv]= VertexOnGeom(*Gh[i].to,Gh[i]);  in GeomToTriangles0()
3734       Gh.UnMarkEdges();  in GeomToTriangles0()
3736       for (i=0;i<Gh.nbe;i++) {  in GeomToTriangles0()
3737 	GeometricalEdge & ei = Gh.edges[i];  in GeomToTriangles0()
3979   assert(Gh.nbe);  in MakeGeometricalEdgeToEdge()
3980   Edge **e= new Edge* [this->Gh.nbe];  in MakeGeometricalEdgeToEdge()
3983   for ( i=0;i<Gh.nbe ; i++)  in MakeGeometricalEdgeToEdge()
3989       e[Gh.Number(on)] = ei;  in MakeGeometricalEdgeToEdge()
4005        if (e[Gh.Number(on)])  break; // optimisation  in MakeGeometricalEdgeToEdge()
4006        e[Gh.Number(on)] = ei;  in MakeGeometricalEdgeToEdge()
4011      for ( i=0;i<Gh.nbe ; i++)  in MakeGeometricalEdgeToEdge()
4014 …cerr << " Bug -- the geometrical edge " << i << " is on no edge curve = " << Gh.edges[i].CurveNumb…  in MakeGeometricalEdgeToEdge()
4015 		<< " s0 " << Gh.Number( Gh.edges[i][0]) << " s1  " << Gh.Number( Gh.edges[i][1]) << endl;  in MakeGeometricalEdgeToEdge()
4042   if (&Gh)  in ~Triangles()
4044       if (Gh.NbRef>0) Gh.NbRef--;  in ~Triangles()
4045       else if (Gh.NbRef==0) delete &Gh;  in ~Triangles()
4448 		   <<  Number(edges[i][j]) << " : "  << " on = " << Gh.Number(edges[i].on) ;  in FillHoleInMesh()
4450 		cerr << " vertex " << Gh.Number(edges[i][j].on->gv);  in FillHoleInMesh()
4452 		cerr << "Edges " << Gh.Number(edges[i][j].on->ge);  in FillHoleInMesh()
4467 : Gh(*(pGh?pGh:&Th.Gh)), BTh(*(pBth?pBth:this))  in Triangles()
4469   Gh.NbRef++;  in Triangles()
4977   : Gh(*(new Geometry())), BTh(*this)  in Triangles()
5064   Gh.AfterRead();  in Triangles()