<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<META NAME="Generator" CONTENT="MS Exchange Server version 5.5.2658.34">
<TITLE>RE: [SciPy-user] interp1d question</TITLE>
</HEAD>
<BODY>

<P><FONT SIZE=2>Thanks Robert. I will update my version. However, x is an array with one dimension(Mx1). 'u' on the other hand is an array with MxN dimensions and I was wanting to interpolate column wise as opposed to row wise (that'll save me a couple of operation after reading the data in from a file)</FONT></P>

<P><FONT SIZE=2>Thanks,</FONT>
<BR><FONT SIZE=2>Jimmy</FONT>
</P>
<BR>

<P><FONT SIZE=2>PS. It works row wise...</FONT>
</P>
<BR>
<BR>

<P><FONT SIZE=2>-----Original Message-----</FONT>
<BR><FONT SIZE=2>From: scipy-user-bounces@scipy.org [<A HREF="mailto:scipy-user-bounces@scipy.org">mailto:scipy-user-bounces@scipy.org</A>] On Behalf Of Robert Kern</FONT>
<BR><FONT SIZE=2>Sent: Wednesday, October 18, 2006 4:35 PM</FONT>
<BR><FONT SIZE=2>To: SciPy Users List</FONT>
<BR><FONT SIZE=2>Subject: Re: [SciPy-user] interp1d question</FONT>
</P>

<P><FONT SIZE=2>Touma Jimmy E CTR USAF AFRL/MNGG wrote:</FONT>
<BR><FONT SIZE=2>&gt; Hi all,</FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt; I am trying to interpolate in 1d along a columns as opposed to rows</FONT>
<BR><FONT SIZE=2>&gt; and it doesn't seem to work. The code below emulates reading data</FONT>
<BR><FONT SIZE=2>&gt; where interpolation should be done column wise. The error when</FONT>
<BR><FONT SIZE=2>&gt; running this code follows. If you could please tell me what I am doing</FONT>
<BR><FONT SIZE=2>&gt; wrong and how to interpolate along columns.</FONT>
</P>

<P><FONT SIZE=2>&gt; /usr/lib/python2.4/site-packages/scipy/interpolate/interpolate.py in </FONT>
<BR><FONT SIZE=2>&gt; __init__(self, x, y, kind, axis, copy, bounds_error, fill_value)</FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 143&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; len_x,len_y = shape(oriented_x)[interp_axis], </FONT>
<BR><FONT SIZE=2>&gt; shape(oriented_y)[interp_axis]</FONT>
<BR><FONT SIZE=2>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 144&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if len_x != len_y:</FONT>
<BR><FONT SIZE=2>&gt; --&gt; 145&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; raise ValueError, &quot;x and y arrays must be equal in </FONT>
<BR><FONT SIZE=2>&gt; length along &quot;\</FONT>
<BR><FONT SIZE=2>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 146&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &quot;interpolation axis.&quot;</FONT>
<BR><FONT SIZE=2>&gt;&nbsp;&nbsp;&nbsp;&nbsp; 147&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if len_x &lt; 2 or len_y &lt; 2:</FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt; ValueError: x and y arrays must be equal in length along interpolation </FONT>
<BR><FONT SIZE=2>&gt; axis.</FONT>
<BR><FONT SIZE=2>&gt; WARNING: Failure executing file: &lt;d1.py&gt;</FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt; In [3]: numpy.__version__</FONT>
<BR><FONT SIZE=2>&gt; Out[3]: '1.0b5'</FONT>
<BR><FONT SIZE=2>&gt; </FONT>
<BR><FONT SIZE=2>&gt; In [4]: scipy.__version__</FONT>
<BR><FONT SIZE=2>&gt; Out[4]: '0.3.2'</FONT>
</P>

<P><FONT SIZE=2>You have a *very* old version of scipy. It uses Numeric, not numpy. Please upgrade.</FONT>
</P>

<P><FONT SIZE=2>In recent versions, interp1d does not accept x arrays with more than one </FONT>
<BR><FONT SIZE=2>dimension. I'm not entirely sure why the old version did; it makes no sense.</FONT>
</P>

<P><FONT SIZE=2>-- </FONT>
<BR><FONT SIZE=2>Robert Kern</FONT>
</P>

<P><FONT SIZE=2>&quot;I have come to believe that the whole world is an enigma, a harmless enigma</FONT>
<BR><FONT SIZE=2>&nbsp; that is made terrible by our own mad attempt to interpret it as though it had</FONT>
<BR><FONT SIZE=2>&nbsp; an underlying truth.&quot;</FONT>
<BR><FONT SIZE=2>&nbsp;&nbsp; -- Umberto Eco</FONT>
<BR><FONT SIZE=2>_______________________________________________</FONT>
<BR><FONT SIZE=2>SciPy-user mailing list</FONT>
<BR><FONT SIZE=2>SciPy-user@scipy.org</FONT>
<BR><FONT SIZE=2><A HREF="http://projects.scipy.org/mailman/listinfo/scipy-user" TARGET="_blank">http://projects.scipy.org/mailman/listinfo/scipy-user</A></FONT>
</P>

</BODY>
</HTML>