Mercurial > repos > public > sbplib

changeset 449:0707a7192bc3 feature/grids

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add first derivative operators to Wave2DCurve
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 09 May 2017 14:51:29 +0200
parents 692bf61385c0
children 8d455e49364f
files +scheme/Wave2dCurve.m
diffstat 1 files changed, 5 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
diff -r 692bf61385c0 -r 0707a7192bc3 +scheme/Wave2dCurve.m
--- a/+scheme/Wave2dCurve.m	Tue May 09 14:47:52 2017 +0200
+++ b/+scheme/Wave2dCurve.m	Tue May 09 14:51:29 2017 +0200
@@ -27,6 +27,8 @@
         gamm_u, gamm_v
         lambda
 
+        Dx, Dy % Physical derivatives
+
         x_u
         x_v
         y_u
@@ -153,6 +155,9 @@
             obj.D = obj.Ji*c^2*(Duu + Duv + Dvu + Dvv);
             obj.lambda = lambda;
 
+            obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
+            obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
+
             obj.gamm_u = h_u*ops_u.borrowing.M.d1;
             obj.gamm_v = h_v*ops_v.borrowing.M.d1;
         end