We should also talk to Ondrej about this at SciPy.  Both sympy (through mpmath) and mpmath have matplotlib based function plotting.  I don't think it is adaptive, but I know mpmath can handle singularities.  Also, Ondrej is doing doing his graduate with with a group that does adaptive finite elements, so he would also be familiar with such algorithms.  I am sure that sympy and mpmath (Sage too as well) would be some of the main users of function plotting and would love to see this happen.

Cheers,

Brian

On Wed, Aug 12, 2009 at 4:28 AM, John Hunter <jdh2358@gmail.com> wrote:
We would like to add function plotting to mpl, but to do this right we
need to be able to adaptively sample a function evaluated over an
interval so that some tolerance condition is satisfied, perhaps with
both a relative and absolute error tolerance condition.  I am a bit
out of my area of competency here, eg I do not know exactly how the
tolerance condition should be specified, but I suspect some of you
here may be experts on this.  Does anyone have some code compatible
with the BSD license, preferably based on numpy but we would consider
an extension code or scipy solution, for doing this?

The functionality we have in mind is provided in matlab with fplot

 http://www.mathworks.com/access/helpdesk/help/techdoc/index.html?/access/helpdesk/help/techdoc/ref/fplot.html

We would like 1D and 2D versions of this ideally.  If anyone has some
suggestions, let me know.

Thanks,
JDH
_______________________________________________
NumPy-Discussion mailing list
NumPy-Discussion@scipy.org
http://mail.scipy.org/mailman/listinfo/numpy-discussion