Thread Previous • Date Previous • Date Next • Thread Next |
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 >
Thread Previous • Date Previous • Date Next • Thread Next |