f_interpolation = scipy.interpolate.interp2d(*zip(*matrix(A)),kind='cubic')
plot_interpolation = contour_plot(lambda x,y:
f_interpolation(x,y)[0], (22.423,22.439),(37.07,37.08) ,cmap='jet', contours=numpy.arange(0,1400,100), colorbar=True)
plot_all = plot_interpolation
plot_all.show(,axes_labels=["m", "m"])
I am running this command to interpolate the data contained in matrix(A). matrix(A) is 1800 x 3 dense matrix over Real Double Field. Sage crashes everytime i do this both in a quad core Amd with 8GB of ram(kubuntu 12.04 64bit, in both the binary and compiled from source versions) and a netbook Lenovo s9e with Intel Atom and 1GB of ram(kubuntu 12.04 32bit, binary version of sage).It returns this message:
/home/user/Sage/sage-5.0.1/local/lib/libcsage.so(print_backtrace+0\
x31)[0x7fc18b464ea9]
/home/user/Sage/sage-5.0.1/local/lib/libcsage.so(sigdie+0x14)[0x7f\
c18b464edb]
/home/user/Sage/sage-5.0.1/local/lib/libcsage.so(sage_signal_handl\
er+0x20c)[0x7fc18b464b06]
/lib/x86_64-linux-gnu/libpthread.so.0(+0xfcb0)[0x7fc18f3f4cb0]
/home/user/Sage/sage-5.0.1/local/lib/python2.7/site-packages/scipy\
/interpolate/_fitpack.so(fprank_+0xdb0)[0x7fc15d29f4e0]
/home/user/Sage/sage-5.0.1/local/lib/python2.7/site-packages/scipy\
/interpolate/_fitpack.so(fpsurf_+0x67b8)[0x7fc15d2990d8]
/home/user/Sage/sage-5.0.1/local/lib/python2.7/site-packages/scipy\
/interpolate/_fitpack.so(surfit_+0xa08)[0x7fc15d289c18]
/home/user/Sage/sage-5.0.1/local/lib/python2.7/site-packages/scipy\
/interpolate/_fitpack.so(+0x4c7b)[0x7fc15d26cc7b]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x53a5)[0x7fc18f6f0e45]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x5420)[0x7fc18f6f0ec0]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0x73843)\
[0x7fc18f675843]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyObject_\
Call+0x53)[0x7fc18f64d7a3]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0x58f1f)\
[0x7fc18f65af1f]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyObject_\
Call+0x53)[0x7fc18f64d7a3]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0xac660)\
[0x7fc18f6ae660]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0xa8458)\
[0x7fc18f6aa458]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyObject_\
Call+0x53)[0x7fc18f64d7a3]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x2b2a)[0x7fc18f6ee5ca]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0x7374c)\
[0x7fc18f67574c]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyObject_\
Call+0x53)[0x7fc18f64d7a3]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x2b2a)[0x7fc18f6ee5ca]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x5420)[0x7fc18f6f0ec0]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x614b)[0x7fc18f6f1beb]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCode+0x32)[0x7fc18f6f2bd2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x515f)[0x7fc18f6f0bff]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCode+0x32)[0x7fc18f6f2bd2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyRun_Fil\
eExFlags+0xb0)[0x7fc18f715340]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0xe6cb8)\
[0x7fc18f6e8cb8]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x53a5)[0x7fc18f6f0e45]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCode+0x32)[0x7fc18f6f2bd2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x515f)[0x7fc18f6f0bff]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCode+0x32)[0x7fc18f6f2bd2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyRun_Fil\
eExFlags+0xb0)[0x7fc18f715340]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(+0xe6cb8)\
[0x7fc18f6e8cb8]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alFrameEx+0x53a5)[0x7fc18f6f0e45]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCodeEx+0x855)[0x7fc18f6f2a95]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyEval_Ev\
alCode+0x32)[0x7fc18f6f2bd2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyRun_Int\
eractiveOneFlags+0x182)[0x7fc18f7162c2]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyRun_Int\
eractiveLoopFlags+0x4e)[0x7fc18f7164fe]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(PyRun_Any\
FileExFlags+0x4c)[0x7fc18f716aac]
/home/user/Sage/sage-5.0.1/local/lib/libpython2.7.so.1.0(Py_Main+0\
xbe5)[0x7fc18f729295]
/lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed)[0x7fc18f04976d]
/home/user/Sage/sage-5.0.1/local/bin/python[0x400661]
------------------------------------------------------------------------\
Unhandled SIGSEGV: A segmentation fault occurred in Sage.
This probably occurred because a *compiled* component of Sage has a bug
in it and is not properly wrapped with sig_on(), sig_off(). You might
want to run Sage under gdb with 'sage -gdb' to debug this.
Sage will now terminate.
------------------------------------------------------------------------\
If matrix(a) is smaller(like 600 x 3) i do get results and a graph, but even then, the interpolation stops soon and produces a graph which makes no sense. It warns about:
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=3,3 nx,ny=11,15 m=72 fp=0.000114 s=0.000000
Any ideas?
