Mercurial > repos > public > sbplib

comparison +sbp/D1Staggered.m @ 638:17dffd56f2fe feature/d1_staggered

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fix and improved comments in D1Staggered wrapper class.
author Martin Almquist <malmquist@stanford.edu>
date Fri, 10 Nov 2017 15:57:21 -0800
parents eead18a8964d
children
comparison
equal deleted inserted replaced
637:eead18a8964d 638:17dffd56f2fe
1 classdef D1Staggered < sbp.OpSet 1 classdef D1Staggered < sbp.OpSet
2 properties 2 properties
3 % x_primal: "primal" grid with m points. Equidistant. Called Plus grid in Ossian's paper. 3 % x_primal: "primal" grid with m points. Equidistant. Called Plus grid in Ossian's paper.
4 % x_dual: "dual" grid with m+1 points. Called Minus grid in Ossian's paper. 4 % x_dual: "dual" grid with m+1 points. Called Minus grid in Ossian's paper.
5
6 % D1_primal takes FROM dual grid TO primal grid
7 % D1_dual takes FROM primal grid TO dual grid
5 8
6 D1_primal % SBP operator approximating first derivative 9 D1_primal % SBP operator approximating first derivative
7 D1_dual % SBP operator approximating first derivative 10 D1_dual % SBP operator approximating first derivative
8 H_primal % Norm matrix 11 H_primal % Norm matrix
9 H_dual % Norm matrix 12 H_dual % Norm matrix
53 D1_primal = Pp\Qp; 56 D1_primal = Pp\Qp;
54 D1_dual = Pm\Qm; 57 D1_dual = Pm\Qm;
55 58
56 obj.D1_primal = D1_primal; 59 obj.D1_primal = D1_primal;
57 obj.D1_dual = D1_dual; 60 obj.D1_dual = D1_dual;
58 pbj.H_primal = Pp; 61 obj.H_primal = Pp;
59 obj.H_dual = Pm; 62 obj.H_dual = Pm;
60 63
61 obj.e_primal_l = sparse(m_primal,1); 64 obj.e_primal_l = sparse(m_primal,1);
62 obj.e_primal_r = sparse(m_primal,1); 65 obj.e_primal_r = sparse(m_primal,1);
63 obj.e_primal_l(1) = 1; 66 obj.e_primal_l(1) = 1;