Path: csiph.com!newsfeed.hal-mli.net!feeder3.hal-mli.net!newsfeed.hal-mli.net!feeder1.hal-mli.net!newsfeed.xs4all.nl!newsfeed3.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail
Return-Path: <oscar.j.benjamin@gmail.com>
X-Original-To: python-list@python.org
Delivered-To: python-list@mail.python.org
X-Spam-Status: OK 0.028
X-Spam-Evidence: '*H*': 0.94; '*S*': 0.00; 'elements.': 0.05; 'indices': 0.07; 'cc:addr:googlegroups.com': 0.09; 'cc:addr :python-list': 0.10; 'subject:error': 0.11; 'benjamin': 0.16; 'cc:name:python list': 0.16; 'numpy': 0.16; 'wrote:': 0.17; 'import': 0.21; 'supposed': 0.21; 'subject:skip:i 10': 0.22; 'wednesday,': 0.22; "i've": 0.23; 'cc:no real name:2**0': 0.24; 'cc:addr:python.org': 0.25; 'header:In-Reply-To:1': 0.25; 'used,': 0.27; 'cc:2**2': 0.27; 'message-id:@mail.gmail.com': 0.27; "doesn't": 0.28; 'arrays': 0.29; 'function': 0.30; 'received:google.com': 0.34; 'received:209.85': 0.35; 'really': 0.36; 'skip:u 20': 0.36; 'but': 0.36; 'subject:with': 0.36; 'received:209': 0.37; 'subject:: ': 0.38; 'where': 0.40; 'skip:u 10': 0.60; 'skip:n 10': 0.63; 'more': 0.63; '2013': 0.84; 'isaac': 0.84; 'oscar': 0.84; 'hundred': 0.95; 'won': 0.96
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:x-received:in-reply-to:references:date:message-id :subject:from:to:cc:content-type; bh=9LokI4fuY/ioDKkVvDZPRY27Goz32mOLYI9okTxECV8=; b=zg9ktYEkMDQHqo9OwoS7tjO0lPGS+DoW0sL4kRMBN6nmSU0XPLX5XLOtuKun0LzcYP +3CLU7DxMHRSLP0infeFYFqRBzS7mAJCVH2kvh2gkVSwgUv/4WYtd4l7z+7NHBV2xM3w aMMzGruSK9F1XSZFdUqwmD/5JzeeWDBN95ng9fDg4IZE2EcZayIiDIRew4ZIhbTXvfe9 14K1k1gYj0sPmgHoPKT0KlgM44CTQaRhQG5X192EU8TqwRk6dEvwiauiPM9tgxFNUbmm TcEBA0YU2gVKDXxD7ZBW979DbWa8JdKwHd5JORAViP2MIl4X0qq97tKjNNjL1BdRcYyH wC3A==
MIME-Version: 1.0
X-Received: by 10.112.37.200 with SMTP id a8mr578296lbk.92.1359023455622; Thu, 24 Jan 2013 02:30:55 -0800 (PST)
In-Reply-To: <94ef456a-af84-48fe-8007-df42c4a625c4@googlegroups.com>
References: <0909ee56-c304-4b9d-8be9-97c7a3261d27@googlegroups.com> <j4l4t9-skt.ln1@satorlaser.homedns.org> <mailman.880.1358935694.2939.python-list@python.org> <8daf86cf-386c-44ae-83d5-25fb04572c00@googlegroups.com> <mailman.898.1358952056.2939.python-list@python.org> <a5662c8c-67c3-47e6-a8b7-2eddbc272657@googlegroups.com> <mailman.910.1358959906.2939.python-list@python.org> <94ef456a-af84-48fe-8007-df42c4a625c4@googlegroups.com>
Date: Thu, 24 Jan 2013 10:30:55 +0000
Subject: Re: Memory error with quadratic interpolation
From: Oscar Benjamin <oscar.j.benjamin@gmail.com>
To: Isaac Won <winefrog@gmail.com>
Content-Type: text/plain; charset=ISO-8859-1
Cc: Python List <python-list@python.org>, Ulrich Eckhardt <ulrich.eckhardt@dominolaser.com>, comp.lang.python@googlegroups.com
X-BeenThere: python-list@python.org
X-Mailman-Version: 2.1.15
Precedence: list
List-Id: General discussion list for the Python programming language <python-list.python.org>
List-Unsubscribe: <http://mail.python.org/mailman/options/python-list>, <mailto:python-list-request@python.org?subject=unsubscribe>
List-Archive: <http://mail.python.org/pipermail/python-list/>
List-Post: <mailto:python-list@python.org>
List-Help: <mailto:python-list-request@python.org?subject=help>
List-Subscribe: <http://mail.python.org/mailman/listinfo/python-list>, <mailto:python-list-request@python.org?subject=subscribe>
Newsgroups: comp.lang.python
Message-ID: <mailman.952.1359023458.2939.python-list@python.org>
Lines: 23
NNTP-Posting-Host: 2001:888:2000:d::a6
X-Trace: 1359023458 news.xs4all.nl 6964 [2001:888:2000:d::a6]:37434
X-Complaints-To: abuse@xs4all.nl
Xref: csiph.com comp.lang.python:37550

On 23 January 2013 17:33, Isaac Won <winefrog@gmail.com> wrote:
> On Wednesday, January 23, 2013 10:51:43 AM UTC-6, Oscar Benjamin wrote:
>> On 23 January 2013 14:57, Isaac Won <winefrog@gmail.com> wrote:
>>
>> > On Wednesday, January 23, 2013 8:40:54 AM UTC-6, Oscar Benjamin wrote:
>>
>> Unless I've misunderstood how this function is supposed to be used, it
>> just doesn't really seem to work for arrays of much more than a few
>> hundred elements.
>>

The solution is to use UnivariateSpline. I don't know what the
difference is but it works where the other fails:

import numpy as np
from scipy.interpolate import UnivariateSpline
x = np.array(10000 * [0.0], float)
indices = np.arange(len(x))
interp = UnivariateSpline(indices, x, k=2)
print(interp(1.5))


Oscar