Re: problem with the octave built-in function "kron"

[Stephen Montgomery-Smith]

Use speye instead of eye.


On 05/03/2013 10:02 AM, scigeek wrote:
> I am trying to run the "kron" function in octave, but as soon as I run it, it
> gives the following line of error.
> 
> error: memory exhausted or requested size is too large for octave range
> type--
> 
> here is the piece of code that I am running:
> =================================
> rx = 800, ry =571
> 
> x = linspace(0,1,rx);
> y = linspace(0,1,ry);
> dx = x(2)-x(1); dy = y(2)-y(1);
> ex = ones(rx,1); ey=ones(ry,1);
> Ix = eye(rx); Iy = eye(ry);
> Dx= (spdiags([ex -2*ex ex],[-1 0 1],rx,rx))/dx^2; 
> Dy= (spdiags([ey -2*ey ey],[-1 0 1],ry,ry))/dy^2;
> L = kron(Iy,Dx)+kron(Dy,Ix);
> ================================
> 
> Can any body help ?
> Thanks
> 
> 
> 
> --
> View this message in context: 
> http://octave.1599824.n4.nabble.com/problem-with-the-octave-built-in-function-kron-tp4652601.html
> Sent from the Octave - General mailing list archive at Nabble.com.
> _______________________________________________
> Help-octave mailing list
> address@hidden
> https://mailman.cae.wisc.edu/listinfo/help-octave
> 
>