retrait des print inutiles

fluid/fluid_problem.c

@@ -1211,9 +1211,9 @@ void fluid_problem_apply_boundary_conditions(FluidProblem *problem)
 static void fluid_problem_compute_node_volume(FluidProblem *problem) {
   free(problem->node_volume);
   Mesh *mesh = problem->mesh;
-  problem->node_volume = (double*)malloc(sizeof(double)*mesh->n_nodes);
+  problem->node_volume = malloc(sizeof(double)*mesh->n_nodes);
   int independent_nodes = mesh->n_nodes - mesh->n_periodic;
-  double* periodic_node_volume = (double*)malloc(sizeof(double)*independent_nodes);
+  double* periodic_node_volume = malloc(sizeof(double)*independent_nodes);
   for (int i = 0; i < independent_nodes; ++i){
     periodic_node_volume[i] = 0;
   }
@@ -1223,9 +1223,10 @@ static void fluid_problem_compute_node_volume(FluidProblem *problem) {
     const int *el = &mesh->elements[iel*N_N];
     const double detj = mesh_dxidx(mesh, iel, dxidx);
     const double vol = node_volume_from_detj(detj);
-    for (int i = 0; i < N_N; ++i)
-      // problem->node_volume[el[i]] += vol;
+    for (int i = 0; i < N_N; ++i){
       periodic_node_volume[periodic_mapping[el[i]]] += vol;
+      // problem->node_volume[el[i]] += vol;
+    }
   }
   for(int i = 0; i < problem->mesh->n_nodes; ++i){
     problem->node_volume[i] = periodic_node_volume[periodic_mapping[i]];
@@ -1236,7 +1237,7 @@ static void fluid_problem_compute_node_volume(FluidProblem *problem) {
 static void compute_porosity(Mesh *mesh, double *node_volume, double *porosity, int n_particles, double *particle_position, double *particle_volume, int *particle_element_id, double *particle_uvw)
 {
   int independent_nodes = mesh->n_nodes - mesh->n_periodic;
-  double* periodic_porosity = (double*)malloc(sizeof(double)*independent_nodes);
+  double* periodic_porosity = malloc(sizeof(double)*independent_nodes);
   for (int i = 0; i < independent_nodes; ++i){
     periodic_porosity[i] = 0;
   }

fluid/mesh.c

@@ -182,7 +182,6 @@ Mesh *mesh_load_msh(const char *filename)
   }
   for (int i = 0; i < n_enodes; ++i) {
     m->elements[i] = get_node_id(m->elements[i], nodes_map, m->n_nodes);
-    printf("m->elements[%d] = %d\n",i,m->elements[i]);
   }
 
 

testcases/periodicTest/mesh.geo

 L = 1;
 H = 1;
 y = 0;
-lc = 1;//0.05;
+lc = 0.05;
 
 Point(1) = {0,0,0,lc};
 Point(2) = {L,0,0,lc};