Mercurial > repos > public > sbplib

comparison +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
comparison
equal deleted inserted replaced
448:692bf61385c0 449:0707a7192bc3
25 du_s, dv_s 25 du_s, dv_s
26 du_n, dv_n 26 du_n, dv_n
27 gamm_u, gamm_v 27 gamm_u, gamm_v
28 lambda 28 lambda
29 29
30 Dx, Dy % Physical derivatives
31
30 x_u 32 x_u
31 x_v 33 x_v
32 y_u 34 y_u
33 y_v 35 y_v
34 end 36 end
151 obj.a12 = a12; 153 obj.a12 = a12;
152 obj.a22 = a22; 154 obj.a22 = a22;
153 obj.D = obj.Ji*c^2*(Duu + Duv + Dvu + Dvv); 155 obj.D = obj.Ji*c^2*(Duu + Duv + Dvu + Dvv);
154 obj.lambda = lambda; 156 obj.lambda = lambda;
155 157
158 obj.Dx = spdiag( y_v./J)*Du + spdiag(-y_u./J)*Dv;
159 obj.Dy = spdiag(-x_v./J)*Du + spdiag( x_u./J)*Dv;
160
156 obj.gamm_u = h_u*ops_u.borrowing.M.d1; 161 obj.gamm_u = h_u*ops_u.borrowing.M.d1;
157 obj.gamm_v = h_v*ops_v.borrowing.M.d1; 162 obj.gamm_v = h_v*ops_v.borrowing.M.d1;
158 end 163 end
159 164
160 165