Introduction to Sparse Matrices In Scilab - Projects

More documents

Recommendations

Info

-->A = scibench_poissonA ( 500 );--> size (A)ans =250000. 250000.-->nnz (A)ans =1248000.11.5 ConclusionScilab 5 can manage sparse matrices and solve partial differential equationssuch as the Poisson equation for example. Indeed, Scilab provides severalsparse linear equation solvers, including a sparse backslash operator,iterative methods (e.g. the pre-conditionned conjugate gradient algorithm)and other direct solvers such as the UMFPACK or TAUCS modules. Withthese tools we can solve huge systems of linear equations, because Scilab onlystores the nonzero entries of these massive matrices. The solvers may fail ifthe memory required to store the matrix is beyond the capacity of Scilab,but this happens only for huge matrices. In this case, the future Scilab 6 mayhelp to overcome this limitation, by removing the use of the stack which isused by Scilab 5 (see [5] for more details on this topic). Another limitation ofScilab is that there is currently no preconditionner for sparse linear systemsof equations. This limitation should be removed once the spilu module isready for a release.Références[1] R. Barrett, M. Berry, T. F. Chan, J. Demmel, J. Donato, J. Dongarra,V. Eijkhout, R. Pozo, C. Romine, and H. Van der Vorst. Templates forthe Solution of Linear Systems : Building Blocks for Iterative Methods,2nd Edition. SIAM, Philadelphia, PA, 1994.[2] François Delebecque, Carlos Klimann, and Serge Steer. Basile, un syst ?de cao pour l’automatique, version 3.0 : guide d’utilisation. 1989.[3] Kenneth Kundert and Alberto Sangiovanni-Vincentelli. Sparse 1.3.http://www.netlib.org/sparse/.34
[4] Kenneth S. Kundert. Sparse matrix techniques. In Albert E. Ruehli, editor,Circuit Analysis, Simulation and Design, volume 3 of pt. 1. North-Holland, 1986.[5] Consortium Scilab DIGITEO Michael Baudin. Programming inscilab. http://forge.scilab.org/index.php/p/docprogscilab/downloads/, 2011.[6] The Scilab Consortium Michael Baudin. Solving poisson pde with sparsematrices. http://wiki.scilab.org, 2011.[7] Bruno Pincon and Karim Ramdani. Scilab tools for pde : Applicationto time-reversal. http://www.iecn.u-nancy.fr/~pincon/scilab/expose_taiwan.pdf, 2004.[8] Héctor E. Rubio Scola. Implementation of lipsol in scilab. 1997. http://hal.inria.fr/inria-00069956/en.[9] Neeraj Sharma and Matthias K. Gobbert. A comparativeevaluation of matlab, octave, freemat, and scilab for researchand teaching. http://userpages.umbc.edu/~gobbert/papers/SharmaGobbertTR2010.pdf, 2010.[10] Wikipedia. Sparse matrix — wikipedia, the free encyclopedia, 2010.[Online ; accessed 11-June-2010].35
Page 2: 4 Cholesky factorizations 155 Funct
Page 6 and 7: sp2adjspeyesponessprandspzerosfulls
Page 8 and 9: -->D= sparse (C)D =( 3, 5) sparse m
Page 10 and 11: lufactlusolvelugetludelchfactchsolv
Page 12 and 13: -->x= lusolve (h,b)x =1.0.50.333333
Page 14 and 15: 0 0 0 0 28 20 71 39 0 00 10 0 0 32
Page 16 and 17: PlotSparseReadHBSparsecond2spcondes
Page 18 and 19: umf licenseumf ludelumf lufactumf l
Page 20 and 21: Read sparse doubles matrices in gat
Page 22 and 23: int read_sparse ( char * fname , un
Page 24 and 25: dnaupddneupddsaupddsaupdznaupdzneup
Page 26 and 27: Information about MatrixMarket file
Page 28 and 29: Sparse matrices can be managed in S
Page 30 and 31: The solver argument is designed so
Page 32 and 33: function u = mysolverTAUCS (N,b)A =

Introduction to Sparse Matrices In Scilab - Projects

Create successful ePaper yourself

Delete template?

Save as template?