dolfin team mailing list archive
-
dolfin team
-
Mailing list archive
-
Message #00116
Re: Expression templates
The same thing will be possible with code generation. In the same way
as you create a PDE today in DOLFIN,
Poisson poisson(f);
you will be to define functions (coefficients) and constants.
/Anders
On Sun, Oct 24, 2004 at 02:54:48PM -0700, Hugo Duncan wrote:
> On Sun, 24 Oct 2004 <dkouroun@xxxxxxxxx> wrote:
>
> > The good thing is
> > that either way, one is finally able to write a variational formulation
> > like:
> >
> > Variable<Simplex, P2, _3D> u,v;
> >
> > f = u*v + k*grad(u)*grad(v); //where k is constant.
>
>
> Drosos,
>
> Good to hear that you have something in progress. What I would like
> to see is something like:
>
> Variable u,v;
> Constant k;
> Function f=u*v + k*grad(u)*grad(v); //where k is constant.
> Form<Simplex, tetrahedron<ThreeD>, order<5> > form(f);
>
> This would allow you to build libraries of functions, and to try out
> different elements.
>
> I don't know how it play's with dolfin's solvers, but being able to
> perturb a system (eg. change k, above) and re-solve using the previous
> solution as a starting point would be of interest. This
> is one reason to prefer an expression template approach over
> an external generator.
>
> Hugo
>
>
>
> __________________________________
> Do you Yahoo!?
> Yahoo! Mail Address AutoComplete - You start. We finish.
> http://promotions.yahoo.com/new_mail
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/cgi-bin/mailman/listinfo/dolfin-dev
>
References