Commit b993a928 authored by Michel Henry's avatar Michel Henry
Browse files

remove Symmetric type in fluid C files

parent 320a778d
Pipeline #9357 passed with stages
in 3 minutes and 54 seconds
......@@ -226,7 +226,7 @@ static void f_boundary(WeakBoundary *wbnd, FluidProblem *problem,const double *n
for (int jD = 0; jD < D; ++jD)
c_du_o_c[iD][jD] = ds[(U+iD)*D+jD] -u[iD]/c*dc[jD];
}
if(wbnd->type != BND_SYMMETRY && wbnd->compute_viscous_term == 1){
if(wbnd->compute_viscous_term == 1){
for (int id = 0; id < D; ++id) {
f0[U+id] += sigma*(u[id]-uext[id]+s_c*n[id]);
f00[(U+id)*n_fields+U+id] += (vid<0?0:sigma);
......
......@@ -44,7 +44,7 @@ typedef struct {
} StrongBoundary;
typedef struct FluidProblem FluidProblem;
typedef enum {BND_WALL=0,BND_OPEN=1,BND_SYMMETRY=2} BoundaryType;
typedef enum {BND_WALL=0,BND_OPEN=1} BoundaryType;
typedef struct {
char *tag;
BoundaryType type;
......
......@@ -240,7 +240,7 @@ class FluidProblem :
pid = 0
bndcb = BNDCB(_Bnd(cb_or_value, self._dim).apply)
self.weak_cb_ref.append(bndcb)
self._lib.fluid_problem_set_weak_boundary(self._fp, tag.encode(), c_int(0), bndcb, c_int(vid), c_int(pid), c_int(-1), c_int(-1), c_int(-1))
self._lib.fluid_problem_set_weak_boundary(self._fp, tag.encode(), c_int(0), bndcb, c_int(vid), c_int(pid), c_int(-1), c_int(-1), c_int(0))
def set_open_boundary(self, tag, velocity=None, pressure=None, porosity=1, concentration=1, compute_viscous_term=1):
"""Sets the weak boundaries (=normal fluxes) for the fluid problem.
......
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