bd output for vector GP
0001 function out=vgpbra(p,u,lam) 0002 % bd output for vector GP 0003 rmax1=max(u(1:p.np)); imax1=max(u(p.np+1:2*p.np)); 0004 rmax2=max(u(2*p.np+1:3*p.np)); imax2=max(u(3*p.np+1:4*p.np)); 0005 out=[rmax1; imax1; imax1/rmax1; rmax2; imax2; imax2/rmax2; ... 0006 triint(u(1:p.np).^2,p.points,p.tria)];