Dear Dan and Raymond, Thanks for your help to both. After your suggestions I review my code. But changing X and Y like
X = mesh.x Y = mesh.y makes the code seriously slow. Since in myFunc I'll make calculations based on t and x,y index and after that I'll have a a matrix of calculated values with size (nx/dx)*(ny*dy) I set matrix directly to source term with setValue function. I hope this approach is true for calculations. Last form of the code can be seen at: http://pastie.org/6636947 Cheers, Serbulent 2013/3/19 Daniel Wheeler <[email protected]>: > On Tue, Mar 19, 2013 at 4:35 PM, Raymond Smith <[email protected]> wrote: >> >> Oh, that's much cleaner. Thanks for that, Dan. > > > No problem. FYI, there is an open ticket on this issue > http://matforge.org/fipy/ticket/447. > > -- > Daniel Wheeler > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
