[Numpy-discussion] Numpy and iterative procedures

Nils Wagner nwagner@iam.uni-stuttgart...
Fri Feb 16 09:18:32 CST 2007


Nadav Horesh wrote:
> At first glance it doesn't look hard to, at least, avoid looping over i, by replacing [i] by [:-2], [i+1] by [1:-1] and [i+2] by [2:]. But I might be wrong. Can you submit the piece of code with at least the most internal loop?
>
>    Nadav.
>   
I guess he is looking for an implementation of

http://en.wikipedia.org/wiki/Successive_over-relaxation

Nils


> -----Original Message-----
> From:	numpy-discussion-bounces@scipy.org on behalf of Geoffrey Zhu
> Sent:	Thu 15-Feb-07 18:32
> To:	Discussion of Numerical Python
> Cc:	
> Subject:	Re: [Numpy-discussion] Numpy and iterative procedures
>
> Thanks Chuck.
>  
> I am trying to use Successive Over-relaxation to solve linear equations
> defined by M*v=q. 
>  
> There are several goals:
>  
> 1. Eventually (in production) I need it to be fast.
> 2. I am playing with the guts of the algorithm for now, to see how it
> works. that means i need some control for now.
> 3. Even in production, there is a chance i'd like to have the ability to
> tinker with the algorithm. 
>  
>
>   _____  
>
> From: numpy-discussion-bounces@scipy.org
> [mailto:numpy-discussion-bounces@scipy.org] On Behalf Of Charles R
> Harris
> Sent: Thursday, February 15, 2007 10:11 AM
> To: Discussion of Numerical Python
> Subject: Re: [Numpy-discussion] Numpy and iterative procedures
>
>
>
>
>
> On 2/15/07, Geoffrey Zhu <gzhu@peak6.com> wrote: 
>
> 	Hi,
> 	
> 	I am new to numpy. I'd like to know if it is possible to code
> efficient
> 	iterative procedures with numpy.
> 	
> 	Specifically, I have the following problem.
> 	
> 	M is an N*N matrix. Q is a N*1 vector. V is an N*1 vector I am
> trying to 
> 	find iteratively from the initial value V_0. The procedure is
> simply to
> 	calculate
> 	
> 	V_{n+1}[i]=3D1/M[I,i]*(q[i]-
> 	(M[i,1]*v_{n+1}[1]+M[I,2]*v_{n+1}[2]+..+M[i,i-1]*v_{n+1}[i-1]) -
> 	(M[I,i+1]*v_{n}[i+1]+M[I,i+2]*v_{n}[i+2]+..+M[I,N]*v_{n}[N])) 
> 	
> 	I do not see that this is something that can esaily be
> vectorized, is
> 	it?
>
>
> I think it would be better if you stated what the actual problem is. Is
> it a differential equation, for instance. That way we can determine what
> the problem class is and what algorithms are available to solve it. 
>
> Chuck
>
>
>
>
> _______________________________________________________
>
> The  information in this email or in any file attached
> hereto is intended only for the personal and confiden-
> tial  use  of  the individual or entity to which it is
> addressed and may contain information that is  propri-
> etary  and  confidential.  If you are not the intended
> recipient of this message you are hereby notified that
> any  review, dissemination, distribution or copying of
> this message is strictly prohibited.  This  communica-
> tion  is  for information purposes only and should not
> be regarded as an offer to sell or as  a  solicitation
> of an offer to buy any financial product. Email trans-
> mission cannot be guaranteed to be  secure  or  error-
> free. P6070214
>
>
>
>
>   
> ------------------------------------------------------------------------
>
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion@scipy.org
> http://projects.scipy.org/mailman/listinfo/numpy-discussion
>   
 


More information about the Numpy-discussion mailing list