Mercurial > repos > public > sbplib

comparison diracDiscrCurve.m @ 1246:25efceb0c392 feature/dirac_discr

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Apply some nitpicking
author Jonatan Werpers <jonatan@werpers.com>
date Wed, 20 Nov 2019 21:02:06 +0100
parents 745dc1f1a069
children 663eb90a4559
comparison
equal deleted inserted replaced
1245:0a1c64d3c717 1246:25efceb0c392
1 % 2-dimensional delta function for single-block curvilinear grid
2 % x_s: source point coordinate vector, e.g. [x; y] or [x; y; z].
3 % g: single-block grid containing the source
4 % m_order: Number of moment conditions
5 % s_order: Number of smoothness conditions
6 % order: Order of SBP derivative approximations
7 % opSet: Cell array of function handle to opSet generator
1 function d = diracDiscrCurve(x_s, g, m_order, s_order, order, opSet) 8 function d = diracDiscrCurve(x_s, g, m_order, s_order, order, opSet)
2 % 2-dimensional delta function for single-block curvilinear grid
3 % x_s: source point coordinate vector, e.g. [x; y] or [x; y; z].
4 % g: single-block grid containing the source
5 % m_order: Number of moment conditions
6 % s_order: Number of smoothness conditions
7 % order: Order of SBP derivative approximations
8 % opSet: Cell array of function handle to opSet generator
9
10 default_arg('order', m_order); 9 default_arg('order', m_order);
11 default_arg('opSet', {@sbp.D2Variable, @sbp.D2Variable}); 10 default_arg('opSet', {@sbp.D2Variable, @sbp.D2Variable});
12 11
13 dim = length(x_s); 12 dim = length(x_s);
14 assert(dim == 2, 'diracDiscrCurve: Only implemented for 2d.'); 13 assert(dim == 2, 'diracDiscrCurve: Only implemented for 2d.');
32 H_u = ops_u.H; 31 H_u = ops_u.H;
33 H_v = ops_v.H; 32 H_v = ops_v.H;
34 33
35 % Get delta function for logical grid and scale by Jacobian 34 % Get delta function for logical grid and scale by Jacobian
36 d = (1./J).*diracDiscr(g, [uS; vS], m_order, s_order, {H_u, H_v}); 35 d = (1./J).*diracDiscr(g, [uS; vS], m_order, s_order, {H_u, H_v});
37
38 end 36 end
39 37
40 function J = jacobian(g, opSet, order) 38 function J = jacobian(g, opSet, order)
41 m = g.size(); 39 m = g.size();
42 m_u = m(1); 40 m_u = m(1);