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()
{

}

BOITE_3D OT_GEOMETRIE::get_boite_3D(MG_VOLUME* mgvolume)
{
  #ifdef ALL_OCC 
  TopoDS_Shape shape1;
  shape1=dynamic_cast<OCC_VOLUME*>(mgvolume)->get_TopoDS_Solid();
  Bnd_Box box_shape1;
  BRepBndLib::AddOptimal(shape1,box_shape1);
  return BOITE_3D(box_shape1.CornerMin().X(),box_shape1.CornerMin().Y(),box_shape1.CornerMin().Z(),box_shape1.CornerMax().X(),box_shape1.CornerMax().Y(),box_shape1.CornerMax().Z());
#endif
}


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;
  
  if(mgeletopo1->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE && mgeletopo2->get_type()==MG_ELEMENT_TOPOLOGIQUE::TYPE_ELEMENT_TOPOLOGIQUE::FACE)
  {
    MG_FACE* face1=(MG_FACE*)mgeletopo1;
    MG_FACE* face2=(MG_FACE*)mgeletopo2;
    get_distance_min_mg_face_mg_face_echantillonnage(face1,face2,distance);
  }
  else
  {
    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();
  }
  
// 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_mg_face_mg_face_echantillonnage(MG_FACE* mgface1, MG_FACE* mgface2, double& distance,int nb_pas)
{
#ifdef ALL_OCC  
  distance=std::numeric_limits<double>::max();
  double umin,umax,vmin,vmax;
  OT_GEOMETRIE::get_param_face(mgface1,umin,umax,vmin,vmax);
  if(mgface1->get_surface()->est_periodique_u())
  {
    umin=0;
    umax=mgface1->get_surface()->get_periode_u();
  }
  if(mgface1->get_surface()->est_periodique_v())
  {
    vmin=0;
    vmax=mgface1->get_surface()->get_periode_v();
  }
  double pas_u = (umax-umin)/nb_pas;
  double pas_v = (vmax-vmin)/nb_pas;
  double uv1[2];
  for(long i=0;i<nb_pas;i++)
  {
    uv1[0]=umin+i*pas_u;
    for(long j=0;j<nb_pas;j++)
    {
      uv1[1] = vmin+j*pas_v;
      double xyz1[3];      
      mgface1->evaluer(uv1,xyz1);
      if(OT_GEOMETRIE::est_dans_mg_face(xyz1,mgface1)==FAIL) continue;
      double xyz2[3];
      if(OT_GEOMETRIE::projection_au_plus_pres_sur_mg_face(xyz1,mgface1,mgface2,xyz2)==OK)
      {
            double dist = 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]));
        if(dist<distance) distance=dist;
      }
    }
  } 
  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:	968
Committed:	Sun Sep 16 15:27:49 2018 UTC (6 years, 7 months ago) by couturad
File size:	32552 byte(s)
Log Message:	Ajout d'une condition de sortie et d'un renvoi d'erreur pour le mailleur FEM. Subdivision des fichiers mstruct_analyse.h/.cpp en sous fichiers pour une meilleure lisibilite. Ajout d'une analyse des modules d'elasticite. Ajout d'une analyse de l'energie. Reconfiguration du main de microstructure.exe (suppression d'actions obsolètes). Reconfiguration des fichiers generer_nb_ves, post_process.