-
Notifications
You must be signed in to change notification settings - Fork 0
/
MeshAdapt.h
49 lines (44 loc) · 1.57 KB
/
MeshAdapt.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
#ifndef _MESHADAPT_H
#include "Solution.h" // En fait on en a besoin qu'à partir de timescheme !
#include "DataFile.h"
#include "Dense"
#include "Sparse"
class Mesh_Adapt
{
private:
Eigen::VectorXd _rho; // inutile ; peut être pas car on fait une référence vers sol
//Vecteur contenant les taille dy de chaques mailles
//_Dy=nouveau maillage dans la boucle Dyk+1(n)
//_Dystar=maillage intermédiaire
//_Dyold=maillage initiale (CI pour l'advection)
//_Dyprevious=maillage precedent le nouveau dans la boucle Dyk(n)
Eigen::VectorXd _Dy,_Dystar,_Dyold, _Dyprevious;
//Vecteur contenant les positions des "noeuds" de chaques mailles selon Y
Eigen::VectorXd _Y,_Yold; // _Y1 pas bon défini sur le mauvais maillage
Eigen::VectorXd _v;
//Pointeur vers le fichier de données
DataFile* _df;
public:
//Constructeur
Mesh_Adapt(DataFile* data_file);
Eigen::VectorXd Maillage_non_uniforme();
//Eigen::VectorXd Maillage_non_uniforme2();
void Update(Solution & sol);
void Update2(Solution & sol);
Eigen::VectorXd Derive_y_2(Eigen::VectorXd rho);
const Eigen::VectorXd Get_Dy() const {return _Dy;};
//Fonction numéro de cellule
const Eigen::VectorXd Get_Dystar() const {return _Dystar;};
const Eigen::VectorXd Get_Dyold() const {return _Dyold;};
void Maillage_Dystar();
double NormLinf();
int cellule (double distance);
void Update_Dyold();
void Update_Dystar_vitesse();
void Update_Dyprevious();
void vitesse();
const Eigen::VectorXd Get_vitesse() const {return _v;}
void Affichage(std::string text,Solution sol);
};
#define _MESHADAPT_H
#endif