[Yt-dev] HaloProfiler Problems

[Britton Smith]

I've just committed a fix for this that should take care of cases like
this.  The problem turned out to be that r_max for the halo was actually
smaller than r_min, so the profiling itself failed.  I've cut the minimum
radius by a factor of 2 to be a little more generous, and also now require
that r_max / r_min be at least 2 as well.  There is also an addition check
to catch halos that make it through to the profiling but create no valid
data.  I've tested both of these and they seem to work.  Let me know if this
comes up again.

Britton

On Mon, Jan 26, 2009 at 10:50 AM, Stephen Skory <stephenskory at yahoo.com>wrote:

> All,
>
> +        virial = {}
> +        if (len(profile[overdensity_
> field]) < 2):
> +            for field in fields:
> +                virial[field] = 0.0
> +            return virial
>
>
> This isn't going to fix it. 'print profile[overdensity_field]' gives this
> for the bad halo:
>
>
> [ NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
>  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
>  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
>  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN
>  NaN  NaN  NaN  NaN]
>
> Again, sorry for not following up more quickly on this.
>
>
>  _______________________________________________________
> sskory at physics.ucsd.edu           o__  Stephen Skory
> http://physics.ucsd.edu/~sskory/ <http://physics.ucsd.edu/%7Esskory/> _.>/
> _Graduate Student
> ________________________________(_)_\(_)_______________
> _______________________________________________
> Yt-dev mailing list
> Yt-dev at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-dev-spacepope.org/attachments/20090126/e416eb5d/attachment.htm>