sbplib: +scheme/Euler1d.m annotate

annotate +scheme/Euler1d.m @ 49:8f0c2dc747dd

Made lots of updates to Euler1d.

author	Jonatan Werpers <jonatan@werpers.com>
date	Thu, 05 Nov 2015 16:33:53 -0800
parents	48b6fb693025
children	75ebf5d3cfe5

rev	line source
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	1 classdef Euler1d < scheme.Scheme
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	2 properties
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	4 N % Number of points total
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	5 h % Grid spacing
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	6 u % Grid values
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	7 x % Values of x and y for each
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	8 order % Order accuracy for the approximation
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	9
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	10 D % non-stabalized scheme operator
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	11 Fx
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	12 M % Derivative norm
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	13 gamma
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	14
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	15 T
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	16 p
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	17
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	18
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	19 H % Discrete norm
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	20 Hi
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	21 e_l, e_r, e_L, e_R;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	22
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	23 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	24
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	25 methods
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	26 function obj = Euler1d(m,xlim,order,gama,opsGen,do_upwind)
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	27 default_arg('opsGen',@sbp.Ordinary);
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	28 default_arg('gama', 1.4);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	29 default_arg('do_upwind', false);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	30 gamma = gama;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	31
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	32 [x, h] = util.get_grid(xlim{:},m);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	33
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	34 if do_upwind
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	35 ops = spb.Upwind(m,h,order);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	36 Dp = ops.derivatives.Dp;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	37 Dm = ops.derivatives.Dm;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	38
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	39 printExpr('issparse(Dp)');
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	40 printExpr('issparse(Dm)');
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	41
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	42 D1 = (Dp + Dm)/2;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	43 else
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	44 ops = opsGen(m,h,order);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	45 D1 = sparse(ops.derivatives.D1);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	46 end
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	47
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	48 H = sparse(ops.norms.H);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	49 Hi = sparse(ops.norms.HI);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	50 e_l = sparse(ops.boundary.e_1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	51 e_r = sparse(ops.boundary.e_m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	52
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	53 I_x = speye(m);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	54 I_3 = speye(3);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	55
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	56
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	57 D1 = kr(D1, I_3);
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	58 if do_upwind
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	59 Ddisp = kr(Ddisp,I_3);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	60 end
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	61
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	62 % Norms
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	63 obj.H = kr(H,I_3);
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	64 obj.Hi = kr(Hi,I_3);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	65
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	66 % Boundary operators
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	67 obj.e_l = e_l;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	68 obj.e_r = e_r;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	69 obj.e_L = kr(e_l,I_3);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	70 obj.e_R = kr(e_r,I_3);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	71
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	72 obj.m = m;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	73 obj.h = h;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	74 obj.order = order;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	75
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	76 % Man har Q_t+F_x=0 i 1D Euler, där
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	77 % q=[rho, rho*u, e]^T
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	78 % F=[rhou, rhou^2+p, (e+p)*u] ^T
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	79 % p=(gamma-1)(e-rhou^2/2);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	80
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	81 %Solving on form q_t + F_x = 0
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	82 function o = F(q)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	83 o = [q(2); q(2).^2/q(1) + p(q); (q(3)+p(q))*q(2)/q(1)];
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	84 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	85
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	86 % Equation of state
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	87 function o = p(q)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	88 o = (gamma-1)*(q(3)-q(2).^2/q(1)/2);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	89 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	90
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	91 % R =
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	92 % [sqrt(2(gamma-1))rho , rho , rho ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	93 % sqrt(2(gamma-1))rhou , rho(u+c) , rho*(u-c) ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	94 % sqrt(2(gamma-1))rhou^2/2, e+(gamma-1)(e-rhou^2/2)+rhouc, e+(gamma-1)(e-rhou^2/2)-rhou*c]);
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	95 function o = T(q)
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	96 rho = q(1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	97 u = q(2)/q(1);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	98 e = q(3);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	99
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	100 c = sqrt(gamma*p(q)/rho);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	101
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	102 sqrt2gamm = sqrt(2*(gamma-1));
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	103
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	104 o = [
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	105 sqrt2gamm*rho , rho , rho ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	106 sqrt2gammrhou , rho(u+c) , rho(u-c) ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	107 sqrt2gammrhou^2/2, e+(gamma-1)(e-rhou^2/2)+rhouc , e+(gamma-1)(e-rhou^2/2)-rhouc
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	108 ];
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	109 % Devide columns by stuff to get rid of extra comp?
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	110 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	111
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	112 function o = Fx(q)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	113 o = zeros(size(q));
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	114 for i = 1:3:3*m
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	115 o(i:i+2) = F(q(i:i+2));
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	116 end
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	117 o = D1*o;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	118 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	119
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	120
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	121 % A=RLambdainv(R), där Lambda=diag(u, u+c, u-c) (c är ljudhastigheten)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	122 % c^2=gamma*p/rho
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	123 % function o = A(rho,u,e)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	124 % end
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	125 if do_upwind
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	126 obj.D = @(q)-Fx(q) + Ddisp(1)u;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	127 else
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	128 obj.D = @(q)-Fx(q);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	129 end
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	130
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	131 obj.Fx = @Fx;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	132 obj.T = @T;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	133 obj.u = x;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	134 obj.x = kr(x,ones(3,1));
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	135 obj.p = @p;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	136 obj.gamma = gamma;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	137 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	138
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	139
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	140 % Enforces the boundary conditions
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	141 % w+ = R*w- + g(t)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	142 function closure = boundary_condition(obj,boundary, type, varargin)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	143 [e_s,e_S,s] = obj.get_boundary_ops(boundary);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	144
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	145 % Boundary condition on form
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	146 % w_in = R*w_out + g, where g is data
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	147
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	148 % How to handle when the number of BC we want to set changes
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	149 % How to handel u = 0 as an example
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	150
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	151 % Måste sätta in s och fixa tecken som de ska vara
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	152
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	153 % Kanske ska man tala om vilka karakteristikor som är in och ut i anropet
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	154 % Man kan sen kolla om det stämmer (men hur blir det med u=0?)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	155 % Och antingen tillåta att man skickar in flera alternativ och väljer automatiskt
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	156 % eller låta koden utanför bestämma vilken penalty som ska appliceras.
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	157
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	158 switch type
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	159 case 'wall'
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	160 closure = obj.boundary_condition_wall(boundary);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	161 otherwise
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	162 error('Unsupported bc type: %s', type);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	163 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	164
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	165 % T = obj.T(e_S*v);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	166
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	167 % c = sqrt(gamma*p/rho);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	168 % l = [u, u+c, u-c];
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	169
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	170 % p_in = find(s*l <= 0);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	171 % p_ot = find(s*l > 0);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	172
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	173 % p = [p_in, p_ot] % Permutation to sort
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	174 % pt(p) = 1:length(p); % Inverse permutation
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	175
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	176 % L_in = diag(abs(l(p_in)));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	177 % L_ot = diag(abs(l(p_ot)));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	178
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	179 % Lambda = diag(u, u+c, u-c);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	180
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	181 % tau = -e_Ssparse(T[L_in; R'*L_in]); % Något med pt
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	182
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	183
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	184 % w_s = T'(e_Sv);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	185 % w_in = w_s(p_in);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	186 % w_ot = w_s(p_ot);
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	187
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	188
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	189 % function closure_fun(q,t)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	190 % obj.Hi * tau * (w_in - R*w_ot - g(t));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	191 % end
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	192
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	193 % function closure_fun_indep(q,t)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	194 % obj.Hi * tau * (w_in - R*w_ot - g;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	195 % end
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	196
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	197
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	198 % switch class(g)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	199 % case 'double'
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	200 % closure = @closure_fun;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	201 % case 'function_handle'
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	202 % closure = @closure_fun_indep;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	203 % otherwise
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	204 % error('Wierd data argument!');
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	205 % end
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	206
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	207 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	208
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	209 % Set wall boundary condition v = 0.
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	210 function closure = boundary_condition_wall(obj,boundary)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	211 [e_s,e_S,s] = obj.get_boundary_ops(boundary);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	212
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	213 % v = 0 corresponds to
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	214 % L = [0 1 0];
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	215 % g = 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	216 %
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	217 % Tp =
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	218 % R = -(u-c)/(u+c)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	219
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	220 % tau = alpha * (u+c)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	221 % (alpha+1)(u+c) + 1/4* alpha^2\|u-c\| <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	222 % 4*(alpha+1)(u+c) + alpha^2\|u-c\| <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	223 % 4 * (alpha+1)(u+c) + alpha^2\|u\| + alpha^2*c <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	224 % \|u\|(sgn(u)4 + sgn(u)4alpha + alpha^2) + c*(4 + 4alpha + alpha^2) <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	225 % \|u\|(alpha^2 + 4sgn(u)alpha + 4sgn(u)) + c*(alpha+2)^2 <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	226 % \|u\|[(alpha + 2sgn(u))^2 - 4(sgn(u)-1)] + c(alpha+2)^2 <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	227
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	228
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	229 % om vi låtsas att u = 0:
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	230 % (alpha+1)c + 1/4 * alpha^2*c <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	231 % alpha^2 + 4*alpha +4 <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	232 % (alpha + 2)^2 <= 0
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	233 % alpha = -2 gives tau = -2*c;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	234
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	235
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	236 % Vill vi sätta penalty på karateristikan som är nära noll också eller vill
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	237 % vi låta den vara fri?
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	238
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	239
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	240 switch s
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	241 case -1
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	242 p_in = 2;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	243 p_zero = 1;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	244 p_ot = 3;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	245 case 1
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	246 p_in = 3;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	247 p_zero = 1;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	248 p_ot = 2;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	249 otherwise
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	250 error();
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	251 end
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	252
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	253 p = [p_in, p_zero, p_ot]; % Permutation to sort
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	254 pt(p) = 1:length(p); % Inverse permutation
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	255
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	256 function o = closure_fun(q)
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	257 p = obj.p(q);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	258
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	259 q_s = e_S'*q;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	260 rho = q_s(1);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	261 u = q_s(2)/rho;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	262 c = sqrt(obj.gamma*p/rho);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	263
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	264 T = obj.T(q_s);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	265 R = -(u-c)/(u+c);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	266 % l = [u, u+c, u-c];
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	267
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	268 % p_in = find(s*l <= 0);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	269 % p_ot = find(s*l > 0);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	270
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	271
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	272 tau1 = -2*c;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	273 tau2 = [0; 0]; % Penalty only on ingoing char.
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	274
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	275 % L_in = diag(abs(l(p_in)));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	276 % L_ot = diag(abs(l(p_ot)));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	277
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	278 tauHat = [tau1; tau2];
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	279 tau = -se_Ssparse(T*tauHat(pt));
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	280
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	281 w_s = inv(T)*q_s;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	282 % w_s = T\q_s;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	283 % w_s = Tinv * q_s; % Med analytisk matris
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	284 w_in = w_s(p_in);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	285 w_ot = w_s(p_ot);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	286
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	287 o = obj.Hi * tau * (w_in - R*w_ot);
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	288 end
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	289
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	290 closure = @closure_fun;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	291 end
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	292
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	293 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	294 error('NOT DONE')
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	295 % u denotes the solution in the own domain
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	296 % v denotes the solution in the neighbour domain
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	297 [e_u,d1_u,d2_u,d3_u,s_u,gamm_u,delt_u, halfnorm_inv] = obj.get_boundary_ops(boundary);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	298 [e_v,d1_v,d2_v,d3_v,s_v,gamm_v,delt_v] = neighbour_scheme.get_boundary_ops(neighbour_boundary);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	299
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	300 tuning = 2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	301
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	302 alpha_u = obj.alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	303 alpha_v = neighbour_scheme.alpha;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	304
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	305 tau1 = ((alpha_u/2)/delt_u + (alpha_v/2)/delt_v)/2*tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	306 % tau1 = (alpha_u/2 + alpha_v/2)/(2delt_u)tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	307 tau4 = s_u*alpha_u/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	308
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	309 sig2 = ((alpha_u/2)/gamm_u + (alpha_v/2)/gamm_v)/2*tuning;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	310 sig3 = -s_u*alpha_u/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	311
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	312 phi2 = s_u*1/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	313
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	314 psi1 = -s_u*1/2;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	315
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	316 tau = tau1e_u + tau4d3_u;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	317 sig = sig2d1_u + sig3d2_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	318 phi = phi2*d1_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	319 psi = psi1*e_u ;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	320
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	321 closure = halfnorm_inv(taue_u' + sigd1_u' + phialpha_ud2_u' + psialpha_u*d3_u');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	322 penalty = -halfnorm_inv(taue_v' + sigd1_v' + phialpha_vd2_v' + psialpha_v*d3_v');
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	323 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	324
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	325 % Ruturns the boundary ops and sign for the boundary specified by the string boundary.
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	326 % The right boundary is considered the positive boundary
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	327 function [e,E,s] = get_boundary_ops(obj,boundary)
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	328 switch boundary
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	329 case 'l'
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	330 e = obj.e_l;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	331 E = obj.e_L;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	332 s = -1;
49 8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	333 case 'r'
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	334 e = obj.e_r;
8f0c2dc747dd Made lots of updates to Euler1d. Jonatan Werpers <jonatan@werpers.com> parents: 0 diff changeset	335 E = obj.e_R;
0 48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	336 s = 1;
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	337 otherwise
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	338 error('No such boundary: boundary = %s',boundary);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	339 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	340 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	341
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	342 function N = size(obj)
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	343 N = prod(obj.m);
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	344 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	345
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	346 end
48b6fb693025 Initial commit. Jonatan Werpers <jonatan@werpers.com> parents: diff changeset	347 end

Mercurial > repos > public > sbplib

annotate +scheme/Euler1d.m @ 49:8f0c2dc747dd