vkjac

 jacobian for von karman-plate

0001 function [c,fu,flam,b]=vkjac(p,u,lam) 
0002 % jacobian for von karman-plate
0003 [c,a,f,b]=vkf(p,u,lam);  % re-use c from vkf
0004 b=0;fu=zeros(p.neq^2,p.nt); v1=ones(1,p.nt); 
0005 fu(11,:)=-v1; fu(33,:)=-v1; fu(45,:)=-v1; fu(56,:)=-v1; 
0006 fu(67,:)=-v1; fu(78,:)=-v1; fu(89,:)=-v1; fu(100,:)=-v1;
0007 u5=pdeintrp(p.points,p.tria,u(4*p.np+1:5*p.np)); 
0008 u6=pdeintrp(p.points,p.tria,u(5*p.np+1:6*p.np)); 
0009 u7=pdeintrp(p.points,p.tria,u(6*p.np+1:7*p.np)); 
0010 u8=pdeintrp(p.points,p.tria,u(7*p.np+1:8*p.np)); 
0011 u9=pdeintrp(p.points,p.tria,u(8*p.np+1:9*p.np)); 
0012 u10=pdeintrp(p.points,p.tria,u(9*p.np+1:10*p.np));
0013 fu(42,:)=-u9; fu(52,:)=-u8; fu(62,:)=2*u10; 
0014 fu(72,:)=-u6; fu(82,:)=-u5; fu(92,:)=2*u7;
0015 fu(44,:)=u6; fu(54,:)=u5; fu(64,:)=-2*u7; 
0016 flam=zeros(p.neq,p.nt);flam(2,:)=u5;