CALL_EIG_SCHWARZ_2D Script file for pre and post processing eig_schwarz_2d.m
0001 % CALL_EIG_SCHWARZ_2D Script file for pre and post processing eig_schwarz_2d.m 0002 % 0003 % 0004 0005 % Written by Paola Gervasio 0006 % $Date: 2007/04/01$ 0007 0008 xa=-1;xb=1; % Omega=(xa,xb) x (ya,yb) 0009 ya=-1;yb=1; 0010 cb='dddd'; % eig_schwarz_2d works only if cb='dddd'; 0011 gammax=[]; gammay=[]; % if SEM decomposition is not uniform: 0012 % they are the arrays with intefaces positions 0013 param=zeros(20,1); 0014 param(1)=2; % 1:P=I, 2: P=P_as 0015 param(2)=2; % number of levels for overlapping elements 0016 fprintf('nx nex lam_max(S) lam_min(S) k(S) \n') 0017 for nx=8 % polynomial degree in each element along x-direction 0018 ny=nx; % polynomial degree in each element along y-direction 0019 for nex=3:10; 0020 ney=nex; % decomposition of Omega in nex x ney rectangles 0021 0022 % call eig_schwarz_2d 0023 0024 [param]=eig_schwarz_2d(xa,xb,ya,yb,gam,... 0025 cb,nex,nx,ney,ny,gammax,gammay,param); 0026 0027 % output 0028 fprintf('nx=%d,nex=%d,l_max(P^(-1)A)=%11.4e, l_min(P^(-1)A)=%11.4e, k(P^(-1)A)=%11.4e \n',... 0029 nx,nex,param(11),param(12),param(13)); 0030 0031 end 0032 end