dolfin team mailing list archive

[dkouroun@xxxxxxxxx]

> This can be fixed (using Triangle or TetGen for the generation of the
> mesh). If you're familiar with any any of these, feel free to
> contribute. I can give you suggestions.
> 

I am familiar with both! I am all ears!

> > FemSpace Vh(Th, P1), Uh(Th, P1); //P1 langrangian elements
> > 
> > Problem WeakFormulation( grad(Vh)*grad(Uh) - f*Vh ); //for a Poisson
> problem
> > or even
> > Problem WeakFormulation( Vh.dx*Uh.dx+Vh.dy*Uh.dy - f*Vh ); //
> 
> This will happen in the future. In fact, this part will be moved outside
> of DOLFIN (but still look as if it were part of DOLFIN). The evaluation
> of variational forms will be moved to other subprojects within the
> FEniCS project.
> 
> Sure, you will have all that and more. But understand that DOLFIN (and
> FEniCS) is a research project which means that our first priority is to
> get things right.
> 
So do you have any suggestions on how to implement this! I had
expression-templates in mind but I would like to use something simpler. Do you 
have any idea how should I design the element classes and that class which reads
the weak formulation? It is very important to me since if implemented 
once it will make my code dimension independed and will save lots of hours of
debugging! 

Best Wishes!
Drosos.