Mercurial > repos > public > sbplib

comparison +scheme/LaplaceCurvilinear.m @ 552:ffaf13533c27 feature/grids/laplace_refactor

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Reorganize properties of the class
author Jonatan Werpers <jonatan@werpers.com>
date Tue, 29 Aug 2017 10:52:21 +0200
parents c5a7a13c03dc
children 63d5c07943dd
comparison
equal deleted inserted replaced
551:c5a7a13c03dc 552:ffaf13533c27
5 5
6 grid 6 grid
7 7
8 order % Order accuracy for the approximation 8 order % Order accuracy for the approximation
9 9
10 D % non-stabalized scheme operator 10 a,b % Parameters of the operator
11 M % Derivative norm 11
12 a,b 12
13 % Inner products and operators for physical coordinates
14 D % Laplace operator
15 H, Hi % Inner product
16 e_w, e_e, e_s, e_n
17 d_w, d_e, d_s, d_n % Normal derivatives at the boundary
18 H_w, H_e, H_s, H_n % Boundary inner products
19 Dx, Dy % Physical derivatives
20 M % Gradient inner product
21
22 % Metric coefficients
13 J, Ji 23 J, Ji
14 a11, a12, a22 24 a11, a12, a22
15 25 x_u
16 H % Discrete norm 26 x_v
17 Hi 27 y_u
28 y_v
29
30 % Inner product and operators for logical coordinates
18 H_u, H_v % Norms in the x and y directions 31 H_u, H_v % Norms in the x and y directions
32 Hi_u, Hi_v
19 Hu,Hv % Kroneckerd norms. 1'*Hx*v corresponds to integration in the x dir. 33 Hu,Hv % Kroneckerd norms. 1'*Hx*v corresponds to integration in the x dir.
20 Hi_u, Hi_v
21 Hiu, Hiv 34 Hiu, Hiv
22 e_w, e_e, e_s, e_n
23 du_w, dv_w 35 du_w, dv_w
24 du_e, dv_e 36 du_e, dv_e
25 du_s, dv_s 37 du_s, dv_s
26 du_n, dv_n 38 du_n, dv_n
27 gamm_u, gamm_v 39 gamm_u, gamm_v
28 lambda 40 lambda
29
30 Dx, Dy % Physical derivatives
31
32 x_u
33 x_v
34 y_u
35 y_v
36 end 41 end
37 42
38 methods 43 methods
39 % Implements a*div(b*grad(u)) as a SBP scheme 44 % Implements a*div(b*grad(u)) as a SBP scheme
40 % TODO: Implement proper H, it should be the real physical quadrature, the logic quadrature may be but in a separate variable (H_logic?) 45 % TODO: Implement proper H, it should be the real physical quadrature, the logic quadrature may be but in a separate variable (H_logic?)