meshq1_ie

 MESHQ1_IE Construction of extended Q1 mesh

 [novl,jacx,jacy]=meshq1_ie(nov,nvl,xy,ipar);

 called by stiffq1.m

 Input: 
        nov = local -global map, previously generated by cosnov_2d
        nvl = number of nodes in each spectral element
        xy = column array with global mesh, length: noe=nov(npdx*npdy,ne)
        ipar = arry woth internal parameters.

 Output:
        novl =  2-indexes array. It maps  Q1 mesh on the Qn extended element 
        jacx, jacy = jacobians of Q1 mesh on the extended element

0001 function [novl,jacx,jacy]=meshq1_ie(nov,nvl,xy,ipar);
0002 % MESHQ1_IE Construction of extended Q1 mesh
0003 %
0004 % [novl,jacx,jacy]=meshq1_ie(nov,nvl,xy,ipar);
0005 %
0006 % called by stiffq1.m
0007 %
0008 % Input:
0009 %        nov = local -global map, previously generated by cosnov_2d
0010 %        nvl = number of nodes in each spectral element
0011 %        xy = column array with global mesh, length: noe=nov(npdx*npdy,ne)
0012 %        ipar = arry woth internal parameters.
0013 %
0014 % Output:
0015 %        novl =  2-indexes array. It maps  Q1 mesh on the Qn extended element
0016 %        jacx, jacy = jacobians of Q1 mesh on the extended element
0017 %
0018 
0019 %   Written by Paola Gervasio
0020 %   $Date: 2007/04/01$
0021 
0022 %
0023 
0024 noe=nvl(1);ne=ipar(6);npxl=ipar(1);npyl=ipar(2);mn=ipar(3);
0025 nexl=nvl(2); neyl=nvl(3);
0026 novl=zeros(4,ne);
0027 npx=nexl+1;
0028 jacx=zeros(ne,1); jacy=zeros(ne,1);
0029 for i=1:nexl
0030 for j=1:neyl
0031     k=(j-1)*npx+i;  iel=(j-1)*nexl+i;
0032     i1=k;i2=k+1;i3=k+npx;i4=k+npx+1;
0033     novl(1:4,iel)=[i1;i2;i3;i4];
0034     jacx(iel)=(xy(i2,1)-xy(i1,1))*0.5;
0035     jacy(iel)=(xy(i3,2)-xy(i1,2))*0.5;
0036 end
0037 end
0038 
0039 
0040             
0041        
0042     
0043