diff --git a/grudge/models/advection.py b/grudge/models/advection.py index 707a729d369898597c2f3b7f236e9b1e82fc1b02..68b1a0bad18c4584974a70ccab45fe5975b92465 100644 --- a/grudge/models/advection.py +++ b/grudge/models/advection.py @@ -172,16 +172,15 @@ class VariableCoefficientAdvectionOperator(HyperbolicOperator): self.flux(pair)) return sym.InverseMassOperator()( - np.dot( - sym.stiffness_t(self.ambient_dim), sym.cse(self.v*u)) - - sym.FaceMassOperator()( - flux(sym.int_tpair(u)) - + flux(sym.bv_tpair(sym.BTAG_ALL, u, bc_in)) + np.dot(sym.stiffness_t(self.ambient_dim), sym.cse(self.v*u)) + - sym.FaceMassOperator()( + flux(sym.int_tpair(u)) + + flux(sym.bv_tpair(sym.BTAG_ALL, u, bc_in)) - # FIXME: Add back support for inflow/outflow tags - #+ flux(sym.bv_tpair(self.inflow_tag, u, bc_in)) - #+ flux(sym.bv_tpair(self.outflow_tag, u, bc_out)) - )) + # FIXME: Add back support for inflow/outflow tags + #+ flux(sym.bv_tpair(self.inflow_tag, u, bc_in)) + #+ flux(sym.bv_tpair(self.outflow_tag, u, bc_out)) + )) # }}} # vim: foldmethod=marker