Mercurial > repos > public > sbplib
changeset 946:706d1c2b4199 feature/utux2D
Raname opts to type in a bunch of interface methods
| author | Jonatan Werpers <jonatan@werpers.com> |
|---|---|
| date | Wed, 05 Dec 2018 15:57:41 +0100 |
| parents | af397cc35239 |
| children | d1c13505182a |
| files | +scheme/Hypsyst2d.m +scheme/Hypsyst2dCurve.m +scheme/Hypsyst3d.m +scheme/Hypsyst3dCurve.m +scheme/Laplace1D.m +scheme/LaplaceCurvilinear.m +scheme/Schrodinger.m +scheme/Utux.m +scheme/Wave2d.m +scheme/Wave2dCurve.m |
| diffstat | 10 files changed, 14 insertions(+), 18 deletions(-) [+] |
line wrap: on
line diff
--- a/+scheme/Hypsyst2d.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Hypsyst2d.m Wed Dec 05 15:57:41 2018 +0100 @@ -93,8 +93,8 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) - error('An interface function does not exist yet'); + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) + error('Not implemented'); end function N = size(obj)
--- a/+scheme/Hypsyst2dCurve.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Hypsyst2dCurve.m Wed Dec 05 15:57:41 2018 +0100 @@ -128,8 +128,8 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) - error('An interface function does not exist yet'); + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) + error('Not implemented'); end function N = size(obj)
--- a/+scheme/Hypsyst3d.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Hypsyst3d.m Wed Dec 05 15:57:41 2018 +0100 @@ -168,8 +168,8 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) - error('An interface function does not exist yet'); + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) + error('Not implemented'); end function N = size(obj)
--- a/+scheme/Hypsyst3dCurve.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Hypsyst3dCurve.m Wed Dec 05 15:57:41 2018 +0100 @@ -264,8 +264,8 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) - error('An interface function does not exist yet'); + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) + error('Not implemented'); end function N = size(obj)
--- a/+scheme/Laplace1D.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Laplace1D.m Wed Dec 05 15:57:41 2018 +0100 @@ -83,7 +83,7 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) % u denotes the solution in the own domain % v denotes the solution in the neighbour domain @@ -98,7 +98,7 @@ gamm_v = neighbour_scheme.gamm; tuning = 1.1; - + tau1 = -(a_u/gamm_u + a_v/gamm_v) * tuning; tau2 = 1/2*a_u; sig1 = -1/2;
--- a/+scheme/LaplaceCurvilinear.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/LaplaceCurvilinear.m Wed Dec 05 15:57:41 2018 +0100 @@ -383,9 +383,6 @@ % % I -- the indices of the boundary points in the grid matrix function [e, d, gamm, H_b, I] = get_boundary_ops(obj, boundary) - - % gridMatrix = zeros(obj.m(2),obj.m(1)); - % gridMatrix(:) = 1:numel(gridMatrix); ind = grid.funcToMatrix(obj.grid, 1:prod(obj.m)); switch boundary @@ -419,7 +416,6 @@ case {'s','n'} gamm = obj.gamm_v; end - end function N = size(obj)
--- a/+scheme/Schrodinger.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Schrodinger.m Wed Dec 05 15:57:41 2018 +0100 @@ -90,7 +90,7 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) % u denotes the solution in the own domain % v denotes the solution in the neighbour domain [e_u,d_u,s_u] = obj.get_boundary_ops(boundary);
--- a/+scheme/Utux.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Utux.m Wed Dec 05 15:57:41 2018 +0100 @@ -68,7 +68,7 @@ end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) switch boundary % Upwind coupling case {'l','left'}
--- a/+scheme/Wave2d.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Wave2d.m Wed Dec 05 15:57:41 2018 +0100 @@ -158,7 +158,7 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) % u denotes the solution in the own domain % v denotes the solution in the neighbour domain [e_u,d_u,s_u,gamm_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
--- a/+scheme/Wave2dCurve.m Wed Dec 05 10:55:50 2018 +0100 +++ b/+scheme/Wave2dCurve.m Wed Dec 05 15:57:41 2018 +0100 @@ -243,7 +243,7 @@ end end - function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,opts) + function [closure, penalty] = interface(obj, boundary, neighbour_scheme, neighbour_boundary, type) % u denotes the solution in the own domain % v denotes the solution in the neighbour domain tuning = 1.2;