[petsc-users] kronecker products

[Jed Brown]

On Thu, Nov 4, 2010 at 10:39, Benjamin Sanderse <B.Sanderse at cwi.nl> wrote:

> I am using a finite volume method on structured grids and things are really
> simple.
> To give you an idea, here is the matrix that calculates du/dx, given a
> vector with u. For simplicity I use the Matlab way to construct such a
> thing:
>
> diagonal       = (1/deltax)*ones(Nx,1);
> Dux_1D         = spdiags([-diagonal diagonal],[0 1],Nx-1,Nx); %
> differencing matrix in 1D
> Dux_2D         = kron(speye(Ny),Dux_1D));
> Dux_3D         = kron(speye(Nz),Dux_2D));
>

Just form the matrix, this isn't a good way to build it in parallel.

Jed
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101104/5048d22b/attachment.htm>