


(local) jacobian for AC with global coupling lam<u> in jacobian ignored here and dealt with in (b)gclss.m


0001 function [c,fu,flam,b]=acgcjac(p,u,lam) 0002 % (local) jacobian for AC with global coupling 0003 % lam<u> in jacobian ignored here and dealt with in (b)gclss.m 0004 global eta; % set in acinit resp. acf after remeshing 0005 um=eta*u; u=pdeintrp(p.points,p.tria,u); 0006 c=p.d; flam=um; b=0; fu=p.c1+2*p.c2*u+3*p.c3*u.^2+5*p.c5*u.^4;