Commit 537ac9bf authored by Stephan Rave's avatar Stephan Rave
Browse files

[discretizers.builtin] fix AdvectionOperatorP1

parent 4fc3893a
Pipeline #65092 passed with stages
in 80 minutes and 9 seconds
......@@ -718,7 +718,7 @@ class AdvectionOperatorP1(NumpyMatrixBasedOperator):
self.logger.info('Calculate all local scalar products beween gradients ...')
D = self.advection_function(self.grid.centers(0), mu=mu)
SF_INTS = - np.einsum('pc,eqi,c,e,ec->eqp', SFQ, SF_GRADS, w, g.integration_elements(0), D).ravel()
SF_INTS = - np.einsum('pc,eqi,c,e,ei->eqp', SFQ, SF_GRADS, w, g.integration_elements(0), D).ravel()
del D
del SF_GRADS
......
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