[yt-users] Error in making light ray for a gadget dataset

[Britton Smith]

Hi Bili,

My apologies for the long delay in replying to this thread.  By now, I
think you may be working with Trident.  If things are working, that's
great.  However, Trident relies on much of the machinery that you were
trying to use, so it's still possible that you see this error.  The problem
is likely related to how the LightRay generator is calculating the
trajectory for the ray and not with the underlying ray machinery which it
relies on.  If it's possible for you to give me access to this dataset, I
would be happy to take a look at this.

Britton

On Mon, Feb 15, 2016 at 9:30 PM, Nathan Goldbaum <nathan12343 at gmail.com>
wrote:

> I think Cameron, Britton, and Devin have done this for their Trident
> project - you might want to contact them.
>
> On Mon, Feb 15, 2016 at 3:23 PM, Bili Dong - Gmail <qobilidop at gmail.com>
> wrote:
>
>> Hi Matt,
>>
>> I've tried to construct the ray and ortho_ray objects, neither of which
>> can reproduce the error. They work just fine. I haven't succeeded in
>> tracing down the error yet, but I'll keep working on that.
>>
>> Besides, I'm wondering if there are examples of successfully making light
>> rays in a particle dataset? That may help figuring out what might go wrong.
>>
>> Thanks,
>> Bili
>>
>> On Fri, Feb 12, 2016 at 12:42 PM, Matthew Turk <matthewturk at gmail.com>
>> wrote:
>>
>>> Hi Bili,
>>>
>>> Hm, I'm honestly not sure.  This looks like an oddity in the way the
>>> octree is being constructed in spatial dimensions.  Any chance you can
>>> reproduce this with a ray or ortho_ray object, which should simplify the
>>> number of moving parts?
>>>
>>> On Tue, Feb 9, 2016 at 2:04 PM, Bili Dong - Gmail <qobilidop at gmail.com>
>>> wrote:
>>>
>>>> Hi yt users,
>>>>
>>>> I'm trying to make a light ray for a gadget dataset using this script:
>>>> http://paste.yt-project.org/show/6234/
>>>>
>>>> But I get the following error:
>>>>
>>>> ```
>>>>
>>>> /Users/qobilidop/repos/yt/yt/geometry/geometry_handler.py in fcoords(self)    320             c = obj.select_fcoords(self.dobj)    321             if c.shape[0] == 0: continue--> 322             ci[ind:ind+c.shape[0], :] = c    323             ind += c.shape[0]    324         return ci
>>>> ValueError: could not broadcast input array from shape (81,3) into shape (40,3)
>>>>
>>>> ```
>>>>
>>>>
>>>> My yt is on the changeset d47150dfbde6.
>>>>
>>>>
>>>> I'm wondering if I'm doing it the correct way? And how could I solve the problem?
>>>>
>>>>
>>>> Thanks,
>>>>
>>>> Bili
>>>>
>>>>
>>>> _______________________________________________
>>>> yt-users mailing list
>>>> yt-users at lists.spacepope.org
>>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>>>
>>>>
>>>
>>> _______________________________________________
>>> yt-users mailing list
>>> yt-users at lists.spacepope.org
>>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>>
>>>
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>>
>
> _______________________________________________
> yt-users mailing list
> yt-users at lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20160224/cc3f3b5c/attachment.htm>