resi

PURPOSE

compute residual for -div(c grad u)+au-f-b grad u=0,

SYNOPSIS

function r=resi(p,u,lam)

DESCRIPTION

 compute residual for -div(c grad u)+au-f-b grad u=0,

CROSS-REFERENCE INFORMATION

This function calls:

assemadv assemble K for advective terms.

This function is called by:

getGlam calc Glam
getder calc Gu, Glam depending on p.jsw;
inistep initial step in cont; yields first 2 points on branch and first tangent
jaccheck compare numerical and assembled Jacobian; useful to check coding, timing,
nloopext newton loop (with damping parameter alpha) for extended system
nlooppde newton-loop (with damping parameter alpha) for pde alone,
resinj residual, called in numerical G_u via numjac, dummy interface
sgndetA compute sign(detA) from data in p
specGu computes and plots neig eigenvalues nearest to origin in figure p.ifig
swibra load data from file pre/fname and attempt branch switch,
tint time integration

SOURCE CODE

0001 function r=resi(p,u,lam)
0002 % compute residual for -div(c grad u)+au-f-b grad u=0,
0003 [c,a,f,b]=p.f(p,u,lam); bc=p.bcf(p,u,lam); 
0004 [K,F]=assempde(bc,p.points,p.edges,p.tria,c,a,f); 
0005 if(any(b)); Kadv=assemadv(p.points,p.tria,b); K=K-Kadv; end 
0006 r=K*u-F;