sbplib: +scheme/Hypsyst2d.m annotate

annotate +scheme/Hypsyst2d.m @ 366:7ada2db63268 feature/hypsyst

Tried to speed up the matrix set-up a little

author	Ylva Rydin <ylva.rydin@telia.com>
date	Tue, 13 Dec 2016 09:28:33 +0100
parents	9b3d7fc61a36
children	9d1fc984f40d

rev	line source
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	1 classdef Hypsyst2d < scheme.Scheme
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	2 properties
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
293 2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	4 n %size of system
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	5 h % Grid spacing
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	6 x,y % Grid
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	7 X,Y % Values of x and y for each grid point
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	8 order % Order accuracy for the approximation
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	9
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	10 D % non-stabalized scheme operator
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	11 A, B, E
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	12
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	13 H % Discrete norm
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	14 % Norms in the x and y directions
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	15 Hxi,Hyi % Kroneckerd norms. 1'Hxv corresponds to integration in the x dir.
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	16 I_x,I_y, I_N
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	17 e_w, e_e, e_s, e_n
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	18 params %parameters for the coeficient matrice
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	19 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	20
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	21
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	22 methods
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	23 function obj = Hypsyst2d(m, lim, order, A, B, E, params)
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	24 default_arg('E', [])
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	25 xlim = lim{1};
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	26 ylim = lim{2};
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	27
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	28 if length(m) == 1
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	29 m = [m m];
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	30 end
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	31
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	32 obj.A=A;
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	33 obj.B=B;
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	34 obj.E=E;
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	35
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	36 m_x = m(1);
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	37 m_y = m(2);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	38 obj.params = params;
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	39
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	40 ops_x = sbp.D2Standard(m_x,xlim,order);
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	41 ops_y = sbp.D2Standard(m_y,ylim,order);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	42
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	43 obj.x = ops_x.x;
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	44 obj.y = ops_y.x;
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	45
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	46 obj.X = kr(obj.x,ones(m_y,1));
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	47 obj.Y = kr(ones(m_x,1),obj.y);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	48
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	49 Aevaluated = obj.evaluateCoefficientMatrix(A, obj.X, obj.Y);
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	50 Bevaluated = obj.evaluateCoefficientMatrix(B, obj.X, obj.Y);
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	51 Eevaluated = obj.evaluateCoefficientMatrix(E, obj.X, obj.Y);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	52
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	53 obj.n = length(A(obj.params,0,0));
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	54
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	55 I_n = eye(obj.n);I_x = speye(m_x);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	56 obj.I_x = I_x;
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	57 I_y = speye(m_y);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	58 obj.I_y = I_y;
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	59
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	60
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	61 D1_x = kr(I_n, ops_x.D1, I_y);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	62 obj.Hxi = kr(I_n, ops_x.HI, I_y);
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	63 D1_y = kr(I_n, I_x, ops_y.D1);
cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	64 obj.Hyi = kr(I_n, I_x, ops_y.HI);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	65
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	66 obj.e_w = kr(I_n, ops_x.e_l, I_y);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	67 obj.e_e = kr(I_n, ops_x.e_r, I_y);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	68 obj.e_s = kr(I_n, I_x, ops_y.e_l);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	69 obj.e_n = kr(I_n, I_x, ops_y.e_r);
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	70
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	71 obj.m=m;
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	72 obj.h=[ops_x.h ops_y.h];
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	73 obj.order=order;
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	74
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	75 obj.D=-AevaluatedD1_x-BevaluatedD1_y-Eevaluated;
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	76
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	77 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	78
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	79 % Closure functions return the opertors applied to the own doamin to close the boundary
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	80 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	81 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	82 % type is a string specifying the type of boundary condition if there are several.
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	83 % data is a function returning the data that should be applied at the boundary.
293 2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	84 function [closure, penalty] = boundary_condition(obj,boundary,type,L)
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	85 default_arg('type','char');
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	86 switch type
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	87 case{'c','char'}
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	88 [closure,penalty]=boundary_condition_char(obj,boundary);
293 2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	89 case{'general'}
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	90 [closure,penalty]=boundary_condition_general(obj,boundary,L);
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	91 otherwise
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	92 error('No such boundary condition')
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	93 end
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	94 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	95
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	96 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	97 error('An interface function does not exist yet');
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	98 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	99
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	100 function N = size(obj)
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	101 N = obj.m;
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	102 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	103
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	104 function [ret] = evaluateCoefficientMatrix(obj, mat, X, Y)
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	105 params=obj.params;
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	106
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	107 if isa(mat,'function_handle')
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	108 [rows,cols]=size(mat(params,0,0));
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	109 matVec=mat(params,X',Y');
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	110 matVec=sparse(matVec);
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	111 side=max(length(X),length(Y));
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	112 else
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	113 matVec=mat;
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	114 [rows,cols]=size(matVec);
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	115 side=max(length(X),length(Y));
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	116 cols=cols/side;
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	117 end
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	118 ret=cell(rows,cols);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	119
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	120 for ii=1:rows
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	121 for jj=1:cols
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	122 ret{ii,jj}=diag(matVec(ii,(jj-1)side+1:jjside));
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	123 end
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	124 end
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	125 ret=cell2mat(ret);
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	126 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	127
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	128
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	129 function [closure, penalty]=boundary_condition_char(obj,boundary)
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	130 params=obj.params;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	131 x=obj.x; y=obj.y;
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	132
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	133 switch boundary
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	134 case {'w','W','west'}
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	135 e_=obj.e_w;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	136 mat=obj.A;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	137 boundPos='l';
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	138 Hi=obj.Hxi;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	139 [V,Vi,D,signVec]=obj.matrixDiag(mat,x(1),y);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	140 side=max(length(y));
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	141 case {'e','E','east'}
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	142 e_=obj.e_e;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	143 mat=obj.A;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	144 boundPos='r';
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	145 Hi=obj.Hxi;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	146 [V,Vi,D,signVec]=obj.matrixDiag(mat,x(end),y);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	147 side=max(length(y));
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	148 case {'s','S','south'}
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	149 e_=obj.e_s;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	150 mat=obj.B;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	151 boundPos='l';
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	152 Hi=obj.Hyi;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	153 [V,Vi,D,signVec]=obj.matrixDiag(mat,x,y(1));
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	154 side=max(length(x));
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	155 case {'n','N','north'}
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	156 e_=obj.e_n;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	157 mat=obj.B;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	158 boundPos='r';
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	159 Hi=obj.Hyi;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	160 [V,Vi,D,signVec]=obj.matrixDiag(mat,x,y(end));
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	161 side=max(length(x));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	162 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	163
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	164 pos=signVec(1); zeroval=signVec(2); neg=signVec(3);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	165
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	166 switch boundPos
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	167 case {'l'}
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	168 tau=sparse(obj.n*side,pos);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	169 Vi_plus=Vi(1:pos,:);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	170 tau(1:pos,:)=-abs(D(1:pos,1:pos));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	171 closure=Hie_VtauVi_plus*e_';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	172 penalty=-Hie_VtauVi_plus;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	173 case {'r'}
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	174 tau=sparse(obj.n*side,neg);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	175 tau((pos+zeroval)+1:obj.nside,:)=-abs(D((pos+zeroval)+1:obj.nside,(pos+zeroval)+1:obj.n*side));
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	176 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	177 closure=Hie_VtauVi_minus*e_';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	178 penalty=-Hie_VtauVi_minus;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	179 end
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	180 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	181
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	182
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	183 function [closure,penalty]=boundary_condition_general(obj,boundary,L)
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	184 params=obj.params;
8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	185 x=obj.x; y=obj.y;
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	186
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	187 switch boundary
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	188 case {'w','W','west'}
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	189 e_=obj.e_w;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	190 mat=obj.A;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	191 boundPos='l';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	192 Hi=obj.Hxi;
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	193 [V,Vi,D,signVec]=obj.matrixDiag(mat,x(1),y);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	194 L=obj.evaluateCoefficientMatrix(L,x(1),y);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	195 side=max(length(y));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	196 case {'e','E','east'}
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	197 e_=obj.e_e;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	198 mat=obj.A;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	199 boundPos='r';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	200 Hi=obj.Hxi;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	201 [V,Vi,D,signVec]=obj.matrixDiag(mat,x(end),y);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	202 L=obj.evaluateCoefficientMatrix(L,x(end),y);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	203 side=max(length(y));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	204 case {'s','S','south'}
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	205 e_=obj.e_s;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	206 mat=obj.B;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	207 boundPos='l';
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	208 Hi=obj.Hyi;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	209 [V,Vi,D,signVec]=obj.matrixDiag(mat,x,y(1));
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	210 L=obj.evaluateCoefficientMatrix(L,x,y(1));
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	211 side=max(length(x));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	212 case {'n','N','north'}
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	213 e_=obj.e_n;
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	214 mat=obj.B;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	215 boundPos='r';
297 cd30b22cee56 Have tried to make a curvelinear sheme for hypsysts. Does not really work yet... Ylva Rydin <ylva.rydin@telia.com> parents: 296 diff changeset	216 Hi=obj.Hyi;
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	217 [V,Vi,D,signVec]=obj.matrixDiag(mat,x,y(end));
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	218 L=obj.evaluateCoefficientMatrix(L,x,y(end));
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	219 side=max(length(x));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	220 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	221
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	222 pos=signVec(1); zeroval=signVec(2); neg=signVec(3);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	223
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	224 switch boundPos
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	225 case {'l'}
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	226 tau=sparse(obj.n*side,pos);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	227 Vi_plus=Vi(1:pos,:);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	228 Vi_minus=Vi(pos+zeroval+1:obj.n*side,:);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	229 V_plus=V(:,1:pos);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	230 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	231
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	232 tau(1:pos,:)=-abs(D(1:pos,1:pos));
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	233 R=-inv(LV_plus)(L*V_minus);
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	234 closure=Hie_Vtau(Vi_plus-RVi_minus)e_';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	235 penalty=-Hie_Vtauinv(LV_plus)L;
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	236 case {'r'}
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	237 tau=sparse(obj.n*side,neg);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	238 tau((pos+zeroval)+1:obj.nside,:)=-abs(D((pos+zeroval)+1:obj.nside,(pos+zeroval)+1:obj.n*side));
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	239 Vi_plus=Vi(1:pos,:);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	240 Vi_minus=Vi((pos+zeroval)+1:obj.n*side,:);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	241
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	242 V_plus=V(:,1:pos);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	243 V_minus=V(:,(pos+zeroval)+1:obj.n*side);
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	244 R=-inv(LV_minus)(L*V_plus);
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	245 closure=Hie_Vtau(Vi_minus-RVi_plus)e_';
3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	246 penalty=-Hie_Vtauinv(LV_minus)L;
294 8ff6ec6249e8 "General" boundary conditions implemented Ylva Rydin <ylva.rydin@telia.com> parents: 293 diff changeset	247 end
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	248 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	249
da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	250
292 3d275c5e45b3 Changed how the matrices are built Ylva Rydin <ylva.rydin@telia.com> parents: 291 diff changeset	251 function [V,Vi, D,signVec]=matrixDiag(obj,mat,x,y)
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	252 params=obj.params;
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	253 syms xs ys
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	254 [V, D]=eig(mat(params,xs,ys));
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	255 Vi=inv(V);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	256 xs=x;
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	257 ys=y;
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	258
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	259 side=max(length(x),length(y));
293 2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	260 Dret=zeros(obj.n,side*obj.n);
2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	261 Vret=zeros(obj.n,side*obj.n);
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	262 Viret=zeros(obj.n,side*obj.n);
293 2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	263 for ii=1:obj.n
2d604d16842c Works with varying coefficients and char boundary condition Ylva Rydin <ylva.rydin@telia.com> parents: 292 diff changeset	264 for jj=1:obj.n
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	265 Dret(jj,(ii-1)side+1:sideii)=eval(D(jj,ii));
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	266 Vret(jj,(ii-1)side+1:sideii)=eval(V(jj,ii));
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	267 Viret(jj,(ii-1)side+1:sideii)=eval(Vi(jj,ii));
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	268 end
807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	269 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	270
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	271 D=sparse(Dret);
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	272 V=sparse(Vret);
9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	273 Vi=sparse(Viret);
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	274 V=obj.evaluateCoefficientMatrix(V,x,y);
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	275 Vi=obj.evaluateCoefficientMatrix(Vi,x,y);
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	276 D=obj.evaluateCoefficientMatrix(D,x,y);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	277 DD=diag(D);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	278
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	279 poseig=(DD>0);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	280 zeroeig=(DD==0);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	281 negeig=(DD<0);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	282
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	283 D=diag([DD(poseig); DD(zeroeig); DD(negeig)]);
d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	284 V=[V(:,poseig) V(:,zeroeig) V(:,negeig)];
352 9b3d7fc61a36 Changed operator in Utux Ylva Rydin <ylva.rydin@telia.com> parents: 301 diff changeset	285 Vi=[Vi(poseig,:); Vi(zeroeig,:); Vi(negeig,:)];
301 d9860ebc3148 HypsystCurve2D Seems to work (Converges with MMS) Ylva Rydin <ylva.rydin@telia.com> parents: 297 diff changeset	286 signVec=[sum(poseig),sum(zeroeig),sum(negeig)];
291 807dfe8be3ec Have made a lot of stupid changes in hypsyst in order to find a stupid bug Ylva Rydin <ylva.rydin@telia.com> parents: 290 diff changeset	287 end
295 da0131655035 Fixed some formatting and naming. Jonatan Werpers <jonatan@werpers.com> parents: 294 diff changeset	288
290 d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	289 end
d32f674bcbe5 A first attempt to make a general scheme fo hyperbolic systems Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	290 end

Mercurial > repos > public > sbplib

annotate +scheme/Hypsyst2d.m @ 366:7ada2db63268 feature/hypsyst