tint

 time integration 
 do nt semi-implicit (Euler)steps for plot each pmod-th step

0001 function p=tint(p,dt,nt,pmod)
0002 % time integration
0003 % do nt semi-implicit (Euler)steps for plot each pmod-th step
0004 fprintf('running tint, see figure %g\n',p.ifig); 
0005 a=1;c=0;f=zeros(p.neq,1); % assemble multi-component mass-matrix bM
0006 bc=p.bcf(p,p.u,p.lam); 
0007 [K,bM,F,Q,G,H,R]=assempde(bc,p.points,p.edges,p.tria,c,a,f); 
0008 nc=0; 
0009 while(nc<nt) 
0010   [cc,aa,ff,b]=p.f(p,p.u,p.lam); p.bc=p.bcf(p,p.u,p.lam); 
0011   [K,F]=assempde(bc,p.points,p.edges,p.tria,cc,aa,ff);
0012   if(any(b)) Kadv=assemadv(p.points,p.tria,b); K=K-Kadv; end 
0013   L=bM+dt*K; p.u=L\(bM*p.u+dt*F); nc=nc+1;
0014   if(mod(nc,pmod)==0); plotsol(p,p.ifig,p.pcmp,p.pstyle); 
0015       r=norm(resi(p,p.u,p.lam),p.normsw); fprintf('res=%g\n',r); 
0016       drawnow; end 
0017 end