Commit ef5d2dd2 authored by Martin Averseng's avatar Martin Averseng
Browse files

spi

parent 51d1458d
......@@ -14,7 +14,7 @@ c = spirale;x1 = -1.5; x2 = 1.5; y1 = -1.5; y2 = 1.8; theta_inc = 0;
% c = Vcurve; x1 = -4; x2 = 4; y1 = -2; y2 = 3; theta_inc = pi/2;
k = 50*pi/2; N = fix(5*k*2)+1;
k = 400*pi/2; N = fix(5*k*2)+1;
m = meshCurve(c,N,'varChange',{@cos,[-pi,0]});
edges = bnd(m);
% Weight definition :
......@@ -80,8 +80,8 @@ figure;
plot(m.vtx(:,1),real(mu));
x = linspace(x1,x2,1e3);
y = linspace(y1,y2,1e3);
x = linspace(x1,x2,4e3);
y = linspace(y1,y2,4e3);
[X,Y] = meshgrid(x,y);
M = [X(:),Y(:),0*X(:)];
NM = size(M,1);
......@@ -92,7 +92,7 @@ K = H0Kernel(k);
Ny = size(Y,1);
Wyomega2 = spdiags(Wy.*omega2(Y),0,Ny,Ny);
Mv = uqm(ntimes(Vh),Gamma);
a = 1.5/sqrt(sqrt(NM*Ny));
a = 1/sqrt(sqrt(NM*Ny));
[Dx,Dy] = offline_dEBD(K,M,Y,a,tol);
A = AbstractMatrix([],...
@(V)(Dx(Wyomega2*(Mv{1}*V))...
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment