call_eig_schur_2d

 CALL_EIG_SCHUR_2D Script for pre and post processing eig_schur_2d.m

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