Unverified Commit 962e6870 authored by Stephan Rave's avatar Stephan Rave Committed by GitHub

Merge pull request #1099 from pymor/fix_advection_operator

Fix AdvectionOperatorP1
parents 778fc6c6 537ac9bf
Pipeline #65119 failed with stages
in 56 minutes and 49 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