Hi Burak,
The getGrad operator doesn't include the boundary condition in its
evaluation. Reading <
http://www.ctcms.nist.gov/fipy/documentation/FAQ.html#i-can-t-get-the-fixedvalue-or-fixedflux-boundary-condition-objects-to-work-right-what-do-i-do-now>
may help. The next version of FiPy will use a
it Burak Atakan writes:
> print T.getGrad().getValue()
> print T.getValue()
> # END
>
> The output is:
> [[ -500. -1000. -1000. -999.9998
> -1000.0006 -999.9983 -1000.0035 -999.9949
> -1000.0049 -500.0161]]
> [ 395
Dear all,
I was wondering if there is an error in getGrad() within FiPy. For testing I
solved a really simple 1D differential equation ending up with a linear slope.
The value which I get from getGrad() are ok for all cells, except the ones at
the boundaries, where they are off by a factor of 2. (S