dolfin team mailing list archive

[Hugo Duncan <hugo_duncan@xxxxxxxxx>]

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