← Back to team overview

dolfin team mailing list archive

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