Re: System of partial, differential and algebraic equations

[Jonathan Guyer]

On Dec 17, 2007, at 12:43 PM, I wrote:

You can also get the gradient at the face between cells, which happens to sit exactly at r = R, with:

    >>> y.getFaceGrad()[rF == 4][0]
    array([ 8.])

Oops. Make that

    >>> rF = mesh.getFaceCenters()[...,0]
    >>> y.getFaceGrad()[rF == R][0]
    array([ 8.])