12.07.2015 Views

DOLFIN User Manual - FEniCS Project

DOLFIN User Manual - FEniCS Project

DOLFIN User Manual - FEniCS Project

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 3Linear algebra◮ Developer’s note: This chapter is currently being written. . .<strong>DOLFIN</strong> does not have its own linear algebra routines, but instead uses theexternal package PETSc [11] for linear algebra functionality.PETSc is a suite of data structures and routines for the scalable (parallel)solution of scientific applications modeled by partial differential equations.It employs the MPI standard for all message-passing communication.For convenience <strong>DOLFIN</strong> provides wrappers for some of the most commonlinear algebra functionality. For more advanced usage, <strong>DOLFIN</strong> providesdirect access to the PETSc pointers to be used with the standard PETScinterface.3.1 Matrices and vectorsThe matrix class Matrix provides wrappers for initializing a sequential sparsematrix, or a sequential sparse matrix in block compressed row format. Forparallel matrices the PETSc interface has to be used directly.The code for initializing a sequential sparse M × N matrix takes the form:25

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

Saved successfully!

Ooh no, something went wrong!