[yt-users] Error: cannot convert float infinity to integer
Jean-Claude Passy
jcpassy at gmail.com
Wed Apr 7 08:58:50 PDT 2010
Hi Matt,
nope, the error does not go away. :-(
JC
Matthew Turk a écrit :
> Hi Jean-Claude,
>
> Does the error go away if you supply x_bounds=(1e-30, 1.0) to the
> add_profile_sphere call? If so, perhaps we just need a better
> mechanism for guessing the center for calculating radius in the plot
> collection object.
>
> -Matt
>
>
>
> On Apr 6, 2010, at 4:13 PM, Jean-Claude Passy <jcpassy at gmail.com> wrote:
>
>
>> Dear all,
>>
>> I get the following error when trying to plot a radial profile of
>> the density. Does anyone have an idea how to solve this ?
>>
>> Thanks a lot,
>>
>> Jean-Claude
>>
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> #####################################################################
>>
>> In [42]: ls
>>
>> CommonEnvelope0001 CommonEnvelope0001.configure
>> CommonEnvelope0001.harrays CommonEnvelope0001.taskmap0000
>>
>> CommonEnvelope0001.boundary CommonEnvelope0001.cpu0000
>> CommonEnvelope0001.hierarchy CommonEnvelope0001.taskmap0001
>>
>> CommonEnvelope0001.boundary.hdf CommonEnvelope0001.cpu0001
>> CommonEnvelope0001.memorymap CommonEnvelope0001.yt
>>
>> In [43]: import yt.raven as R
>>
>> In [44]: import yt.lagos as L
>>
>> In [45]: filen = 'CommonEnvelope0001'
>>
>> In [46]: data = L.EnzoStaticOutput(filen)
>>
>> In [47]: plot_data = R.PlotCollection(data)
>>
>> yt INFO 2010-04-06 15:39:02,963 Created plot
>> collection with default plot-center = [0.4921875, 0.4921875,
>> 0.4921875]
>>
>> In [48]: plot_data.add_profile_sphere(0.5,'1',["Radius", "Density"])
>>
>> yt INFO 2010-04-06 15:39:56,695 Getting field Radius
>> from 1
>>
>> yt INFO 2010-04-06 15:39:56,696 Getting field x from 1
>>
>> yt INFO 2010-04-06 15:39:56,857 Getting field y from 1
>>
>> yt INFO 2010-04-06 15:39:56,919 Getting field z from 1
>>
>> yt INFO 2010-04-06 15:39:56,987 Getting field Density
>> from 1
>>
>> Out[48]: <yt.raven.PlotTypes.Profile1DPlot object at 0x3724d90>
>>
>> In [49]: plot_data.save('DensityProfile', format = 'png')
>>
>> ---
>> ---
>> ---------------------------------------------------------------------
>>
>> OverflowError Traceback (most recent
>> call last)
>>
>> /rpod2/jcpassy/Enzo/enzo/src/yt/scripts/iyt in <module>()
>>
>> ----> 1
>> 2
>> 3
>> 4
>> 5
>> /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotCollection.pyc in save
>> (self, basename, format, override, force_save)
>>
>> 79 fn.append(plot.save_image(basename, \
>>
>> 80 format=format, submit=self._run_id,
>>
>> ---> 81 override=override,
>> force_save=force_save))
>>
>> 82 if self.submit:
>>
>> 83 im = plot.im.copy()
>>
>> /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotTypes.pyc in save_image
>> (self, prefix, format, submit, override, force_save)
>>
>> 110 *override* will force no filename generation beyond
>> the prefix.
>>
>> 111 """
>>
>> --> 112 self._redraw_image()
>>
>> 113 if not override:
>>
>> 114 self._generate_prefix(prefix)
>>
>> /rpod2/jcpassy/Enzo/enzo/src/yt/yt/raven/PlotTypes.pyc in
>> _redraw_image(self)
>>
>> 711 func(self.data[self.fields[0]][indices],
>>
>> 712 self.data[self.fields[1]][indices],
>>
>> --> 713 **self.plot_options)
>>
>> 714 self._autoset_label(self.fields[0], self.set_x_label,
>> 'x')
>>
>> 715 self._autoset_label(self.fields[1], self.set_y_label,
>> 'y')
>>
>> /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/
>> axes.pyc in loglog(self, *args, **kwargs)
>>
>> 3566 b = self._hold
>>
>> 3567 self._hold = True # we've already processed the hold
>>
>> -> 3568 l = self.plot(*args, **kwargs)
>>
>> 3569 self._hold = b # restore the hold
>>
>> 3570
>> /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/
>> axes.pyc in plot(self, *args, **kwargs)
>>
>> 3435
>> 3436
>> -> 3437 self.autoscale_view(scalex=scalex, scaley=scaley)
>>
>> 3438 return lines
>>
>> 3439
>> /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/
>> axes.pyc in autoscale_view(self, tight, scalex, scaley)
>>
>> 1625 self.set_xbound(XL)
>>
>> 1626 if scaley and self._autoscaleYon:
>>
>> -> 1627 YL = self.yaxis.get_major_locator().view_limits
>> (y0, y1)
>>
>> 1628 self.set_ybound(YL)
>>
>> 1629
>> /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/
>> ticker.pyc in view_limits(self, vmin, vmax)
>>
>> 1117 vmin = minpos
>>
>> 1118
>> -> 1119 if not is_decade(vmin,self._base): vmin = decade_down
>> (vmin,self._base)
>>
>> 1120 if not is_decade(vmax,self._base): vmax = decade_up
>> (vmax,self._base)
>>
>> 1121
>> /rpod2/jcpassy/yt-x86_64/lib/python2.6/site-packages/matplotlib/
>> ticker.pyc in is_decade(x, base)
>>
>> 1028 def is_decade(x,base=10):
>>
>> 1029 lx = math.log(x)/math.log(base)
>>
>> -> 1030 return lx==int(lx)
>>
>> 1031
>> 1032 class LogLocator(Locator):
>>
>> OverflowError: cannot convert float infinity to integer
>>
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> ###
>> #####################################################################
>>
>> _______________________________________________
>> yt-users mailing list
>> yt-users at lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
More information about the yt-users
mailing list