patch_se

 PATCH_SE Imposes strong continuity of normal derivatives across interfaces

 called by lap_2d

0001 function [A,f]=patch_se(A,f,ifro,nov,dx,jacx,dy,jacy);
0002 % PATCH_SE Imposes strong continuity of normal derivatives across interfaces
0003 %
0004 % called by lap_2d
0005 %
0006 
0007 %   Written by Paola Gervasio
0008 %   $Date: 2007/04/01$
0009 
0010 [ldnov,ne]=size(nov);
0011 noe=length(f); 
0012 npdx=length(dx);
0013 npdy=length(dy);
0014 
0015 
0016 %  clear the rows of A related to interface unknowns
0017 
0018 for i=1:noe
0019 if (ifro(i)==-1)
0020 A(i,:)=zeros(1,noe);
0021 f(i)=0;
0022 end
0023 end
0024 
0025 % derivative
0026 mn=npdx*npdy;
0027 for ie=1:ne
0028 Al=patch_sp(dx,jacx(ie),dy,jacy(ie));
0029 for j=1:npdy
0030 for i=1:npdx
0031 k=(j-1)*npdx+i;
0032 ii=nov(k,ie);
0033 if(ifro(ii)==-1)
0034 A(ii,nov(1:mn,ie))=A(ii,nov(1:mn,ie))+Al(k,:);
0035 end
0036 end
0037 end
0038 end
0039 return