Commit 02db5e4f authored by Jonathan Lambrechts's avatar Jonathan Lambrechts
Browse files

clean

parent caefa10b
Pipeline #4539 passed with stage
in 2 minutes and 11 seconds
......@@ -636,7 +636,6 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, double dt, d
dx[0][0]*dx[1][1] - dx[0][1]*dx[1][0]};
const double detbnd = sqrt(nn[0]*nn[0] + nn[1]*nn[1] + nn[2]*nn[2]);
const double n[3] = {-nn[0]/detbnd,-nn[1]/detbnd,-nn[2]/detbnd};
//printf("%g,%g,%g,%s\n",n[0],n[1],n[2],mesh->phys_name[iphys]);
#endif
for (int iqp = 0; iqp < N_LQP; ++iqp) {
double xi[D];
......@@ -880,7 +879,7 @@ static void fluid_problem_assemble_system(FluidProblem *problem, double *rhs, co
}
}
node_force_volume(problem, solution_old, dt, all_local_vector, all_local_matrix);
//compute_weak_boundary_conditions(problem, dt, all_local_vector, all_local_matrix);
compute_weak_boundary_conditions(problem, dt, all_local_vector, all_local_matrix);
fluid_problem_volume(problem, solution_old, dt, all_local_vector, all_local_matrix);
for (int iel=0; iel < mesh->n_elements; ++iel) {
double *local_vector = &all_local_vector[local_size*iel];
......
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