
Release 0.1 
  Appetizer
  Lagrange elements of general order on triangles 
  Lagrange elements of general order on tetrahedra 
  Integration on general triangles 
  Integration on general tetrahedra 
  Poisson element matrices on general triangles 
  Poisson element matrices on general tetrahedra 
  Simple regression testing 
  Barycentric coordinates  
  Bezier coordinates  
  Degrees of freedom 
  Document the current features  
  Clean up

Release 0.2 
  Bernstein polynomials 
  RaviartThomas and Nedelec elements
  Stokes element matrices on general triangles/tetrahedra
  autoconf regression testing 
  Python support by using SWIG and SWIGiNaC   
  Simple unit square/box mesh for academic testing of new methods
  Global matrices and dofs on the academic mesh 
  Make the example with the nonlinear PDE and the computed Jacobian clean
  Implement classes that for vectors and sums of expressions  
  Document the current features  

Release 0.3 
  Code generation of element basis functions and  element matrices 
  Optimization of code generation based on tree traversal
  More complicated element methods, like the newly introduced 
    elements for the mixed formulation of elasticity problems
  Interface to mesh packages, e.g., NetGen. 
  Interface to linear algebra packages, e.g., Trilinos. 
  Document the current features  



Random comments: 

Slicing of lst ? 

Symbol factory is needed. Have just wasted several hours!  

Should be possible to construct polygons without a prefix!  (ok)   

Should the order of the nodes in the Polygon constructor matter or
should the Polygon's determine the order themselves. In the latter
case a mapping must be stored.


The Poisson matrix maker should probably be a class ? 
This will make sure that doxygen generates documentation

Need to check the integration on the reference element in 3D. 

Check the code that is supposed to fix the comparison of ex. 

Check if it is possible to take integrate*2 and then eval_integ*2 (it is with tree traversal)  

The test bed is now a mix between a sandbox and actual tests. 
This needs to be dealt with. 

What about Python


Design:  
      1. What should be the basic vector structure ? 
         lst and exvector does not provide +-/* etc. 
	 matrix does and is more general, but slightly
	 inelegant 

      2. Error handling







