On Wed, Feb 25, 2009 at 09:34:44PM +0000, A Navaei wrote: > Is the linear algebra support in dolfin capable of automatically > generating the sparse linear system of equations? (I know that ffc > effectively does this, but note that this question does not involve > solving a PDE using FEM, but a linear system of equations). > > For example, assume we have some busy sets of linear equations and we > wish to convert them to Ax = b so that we can feed them to the linear > solvers. Is there any way of automatically obtaining A and b avoiding > manual calculations and hard-coding?
Yes! Isn't that what you are using already? Define a bilinear form a, a linear form L, and then just assemble the matrix and vector by A = assemble(a) b = assemble(L) -- Anders
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
