Flux of PDE solution
[cgxu,cgyu] = pdecgrad(p,t,c,u) [cgxu,cgyu] = pdecgrad(p,t,c,u,time) [cgxu,cgyu] = pdecgrad(p,t,c,u,time,sdl)
Row i of cgxu contains
Row i of cgyu contains
There is one column for each triangle in t in both cgxu and cgyu.
The geometry of the PDE problem is given by the mesh data p and t. Details on the mesh data representation can be found in the entry on initmesh.
The format for the solution vector u is described in assempde.
The scalar optional argument time is used for parabolic and hyperbolic problems, if c depends on t, the time.
The optional argument sdl restricts the computation to the subdomains in the list sdl.