0001 function p=rbconvinit(p)
0002
0003 p=stanparam(p); p.pstyle=2;
0004 pre=sprintf('%s',inputname(1)); p=setfn(p,pre);
0005 p.neq=3; p.f=@rbconvf; p.jac=@rbconvjac;
0006 p.bcf=@rbconvbc_noslip; p.outfu=@stanbra;
0007 lx=2; lz=0.5; p.geo=rec(lx,lz);
0008 p=stanmesh(p,100,25); p=setbmesh(p);
0009 p.xi=1/p.np;p.tau=zeros(p.neq*p.np+1,1); p.tau(1)=1;
0010 p.ds=0.1; p.dsmin=0.00001; p.dlammax=5; p.nsteps=10; p.neig=20;
0011 p.bisecmax=5; p.amod=0; p.pcomp=2; p.tol=1e-10;
0012 p.lam=740; p.lammax=900; p.u=0*ones(p.neq*p.np,1); p.errbound=1e-1;
0013