[SciPy-User] [SciPy-user] ValueError: arrays must have same number of dimensions

[Jorge E. ´Sanchez Sanchez]

Hi,

it seems to me that you don't know that loops in python end before the last index value of the loop,
you just need to do something like:
 Np1 = N+1
for j in range (1, Np1):

HTH
Best Regards
Jorge


> Date: Thu, 28 Jun 2012 08:32:47 -0700
> From: surfcast23 at gmail.com
> To: scipy-user at scipy.org
> Subject: [SciPy-User] [SciPy-user] ValueError: arrays must have same number	of dimensions
> 
> 
> Hi I am trying to translate a Matlab code from Trefethen's Spectral Methods
> in MATLAB to Python. I am running into a problem with numpy.concatenate
> wanting the arrays to have the same number of dimensions.
> Here is part of the Matlab code that I am rewriting 
> 
> 
> N = 512; h = 2*pi/N; x = h*(1:N); t = 0; dt = h/4;
> a = .1;
> c = a + sin (x-1).^2;
> v = exp(-100*(x-1).^2); vold = exp(-100*(x-a*dt-1).^2);
> 
> column = [0 .5*(-1).^(1:N-1).*cot((1:N-1)*h/2)];
> D = toeplitz(column,-column);
> 
> % Time-stepping by leap frog formula:
> tmax = 15; tplot = .15; clf, drawnow, set(gcf,'renderer','zbuffer')
> plotgap = round(tplot/dt); dt = tplot/plotgap;
> nplots = round(tmax/tplot);
> data = [v; zeros(nplots,N)]; tdata = t;
> 
> 
> What I am trying in Python
> 
> for j in range(1,N):
>  tmax = 15;   
>  tplot = .15;
>  h = 2*pi/N;
>  dt = h/4; 
>  nplots = around(tmax/tplot);
>  plotgap = around(tplot/dt);
>  dt = tplot/plotgap;  
>  x = h*arange(1,N);
>  c = a+np.sin(x-1)**2;
>  v = exp(-100*(x-1)**2);
>  vold = exp(-100*(x-a*dt-1)**2);
>  
>  column = ((0.5*(-1)**arange(1,N-1))*cotdg(arange(1,N-1))*(h/2));
>  D = toeplitz(column,-column);
>  
>  k = np.zeros(((nplots,N)));
>  data = np.concatenate((v, k),axis = 0);
>  tdata = t;
> 
> 
> I have highlighted the vectors in question.
> v.shape gives (511,)
> 
> k.shape gives (100, 512)
> 
> I would appreciate any help with this. Thank you!
> 
> -- 
> View this message in context: http://old.nabble.com/ValueError%3A-arrays-must-have-same-number-of-dimensions-tp34086886p34086886.html
> Sent from the Scipy-User mailing list archive at Nabble.com.
> 
> _______________________________________________
> SciPy-User mailing list
> SciPy-User at scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20120628/3cd383bf/attachment.html>