Mercurial > repos > public > sbplib

changeset 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 8aa0909125a4
children bc2be302df90
files +scheme/Elastic2dCurvilinearAnisotropicUpwind.m
diffstat 1 files changed, 3 insertions(+), 2 deletions(-) [+]
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;