[yt-users] Treecode question

[Matthew Turk]

Hi Andrew,

It should be changed.  In theory, they'll usually be the same in most
cases, but it is a typo and it's fixed.

(Also, the treecode was all Stephen Skory -- I agree with your
sentiment, it's a huge benefit, and I'm very appreciative of it too!)

-Matt

On Fri, May 27, 2011 at 10:16 AM, Andrew Davis <andrew.davis at yale.edu> wrote:
> Hi folks,
>
> I just updated my yt yesterday to use the new treecode algorithm in calculating
> binding energies.  Many many thanks to those who worked on it and tested it.
> However,  I'm looking at the code, and I see this in the function (on line 388
> of derived_quantities.py):
>
>        octree = Octree(cover_ActiveDimensions, 1, True)
>        #print 'here', cover_ActiveDimensions
>        # Now discover what levels this data comes from, not assuming
>        # symmetry.
>        dxes = na.unique(data['dx']) # unique returns a sorted array,
>        dyes = na.unique(data['dy']) # so these will all have the same
> --->>>>  dzes = na.unique(data['dx']) # order.
>        # We only need one dim to figure out levels, we'll use x.
>        dx = 1./data.pf.domain_dimensions[0]
>
> I've marked the line where I would think it should be getting data['dz'] instead
> of 'dx'.  Am I just not getting this or should this be changed?
>
> Thanks!
>
> Andrew
> *************************************************************
> **  Andrew J. Davis             andrew.davis at yale.edu      **
> **  Dept. of Astronomy          203-432-5119               **
> **  Yale University             www.astro.yale.edu/adavis  **
> *************************************************************
>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>