Mercurial > repos > public > sbplib

changeset 650:8e55298657b9 feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add characteristic BC
author Martin Almquist <malmquist@stanford.edu>
date Wed, 15 Nov 2017 14:56:52 -0800
parents 1bdbe026abbc
children 4ee7d15bd8e6
files +scheme/Staggered1DAcoustics.m
diffstat 1 files changed, 11 insertions(+), 8 deletions(-) [+]
line wrap: on
line diff
diff -r 1bdbe026abbc -r 8e55298657b9 +scheme/Staggered1DAcoustics.m
--- a/+scheme/Staggered1DAcoustics.m	Tue Nov 14 15:40:06 2017 -0800
+++ b/+scheme/Staggered1DAcoustics.m	Wed Nov 15 14:56:52 2017 -0800
@@ -116,14 +116,6 @@
             % No other types implemented yet
 
             % BC on the form Lu - g = 0;
-            switch type
-            case 'p'
-                L = [1, 0];
-            case 'v'
-                L = [0, 1];
-            otherwise
-                error('Boundary condition not implemented.');
-            end
 
             % Diagonalize B
             B = obj.B;
@@ -146,6 +138,17 @@
             Tin = T(:,Iin);
             Tout = T(:,Iout);
 
+            switch type
+            case 'p'
+                L = [1, 0];
+            case 'v'
+                L = [0, 1];
+            case 'characteristic'
+                L = Tin';
+            otherwise
+                error('Boundary condition not implemented.');
+            end
+
             % Penalty parameters
             A = obj.A;
             sigma = [0; 0];