geometrie/src/ot_geometrie.cpp

#include "ot_geometrie.h"
#include "ot_mathematique.h"
#ifdef ALL_OCC
#include <BRepExtrema_DistShapeShape.hxx>
#include <BRepExtrema_DistanceSS.hxx>
#include <GeomAPI_ProjectPointOnCurve.hxx>
#include <GeomAPI_ProjectPointOnSurf.hxx>
#include <ShapeAnalysis.hxx>
#include <BRepBuilderAPI_MakeVertex.hxx>
#include <BRepAlgoAPI_Common.hxx>
#include <BRepGProp.hxx>
#include <GProp_GProps.hxx>
#include <BRepExtrema_ShapeProximity.hxx>
#include <limits>

#endif

OT_GEOMETRIE::OT_GEOMETRIE()
{

}

OT_GEOMETRIE::~OT_GEOMETRIE()
{

}

int OT_GEOMETRIE::get_lien_topologique(MG_ELEMENT_TOPOLOGIQUE* mgeletopo1, MG_ELEMENT_TOPOLOGIQUE* mgeletopo2)
{
  TPL_MAP_ENTITE<MG_SOMMET*> tpl_map_sommet1;
  TPL_MAP_ENTITE<MG_SOMMET*> tpl_map_sommet2;
  TPL_MAP_ENTITE<MG_ARETE*> tpl_map_arete1;
  TPL_MAP_ENTITE<MG_ARETE*> tpl_map_arete2;
  TPL_MAP_ENTITE<MG_FACE*> tpl_map_face1;
  TPL_MAP_ENTITE<MG_FACE*> tpl_map_face2;
  get_map_mg_ele_topo_sous_jacent(mgeletopo1,tpl_map_sommet1,tpl_map_arete1,tpl_map_face1);
  get_map_mg_ele_topo_sous_jacent(mgeletopo2,tpl_map_sommet2,tpl_map_arete2,tpl_map_face2);
  if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE) tpl_map_face1.ajouter(dynamic_cast<MG_FACE*>(mgeletopo1));
  else if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE) tpl_map_arete1.ajouter(dynamic_cast<MG_ARETE*>(mgeletopo1));
  else if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET) tpl_map_sommet1.ajouter(dynamic_cast<MG_SOMMET*>(mgeletopo1));
  if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE) tpl_map_face2.ajouter(dynamic_cast<MG_FACE*>(mgeletopo2));
  else if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE) tpl_map_arete2.ajouter(dynamic_cast<MG_ARETE*>(mgeletopo2));
  else if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET) tpl_map_sommet2.ajouter(dynamic_cast<MG_SOMMET*>(mgeletopo2));
  bool trouve=false;
  TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face1;
  MG_FACE *face1=tpl_map_face1.get_premier(it_face1);
  while(trouve==false && face1!=NULL)
  {
    if(tpl_map_face2.existe(face1))
    {
      trouve=true;
      return TYPE_LIEN_TOPOLOGIQUE::FACE;
    }
    else face1=tpl_map_face1.get_suivant(it_face1);
  }
  TPL_MAP_ENTITE<MG_ARETE*>::ITERATEUR it_arete1;
  MG_ARETE *arete1=tpl_map_arete1.get_premier(it_arete1);
  while(trouve==false && arete1!=NULL)
  {
    if(tpl_map_arete2.existe(arete1))
    {
      trouve=true;
      return TYPE_LIEN_TOPOLOGIQUE::ARETE;
    }
    else arete1=tpl_map_arete1.get_suivant(it_arete1);
  }
  TPL_MAP_ENTITE<MG_SOMMET*>::ITERATEUR it_sommet1;
  MG_SOMMET *sommet1=tpl_map_sommet1.get_premier(it_sommet1);
  while(trouve==false && sommet1!=NULL)
  {
    if(tpl_map_sommet2.existe(sommet1))
    {
      trouve=true;
      return TYPE_LIEN_TOPOLOGIQUE::SOMMET;
    }
    else sommet1=tpl_map_sommet1.get_suivant(it_sommet1);
  }
  return TYPE_LIEN_TOPOLOGIQUE::AUCUN;
}


int OT_GEOMETRIE::declage_mg_face_mg_face(MG_FACE* mgface1, MG_FACE* mgface2, double& moyenne_distance, double& ecart_type_distance, int nb_pas, double eps)
{
  moyenne_distance=0;
  ecart_type_distance=0;
  std::vector<double> vector_distance;
  double umin1,umax1,vmin1,vmax1;
  get_param_face(mgface1,umin1,umax1,vmin1,vmax1);
  double pas_u = (umax1-umin1)/nb_pas;
  double pas_v = (vmax1-vmin1)/nb_pas;
  double uv[2];
  for(int i=0;i<nb_pas;i++)
  {
    uv[0]=umin1+i*pas_u;
    for(int j=0;j<nb_pas;j++)
    {
      uv[1]=vmin1+j*pas_v;
      double xyz1[3];
      mgface1->evaluer(uv,xyz1);
      double xyz2[3];
//       int ret=projection_au_plus_pres_sur_mg_face(xyz1,mgface1,mgface2,xyz2);
      int ret = projection_orthogonale_sur_mg_face(xyz1,mgface2,xyz2);
      if(ret==OK)
      {
        double distance = sqrt((xyz2[0]-xyz1[0])*(xyz2[0]-xyz1[0])+(xyz2[1]-xyz1[1])*(xyz2[1]-xyz1[1])+(xyz2[2]-xyz1[2])*(xyz2[2]-xyz1[2]));
        vector_distance.push_back(distance);
        moyenne_distance+=distance;
      }
    }
  }
  if(vector_distance.size()==0) return FAIL;
  moyenne_distance=moyenne_distance/vector_distance.size();
  std::vector<double>::iterator it;
  for(it=vector_distance.begin();it!=vector_distance.end();it++)
  {
    ecart_type_distance+=(*it-moyenne_distance)*(*it-moyenne_distance);
  }
  ecart_type_distance=sqrt(ecart_type_distance*(1.0/(vector_distance.size()-1.0)));
  return OK;
}


int OT_GEOMETRIE::projection_orthogonale_sur_mg_arete(double* xyz, MG_ARETE* mgarete, double& t)
{
#ifdef ALL_OCC
  gp_Pnt point(xyz[0],xyz[1],xyz[2]);
  TopoDS_Edge edge = ((OCC_ARETE*)mgarete)->get_TopoDS_Edge();
  GeomAPI_ProjectPointOnCurve projecteur;
  double t0,t1;
  Handle(Geom_Curve) courbe_edge = BRep_Tool::Curve(edge,t0,t1);
  projecteur.Init(point,courbe_edge,t0,t1);
  if(projecteur.NbPoints()==0) return FAIL;
  t = projecteur.LowerDistanceParameter();
  return OK;
#endif
}

int OT_GEOMETRIE::projection_orthogonale_sur_mg_face(double* xyz, MG_FACE* mgface,double* xyz2)
{
#ifdef ALL_OCC
  gp_Pnt point(xyz[0],xyz[1],xyz[2]);
  TopoDS_Face face = ((OCC_FACE*)mgface)->get_TopoDS_Face();
  GeomAPI_ProjectPointOnSurf projecteur;
  double u0,u1,v0,v1;
  ShapeAnalysis::GetFaceUVBounds(face,u0,u1,v0,v1);
  Handle(Geom_Surface) surface = BRep_Tool::Surface(face);
//   projecteur.Init(point,surface,u0,u1,v0,v1);
  projecteur.Init(point,surface);
  if(!projecteur.IsDone()) return FAIL;
  if(projecteur.NbPoints()==0) return FAIL;
  gp_Pnt point2 = projecteur.NearestPoint();
  xyz2[0]= point2.X();
  xyz2[1]= point2.Y();
  xyz2[2]= point2.Z();
  return OK;
#endif
}

int OT_GEOMETRIE::projection_au_plus_pres_sur_mg_face(double* xyz1,MG_FACE* mgface1, MG_FACE* mgface2, double* xyz2)
{
#ifdef ALL_OCC
  gp_Pnt point(xyz1[0],xyz1[1],xyz1[2]);  
  BRepBuilderAPI_MakeVertex makevertex(point);
  TopoDS_Vertex vertex=makevertex.Vertex();
  TopoDS_Face face1 = ((OCC_FACE*)mgface1)->get_TopoDS_Face();
  TopoDS_Face face2 = ((OCC_FACE*)mgface2)->get_TopoDS_Face();
  BRepClass_FaceClassifier face1_classifier(face1,point,Precision::Confusion());
  if(face1_classifier.State()==TopAbs_OUT) return FAIL;
  BRepExtrema_DistShapeShape distshapeshape(vertex,face2);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  double distance_min=std::numeric_limits< double >::max();
  double tmp[3];
  for(int i=1;i==distshapeshape.NbSolution();i++)
  {
    gp_Pnt point2= distshapeshape.PointOnShape2(i);
    tmp[0]=point2.X();
    tmp[1]=point2.Y();
    tmp[2]=point2.Z();
    double distance=sqrt((tmp[0]-xyz1[0])*(tmp[0]-xyz1[0])+(tmp[1]-xyz1[1])*(tmp[1]-xyz1[1])+(tmp[2]-xyz1[2])*(tmp[2]-xyz1[2]));
    if(distance<distance_min)
    {
      distance_min=distance;
      xyz2[0]=tmp[0];
      xyz2[1]=tmp[1];
      xyz2[2]=tmp[2];
    }
  }
  return OK;
#endif
}

int OT_GEOMETRIE::projection_au_plus_pres_sur_mg_arete(double* xyz, MG_ARETE* mgarete, double* xyz2)
{
#ifdef ALL_OCC
  gp_Pnt point(xyz[0],xyz[1],xyz[2]);  
  BRepBuilderAPI_MakeVertex makevertex(point);
  TopoDS_Vertex vertex=makevertex.Vertex();
  TopoDS_Edge edge = ((OCC_ARETE*)mgarete)->get_TopoDS_Edge();
  BRepExtrema_DistShapeShape distshapeshape(vertex,edge);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  double distance_min=std::numeric_limits< double >::max();
  double tmp[3];
  for(int i=1;i==distshapeshape.NbSolution();i++)
  {
    gp_Pnt point2= distshapeshape.PointOnShape2(i);
    tmp[0]=point2.X();
    tmp[1]=point2.Y();
    tmp[2]=point2.Z();
    double distance=sqrt((tmp[0]-xyz[0])*(tmp[0]-xyz[0])+(tmp[1]-xyz[1])*(tmp[1]-xyz[1])+(tmp[2]-xyz[2])*(tmp[2]-xyz[2]));
    if(distance<distance_min)
    {
      distance_min=distance;
      xyz2[0]=tmp[0];
      xyz2[1]=tmp[1];
      xyz2[2]=tmp[2];
    }
  }
  return OK;
#endif
}


void OT_GEOMETRIE::get_param_face(MG_FACE* mgface, double &umin, double &umax, double &vmin, double &vmax)
{
#ifdef ALL_OCC  
  TopoDS_Face face = ((OCC_FACE*)mgface)->get_TopoDS_Face();
  ShapeAnalysis::GetFaceUVBounds(face,umin,umax,vmin,vmax);
#endif
}


int OT_GEOMETRIE::get_distance_min_mg_eletopo_mg_eletopo(MG_ELEMENT_TOPOLOGIQUE* mgeletopo1, MG_ELEMENT_TOPOLOGIQUE* mgeletopo2, double& distance)
{
#ifdef ALL_OCC 
  TopoDS_Shape shape1;
  TopoDS_Shape shape2;
  if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
  {
    shape1=dynamic_cast<OCC_SOMMET*>(mgeletopo1)->get_TopoDS_Vertex();
  }
  else if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
  {
    shape1=dynamic_cast<OCC_ARETE*>(mgeletopo1)->get_TopoDS_Edge();
  }
  else if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    shape1=dynamic_cast<OCC_FACE*>(mgeletopo1)->get_TopoDS_Face();
  }
  else if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::VOLUME)
  {
    shape1=dynamic_cast<OCC_VOLUME*>(mgeletopo1)->get_TopoDS_Solid();
  }
  else return FAIL;
  if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
  {
    shape2=dynamic_cast<OCC_SOMMET*>(mgeletopo2)->get_TopoDS_Vertex();
  }
  else if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
  {
    shape2=dynamic_cast<OCC_ARETE*>(mgeletopo2)->get_TopoDS_Edge();
  }
  else if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    shape2=dynamic_cast<OCC_FACE*>(mgeletopo2)->get_TopoDS_Face();
  }
  else if(mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::VOLUME)
  {
    shape2=dynamic_cast<OCC_VOLUME*>(mgeletopo2)->get_TopoDS_Solid();
  }
  else return FAIL;

  BRepExtrema_DistShapeShape distshapeshape(shape2,shape1,Extrema_ExtFlag_MINMAX);
  distshapeshape.SetDeflection(Precision::Angular());
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  distance=std::numeric_limits< double >::max();
  for(int i=1;i<distshapeshape.NbSolution()+1;i++)
  {
    gp_Pnt p1=distshapeshape.PointOnShape1(i);
    gp_Pnt p2=distshapeshape.PointOnShape2(i);
    double dist=sqrt((p2.X()-p1.X())*(p2.X()-p1.X())+(p2.Y()-p1.Y())*(p2.Y()-p1.Y())+(p2.Z()-p1.Z())*(p2.Z()-p1.Z()));
    if(dist<distance)distance=dist;
  }
  distance=distshapeshape.Value();
  
  
//   Bnd_Box box_shape1;
//   BRepBndLib::Add(shape1, box_shape1);
//   Bnd_Box box_shape2;
//   BRepBndLib::Add(shape2, box_shape2);
//   Bnd_Box box_all;
//   box_all.Add(box_shape1);
//   box_all.Add(box_shape2);
//   double dist_ref=box_all.CornerMin().Distance(box_all.CornerMax());
//   BRepExtrema_DistanceSS distss(shape1,shape2,box_shape1,box_shape2,dist_ref,Extrema_ExtFlag_MINMAX);
//   if(!distss.IsDone()) return FAIL;
//   distance=distss.DistValue();
  
  return OK;  
#endif
}


int OT_GEOMETRIE::get_distance_min_mg_volume_mg_volume(MG_VOLUME* mgvolume1, MG_VOLUME* mgvolume2,double &distance)
{
#ifdef ALL_OCC  
  TopoDS_Solid solid1 = ((OCC_VOLUME*)mgvolume1)->get_TopoDS_Solid();
  TopoDS_Solid solid2 = ((OCC_VOLUME*)mgvolume2)->get_TopoDS_Solid();
  BRepExtrema_DistShapeShape distshapeshape(solid1,solid2,Extrema_ExtFlag_MIN);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  distance=distshapeshape.Value();  
  return OK;
#endif
}

int OT_GEOMETRIE::get_distance_min_mg_face_mg_volume(MG_FACE* mgface, MG_VOLUME* mgvolume, double& distance)
{
#ifdef ALL_OCC  
  TopoDS_Face face = ((OCC_FACE*)mgface)->get_TopoDS_Face();
  TopoDS_Solid solid = ((OCC_VOLUME*)mgvolume)->get_TopoDS_Solid();
  BRepExtrema_DistShapeShape distshapeshape(face,solid);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  distance=distshapeshape.Value();  
  return OK;
#endif
}

int OT_GEOMETRIE::get_distance_min_mg_face_mg_face(MG_FACE* mgface1, MG_FACE* mgface2, double& distance)
{
#ifdef ALL_OCC  
  TopoDS_Face face1 = ((OCC_FACE*)mgface1)->get_TopoDS_Face();
  TopoDS_Face face2 = ((OCC_FACE*)mgface2)->get_TopoDS_Face();
  BRepExtrema_DistShapeShape distshapeshape(face1,face2);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  distance=distshapeshape.Value();  
  return OK;
#endif
}


int OT_GEOMETRIE::get_distance_min_liste_mg_face_mg_volume(TPL_MAP_ENTITE< MG_FACE* >& map_face, MG_VOLUME* mgvolume,double &distance)
{
  distance=std::numeric_limits<double>::max();
  TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it_face;
  for(MG_FACE* face=map_face.get_premier(it_face);face!=NULL;face=map_face.get_suivant(it_face))
  {
    double distance_i;
    if(OT_GEOMETRIE::get_distance_min_mg_face_mg_volume(face,mgvolume,distance_i)==FAIL) return FAIL;
    if(distance_i<distance) distance=distance_i;
  }
  return OK;
}

int OT_GEOMETRIE::get_distance_min_mg_arete_mg_arete(MG_ARETE* mgarete1, MG_ARETE* mgarete2, double& distance)
{
  #ifdef ALL_OCC  
  TopoDS_Edge edge1 = ((OCC_ARETE*)mgarete1)->get_TopoDS_Edge();
  TopoDS_Edge edge2 = ((OCC_ARETE*)mgarete2)->get_TopoDS_Edge();
  BRepExtrema_DistShapeShape distshapeshape(edge1,edge2);
  distshapeshape.Perform();
  if(!distshapeshape.IsDone()) return FAIL;
  distance=distshapeshape.Value();  
  return OK;
#endif
}


double OT_GEOMETRIE::get_volume_intersection(MG_VOLUME* mgvolume1, MG_VOLUME* mgvolume2, double eps)
{
#ifdef ALL_OCC 
  TopoDS_Solid solid1 = ((OCC_VOLUME*)mgvolume1)->get_TopoDS_Solid();
  TopoDS_Solid solid2 = ((OCC_VOLUME*)mgvolume2)->get_TopoDS_Solid();
  BRepAlgoAPI_Common brep_common(solid1,solid2);
  if(!brep_common.IsDone()) return 0.0;
  if(brep_common.Shape().IsNull()) std::cerr << "NULL SHAPE !!!!!!!!!!!" << std::endl;
  GProp_GProps props;
  BRepGProp::VolumeProperties(brep_common.Shape(),props,eps);
  return props.Mass();
#endif
}

double OT_GEOMETRIE::get_volume(MG_VOLUME* mgvolume,double eps)
{
#ifdef ALL_OCC 
  TopoDS_Solid solid = ((OCC_VOLUME*)mgvolume)->get_TopoDS_Solid();
  GProp_GProps props;
  BRepGProp::VolumeProperties(solid,props,eps);
  return props.Mass();
#endif
}

double OT_GEOMETRIE::get_longueur(MG_ARETE* mgarete)
{
#ifdef ALL_OCC
  TopoDS_Edge edge = ((OCC_ARETE*)mgarete)->get_TopoDS_Edge();
  GProp_GProps props;
  BRepGProp::LinearProperties(edge,props);
  return props.Mass();
#endif
}

double OT_GEOMETRIE::get_aire(MG_FACE* mgface,double eps)
{
  #ifdef ALL_OCC
  TopoDS_Face face = ((OCC_FACE*)mgface)->get_TopoDS_Face();
  GProp_GProps props;
  BRepGProp::SurfaceProperties(face,props,eps);
  return props.Mass();
#endif
}

int OT_GEOMETRIE::est_dans_mg_face(double* xyz, MG_FACE* mgface)
{
#ifdef ALL_OCC
  TopoDS_Face face = ((OCC_FACE*)mgface)->get_TopoDS_Face();
  gp_Pnt point(xyz[0],xyz[1],xyz[2]);  
  BRepClass_FaceClassifier face1_classifier(face,point,Precision::Confusion());
  if(face1_classifier.State()==TopAbs_OUT) return FAIL;
  return OK;
#endif
}


void OT_GEOMETRIE::get_map_mg_ele_topo_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo,
                                                   TPL_MAP_ENTITE< MG_SOMMET* >& map_sommet,
                                                   TPL_MAP_ENTITE< MG_ARETE* >& map_arete,
                                                   TPL_MAP_ENTITE< MG_FACE* >& map_face)
{
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> map_ele_topo;
  mgeletopo->get_topologie_sousjacente(&map_ele_topo);
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
  for(MG_ELEMENT_TOPOLOGIQUE *ele=map_ele_topo.get_premier(it);ele!=NULL;ele=map_ele_topo.get_suivant(it))
  {
    if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
    {
      map_sommet.ajouter(dynamic_cast<MG_SOMMET*>(ele));
    }
    else if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
    {
      map_arete.ajouter(dynamic_cast<MG_ARETE*>(ele));
    }
    else if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
    {
      map_face.ajouter(dynamic_cast<MG_FACE*>(ele));
    }
  }
}


void OT_GEOMETRIE::get_map_mg_sommet_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_SOMMET* >& map_sommet)
{
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> map_ele_topo;
  mgeletopo->get_topologie_sousjacente(&map_ele_topo);
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
  for(MG_ELEMENT_TOPOLOGIQUE *ele=map_ele_topo.get_premier(it);ele!=NULL;ele=map_ele_topo.get_suivant(it))
  {
    if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET)
    {
      map_sommet.ajouter(dynamic_cast<MG_SOMMET*>(ele));
    }
  }
}

void OT_GEOMETRIE::get_map_mg_arete_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_ARETE* >& map_arete)
{
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> map_ele_topo;
  mgeletopo->get_topologie_sousjacente(&map_ele_topo);
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
  for(MG_ELEMENT_TOPOLOGIQUE *ele=map_ele_topo.get_premier(it);ele!=NULL;ele=map_ele_topo.get_suivant(it))
  {
    if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE)
    {
      map_arete.ajouter(dynamic_cast<MG_ARETE*>(ele));
    }
  }
}

void OT_GEOMETRIE::get_map_mg_face_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_FACE* >& map_face)
{
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*> map_ele_topo;
  mgeletopo->get_topologie_sousjacente(&map_ele_topo);
  TPL_MAP_ENTITE<MG_ELEMENT_TOPOLOGIQUE*>::ITERATEUR it;
  for(MG_ELEMENT_TOPOLOGIQUE *ele=map_ele_topo.get_premier(it);ele!=NULL;ele=map_ele_topo.get_suivant(it))
  {
    if(ele->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
    {
      map_face.ajouter(dynamic_cast<MG_FACE*>(ele));
    }
  }
}

void OT_GEOMETRIE::get_map_mg_point_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_POINT* >& map_point)
{
  TPL_MAP_ENTITE<MG_SOMMET*> map_sommet;
  get_map_mg_sommet_sous_jacent(mgeletopo,map_sommet);
  TPL_MAP_ENTITE<MG_SOMMET*>::ITERATEUR it;
  for(MG_SOMMET* sommet=map_sommet.get_premier(it);sommet!=NULL;sommet=map_sommet.get_suivant(it))
    map_point.ajouter(sommet->get_point());
}

void OT_GEOMETRIE::get_map_mg_courbe_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_COURBE* >& map_courbe)
{
  TPL_MAP_ENTITE<MG_ARETE*> map_arete;
  get_map_mg_arete_sous_jacent(mgeletopo,map_arete);
  TPL_MAP_ENTITE<MG_ARETE*>::ITERATEUR it;
  for(MG_ARETE* arete=map_arete.get_premier(it);arete!=NULL;arete=map_arete.get_suivant(it))
    map_courbe.ajouter(arete->get_courbe());
}

void OT_GEOMETRIE::get_map_mg_surface_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mgeletopo, TPL_MAP_ENTITE< MG_SURFACE* >& map_surface)
{
  TPL_MAP_ENTITE<MG_FACE*> map_face;
  get_map_mg_face_sous_jacent(mgeletopo,map_face);
  TPL_MAP_ENTITE<MG_FACE*>::ITERATEUR it;
  for(MG_FACE* face=map_face.get_premier(it);face!=NULL;face=map_face.get_suivant(it))
    map_surface.ajouter(face->get_surface());
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_SOMMET* sommet,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  if(map_mg_element_geometrique!=NULL) map_mg_element_geometrique->ajouter(sommet->get_point());
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_COSOMMET* cosommet,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  if(map_mg_element_topologique!=NULL) map_mg_element_topologique->ajouter(cosommet->get_sommet());
  get_map_mg_element_sous_jacent(cosommet->get_sommet(),
                                 map_mg_element_topologique,
                                 map_mg_element_cotopologique,
                                 map_mg_element_geometrique);
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_ARETE* arete,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  if(map_mg_element_cotopologique!=NULL)
  {
    map_mg_element_cotopologique->ajouter(arete->get_cosommet1());
    map_mg_element_cotopologique->ajouter(arete->get_cosommet2());
  }
  if(map_mg_element_geometrique!=NULL) map_mg_element_geometrique->ajouter(arete->get_courbe());
  get_map_mg_element_sous_jacent(arete->get_cosommet1(),
                                 map_mg_element_topologique,
                                 map_mg_element_cotopologique,
                                 map_mg_element_geometrique);
  get_map_mg_element_sous_jacent(arete->get_cosommet2(),
                                 map_mg_element_topologique,
                                 map_mg_element_cotopologique,
                                 map_mg_element_geometrique);
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_COARETE* coarete,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  if(map_mg_element_topologique!=NULL) map_mg_element_topologique->ajouter(coarete->get_arete());
  get_map_mg_element_sous_jacent(coarete->get_arete(),
                                 map_mg_element_topologique,
                                 map_mg_element_cotopologique,
                                 map_mg_element_geometrique);
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_BOUCLE* boucle,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<boucle->get_nb_mg_coarete();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(boucle->get_mg_coarete(i));
    get_map_mg_element_sous_jacent(boucle->get_mg_coarete(i),
                                   map_mg_element_topologique,
                                   map_mg_element_cotopologique,
                                   map_mg_element_geometrique);
  }
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_POUTRE* poutre,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<poutre->get_nb_mg_boucle();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(poutre->get_mg_boucle(i));
    get_map_mg_element_sous_jacent(poutre->get_mg_boucle(i),
                                   map_mg_element_topologique,
                                   map_mg_element_cotopologique,
                                   map_mg_element_geometrique);
  }
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_FACE* face,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<face->get_nb_mg_boucle();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(face->get_mg_boucle(i));
    get_map_mg_element_sous_jacent(face->get_mg_boucle(i),
                                   map_mg_element_topologique,
                                   map_mg_element_cotopologique,
                                   map_mg_element_geometrique);
  }
  if(map_mg_element_geometrique!=NULL) map_mg_element_geometrique->ajouter(face->get_surface());
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_COFACE* coface,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  if(map_mg_element_topologique!=NULL) map_mg_element_topologique->ajouter(coface->get_face());
  get_map_mg_element_sous_jacent(coface->get_face(),
                                 map_mg_element_topologique,
                                 map_mg_element_cotopologique,
                                 map_mg_element_geometrique);
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_COQUILLE* coquille,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<coquille->get_nb_mg_coface();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(coquille->get_mg_coface(i));
    get_map_mg_element_sous_jacent(coquille->get_mg_coface(i),
                                   map_mg_element_topologique,
                                   map_mg_element_cotopologique,
                                   map_mg_element_geometrique);
  }
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_VOLUME* volume,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<volume->get_nb_mg_coquille();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(volume->get_mg_coquille(i));
    get_map_mg_element_sous_jacent(volume->get_mg_coquille(i),
                                  map_mg_element_topologique,
                                  map_mg_element_cotopologique,
                                  map_mg_element_geometrique);
  }  
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_COQUE* coque,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  for(int i=0;i<coque->get_nb_mg_coquille();i++)
  {
    if(map_mg_element_cotopologique!=NULL) map_mg_element_cotopologique->ajouter(coque->get_mg_coquille(i));
    get_map_mg_element_sous_jacent(coque->get_mg_coquille(i),
                                  map_mg_element_topologique,
                                  map_mg_element_cotopologique,
                                  map_mg_element_geometrique);
  }  
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_ELEMENT_TOPOLOGIQUE* mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  switch(mg_element_topologique->get_type())
  {
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::SOMMET:
    {
      MG_SOMMET* sommet = (MG_SOMMET*)mg_element_topologique;
      get_map_mg_element_sous_jacent(sommet,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::ARETE:
    {
      MG_ARETE* arete = (MG_ARETE*)mg_element_topologique;
      get_map_mg_element_sous_jacent(arete,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE:
    {
      MG_FACE* face = (MG_FACE*)mg_element_topologique;
      get_map_mg_element_sous_jacent(face,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::VOLUME:
    {
      MG_VOLUME* volume = (MG_VOLUME*)mg_element_topologique;
      get_map_mg_element_sous_jacent(volume,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::COQUE:
    {
      MG_COQUE* coque = (MG_COQUE*)mg_element_topologique;
      get_map_mg_element_sous_jacent(coque,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::POUTRE:
    {
      MG_POUTRE* poutre = (MG_POUTRE*)mg_element_topologique;
      get_map_mg_element_sous_jacent(poutre,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
  }
}

void OT_GEOMETRIE::get_map_mg_element_sous_jacent(MG_ELEMENT_COTOPOLOGIQUE* mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_TOPOLOGIQUE* >* map_mg_element_topologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_COTOPOLOGIQUE* >* map_mg_element_cotopologique,
                                                  TPL_MAP_ENTITE< MG_ELEMENT_GEOMETRIQUE* >* map_mg_element_geometrique)
{
  switch(mg_element_cotopologique->get_type())
  {
    case MG_ELEMENT_COTOPOLOGIQUE::TYPE_ELEMENT_COTOPOLOGIQUE::COSOMMET:
    {
      MG_COSOMMET* cosommet = (MG_COSOMMET*)mg_element_cotopologique;
      get_map_mg_element_sous_jacent(cosommet,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_COTOPOLOGIQUE::TYPE_ELEMENT_COTOPOLOGIQUE::COARETE:
    {
      MG_COARETE* coarete = (MG_COARETE*)mg_element_cotopologique;
      get_map_mg_element_sous_jacent(coarete,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_COTOPOLOGIQUE::TYPE_ELEMENT_COTOPOLOGIQUE::BOUCLE:
    {
      MG_BOUCLE* boucle = (MG_BOUCLE*)mg_element_cotopologique;
      get_map_mg_element_sous_jacent(boucle,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_COTOPOLOGIQUE::TYPE_ELEMENT_COTOPOLOGIQUE::COFACE:
    {
      MG_COFACE* coface = (MG_COFACE*)mg_element_cotopologique;
      get_map_mg_element_sous_jacent(coface,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
    case MG_ELEMENT_COTOPOLOGIQUE::TYPE_ELEMENT_COTOPOLOGIQUE::COQUILLE:
    {
      MG_COQUILLE* coquille = (MG_COQUILLE*)mg_element_cotopologique;
      get_map_mg_element_sous_jacent(coquille,
                                     map_mg_element_topologique,
                                     map_mg_element_cotopologique,
                                     map_mg_element_geometrique);
      break;
    }
  }
}


Revision:	951
Committed:	Fri Aug 10 15:17:17 2018 UTC (6 years, 9 months ago) by couturad
File size:	29909 byte(s)
Log Message:	-> Ajout de Project Chrono (voir CMakeLists.txt). -> Ajout d'un générateur de microstructure basé sur la dynamique des corps rigides (MSTRUCT_GENERATEUR_DCR). -> Ajout d'un opérateur de décallage de la topologie (MG_CG_OP_TRANSF_DECALLAGE). -> Retrait de «using namespace std» (conflit avec namespace chrono) et modification des fichiers affectés. -> Modification de mailleur2d.cpp afin d'enregistrer un fichier MAGiC (void.magic) lorsque le nombre d'itération dépasse la valeur maximale.