Domain decomposition methods: Difference between revisions

Revision as of 18:58, 15 February 2020

Revision as of 09:07, 23 July 2020

Line 38:

<source lang="matlab">

clear all

N = 4;

a1 = 0; b1 = 1/2;

[T D1 D2 E1 E2 x xsub] = cheb(N,a1,b1); % the diff matrices on [0,1/2] are the same

%as those on [1/2 1].

I = eye(N+1);

H = D2-I;

H1 = [[1 zeros(1,N)]; H(2:end-1,:); [zeros(1,N) 1]];

H1 = [H1 [zeros(N,N+1); -[1 zeros(1,N)]]];

H2 = [D1(1,:); H(2:end-1,:); [zeros(1,N) 1]];

H2 = [[-D1(N+1,:); zeros(N,N+1)] H2];

K = [H1; H2];

F = [zeros(2*N+1,1); 1];

u = K\F;

xx = -cos(pi*(0:N)'/N);

x1 = 1/4*(xx+1); x2 = 1/4*(xx+3);

x = [x1; x2];

uex = (exp(x)-exp(-x))./(exp(1)-exp(-1));

</source>