197 Fem2D::Mesh *bamg2msh(const bamg::Geometry &Gh)  in bamg2msh()  argument
200   int nv=  Gh.nbv;  in bamg2msh()
201   int neb=Gh.nbe;  in bamg2msh()
207       const bamg::GeometricalVertex & vg( Gh[i]);  in bamg2msh()
215       const bamg::GeometricalEdge & e= Gh(ie);  in bamg2msh()
216       int i0=Gh.Number(e[0]),i1=Gh.Number(e[0]);  in bamg2msh()
290   Tn->Gh.AfterRead();  in msh2bamg()
387   Tn->Gh.AfterRead();  in msh2bamg()
433   Geometry * Gh =  new Geometry;  in BuildMesh()  local
438   Gh->name=new char [strlen(filename)+1];  in BuildMesh()
439   strcpy(Gh->name,filename);  in BuildMesh()
443   Gh->MaximalAngleOfCorner =30.00*Pi/180.0;  in BuildMesh()
444   Gh->nbv = 0;  in BuildMesh()
445   Gh->nbvx = nbvx;  in BuildMesh()
447   Gh->nbe = nbe;  in BuildMesh()
448   Gh->edges = new GeometricalEdge[Gh->nbe];  in BuildMesh()
449   bamg::Vertex *vertices =  new Vertex[Gh->nbvx];// correction 2009/07/03  in BuildMesh()
501   Gh->pmin =  vertices[0].r;  in BuildMesh()
502   Gh->pmax =  vertices[0].r;  in BuildMesh()
506       Gh->pmin.x = Min(Gh->pmin.x,vertices[i].r.x);  in BuildMesh()
507       Gh->pmin.y = Min(Gh->pmin.y,vertices[i].r.y);  in BuildMesh()
508       Gh->pmax.x = Max(Gh->pmax.x,vertices[i].r.x);  in BuildMesh()
509       Gh->pmax.y = Max(Gh->pmax.y,vertices[i].r.y);  in BuildMesh()
512   double diameter=Max(Gh->pmax.x-Gh->pmin.x,Gh->pmax.y-Gh->pmin.y);  in BuildMesh()
513   Gh->coefIcoor= (MaxICoor)/diameter;  in BuildMesh()
514   Icoor1 epsI = (Icoor1) (Gh->coefIcoor*eps);  in BuildMesh()
515   ffassert(Gh->coefIcoor >0);  in BuildMesh()
522       cout <<"\t\t"  << "     Geom: min="<< Gh->pmin << "max ="<< Gh->pmax  in BuildMesh()
523            << " hmin = " << Gh->MinimalHmin() <<  endl;  in BuildMesh()
531         vertices[i].i = Gh->toI2(vertices[i].r);  in BuildMesh()
544   Gh->nbvx = nbv;  in BuildMesh()
545   Gh->nbv = nbv;  in BuildMesh()
547   Gh->vertices = new GeometricalVertex[nbv];  in BuildMesh()
548   throwassert(Gh->nbvx >= Gh->nbv);  in BuildMesh()
549   Gh->nbiv = Gh->nbv;  in BuildMesh()
557           Gh->vertices[i].r.x = vertices[kkk].r.x ;  in BuildMesh()
558           Gh->vertices[i].r.y = vertices[kkk].r.y;  in BuildMesh()
559           throwassert(Gh->vertices[i].IsThe());  in BuildMesh()
560           Gh->vertices[i].ReferenceNumber = vertices[kkk].ReferenceNumber  ;  in BuildMesh()
561           Gh->vertices[i].DirOfSearch = NoDirOfSearch;  in BuildMesh()
562           Gh->vertices[i].color =0;  in BuildMesh()
563           Gh->vertices[i].Set();  in BuildMesh()
565            Gh->vertices[i].SetRequired();  in BuildMesh()
567           if(Gh->vertices[i].ReferenceNumber < 0)  in BuildMesh()
568             Gh->vertices[i].SetRequired();  in BuildMesh()
576     cout <<"\t\t"  << "     Record Edges: Nb of Edge " << Gh->nbe <<endl;  in BuildMesh()
577   throwassert(Gh->edges);  in BuildMesh()
578   throwassert (Gh->nbv >0);  in BuildMesh()
582       len = new Real4[Gh->nbv];  in BuildMesh()
583       for(i=0;i<Gh->nbv;i++)  in BuildMesh()
610           Gh->edges[i].ref = mp.label;  in BuildMesh()
611           Gh->edges[i].v[0]=  Gh->vertices + i1;  in BuildMesh()
612           Gh->edges[i].v[1]=  Gh->vertices + i2;  in BuildMesh()
613           R2 x12 = Gh->vertices[i2].r-Gh->vertices[i1].r;  in BuildMesh()
615           Gh->edges[i].tg[0]=zero2;  in BuildMesh()
616           Gh->edges[i].tg[1]=zero2;  in BuildMesh()
617           Gh->edges[i].SensAdj[0] = Gh->edges[i].SensAdj[1] = -1;  in BuildMesh()
618           Gh->edges[i].Adj[0] = Gh->edges[i].Adj[1] = 0;  in BuildMesh()
619           Gh->edges[i].flag = 0;  in BuildMesh()
620           Gh->edges[i].link=0;  in BuildMesh()
622 	  Gh->edges[i].SetRequired();  in BuildMesh()
626               Gh->vertices[i1].color++;  in BuildMesh()
627               Gh->vertices[i2].color++;  in BuildMesh()
640   throwassert(i==Gh->nbe);  in BuildMesh()
645       for (i=0;i<Gh->nbv;i++)  in BuildMesh()
647         if(hvint &&Gh->vertices[i].ReferenceNumber <brefintp +nbvinter)  in BuildMesh()
649             long k =Gh->vertices[i].ReferenceNumber-brefintp;  in BuildMesh()
650             Gh->vertices[i].m=Metric( (*pintern)(k ,2L));  in BuildMesh()
651             Gh->vertices[i].ReferenceNumber = -1; //++ bof bof FH ..  in BuildMesh()
654         else if (Gh->vertices[i].color > 0)  in BuildMesh()
655           Gh->vertices[i].m=  Metric(len[i] /(Real4) Gh->vertices[i].color);  in BuildMesh()
657           Gh->vertices[i].m=  Metric(Hmin);  in BuildMesh()
665   Gh->NbSubDomains=nbsd;  in BuildMesh()
666   if (Gh->NbSubDomains>0)  in BuildMesh()
668         Gh->subdomains = new GeometricalSubDomain[  Gh->NbSubDomains];  in BuildMesh()
678                 Gh->subdomains[i].sens = Nbseg >0 ? 1 : -1;  in BuildMesh()
679                 Gh->subdomains[i].edge=Gh->edges + i1;  in BuildMesh()
680                 Gh->subdomains[i].ref = i;  in BuildMesh()
684   Gh->NbEquiEdges=0;  in BuildMesh()
685   Gh->NbCrackedEdges=0;  in BuildMesh()
688     m=bamg2msh(*Gh);  in BuildMesh()
691       Gh->AfterRead();  in BuildMesh()
692       int nbtx= nbvmax ? nbvmax :  (Gh->nbv*Gh->nbv)/9 +1000;  in BuildMesh()
693       if(verbosity> 99) cout << " ** Gh = " << endl << *Gh << endl << " *** " <<endl; ;  in BuildMesh()
696 	  Th =new Triangles( nbtx ,*Gh);  in BuildMesh()
737 	  Gh->NbRef=0;  in BuildMesh()
738 	  delete Gh;  in BuildMesh()
747   delete Gh;  in BuildMesh()
968   Geometry Gh(s->c_str());  in buildmeshbamg()  local
969   int nbvx = nbvxin ? nbvxin : ((Gh.nbv*Gh.nbv)/9 +1000);  in buildmeshbamg()
970   Triangles * bTh=  new Triangles(nbvx,Gh);  in buildmeshbamg()