On Thu, Sep 15, 2011 at 7:10 PM, kcrisman <kcris...@gmail.com> wrote: > > > On Sep 15, 9:08 am, Dan Drake <dr...@kaist.edu> wrote: >> This is strange: >> >> x, y =var('x y') >> plot3d(sqrt(x^2+y^2)*sin(1/sqrt(x^2+y^2)), (x,-1/2, 1/2), (y, -1/2, 1/2), >> adaptive=True) >> >> fails with "ValueError: cannot convert float NaN to integer". Something >> goes wrong when it partitions up the domain, probably when it looks at >> the origin. Is there a way to avoid this? Is this a bug? > > Or rather, very very close to the origin. I don't have time right > now, but can you confirm that there is not a depth of recursion > keyword like there is in 2d plotting? I would call this a bug. > > -- > To post to this group, send email to sage-support@googlegroups.com > To unsubscribe from this group, send email to > sage-support+unsubscr...@googlegroups.com > For more options, visit this group at > http://groups.google.com/group/sage-support > URL: http://www.sagemath.org >
For the time being there is a hack - x, y, x1, y1 =var('x y x1 y1') plot3d(lambda x,y: limit( limit(sqrt(x1^2+y1^2)*sin(1/sqrt(x1^2+y1^2)), x1=x), y1=y), (x,-1/2, 1/2), (y, -1/2, 1/2), adaptive=True) #slow or plot3d(lambda x,y: limit( limit(sqrt(x1^2+y1^2)*sin(1/sqrt(x1^2+y1^2)), x1=x), y1=y) if x^2+y^2<1e-4 else sqrt(x^2+y^2)*sin(1/sqrt(x^2+y^2)), (x,-1/2, 1/2), (y, -1/2, 1/2), adaptive=True) # dirty but fast -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org