Mercurial > repos > public > sbplib

changeset 111:0e66299592cc

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed some bugs regarding flowState determination.
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 09 Dec 2015 15:48:58 +0100
parents f5ed7ff58115
children b7f40a0c6476
files +scheme/Euler1d.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Euler1d.m	Tue Dec 08 13:03:04 2015 +0100
+++ b/+scheme/Euler1d.m	Wed Dec 09 15:48:58 2015 +0100
@@ -165,7 +165,7 @@
 
         function fs = flowStateL(obj, q)
             q_l = obj.e_L'*q;
-            c = obj.c(q);
+            c = obj.c(q_l);
             v = q_l(2,:)/q_l(1,:);
 
             if v > c
@@ -175,7 +175,7 @@
             elseif v > -c
                 fs = scheme.Euler1d.SUBSONIC_OUTFLOW;
             else
-                fs = scheme.Euler1d.SUPERSONIC_INFLOW;
+                fs = scheme.Euler1d.SUPERSONIC_OUTFLOW;
             end
         end
 
@@ -183,7 +183,7 @@
         %  +-1 for subsonic
         function fs = flowStateR(obj, q)
             q_r = obj.e_R'*q;
-            c = obj.c(q);
+            c = obj.c(q_r);
             v = q_r(2,:)/q_r(1,:);
 
             if v < -c
@@ -193,7 +193,7 @@
             elseif v < c
                 fs = scheme.Euler1d.SUBSONIC_OUTFLOW;
             else
-                fs = scheme.Euler1d.SUPERSONIC_INFLOW;
+                fs = scheme.Euler1d.SUPERSONIC_OUTFLOW;
             end
         end