Mercurial > repos > public > sbplib

diff +scheme/Wave2dCurve.m @ 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 359861563866
children 8d455e49364f
line wrap: on
line diff
--- 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