sbplib: +scheme/Hypsyst2dCurve.m annotate

author	Ylva Rydin <ylva.rydin@telia.com>
date	Tue, 04 Oct 2016 08:40:42 +0200
parents
children	4d8d6eb0c116

rev	line source
298 861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	1 classdef Hypsyst2dCurve < scheme.Scheme
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	2 properties
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	3 m % Number of points in each direction, possibly a vector
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	4 n %size of system
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	5 h % Grid spacing
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	6 X,Y % Values of x and y for each grid point
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	7
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	8 J, Ji
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	9 xi,eta
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	10 Xi,Eta
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	11
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	12 A,B
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	13 X_eta, Y_eta
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	14 X_xi,Y_xi
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	15 order % Order accuracy for the approximation
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	16
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	17 D % non-stabalized scheme operator
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	18 Ahat, Bhat, E
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	19
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	20 H % Discrete norm
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	21 % Norms in the x and y directions
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	22 Hxii,Hetai % Kroneckerd norms. 1'Hxv corresponds to integration in the x dir.
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	23 I_xi,I_eta, I_N, onesN
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	24 e_w, e_e, e_s, e_n
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	25 index_w, index_e,index_s,index_n
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	26 params %parameters for the coeficient matrice
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	27 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	28
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	29
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	30 methods
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	31 function obj = Hypsyst2dCurve(m, order, A, B, E, params, ti)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	32 default_arg('E', [])
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	33 xilim = {0 1};
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	34 etalim = {0 1};
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	35
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	36 if length(m) == 1
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	37 m = [m m];
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	38 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	39 obj.params = params;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	40 obj.A=A;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	41 obj.B=B;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	42
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	43
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	44 obj.Ahat=@(params,x,y,x_eta,y_eta)(A(params,x,y).y_eta-B(params,x,y).x_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	45 obj.Bhat=@(params,x,y,x_xi,y_xi)(B(params,x,y).x_xi-A(params,x,y).y_xi);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	46 obj.E=@(params,x,y,~,~)E(params,x,y);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	47
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	48 m_xi = m(1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	49 m_eta = m(2);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	50 m_tot=m_xi*m_eta;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	51
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	52 ops_xi = sbp.D2Standard(m_xi,xilim,order);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	53 ops_eta = sbp.D2Standard(m_eta,etalim,order);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	54
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	55 obj.xi = ops_xi.x;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	56 obj.eta = ops_eta.x;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	57
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	58 obj.Xi = kr(obj.xi,ones(m_eta,1));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	59 obj.Eta = kr(ones(m_xi,1),obj.eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	60
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	61 obj.n = length(A(obj.params,0,0));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	62 obj.onesN=ones(obj.n);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	63
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	64 obj.index_w=1:m_xi;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	65 obj.index_e=(m_tot-m_xi+1):m_tot;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	66 obj.index_s=1:m_xi:(m_tot-m_xi+1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	67 obj.index_n=(m_xi):m_xi:m_tot;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	68
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	69 I_n = eye(obj.n);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	70 I_xi = speye(m_xi);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	71 obj.I_xi = I_xi;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	72 I_eta = speye(m_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	73 obj.I_eta = I_eta;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	74
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	75 D1_xi = kr(I_n, ops_xi.D1, I_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	76 obj.Hxii = kr(I_n, ops_xi.HI, I_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	77 D1_eta = kr(I_n, I_xi, ops_eta.D1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	78 obj.Hetai = kr(I_n, I_xi, ops_eta.HI);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	79
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	80 obj.e_w = kr(I_n, ops_xi.e_l, I_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	81 obj.e_e = kr(I_n, ops_xi.e_r, I_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	82 obj.e_s = kr(I_n, I_xi, ops_eta.e_l);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	83 obj.e_n = kr(I_n, I_xi, ops_eta.e_r);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	84
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	85 [X,Y] = ti.map(obj.xi,obj.eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	86
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	87 [x_xi,x_eta] = gridDerivatives(X,ops_xi.D1,ops_eta.D1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	88 [y_xi,y_eta] = gridDerivatives(Y,ops_xi.D1, ops_eta.D1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	89
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	90 obj.X=reshape(X,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	91 obj.Y=reshape(Y,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	92 obj.X_xi=reshape(x_xi,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	93 obj.Y_xi=reshape(y_xi,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	94 obj.X_eta=reshape(x_eta,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	95 obj.Y_eta=reshape(y_eta,m_xi*m_eta,1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	96
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	97 Ahat_evaluated = obj.evaluateCoefficientMatrix(obj.Ahat, obj.X, obj.Y,obj.X_eta,obj.Y_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	98 Bhat_evaluated = obj.evaluateCoefficientMatrix(obj.Bhat, obj.X, obj.Y,obj.X_xi,obj.Y_xi);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	99 E_evaluated = obj.evaluateCoefficientMatrix(obj.E, obj.X, obj.Y,[],[]);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	100
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	101 obj.m=m;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	102 obj.h=[ops_xi.h ops_eta.h];
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	103 obj.order=order;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	104 obj.J=x_xi.y_eta - x_eta.y_xi;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	105 obj.Ji =kr(I_n,spdiags(1./obj.J(:),0,m_xim_eta,m_xim_eta));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	106
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	107 obj.D=obj.Ji(-Ahat_evaluatedD1_xi-Bhat_evaluated*D1_eta)-E_evaluated;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	108
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	109 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	110
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	111 % Closure functions return the opertors applied to the own doamin to close the boundary
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	112 % Penalty functions return the opertors to force the solution. In the case of an interface it returns the operator applied to the other doamin.
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	113 % boundary is a string specifying the boundary e.g. 'l','r' or 'e','w','n','s'.
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	114 % type is a string specifying the type of boundary condition if there are several.
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	115 % data is a function returning the data that should be applied at the boundary.
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	116 function [closure, penalty] = boundary_condition(obj,boundary,type,L)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	117 default_arg('type','char');
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	118 switch type
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	119 case{'c','char'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	120 [closure,penalty]=boundary_condition_char(obj,boundary);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	121 case{'general'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	122 [closure,penalty]=boundary_condition_general(obj,boundary,L);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	123 otherwise
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	124 error('No such boundary condition')
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	125 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	126 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	127
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	128 function [closure, penalty] = interface(obj,boundary,neighbour_scheme,neighbour_boundary)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	129 error('An interface function does not exist yet');
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	130 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	131
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	132 function N = size(obj)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	133 N = obj.m;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	134 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	135
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	136 function [ret] = evaluateCoefficientMatrix(obj, mat, X, Y,x_,y_)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	137 params=obj.params;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	138
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	139 if isa(mat,'function_handle')
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	140 [rows,cols]=size(mat(params,0,0,0,0));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	141 x_=kr(x_,obj.onesN);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	142 y_=kr(y_,obj.onesN);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	143 matVec=mat(params,X',Y',x_',y_');
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	144 matVec=sparse(matVec);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	145 side=max(length(X),length(Y));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	146 else
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	147 matVec=mat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	148 [rows,cols]=size(matVec);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	149 side=max(length(X),length(Y));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	150 cols=cols/side;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	151 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	152 ret=kron(ones(rows,cols),speye(side));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	153
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	154 for ii=1:rows
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	155 for jj=1:cols
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	156 ret((ii-1)side+1:iiside,(jj-1)side+1:jjside)=diag(matVec(ii,(jj-1)side+1:jjside));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	157 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	158 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	159 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	160
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	161
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	162 function [closure, penalty]=boundary_condition_char(obj,boundary)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	163 params=obj.params;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	164 xi=obj.xi; eta=obj.eta;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	165 side=max(length(xi),length(eta));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	166
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	167 switch boundary
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	168 case {'w','W','west'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	169 e_=obj.e_w;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	170 mat=obj.Ahat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	171 boundPos='l';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	172 Hi=obj.Hxii;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	173 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi(1),eta,obj.X_eta(obj.index_w),obj.Y_eta(obj.index_w));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	174 case {'e','E','east'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	175 e_=obj.e_e;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	176 mat=obj.Ahat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	177 boundPos='r';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	178 Hi=obj.Hxii;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	179 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi(end),eta,obj.X_eta(obj.index_e),obj.Y_eta(obj.index_e));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	180 case {'s','S','south'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	181 e_=obj.e_s;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	182 mat=obj.Bhat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	183 boundPos='l';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	184 Hi=obj.Hetai;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	185 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi,eta(1),obj.X_xi(obj.index_s),obj.Y_xi(obj.index_s));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	186 case {'n','N','north'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	187 e_=obj.e_n;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	188 mat=obj.Bhat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	189 boundPos='r';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	190 Hi=obj.Hetai;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	191 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi,eta(end),obj.X_xi(obj.index_n),obj.Y_xi(obj.index_n));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	192 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	193
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	194 pos=signVec(1); zeroval=signVec(2); neg=signVec(3);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	195
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	196 switch boundPos
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	197 case {'l'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	198 tau=sparse(obj.nside,posside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	199 Vi_plus=Vi(1:pos*side,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	200 tau(1:posside,:)=-abs(D(1:posside,1:pos*side));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	201 closure=Hie_VtauVi_plus*e_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	202 penalty=-Hie_VtauVi_plus;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	203 case {'r'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	204 tau=sparse(obj.nside,negside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	205 tau((pos+zeroval)side+1:obj.nside,:)=-abs(D((pos+zeroval)side+1:obj.nside,(pos+zeroval)side+1:obj.nside));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	206 Vi_minus=Vi((pos+zeroval)side+1:obj.nside,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	207 closure=Hie_VtauVi_minus*e_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	208 penalty=-Hie_VtauVi_minus;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	209 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	210 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	211
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	212
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	213 function [closure,penalty]=boundary_condition_general(obj,boundary,L)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	214 params=obj.params;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	215 xi=obj.xi; eta=obj.eta;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	216 side=max(length(xi),length(eta));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	217
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	218 switch boundary
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	219 case {'w','W','west'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	220 e_=obj.e_w;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	221 mat=obj.Ahat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	222 boundPos='l';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	223 Hi=obj.Hxii;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	224 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi(1),eta,obj.x_eta,obj.y_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	225 L=obj.evaluateCoefficientMatrix(L,xi(1),eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	226 case {'e','E','east'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	227 e_=obj.e_e;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	228 mat=obj.Ahat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	229 boundPos='r';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	230 Hi=obj.Hxii;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	231 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi(end),eta,obj.x_eta,obj.y_eta);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	232 L=obj.evaluateCoefficientMatrix(L,xi(end),eta,[],[]);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	233 case {'s','S','south'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	234 e_=obj.e_s;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	235 mat=obj.Bhat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	236 boundPos='l';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	237 Hi=obj.Hetai;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	238 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi,eta(1),obj.x_xi,obj.y_xi);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	239 L=obj.evaluateCoefficientMatrix(L,xi,eta(1));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	240 case {'n','N','north'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	241 e_=obj.e_n;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	242 mat=obj.Bhat;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	243 boundPos='r';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	244 Hi=obj.Hetai;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	245 [V,Vi,D,signVec]=obj.matrixDiag(mat,xi,eta(end),obj.x_xi,obj.y_xi);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	246 L=obj.evaluateCoefficientMatrix(L,xi,eta(end));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	247 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	248
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	249 pos=signVec(1); zeroval=signVec(2); neg=signVec(3);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	250
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	251 switch boundPos
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	252 case {'l'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	253 tau=sparse(obj.nside,posside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	254 Vi_plus=Vi(1:pos*side,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	255 Vi_minus=Vi(posside+1:obj.nside,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	256 V_plus=V(:,1:pos*side);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	257 V_minus=V(:,(pos+zeroval)side+1:obj.nside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	258
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	259 tau(1:posside,:)=-abs(D(1:posside,1:pos*side));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	260 R=-inv(LV_plus)(L*V_minus);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	261 closure=Hie_Vtau(Vi_plus-RVi_minus)e_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	262 penalty=-Hie_Vtauinv(LV_plus)L;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	263 case {'r'}
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	264 tau=sparse(obj.nside,negside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	265 tau((pos+zeroval)side+1:obj.nside,:)=-abs(D((pos+zeroval)side+1:obj.nside,(pos+zeroval)side+1:obj.nside));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	266 Vi_plus=Vi(1:pos*side,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	267 Vi_minus=Vi((pos+zeroval)side+1:obj.nside,:);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	268
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	269 V_plus=V(:,1:pos*side);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	270 V_minus=V(:,(pos+zeroval)side+1:obj.nside);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	271 R=-inv(LV_minus)(L*V_plus);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	272 closure=Hie_Vtau(Vi_minus-RVi_plus)e_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	273 penalty=-Hie_Vtauinv(LV_minus)L;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	274 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	275 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	276
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	277
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	278 function [V,Vi, D,signVec]=matrixDiag(obj,mat,x,y,x_,y_)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	279 params=obj.params;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	280 syms xs ys
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	281 if(sum(abs(x_))~=0)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	282 syms xs_
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	283 else
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	284 xs_=0;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	285 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	286
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	287 if(sum(abs(y_))~=0)
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	288 syms ys_;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	289 else
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	290 ys_=0;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	291 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	292
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	293 [V, D]=eig(mat(params,xs,ys,xs_,ys_));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	294 xs=1;ys=1; xs_=x_(1); ys_=y_(1);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	295 DD=eval(diag(D));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	296
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	297 poseig=find(DD>0);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	298 zeroeig=find(DD==0);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	299 negeig=find(DD<0);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	300 syms xs ys xs_ ys_
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	301 DD=diag(D);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	302
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	303 D=diag([DD(poseig);DD(zeroeig); DD(negeig)]);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	304 V=[V(:,poseig) V(:,zeroeig) V(:,negeig)];
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	305 Vi=inv(V);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	306 xs=x;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	307 ys=y;
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	308 xs_=x_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	309 ys_=y_';
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	310
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	311 side=max(length(x),length(y));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	312 Dret=zeros(obj.n,side*obj.n);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	313 Vret=zeros(obj.n,side*obj.n);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	314 Viret=zeros(obj.n,side*obj.n);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	315 for ii=1:obj.n
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	316 for jj=1:obj.n
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	317 Dret(jj,(ii-1)side+1:sideii)=eval(D(jj,ii));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	318 Vret(jj,(ii-1)side+1:sideii)=eval(V(jj,ii));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	319 Viret(jj,(ii-1)side+1:sideii)=eval(Vi(jj,ii));
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	320 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	321 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	322
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	323 D=sparse(Dret);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	324 V=sparse(Vret);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	325 Vi=sparse(Viret);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	326 V=obj.evaluateCoefficientMatrix(V,x,y,x_,y_);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	327 D=obj.evaluateCoefficientMatrix(D,x,y,x_,y_);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	328 Vi=obj.evaluateCoefficientMatrix(Vi,x,y,x_,y_);
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	329 signVec=[length(poseig),length(zeroeig),length(negeig)];
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	330 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	331
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	332 end
861972361f75 The curvelinear works for quads Ylva Rydin <ylva.rydin@telia.com> parents: diff changeset	333 end

298

861972361f75 The curvelinear works for quads