Path: csiph.com!newsfeed.hal-mli.net!feeder3.hal-mli.net!newsfeed.hal-mli.net!feeder1.hal-mli.net!newsreader4.netcologne.de!news.netcologne.de!newsfeed.freenet.ag!news2.euro.net!newsgate.cistron.nl!newsgate.news.xs4all.nl!post.news.xs4all.nl!not-for-mail
Return-Path: <winefrog@gmail.com>
X-Original-To: python-list@python.org
Delivered-To: python-list@mail.python.org
X-Spam-Status: OK 0.011
X-Spam-Evidence: '*H*': 0.98; '*S*': 0.00; 'scipy': 0.05; '__init__': 0.09; 'modules.': 0.09; 'skip:" 80': 0.09; 'to:addr:comp.lang.python': 0.09; 'cc:addr:python-list': 0.10; 'subject:error': 0.11; 'adjusted': 0.16; 'advice:': 0.16; 'benjamin': 0.16; 'cc:name:python list': 0.16; 'wrote:': 0.17; 'refers': 0.17; 'error.': 0.21; 'subject:skip:i 10': 0.22; 'wednesday,': 0.22; 'cc:addr:python.org': 0.25; 'header:In-Reply- To:1': 0.25; 'header:User-Agent:1': 0.26; 'cc:addr:gmail.com': 0.27; 'cc:2**2': 0.27; 'actual': 0.28; 'skip:_ 10': 0.29; 'error': 0.30; 'code': 0.31; 'file': 0.32; 'skip:s 30': 0.33; 'traceback': 0.33; "can't": 0.34; 'received:google.com': 0.34; 'received:209.85': 0.35; 'subject:with': 0.36; 'thank': 0.36; 'received:209': 0.37; 'subject:: ': 0.38; 'sure': 0.38; 'your': 0.60; 'kind': 0.61; 'skip:n 10': 0.63; 'dear': 0.66; '2013': 0.84; 'isaac': 0.84; 'oscar': 0.84; 'won': 0.96
X-Received: by 10.50.13.130 with SMTP id h2mr428993igc.16.1358953059856; Wed, 23 Jan 2013 06:57:39 -0800 (PST)
Newsgroups: comp.lang.python
Date: Wed, 23 Jan 2013 06:57:39 -0800 (PST)
In-Reply-To: <mailman.898.1358952056.2939.python-list@python.org>
Complaints-To: groups-abuse@google.com
Injection-Info: glegroupsg2000goo.googlegroups.com; posting-host=129.118.223.43; posting-account=DOyQBwoAAADIEpwWr89vk1x7RGFM33oc
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>
User-Agent: G2/1.0
X-Google-Web-Client: true
X-Google-IP: 129.118.223.43
MIME-Version: 1.0
Subject: Re: Memory error with quadratic interpolation
From: Isaac Won <winefrog@gmail.com>
To: comp.lang.python@googlegroups.com
Content-Type: text/plain; charset=ISO-8859-1
Cc: Isaac Won <winefrog@gmail.com>, Python List <python-list@python.org>, Ulrich Eckhardt <ulrich.eckhardt@dominolaser.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>
Message-ID: <mailman.901.1358953068.2939.python-list@python.org>
Lines: 61
NNTP-Posting-Host: 2001:888:2000:d::a6
X-Trace: 1358953068 news.xs4all.nl 6877 [2001:888:2000:d::a6]:40780
X-Complaints-To: abuse@xs4all.nl
Xref: csiph.com comp.lang.python:37479

On Wednesday, January 23, 2013 8:40:54 AM UTC-6, Oscar Benjamin wrote:
> On 23 January 2013 14:28, Isaac Won <winefrog@gmail.com> wrote:
> 
> > On Wednesday, January 23, 2013 4:08:13 AM UTC-6, Oscar Benjamin wrote:
> 
> >
> 
> > To Oscar
> 
> > My actual error message is:
> 
> > File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 311, in __init__
> 
> >      self._spline = splmake(x,oriented_y,order=order)
> 
> >   File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 809, in splmake
> 
> >      coefs = func(xk, yk, order, conds, B)
> 
> >   File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 530, in _find_smoothest
> 
> >      u,s,vh = np.dual.svd(B)
> 
> >   File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/linalg/decomp_svd.py", line 91, in svd
> 
> >      full_matrices=full_matrices, overwrite_a = overwrite_a)
> 
> > MemoryError
> 
> 
> 
> Are you sure that's the *whole* error message? The traceback only
> 
> refers to the scipy modules. I can't see the line from your code that
> 
> is generating the error.
> 
> 
> 
> 
> 
> Oscar

Dear Oscar,

Following is full error message after I adjusted following Ulich's advice:

interp = interp1d(indices[not_nan], x[not_nan], kind = 'quadratic') 
File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 311, in __init__
     self._spline = splmake(x,oriented_y,order=order) 
  File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 809, in splmake
     coefs = func(xk, yk, order, conds, B) 
  File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/interpolate/interpolate.py", line 530, in _find_smoothest
     u,s,vh = np.dual.svd(B) 
  File "/lustre/work/apps/python-2.7.1/lib/python2.7/site-packages/scipy/linalg/decomp_svd.py", line 91, in svd
     full_matrices=full_matrices, overwrite_a = overwrite_a) 
MemoryError 
----------------------------------------------------------------------
Thank you,

Hoonill