Hi folks, Plot3d() look amazing! I really like the way lighting is used to because it help understand the shape of the surface.
I have a data file that I want to plot using plot3d, so I import the data using numpy.loadtxt. plot3d() only plots functions so I try to interpolate the dataset using scipy.interpolate.interp2d() and then pass the interpolator as a lambda function to plot3d(). import numpy as np from scipy.interpolate import interp2d x = np.loadtxt("x.txt") y = np.loadtxt("y.txt") z = np.loadtxt("z.txt") data_i = interp_surface(x,y,z,kind='linear',copy=False) # Fails here with the error below plot3d(lambda X,Y: data_i(X,Y), x, y) Unfortunately, the scipy cannot create the interpolation object even though the data is really smooth. I give the error at the bottom of this message. Sage is using scipy version 0.7. I have Scipy version 0.8.0 installed and the interp2d() call works perfectly! So I guess this is a problem with Scipy 0.7. So my question is how can I plot my data!? Is there away of using plot3d 'outside' of Sage? I noticed that it is some how making a call to JMol (or PyMol), would you know how I could do this? Regards, Daniel # # Error message # Warning: No more knots can be added because the number of B-spline coefficients already exceeds the number of data points m. Probably causes: either s or m too small. (fp>s) kx,ky=1,1 nx,ny=21,21 m=360 fp=0.000426 s=0.000000 -- 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