Re: error in getGrad()?

2011-12-20 Thread Daniel Wheeler
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

Re: error in getGrad()?

2011-12-20 Thread Burak Atakan
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

error in getGrad()?

2011-12-20 Thread Burak Atakan
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