dolfin team mailing list archive

[Matthew Knepley <knepley@xxxxxxxxx>]

On Mon, Apr 6, 2009 at 1:09 AM, Anders Logg <logg@xxxxxxxxx> wrote:

> On Sun, Apr 05, 2009 at 05:11:51PM -0500, Matthew Knepley wrote:
> > On Sun, Apr 5, 2009 at 2:54 PM, Anders Logg <logg@xxxxxxxxx> wrote:
> >
> >     Dirichlet BC would need to be added in/after each multiplication and
> >     it should be possible to build it into the mult() operator and make
> it
> >     efficient.
> >
> >
> > I still think the best way to handle this is to eliminate them, as I
> talked
> > about
> > last time at Simula.
> >
> >   Matt
>
> I think we're in position to do that now with the new FunctionSpace
> class. One could allow restrictions to be imposed, for example
>
>  V.restrict(bc);
>
> But I haven't seen any good numbers to indicate it's worth the effort.
> How much better is it to eliminate, condition numbers etc for letting
> the constrained dofs just sit along?


Its not about performance, so much as

  1) Programming ease, ESPECIALLY with nonlinear problems

  2) Separation of storage organization from traversal organization

  3) Cache performance

  4) Interaction with external packages

I think I must have done a terrible job explaining this. I have always coded
both ways and always found that elimination is better.

   Matt


>
> --
> Anders
>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.9 (GNU/Linux)
>
> iEYEARECAAYFAknZnKMACgkQTuwUCDsYZdGTUACcDkA5TX7WO8pJDCDX1LHLPrUD
> UEIAoJgYQ9J9dkrFbhdgPWPQy1uKxkVM
> =Lwul
> -----END PGP SIGNATURE-----
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>
>


-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener