Mercurial > repos > public > sbplib

diff +scheme/Elastic2dCurvilinearAnisotropicUpwind.m @ 1275:8ff3a95ad7cc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change metric computation in CurvilinearUpwind so that one can use odd orders for the diffOp
author Martin Almquist <malmquist@stanford.edu>
date Wed, 03 Jun 2020 10:29:53 -0700
parents 15865fbda16e
children
line wrap: on
line diff
--- a/+scheme/Elastic2dCurvilinearAnisotropicUpwind.m	Tue Jun 02 11:44:38 2020 -0700
+++ b/+scheme/Elastic2dCurvilinearAnisotropicUpwind.m	Wed Jun 03 10:29:53 2020 -0700
@@ -106,10 +106,11 @@
 
             % 1D operators
             opSetMetric = {@sbp.D2VariableCompatible, @sbp.D2VariableCompatible};
+            orderMetric = ceil(order/2)*2;
             m_u = m(1);
             m_v = m(2);
-            ops_u = opSetMetric{1}(m_u, {0, 1}, order);
-            ops_v = opSetMetric{2}(m_v, {0, 1}, order);
+            ops_u = opSetMetric{1}(m_u, {0, 1}, orderMetric);
+            ops_v = opSetMetric{2}(m_v, {0, 1}, orderMetric);
 
             h_u = ops_u.h;
             h_v = ops_v.h;