Update tab.cpp

src/tab.cpp

@@ -6,24 +6,69 @@
 
 using namespace std;
 
-Materiaux::Materiaux(const char* materiaux,double lambda, double kho, double c)
+Materiaux::Materiaux(const char* materiaux,double lambda, double rho, double c)
 {
     materiaux_=materiaux;
     lambda_= lambda;
-    kho_ = kho;
+    rho_ = rho;
     c_ = c;
 }
 
+/**
+Résolution de l'équation de la chaleur à 1 dimension implicite, avec la décomposition LU
+
+@param lambda: coefficient de diffusivité thermique
+@param T: vecteur qui contient les 101 valeures de la température en Kelvin à l'instant n, le long de la barre
+@return: un vecteur qui contient les 101 valeures de la température en Kelvin à l'instant n+1, le long de la barre
+*/
+vector<double> resolution(double lambda, vector<double> T){
+    int n= T.size();
+    vector<double> d(n,0);
+    vector<double> l(n-1,0);
+    d[0]=1+2*lambda;
+
+    for (int i =1;i<n;i++){
+        l[i-1]=(-1)*lambda/d[i-1];
+        d[i]=1+2*lambda+lambda*l[i-1];
+    }
+    vector<double> y(n,0);
+    y[0]=T[0];
+    for (int i =1;i<n;i++){
+        y[i]=T[i]-l[i-1]*y[i-1];
+    }
+   vector<double> x(n,0);
+   x[n-1]=y[n-1]/d[n-1];
+   for (int i =n-2;i>0;i--){
+        x[i]=(y[i]+lambda*x[i+1])/d[i];
+    }
+
+    return x;
+}
 
 /**
-Résolution de l'équation de la chaleur à 1 dimension
+Fait la somme de 2 vecteurs de doubles
 
-La fonction prend en compte les conditions initiales
+@param F: vecteur qui contient les 101 valeures d'apport extérieur de chaleur à l'instant n, le long de la barre
+@param T: vecteur qui contient les 101 valeures de la température en Kelvin à l'instant n, le long de la barre
+@return: la somme des 2 vecteurs
+*/
+vector<double> sum_vect(vector <double> T,vector <double> F){
 
-@param  mat: un matériau de type Materiaux
-@param   u0: la valeur initial de U 
-@param f valeur en Kelvin qui simule un ajout de chaleur
-@return ecris sur le flux de sortie une matrice U[i][j] qui contient la chaleur en Kelvin à l'instant i et à la position x=j
+    vector <double> res(T.size(),0);
+    for (size_t i=0;i<T.size();i++){
+        res[i]= T[i]+F[i];
+    }
+    return res;
+
+}
+
+/**
+Calcule les valeurs de la températures pour chaque chaque couple (x,t), où x,la position, prend 100 valeurs de 0 à 1, et T prend 101 valeurs de 0 à 16
+
+@param mat: materiau de la barre
+@param f: valeur de la température de la source exterieure de chaleur (en K)
+@param u0: valeur de la température ambiante (en K)
+@return: créé les fichiers contenant les tableaux de valeur, ne retourne rien
 */
 void settab(Materiaux mat, int u0, double f)
 {
@@ -32,51 +77,33 @@ void settab(Materiaux mat, int u0, double f)
     double dt, dx;
     dt= T/100.;
     dx= L/100.;
-    double** temp;
-    double coef_1= mat.getlambda()/(dx*dx);
-    double coef_2= dt/(mat.getc()*mat.getkho());
-    double F=16*f*f;
-
-    temp = new double*[100];
-    for (int i = 0; i < 100; i++)
-        temp[i] = new double[101];
-
-      for (int j = 0; j < 101; j++){
-        temp[0][j]=u0;
+    double coef= (mat.getlambda()*dt)/(mat.getc()*mat.getrho()*(dx*dx));
+    double F2=16*f*f;
+    vector<double> F(101,0);
+    F[10]=F2/(mat.getc()*mat.getrho());
+    for(int j=50;j<61;j++){
+        F[j]=F2*0.75/(mat.getc()*mat.getrho());
     }
 
-    for (int i = 0; i < 100; i++){
-        temp[i][100]=u0;
-        temp[i][0]= u0;
-    }
+    vector<vector<double> > temp;
+    vector<double> tp;
+    temp.push_back(vector<double>(101,u0));
 
 
     for (int i = 1; i < 100; i++){
-        for (int j = 1; j < 10; j++){
-            temp[i][j]= temp[i-1][j] + coef_1*coef_2*(temp[i-1][j-1] -2*temp[i-1][j] + temp[i-1][j+1]);
-        }
-
-        temp[i][10]= temp[i-1][10] + coef_1*coef_2*(temp[i-1][9] -2*temp[i-1][10] + temp[i-1][11])+coef_2*F;
-
-        for (int j = 11; j < 50; j++){
-            temp[i][j]= temp[i-1][j] + coef_1*coef_2*(temp[i-1][j-1] -2*temp[i-1][j] + temp[i-1][j+1]);
-        }
-
-        for(int j=50;j<61;j++){
-            temp[i][j]= temp[i-1][j] + coef_1*coef_2*(temp[i-1][j-1] -2*temp[i-1][j] + temp[i-1][j+1])+coef_2*F*0.75;
-        }
-
-        for(int j=61;j<100;j++){
-            temp[i][j]= temp[i-1][j] + coef_1*coef_2*(temp[i-1][j-1] -2*temp[i-1][j] + temp[i-1][j+1]);
-        }
+        tp=temp[i-1];
+        tp[0]=tp[0]+coef*u0;
+        tp[100]=tp[100]+coef*u0;
+        temp.push_back(resolution(coef,sum_vect(tp,F)));
+        temp[i][0]=u0;
+        temp[i][100]=u0;
     }
-
     const char separateur(' ');
     std::ofstream sortie(mat.get_mat(),std::ios::out);
 
 
      for (int i = 0; i < 100; i++){
-        for (int j = 0; j < 100; j++){
+        for (int j = 0; j < 101; j++){
             sortie << temp[i][j] << separateur;
         }
         sortie<< std::endl;
@@ -84,10 +111,6 @@ void settab(Materiaux mat, int u0, double f)
 
      sortie.close();
 
-    for (int i = 0; i < 100; i++)
-         delete[] temp[i];
-    delete[] temp;
-
 }
 /**
 Résolution de l'équation de la chaleur à 2 dimension