# [SciPy-user] Getting coordinates of a level (contour) curve

Johann Cohen-Tanugi cohen@slac.stanford....
Thu Aug 28 17:57:33 CDT 2008

Hi Rob,
I just tried your code after installing PyDSTool from svn. But I get the
following error :
Velocity around curve is always 1, e.g. look at 100th point
norm(Point(sol[100].labels['EP']['data'].V))
=---------------------------------------------------------------------------
NotImplementedError                       Traceback (most recent call last)

/home/cohen/PyCont_LevelCurve.py in <module>()
59 print "\nVelocity around curve is always 1, e.g. look at 100th
point"
60 print "norm(Point(sol[100].labels['EP']['data'].V)) =", \
---> 61       norm(Point(sol[100].labels['EP']['data'].V))
62
63 print "... at which we have travelled distance ds =", \

/usr/lib/python2.5/site-packages/matplotlib/mlab.pyc in norm(x, y)
1783     Deprecated - see numpy.linalg.norm
1784     """
-> 1785     raise NotImplementedError('Deprecated - see numpy.linalg.norm')
1786
1787

NotImplementedError: Deprecated - see numpy.linalg.norm
WARNING: Failure executing file: <PyCont_LevelCurve.py>

My version of numpy is  '1.2.0.dev5694'.
htanks,
Johann

Rob Clewley wrote:
> Attached is a commented example in PyCont for a 2D zero level set that
> defines an ellipse. It's very easy! I'll add this to PyDSTool's
> examples. Thanks for giving me the impetus to do this. Let me know
> what you think.
>
> -Rob
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> SciPy-user mailing list
> SciPy-user@scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-user
>