11.07.2015 Views

Introduction to Sparse Matrices In Scilab - Projects

Introduction to Sparse Matrices In Scilab - Projects

Introduction to Sparse Matrices In Scilab - Projects

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

}{sciprint (" %+fi", pdblImg [k ]);}sciprint ("\n");k ++;}}LhsVar (1) = 0;return 0;9 The ARnoldi PACKageARPACK is a collection of Fortran77 subroutines designed <strong>to</strong> solve largescale eigenvalue problems. The functions available in <strong>Scilab</strong> are presented inthe figure 10.The package is designed <strong>to</strong> compute a few eigenvalues and correspondingeigenvec<strong>to</strong>rs of a general n by n matrix A. It is most appropriate for largesparse or structured matrices A where structured means that a matrix-vec<strong>to</strong>rproduct w = Av requires order n rather than the usual order n 2 floating poin<strong>to</strong>perations. This software is based upon an algorithmic variant of the Arnoldiprocess called the Implicitly Restarted Arnoldi Method (IRAM). When thematrix A is symmetric it reduces <strong>to</strong> a variant of the Lanczos process calledthe Implicitly Restarted Lanczos Method (IRLM). These variants may beviewed as a synthesis of the Arnoldi/Lanczos process with the ImplicitlyShifted QR technique that is suitable for large scale problems. For manystandard problems, a matrix fac<strong>to</strong>rization is not required. Only the action ofthe matrix on a vec<strong>to</strong>r is needed.ARPACK software is capable of solving large scale symmetric, nonsymmetric,and generalized eigenproblems from significant application areas. Thesoftware is designed <strong>to</strong> compute a few (k) eigenvalues with user specified featuressuch as those of largest real part or largest magnitude. S<strong>to</strong>rage requirementsare on the order of nk locations. No auxiliary s<strong>to</strong>rage is required. A se<strong>to</strong>f Schur basis vec<strong>to</strong>rs for the desired k-dimensional eigen-space is computedwhich is numerically orthogonal <strong>to</strong> working precision. Numerically accurateeigenvec<strong>to</strong>rs are available on request.The following is a list of the main features of the library :– Reverse Communication <strong>In</strong>terface.23

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!