sbplib: +scheme/Utux2d.m annotate

annotate +scheme/Utux2d.m @ 1001:514a98f9f90d feature/getBoundaryOp

Add getBoundaryOperator and getBoundaryQuadrature as abstract methods in scheme.Scheme

author	Jonatan Werpers <jonatan@werpers.com>
date	Wed, 16 Jan 2019 16:35:01 +0100
parents	78db023a7fe3
children	8d73fcdb07a5

rev	line source
950 cab047de7f5d Rename 2D schemes to 2d Jonatan Werpers <jonatan@werpers.com> parents: 948 diff changeset	1 classdef Utux2d < scheme.Scheme
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	2 properties
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	4 h % Grid spacing
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	5 grid % Grid
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	6 order % Order accuracy for the approximation
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	7 v0 % Initial data
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	8
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	9 a % Wave speed a = [a1, a2];
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	10 % Can either be a constant vector or a cell array of function handles.
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	11
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	12 H % Discrete norm
595 2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	13 H_x, H_y % Norms in the x and y directions
2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	14 Hi, Hx, Hy, Hxi, Hyi % Kroneckered norms
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	15
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	16 % Derivatives
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	17 Dx, Dy
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	18
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	19 % Boundary operators
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	20 e_w, e_e, e_s, e_n
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	21
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	22 D % Total discrete operator
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	23 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	24
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	25
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	26 methods
950 cab047de7f5d Rename 2D schemes to 2d Jonatan Werpers <jonatan@werpers.com> parents: 948 diff changeset	27 function obj = Utux2d(g ,order, opSet, a)
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	28
459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	29 default_arg('a',1/sqrt(2)*[1, 1]);
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	30 default_arg('opSet',@sbp.D2Standard);
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	31
948 3dd7f87c9a1b Use assertType(...) instead of assert(isa(...)) Jonatan Werpers <jonatan@werpers.com> parents: 942 diff changeset	32 assertType(g, 'grid.Cartesian');
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	33 if iscell(a)
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	34 a1 = grid.evalOn(g, a{1});
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	35 a2 = grid.evalOn(g, a{2});
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	36 a = {spdiag(a1), spdiag(a2)};
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	37 else
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	38 a = {a(1), a(2)};
f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	39 end
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	40
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	41 m = g.size();
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	42 m_x = m(1);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	43 m_y = m(2);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	44 m_tot = g.N();
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	45
595 2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	46 xlim = {g.x{1}(1), g.x{1}(end)};
2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	47 ylim = {g.x{2}(1), g.x{2}(end)};
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	48 obj.grid = g;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	49
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	50 % Operator sets
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	51 ops_x = opSet(m_x, xlim, order);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	52 ops_y = opSet(m_y, ylim, order);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	53 Ix = speye(m_x);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	54 Iy = speye(m_y);
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	55
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	56 % Norms
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	57 Hx = ops_x.H;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	58 Hy = ops_y.H;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	59 Hxi = ops_x.HI;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	60 Hyi = ops_y.HI;
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	61
595 2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	62 obj.H_x = Hx;
2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	63 obj.H_y = Hy;
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	64 obj.H = kron(Hx,Hy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	65 obj.Hi = kron(Hxi,Hyi);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	66 obj.Hx = kron(Hx,Iy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	67 obj.Hy = kron(Ix,Hy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	68 obj.Hxi = kron(Hxi,Iy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	69 obj.Hyi = kron(Ix,Hyi);
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	70
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	71 % Derivatives
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	72 Dx = ops_x.D1;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	73 Dy = ops_y.D1;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	74 obj.Dx = kron(Dx,Iy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	75 obj.Dy = kron(Ix,Dy);
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	76
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	77 % Boundary operators
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	78 obj.e_w = kr(ops_x.e_l, Iy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	79 obj.e_e = kr(ops_x.e_r, Iy);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	80 obj.e_s = kr(Ix, ops_y.e_l);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	81 obj.e_n = kr(Ix, ops_y.e_r);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	82
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	83 obj.m = m;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	84 obj.h = [ops_x.h ops_y.h];
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	85 obj.order = order;
595 2a2f34778ded Make Utux2D work Martin Almquist <malmquist@stanford.edu> parents: 591 diff changeset	86 obj.a = a;
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	87 obj.D = -(a{1}obj.Dx + a{2}obj.Dy);
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	88
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	89 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	90 % Closure functions return the opertors applied to the own domain to close the boundary
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	91 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	92 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	93 % type is a string specifying the type of boundary condition if there are several.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	94 % data is a function returning the data that should be applied at the boundary.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	95 % neighbour_scheme is an instance of Scheme that should be interfaced to.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	96 % neighbour_boundary is a string specifying which boundary to interface to.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	97 function [closure, penalty] = boundary_condition(obj,boundary,type)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	98 default_arg('type','dirichlet');
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	99
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	100 sigma = -1; % Scalar penalty parameter
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	101 switch boundary
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	102 case {'w','W','west','West'}
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	103 tau = sigmaobj.a{1}obj.e_w*obj.H_y;
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	104 closure = obj.Hitauobj.e_w';
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	105
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	106 case {'s','S','south','South'}
743 f4595f14d696 Change schemes to work for special coefficients. Martin Almquist <malmquist@stanford.edu> parents: 666 diff changeset	107 tau = sigmaobj.a{2}obj.e_s*obj.H_x;
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	108 closure = obj.Hitauobj.e_s';
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	109 end
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	110 penalty = -obj.Hi*tau;
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	111
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	112 end
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	113
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	114 % type Struct that specifies the interface coupling.
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	115 % Fields:
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	116 % -- couplingType String, type of interface coupling
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	117 % % Default: 'upwind'. Other: 'centered'
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	118 % -- interpolation: type of interpolation, default 'none'
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	119 % -- interpolationDamping: damping on upstream and downstream sides, when using interpolation.
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	120 % Default {0,0} gives zero damping.
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	121 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary,type)
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	122
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	123 defaultType.couplingType = 'upwind';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	124 defaultType.interpolation = 'none';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	125 defaultType.interpolationDamping = {0,0};
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	126 default_struct('type', defaultType);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	127
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	128 switch type.interpolation
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	129 case {'none', ''}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	130 [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	131 case {'op','OP'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	132 [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	133 otherwise
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	134 error('Unknown type of interpolation: %s ', type.interpolation);
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	135 end
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	136 end
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	137
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	138 function [closure, penalty] = interfaceStandard(obj,boundary,neighbour_scheme,neighbour_boundary,type)
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	139 couplingType = type.couplingType;
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	140
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	141 % Get neighbour boundary operator
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	142 e_neighbour = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	143
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	144 switch couplingType
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	145
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	146 % Upwind coupling (energy dissipation)
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	147 case 'upwind'
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	148 sigma_ds = -1; %"Downstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	149 sigma_us = 0; %"Upstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	150
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	151 % Energy-preserving coupling (no energy dissipation)
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	152 case 'centered'
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	153 sigma_ds = -1/2; %"Downstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	154 sigma_us = 1/2; %"Upstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	155
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	156 otherwise
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	157 error(['Interface coupling type ' couplingType ' is not available.'])
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	158 end
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	159
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	160 switch boundary
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	161 case {'w','W','west','West'}
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	162 tau = sigma_dsobj.a{1}obj.e_w*obj.H_y;
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	163 closure = obj.Hitauobj.e_w';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	164 penalty = -obj.Hitaue_neighbour';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	165 case {'e','E','east','East'}
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	166 tau = sigma_usobj.a{1}obj.e_e*obj.H_y;
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	167 closure = obj.Hitauobj.e_e';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	168 penalty = -obj.Hitaue_neighbour';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	169 case {'s','S','south','South'}
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	170 tau = sigma_dsobj.a{2}obj.e_s*obj.H_x;
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	171 closure = obj.Hitauobj.e_s';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	172 penalty = -obj.Hitaue_neighbour';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	173 case {'n','N','north','North'}
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	174 tau = sigma_usobj.a{2}obj.e_n*obj.H_x;
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	175 closure = obj.Hitauobj.e_n';
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	176 penalty = -obj.Hitaue_neighbour';
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	177 end
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	178
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	179 end
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	180
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	181 function [closure, penalty] = interfaceNonConforming(obj,boundary,neighbour_scheme,neighbour_boundary,type)
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	182
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	183 % User can request special interpolation operators by specifying type.interpOpSet
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	184 default_field(type, 'interpOpSet', @sbp.InterpOpsOP);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	185
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	186 interpOpSet = type.interpOpSet;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	187 couplingType = type.couplingType;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	188 interpolationDamping = type.interpolationDamping;
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	189
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	190 % Get neighbour boundary operator
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	191 e_neighbour = neighbour_scheme.getBoundaryOperator('e', neighbour_boundary);
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	192
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	193 switch couplingType
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	194
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	195 % Upwind coupling (energy dissipation)
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	196 case 'upwind'
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	197 sigma_ds = -1; %"Downstream" penalty
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	198 sigma_us = 0; %"Upstream" penalty
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	199
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	200 % Energy-preserving coupling (no energy dissipation)
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	201 case 'centered'
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	202 sigma_ds = -1/2; %"Downstream" penalty
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	203 sigma_us = 1/2; %"Upstream" penalty
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	204
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	205 otherwise
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	206 error(['Interface coupling type ' couplingType ' is not available.'])
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	207 end
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	208
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	209 int_damp_us = interpolationDamping{1};
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	210 int_damp_ds = interpolationDamping{2};
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	211
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	212 % u denotes the solution in the own domain
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	213 % v denotes the solution in the neighbour domain
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	214 % Find the number of grid points along the interface
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	215 switch boundary
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	216 case {'w','e'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	217 m_u = obj.m(2);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	218 case {'s','n'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	219 m_u = obj.m(1);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	220 end
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	221 m_v = size(e_neighbour, 2);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	222
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	223 % Build interpolation operators
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	224 intOps = interpOpSet(m_u, m_v, obj.order, neighbour_scheme.order);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	225 Iu2v = intOps.Iu2v;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	226 Iv2u = intOps.Iv2u;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	227
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	228 I_local2neighbour_ds = intOps.Iu2v.bad;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	229 I_local2neighbour_us = intOps.Iu2v.good;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	230 I_neighbour2local_ds = intOps.Iv2u.good;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	231 I_neighbour2local_us = intOps.Iv2u.bad;
914 50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	232
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	233 I_back_forth_us = I_neighbour2local_us*I_local2neighbour_us;
50cafc4b9e40 Make default_field overwrite if field exists but is empty. Refactor interface method in Utux2d. Martin Almquist <malmquist@stanford.edu> parents: 910 diff changeset	234 I_back_forth_ds = I_neighbour2local_ds*I_local2neighbour_ds;
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	235
2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	236
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	237 switch boundary
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	238 case {'w','W','west','West'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	239 tau = sigma_dsobj.a{1}obj.e_w*obj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	240 closure = obj.Hitauobj.e_w';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	241 penalty = -obj.HitauI_neighbour2local_ds*e_neighbour';
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	242
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	243 beta = int_damp_ds*obj.a{1}...
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	244 obj.e_wobj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	245 closure = closure + obj.HibetaI_back_forth_dsobj.e_w' - obj.Hibeta*obj.e_w';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	246 case {'e','E','east','East'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	247 tau = sigma_usobj.a{1}obj.e_e*obj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	248 closure = obj.Hitauobj.e_e';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	249 penalty = -obj.HitauI_neighbour2local_us*e_neighbour';
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	250
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	251 beta = int_damp_us*obj.a{1}...
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	252 obj.e_eobj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	253 closure = closure + obj.HibetaI_back_forth_usobj.e_e' - obj.Hibeta*obj.e_e';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	254 case {'s','S','south','South'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	255 tau = sigma_dsobj.a{2}obj.e_s*obj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	256 closure = obj.Hitauobj.e_s';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	257 penalty = -obj.HitauI_neighbour2local_ds*e_neighbour';
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	258
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	259 beta = int_damp_ds*obj.a{2}...
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	260 obj.e_sobj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	261 closure = closure + obj.HibetaI_back_forth_dsobj.e_s' - obj.Hibeta*obj.e_s';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	262 case {'n','N','north','North'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	263 tau = sigma_usobj.a{2}obj.e_n*obj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	264 closure = obj.Hitauobj.e_n';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	265 penalty = -obj.HitauI_neighbour2local_us*e_neighbour';
666 2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	266
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	267 beta = int_damp_us*obj.a{2}...
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	268 obj.e_nobj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	269 closure = closure + obj.HibetaI_back_forth_usobj.e_n' - obj.Hibeta*obj.e_n';
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	270 end
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	271
459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	272
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	273 end
905 459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	274
997 78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	275 % Returns the boundary operator op for the boundary specified by the string boundary.
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	276 % op -- string or a cell array of strings
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	277 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	278 function varargout = getBoundaryOperator(obj, op, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	279
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	280 if ~iscell(op)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	281 op = {op};
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	282 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	283
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	284 for i = 1:numel(op)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	285 switch op{i}
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	286 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	287 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	288 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	289 e = obj.e_w;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	290 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	291 e = obj.e_e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	292 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	293 e = obj.e_s;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	294 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	295 e = obj.e_n;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	296 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	297 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	298 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	299 varargout{i} = e;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	300 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	301 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	302
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	303 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	304
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	305 % Returns square boundary quadrature matrix, of dimension
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	306 % corresponding to the number of boundary points
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	307 %
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	308 % boundary -- string
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	309 function H_b = getBoundaryQuadrature(obj, boundary)
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	310
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	311 switch boundary
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	312 case 'w'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	313 H_b = obj.H_y;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	314 case 'e'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	315 H_b = obj.H_y;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	316 case 's'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	317 H_b = obj.H_x;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	318 case 'n'
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	319 H_b = obj.H_x;
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	320 otherwise
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	321 error('No such boundary: boundary = %s',boundary);
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	322 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	323 end
78db023a7fe3 Add getBoundaryOperator method to all 2d schemes, except obsolete Wave2dCurve and ElastiCurve, which needs a big makeover. Martin Almquist <malmquist@stanford.edu> parents: 950 diff changeset	324
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	325 function N = size(obj)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	326 N = obj.m;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	327 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	328
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	329 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	330
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	331 methods(Static)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	332 % Calculates the matrices needed for the inteface coupling between boundary bound_u of scheme schm_u
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	333 % and bound_v of scheme schm_v.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	334 % [uu, uv, vv, vu] = inteface_coupling(A,'r',B,'l')
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	335 function [uu, uv, vv, vu] = interface_coupling(schm_u,bound_u,schm_v,bound_v)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	336 [uu,uv] = schm_u.interface(bound_u,schm_v,bound_v);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	337 [vv,vu] = schm_v.interface(bound_v,schm_u,bound_u);
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	338 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	339 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	340 end

Mercurial > repos > public > sbplib

annotate +scheme/Utux2d.m @ 1001:514a98f9f90d feature/getBoundaryOp