Mercurial > repos > public > sbplib

changeset 58:24103284e09d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fixed coding errors in Euler1d.
author Jonatan Werpers <jonatan@werpers.com>
date Fri, 13 Nov 2015 16:39:45 -0800
parents 9a647dcccbdd
children e431c1260f52
files +scheme/Euler1d.m
diffstat 1 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Euler1d.m	Thu Nov 12 17:31:33 2015 -0800
+++ b/+scheme/Euler1d.m	Fri Nov 13 16:39:45 2015 -0800
@@ -79,13 +79,13 @@
             %Solving on form q_t + F_x = 0
             function o = F(Q)
                 % Flux: f = [q2; q2.^2/q1 + p(q); (q3+p(q))*q2/q1];
-                o = [Q(2,:); Q(2,:).^2/Q(1,:) + p(Q); (Q(3,:)+p(Q)).*Q(2,:)./Q(1,:)];
+                o = [Q(2,:); Q(2,:).^2./Q(1,:) + p(Q); (Q(3,:)+p(Q)).*Q(2,:)./Q(1,:)];
             end
 
             % Equation of state
             function o = p(Q)
                 % Pressure p = (gamma-1)*(q3-q2.^2/q1/2)
-                o = (gamma-1)*(Q(3,:)-Q(2,:).^2/Q(1,:)/2);
+                o = (gamma-1)*( Q(3,:)-1/2*Q(2,:).^2./Q(1,:) );
             end
 
             function o = c(Q)