Mercurial > repos > public > sbplib
changeset 539:08b6281ba2a9 feature/grids
Add some todos
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Tue, 08 Aug 2017 13:25:19 +0200 |
| parents | 95e41f7a6f1a |
| children | ef8003bdbc0d |
| files | +multiblock/Grid.m +scheme/LaplaceCurvilinear.m +scheme/Scheme.m |
| diffstat | 3 files changed, 7 insertions(+), 1 deletions(-) [+] |
line wrap: on
line diff
--- a/+multiblock/Grid.m Tue Aug 08 08:26:49 2017 +0200 +++ b/+multiblock/Grid.m Tue Aug 08 13:25:19 2017 +0200 @@ -82,6 +82,9 @@ gfs = mat2cell(gf, N, 1); end + % TODO: Split op? + % Should the method to split an operator be moved here instead of being in multiblock.DiffOp? + % Converts a gridfunction to a set of plot matrices % Takes a grid function and and a structured grid. function F = funcToPlotMatrices(obj, gf)
--- a/+scheme/LaplaceCurvilinear.m Tue Aug 08 08:26:49 2017 +0200 +++ b/+scheme/LaplaceCurvilinear.m Tue Aug 08 13:25:19 2017 +0200 @@ -37,6 +37,8 @@ methods % Implements a*div(b*grad(u)) as a SBP scheme + % TODO: Implement proper H, it should be the real physical quadrature, the logic quadrature may be but in a separate variable (H_logic?) + function obj = LaplaceCurvilinear(g ,order, a, b, opSet) default_arg('opSet',@sbp.D2Variable); default_arg('a', 1);
--- a/+scheme/Scheme.m Tue Aug 08 08:26:49 2017 +0200 +++ b/+scheme/Scheme.m Tue Aug 08 13:25:19 2017 +0200 @@ -28,7 +28,8 @@ [closure, penalty] = boundary_condition(obj,boundary,type) [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary) - % op = getBoundaryOperator()?? + % TODO: op = getBoundaryOperator()?? + % makes sense to have it available through a method instead of random properties % Returns the number of degrees of freedom. N = size(obj)