create sparse matrix from diagonal arrays

17 visualizzazioni (ultimi 30 giorni)
I have the 7 diagonal that come from a 3D finite volume discretization in separate arrays. These 7 arrays have 3D shape and correspond to each elements of the stencil for all points: top, bottom, east, west, north south and central point. I would like to create a sparse matrix object from these arrays, but I feel a little lost with the spdiag documentation. Help would be gratly apreaciated.

Risposta accettata

Bill Greene
Bill Greene il 14 Apr 2014
This doesn't exactly answer your question but still may be useful. This example:
shows how to create sparse finite difference matrices for the 1D, 2D, and 3D Poisson equation. For the 3D case, the 7-point stencil is used. For the 1D case, the use of spdiags is straightforward:
K1D = spdiags(ones(n,1)*[-1 2 -1],-1:1,n,n); % 1d Poisson matrix
Possibly you've already seen this simple example. But you'll notice that the kron function is used to create the 2D and 3D matrices. Use of spciags to create these is not nearly as straightforward as in the 1D case.
You might also find this blog useful:
As it explains, the key to efficiently creating sparse matrices in MATLAB is to pre-compute an array of row, column, value triplets and then use a single call to the sparse function to convert these to a sparse matrix.
Bill

Più risposte (2)

Andrew Newell
Andrew Newell il 10 Apr 2014
MATLAB supports sparse matrices, but for 3D arrays you'll need to visit the File Exchange and get N-dimensional sparse arrays by Matt J.
  2 Commenti
Jack Silver
Jack Silver il 11 Apr 2014
Sorry, but I think this is not exactly what I need. Maybe I was not clear enough (english is my third language), so lets try to be more precise :
I would like to solve a PDE using an implicit finite volume discretization. Since the problem is 3D, the matrix should be septadiagonal. I have computed the 7 diagonals (for top, bottom, east, west, north south and central point) in 7 different arrays.
Now, I would like to create a sparse matrix object from these arrays. I think that spdiag is the way to do, but I do not understand how it work
Andrew Newell
Andrew Newell il 11 Apr 2014
I'm still having trouble understanding what you mean. Could you attach a file with some data?

Accedi per commentare.


Jack Silver
Jack Silver il 24 Apr 2014
Sorry for being long to get back and for the lack of clarity. I need something equivalent to this python code

Categorie

Scopri di più su Sparse Matrices in Help Center e File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by