


ICEREDUCEMATRIX - reduce stiffness matrix from g to n
Usage:
Knn= IceReducematrix( Kgg, Gmn, pv_m,pv_n,msize,flag)

0001 function Knn= IceReducematrix( Kgg, Gmn, pv_m,pv_n,msize,flag) 0002 %ICEREDUCEMATRIX - reduce stiffness matrix from g to n 0003 % 0004 % Usage: 0005 % Knn= IceReducematrix( Kgg, Gmn, pv_m,pv_n,msize,flag) 0006 0007 if msize>0, 0008 Kmm=Kgg(pv_m,pv_m); 0009 Kmn=Kgg(pv_m,pv_n); 0010 Knm=Kgg(pv_n,pv_m); 0011 Knn=Kgg(pv_n,pv_n); 0012 0013 Knn= Knn + Knm * Gmn; 0014 0015 if flag~=1, 0016 Kmn = Kmn + Kmm * Gmn; 0017 Knn = Knn + Gmn' * Kmn; 0018 end 0019 0020 else 0021 Knn=Kgg; 0022 end