Commit eae35a08 by Jonathan Lambrechts

### Revert "separation Ã2 fluides"

`This reverts commit 99f25a0e.`
parent 99f25a0e
 ... ... @@ -248,7 +248,6 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, double dt, d const int iel = bnd[0]; const int i0 = bnd[1]; const int i1 = (bnd[1]+1) % 3; const int i2 = (bnd[1]+2) % 3; const uint32_t *el = &mesh->elements[iel*N_N]; const int nodes[2] = {el[i0],el[i1]}; const double *x[2] = {&mesh->x[nodes[0]*3], &mesh->x[nodes[1]*3]}; ... ... @@ -288,10 +287,8 @@ static void compute_weak_boundary_conditions(FluidProblem *problem, double dt, d local_matrix[(i0*n_fields+Q)*local_size + (i0*n_fields+P)] += dphi[i0][k]*n[k]*problem->epsilon/2.*l; local_matrix[(i0*n_fields+Q)*local_size + (i1*n_fields+P)] += dphi[i1][k]*n[k]*problem->epsilon/2.*l; local_matrix[(i0*n_fields+Q)*local_size + (i2*n_fields+P)] += dphi[i2][k]*n[k]*problem->epsilon/2.*l; local_matrix[(i1*n_fields+Q)*local_size + (i0*n_fields+P)] += dphi[i0][k]*n[k]*problem->epsilon/2.*l; local_matrix[(i1*n_fields+Q)*local_size + (i1*n_fields+P)] += dphi[i1][k]*n[k]*problem->epsilon/2.*l; local_matrix[(i1*n_fields+Q)*local_size + (i2*n_fields+P)] += dphi[i2][k]*n[k]*problem->epsilon/2.*l; local_vector[i0+Q*N_SF] += dp[k]*n[k]*problem->epsilon/2.*l; local_vector[i1+Q*N_SF] += dp[k]*n[k]*problem->epsilon/2.*l; ... ...