sbplib: +scheme/Utux2d.m annotate

annotate +scheme/Utux2d.m @ 989:e41c93d7ab08 feature/timesteppers

Merge with default

author	Jonatan Werpers <jonatan@werpers.com>
date	Wed, 09 Jan 2019 08:56:42 +0100
parents	cab047de7f5d
children	78db023a7fe3 f029b97dbc72 c12b84fe9b00

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
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	142 switch neighbour_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	143 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	144 e_neighbour = neighbour_scheme.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	145 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	146 e_neighbour = neighbour_scheme.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	147 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	148 e_neighbour = neighbour_scheme.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	149 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	150 e_neighbour = neighbour_scheme.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	151 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	152
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	153 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	154
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	155 % 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	156 case 'upwind'
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	157 sigma_ds = -1; %"Downstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	158 sigma_us = 0; %"Upstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	159
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	160 % 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	161 case 'centered'
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	162 sigma_ds = -1/2; %"Downstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	163 sigma_us = 1/2; %"Upstream" penalty
0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	164
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	165 otherwise
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	166 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	167 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	168
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 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	170 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	171 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	172 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	173 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	174 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	175 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	176 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	177 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	178 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	179 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	180 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	181 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	182 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	183 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	184 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	185 penalty = -obj.Hitaue_neighbour';
605 0f9d20dbb7ce Add centered interface coupling in addition to upwind Martin Almquist <malmquist@stanford.edu> parents: 595 diff changeset	186 end
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	187
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	188 end
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	189
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	190 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	191
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	192 % 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	193 default_field(type, 'interpOpSet', @sbp.InterpOpsOP);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	194
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	195 interpOpSet = type.interpOpSet;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	196 couplingType = type.couplingType;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	197 interpolationDamping = type.interpolationDamping;
610 b7b3c11fab4d Add interpolation to scheme Martin Almquist <malmquist@stanford.edu> parents: 605 diff changeset	198
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	199 % Get neighbour boundary operator
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 switch neighbour_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	201 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	202 e_neighbour = neighbour_scheme.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	203 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	204 e_neighbour = neighbour_scheme.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	205 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	206 e_neighbour = neighbour_scheme.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	207 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	208 e_neighbour = neighbour_scheme.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	209 end
459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	210
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	211 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	212
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 % 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	214 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	215 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	216 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	217
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	218 % 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	219 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	220 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	221 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	222
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	223 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	224 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	225 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	226
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	227 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	228 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	229
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	230 % 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	231 % 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	232 % 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	233 switch boundary
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	234 case {'w','e'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	235 m_u = obj.m(2);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	236 case {'s','n'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	237 m_u = obj.m(1);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	238 end
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	239 m_v = size(e_neighbour, 2);
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	240
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	241 % Build interpolation operators
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	242 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	243 Iu2v = intOps.Iu2v;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	244 Iv2u = intOps.Iv2u;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	245
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	246 I_local2neighbour_ds = intOps.Iu2v.bad;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	247 I_local2neighbour_us = intOps.Iu2v.good;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	248 I_neighbour2local_ds = intOps.Iv2u.good;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	249 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	250
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	251 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	252 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	253
2d85f17a8aec Add possibility for damping terms at interpolation interface. Martin Almquist <malmquist@stanford.edu> parents: 610 diff changeset	254
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	255 switch boundary
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	256 case {'w','W','west','West'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	257 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	258 closure = obj.Hitauobj.e_w';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	259 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	260
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	261 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	262 obj.e_wobj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	263 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	264 case {'e','E','east','East'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	265 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	266 closure = obj.Hitauobj.e_e';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	267 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	268
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	269 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	270 obj.e_eobj.H_y;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	271 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	272 case {'s','S','south','South'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	273 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	274 closure = obj.Hitauobj.e_s';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	275 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	276
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	277 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	278 obj.e_sobj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	279 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	280 case {'n','N','north','North'}
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	281 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	282 closure = obj.Hitauobj.e_n';
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	283 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	284
942 35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	285 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	286 obj.e_nobj.H_x;
35701c85e356 Make Utux2D work with new interface type etc. Martin Almquist <malmquist@stanford.edu> parents: 914 diff changeset	287 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	288 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	289
459eeb99130f Include type as (optional) input parameter in the interface method of all schemes. Martin Almquist <malmquist@stanford.edu> parents: 743 diff changeset	290
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	291 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	292
591 39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	293 function N = size(obj)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	294 N = obj.m;
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	295 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	296
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	297 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	298
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	299 methods(Static)
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	300 % 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	301 % and bound_v of scheme schm_v.
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	302 % [uu, uv, vv, vu] = inteface_coupling(A,'r',B,'l')
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	303 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	304 [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	305 [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	306 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	307 end
39554f2de783 Add Utux2D scheme Martin Almquist <martin.almquist@it.uu.se> parents: diff changeset	308 end

Mercurial > repos > public > sbplib

annotate +scheme/Utux2d.m @ 989:e41c93d7ab08 feature/timesteppers