polycristal.cpp

Aller à la documentation de ce fichier.

//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//  MAGiC
//####//  Jean Christophe Cuilliere et Vincent FRANCOIS
//####//  Departement de Genie Mecanique - UQTR
//####//------------------------------------------------------------
//####//  MAGIC est un projet de recherche de l equipe ERICCA
//####//  du departement de genie mecanique de l Universite du Quebec a Trois Rivieres
//####//  http://www.uqtr.ca/ericca
//####//  http://www.uqtr.ca/
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//
//####//       polycristal.cpp
//####//
//####//------------------------------------------------------------
//####//------------------------------------------------------------
//####//    COPYRIGHT 2000-2024
//####//  jeu 13 jun 2024 11:54:00 EDT
//####//------------------------------------------------------------
//####//------------------------------------------------------------
#include "polycristal.h"
#include "poly_voro.h"
#include "poly_build_occ.h"
#include "poly_point.h"
 
#include <stdio.h>
#include <random>
#include <vector>
 
Polycristal::Polycristal(char* fichier,double dgval):POLY_AFFICHE(),nomfichier("resultat.brep"),dg(dgval)
{
    FILE *input = fopen(fichier, "rt");
    char chaine[500];
    fgets(chaine, 500, input);
    int nb_point, nTetra;
    sscanf(chaine,"%d %d", &nb_point, &nTetra);
 
    
    for (int i=0; i<nb_point; i++)
    {
      int num;
      double x, y, z;
      fgets(chaine,500,input);
      sscanf(chaine, "%d %lf %lf %lf", &num, &x, &y, &z);
      Poly_Point* pnt_xyz = new Poly_Point(x, y, z);
      list_points.push_back(pnt_xyz);
    }
    fclose(input);
 
  
  
}
 
 
Polycristal::Polycristal(std::vector<double> &listepoint,std::string nom,double dgval):POLY_AFFICHE(),nomfichier(nom),dg(dgval)
{
    int nb_point=listepoint.size()/3.;
    for (int i=0; i<nb_point; i++)
    {
      Poly_Point* pnt_xyz = new Poly_Point(listepoint[3*i],listepoint[3*i+1],listepoint[3*i+2]);
      list_points.push_back(pnt_xyz);
    }
    
 
}
 
Polycristal::Polycristal(int nbParticules,double dgval):POLY_AFFICHE(),nomfichier("resultat.brep"),dg(dgval)
{
  list_points = random_particules(nbParticules);
 
}
 
 
void Polycristal::construit(bool avecstep)
{
  Poly_Voro voro(list_points,dg);
  if (affichageactif) voro.active_affichage(affiche2);
  voro.construit();
  voro.fusion_noeuds();
  Poly_Build_OCC build_occ(&voro,nomfichier);
 if (affichageactif)  build_occ.active_affichage(affiche2);
 build_occ.construit(avecstep);
}
 
 
 
Polycristal::Polycristal(Polycristal& mdd)
{
}
 
 
Polycristal::~Polycristal()
{
}
 
std::vector<Poly_Point*> Polycristal::random_particules(int nbParticules)
{ 
  std::vector<Poly_Point*> list_points;
  
  double x_min= 0.0, x_max= 1.0;
  double y_min= 0.0, y_max= 1.0;
  double z_min= 0.0, z_max= 1.0;
 
  Poly_Point* pnt_1 = new Poly_Point(x_min, y_min, z_min);
  Poly_Point* pnt_2 = new Poly_Point(x_min, y_max, z_min);
  Poly_Point* pnt_3 = new Poly_Point(x_max, y_max, z_min); 
  Poly_Point* pnt_4 = new Poly_Point(x_max, y_min, z_min);   
  Poly_Point* pnt_5 = new Poly_Point(x_min, y_min, z_max);
  Poly_Point* pnt_6 = new Poly_Point(x_min, y_max, z_max);
  Poly_Point* pnt_7 = new Poly_Point(x_max, y_max, z_max);
  Poly_Point* pnt_8 = new Poly_Point(x_max, y_min, z_max);
  
  list_points.push_back(pnt_1);
  list_points.push_back(pnt_2);
  list_points.push_back(pnt_3);
  list_points.push_back(pnt_4);
  list_points.push_back(pnt_5);
  list_points.push_back(pnt_6);
  list_points.push_back(pnt_7);
  list_points.push_back(pnt_8);
  
  std::random_device rd;  // random number engine
  std::mt19937 gen(rd());
  std::uniform_real_distribution<> dis_x(x_min, x_max);
  std::uniform_real_distribution<> dis_y(y_min, y_max);
  std::uniform_real_distribution<> dis_z(z_min, z_max);
  
  for(int i=0; i < nbParticules; i++) {
      double x = dis_x(gen) ;
      double y = dis_y(gen) ;
      double z = dis_z(gen) ;
      Poly_Point* pnt_xyz = new Poly_Point(x, y, z);
      list_points.push_back(pnt_xyz);
  }
  return list_points;
}