Commit 598361bb authored by Matthieu Constant's avatar Matthieu Constant
Browse files

good adapt

parent 92e0df31
Pipeline #4532 failed with stage
in 18 seconds
......@@ -1138,7 +1138,7 @@ void fluid_problem_adapt_gen_mesh(FluidProblem *problem, double lcmax, double lc
double C[N_N],Un[D][N_N],An[N_N],Pn[N_N];
for (int i = 0; i < N_N; ++i){
C[i] = porosity[el[i]];
Pn[i] = solution[el[i]*n_fields+D]-problem->g *problem->rho[0]*mesh->x[el[i]*3+1];
Pn[i] = solution[el[i]*n_fields+D]-0*problem->g *problem->rho[0]*mesh->x[el[i]*3+1];
if(problem->n_fluids ==2) An[i] = solution[el[i]*n_fields+D+1];
for (int j=0; j<D; ++j) {
Un[j][i] = solution[el[i]*n_fields+j];
......@@ -1212,7 +1212,7 @@ void fluid_problem_adapt_gen_mesh(FluidProblem *problem, double lcmax, double lc
double C[N_N],Un[D][N_N],An[N_N],Pn[N_N];
for (int i = 0; i < N_N; ++i){
C[i] = porosity[el[i]];
Pn[i] = solution[el[i]*n_fields+D]-problem->g *problem->rho[0]*mesh->x[el[i]*3+1];
Pn[i] = solution[el[i]*n_fields+D]-0*problem->g *problem->rho[0]*mesh->x[el[i]*3+1];
if(problem->n_fluids ==2) An[i] = solution[el[i]*n_fields+D+1];
for (int j=0; j<D; ++j) {
Un[j][i] = solution[el[i]*n_fields+j];
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment