[SciPy-User] Roots of a signal

[Nils Wagner]

Strange.

python -i test_spline.py
/home/nwagner/local/lib64/python2.6/site-packages/scikits/__init__.py:1: 
UserWarning: Module pkg_resources was already imported 
from 
/home/nwagner/local/lib64/python2.6/site-packages/setuptools-0.6c9-py2.6.egg/pkg_resources.py, 
but 
/home/nwagner/local/lib64/python2.6/site-packages/distribute-0.6.4-py2.6.egg 
is being added to sys.path
   __import__('pkg_resources').declare_namespace(__name__)
/home/nwagner/local/lib64/python2.6/site-packages/scikits/__init__.py:1: 
UserWarning: Module site was already imported from 
/usr/lib64/python2.6/site.pyc, but 
/home/nwagner/local/lib64/python2.6/site-packages/distribute-0.6.4-py2.6.egg 
is being added to sys.path
   __import__('pkg_resources').declare_namespace(__name__)
/home/nwagner/local/lib64/python2.6/site-packages/matplotlib/backends/backend_gtk.py:621: 
DeprecationWarning: Use the new widget gtk.Tooltip
   self.tooltips = gtk.Tooltips()
Traceback (most recent call last):
   File "test_spline.py", line 11, in <module>
     roots = sproot(tck,10)
   File 
"/home/nwagner/local/lib64/python2.6/site-packages/scipy/interpolate/fitpack.py", 
line 576, in sproot
     z,ier=_fitpack._sproot(t,c,k,mest)
MemoryError

>>> scipy.__version__
'0.8.0.dev6212'
>>> numpy.__version__
'1.5.0.dev8087'
>>> show_config()                                                                                                  
atlas_threads_info:
     libraries = ['lapack', 'ptf77blas', 'ptcblas', 
'atlas']
     library_dirs = 
['/home/nwagner/src/ATLAS3.8.2/mybuild/lib']
     language = f77
     include_dirs = 
['/home/nwagner/src/ATLAS3.8.2/include']
blas_opt_info:
     libraries = ['ptf77blas', 'ptcblas', 'atlas']
     library_dirs = 
['/home/nwagner/src/ATLAS3.8.2/mybuild/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.8.2\\""')]
     language = c
     include_dirs = 
['/home/nwagner/src/ATLAS3.8.2/include']
atlas_blas_threads_info:
     libraries = ['ptf77blas', 'ptcblas', 'atlas']
     library_dirs = 
['/home/nwagner/src/ATLAS3.8.2/mybuild/lib']
     language = c
     include_dirs = 
['/home/nwagner/src/ATLAS3.8.2/include']
lapack_opt_info:
     libraries = ['lapack', 'ptf77blas', 'ptcblas', 
'atlas']
     library_dirs = 
['/home/nwagner/src/ATLAS3.8.2/mybuild/lib']
     define_macros = [('ATLAS_INFO', '"\\"3.8.2\\""')]
     language = f77
     include_dirs = 
['/home/nwagner/src/ATLAS3.8.2/include']
lapack_mkl_info:
   NOT AVAILABLE
blas_mkl_info:
   NOT AVAILABLE
mkl_info:
   NOT AVAILABLE


cat /proc/meminfo
MemTotal:        1793140 kB

> 
> your script finishes without problem for me
> 
>>>> roots
> array([  1.32135351e-18,   5.01642799e-01, 
>  9.99383967e-01,
>         1.50051340e+00,   2.00014353e+00, 
>  2.49953905e+00,
>         3.00076843e+00,   3.49870336e+00])
> 
> from the plot (of only the first half) and
> 
>>>> np.diff(y2>0).sum()
> 9
> 
> it looks like there is one root missing
> 
> Josef
>