Mercurial > repos > public > sbplib

changeset 1271:b5025bd67be1 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Allow using D1 D1 in VariableAnisotropicUpwind
author Martin Almquist <malmquist@stanford.edu>
date Sun, 31 May 2020 21:58:40 -0700
parents 9a7cac4202d0
children 15865fbda16e
files +scheme/Elastic2dVariableAnisotropicUpwind.m
diffstat 1 files changed, 11 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
diff -r 9a7cac4202d0 -r b5025bd67be1 +scheme/Elastic2dVariableAnisotropicUpwind.m
--- a/+scheme/Elastic2dVariableAnisotropicUpwind.m	Sun May 31 20:18:45 2020 -0700
+++ b/+scheme/Elastic2dVariableAnisotropicUpwind.m	Sun May 31 21:58:40 2020 -0700
@@ -128,8 +128,17 @@
 
             for i = 1:dim
                 I{i} = speye(m(i));
-                Dp{i} = ops{i}.Dp;
-                Dm{i} = ops{i}.Dm;
+
+                if isprop(ops{i}, 'Dp') && isprop(ops{i}, 'Dm')
+                    Dp{i} = ops{i}.Dp;
+                    Dm{i} = ops{i}.Dm;
+                elseif isprop(ops{i}, 'D1')
+                    Dp{i} = ops{i}.D1;
+                    Dm{i} = ops{i}.D1;
+                else
+                    error('opSet does not have Dp and Dm or D1');
+                end
+
                 H{i} =  ops{i}.H;
                 Hi{i} = ops{i}.HI;
                 e_0{i} = ops{i}.e_l;