> Kent-Andre Mardal wrote: >> BTW: the Assembler code that computes A, b and enforces BC >> simultaneously in a symmetric way, > > Very nice. We can now implement some symmetric linear solvers.
I'm interested in a cholesky solver so I'll see what I can do with CHOLMOD. It is part of the same suite as UMFPACK so it seems a sensible choice. /Dag > > Do you keep the 'Dirichlet' dofs and zero the row and column, or do you > eliminate them from the global system? > > is now faster than the standard >> non-symmetric way of doing it: >> >> timings: >> new : 1.55 >> standard: 1.94 >> > > Which backend is this using? It could be saying more about the > efficiency of the function ident(..), which I suspect is slow for > uBlasMatrix, than the assembly process. > > Garth > >> check demo/pde/sym-dirichlet-bc/speed-up-test.py >> >> Currently only cell integrals and exterior facet integrals are >> implemented. >> >> Kent >> >> On on., 2008-08-13 at 11:46 +0200, DOLFIN wrote: >>> One or more new changesets pushed to the primary dolfin repository. >>> A short summary of the last three changesets is included below. >>> >>> changeset: 4563:3fbd8bd0822d1f363521ccd29afc1f694205ecc5 >>> tag: tip >>> parent: 4562:d912b2383bc98f98be1fa61e25cac3dec72df8fb >>> parent: 4554:b3a0c04c85c2231e87bf187b35bcf6e2e27e8247 >>> user: Kent-Andre Mardal <[EMAIL PROTECTED]> >>> date: Wed Aug 13 11:41:52 2008 +0200 >>> files: dolfin/fem/Assembler.cpp dolfin/fem/Assembler.h >>> description: >>> merge >>> >>> >>> changeset: 4562:d912b2383bc98f98be1fa61e25cac3dec72df8fb >>> user: Kent-Andre Mardal <[EMAIL PROTECTED]> >>> date: Wed Aug 13 11:41:23 2008 +0200 >>> files: demo/pde/sym-dirichlet-bc/demo.py >>> demo/pde/sym-dirichlet-bc/speed-up-test.py >>> description: >>> update text >>> >>> >>> changeset: 4561:4749e1c883b9e21767a2c1142541badb6461866c >>> parent: 4560:5daf9606ed74d52cb52ccf4e2f3cf20ab7152277 >>> parent: 4549:463016a27aadf2d3994d2249b4059feaf264dd5d >>> user: Kent-Andre Mardal <[EMAIL PROTECTED]> >>> date: Wed Aug 13 11:34:33 2008 +0200 >>> files: dolfin/fem/Assembler.cpp dolfin/la/MTL4SparsityPattern.cpp >>> dolfin/la/MTL4SparsityPattern.h dolfin/mesh/ALE.cpp dolfin/mesh/ALE.h >>> dolfin/mesh/ALEType.h dolfin/swig/dolfin_headers.i >>> description: >>> merge >>> >>> ---------------------------------------------------------------------- >>> For more details, visit http://www.fenics.org/hg/dolfin >>> _______________________________________________ >>> DOLFIN-dev mailing list >>> [email protected] >>> http://www.fenics.org/mailman/listinfo/dolfin-dev >> >> _______________________________________________ >> DOLFIN-dev mailing list >> [email protected] >> http://www.fenics.org/mailman/listinfo/dolfin-dev > > _______________________________________________ > DOLFIN-dev mailing list > [email protected] > http://www.fenics.org/mailman/listinfo/dolfin-dev > _______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
