Mercurial > repos > public > sbplib
comparison +scheme/Laplace1d.m @ 1049:0c504a21432d feature/getBoundaryOp
Add getBoundaryQuadrature to all 1d diffOps
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Tue, 22 Jan 2019 11:17:28 -0800 |
parents | 6bc55a773e7c |
children | ae4b090b5299 |
comparison
equal
deleted
inserted
replaced
1048:adbb80e60b10 | 1049:0c504a21432d |
---|---|
124 assertIsMember(boundary, {'l', 'r'}) | 124 assertIsMember(boundary, {'l', 'r'}) |
125 | 125 |
126 o = obj.([op, '_', boundary]) | 126 o = obj.([op, '_', boundary]) |
127 end | 127 end |
128 | 128 |
129 % Returns square boundary quadrature matrix, of dimension | |
130 % corresponding to the number of boundary points | |
131 % | |
132 % boundary -- string | |
133 % Note: for 1d diffOps, the boundary quadrature is the scalar 1. | |
134 function H_b = getBoundaryQuadrature(obj, boundary) | |
135 assertIsMember(boundary, {'l', 'r'}) | |
136 | |
137 H_b = 1; | |
138 end | |
139 | |
129 % Returns the boundary sign. The right boundary is considered the positive boundary | 140 % Returns the boundary sign. The right boundary is considered the positive boundary |
130 % boundary -- string | 141 % boundary -- string |
131 function s = getBoundarySign(obj, boundary) | 142 function s = getBoundarySign(obj, boundary) |
132 assertIsMember(boundary, {'l', 'r'}) | 143 assertIsMember(boundary, {'l', 'r'}) |
133 | 144 |