high order smoother class added : high order meshes are guaranteed to be
valid : some elements are "un-curved" in order to ensure good value of
the distorsion mapping. P2 triangle quality estimate is robust

conservation laws can now take a function rho for a factor in the mass matrix term M = int_ rho Ni Nj
simply add a function * massFactor in your conservation law to use it :-)

- fixed MED node ordering
- polished 64 bit Cocoa Mac OS X FLTK version
- updated copyright

dgshell v 1.0
lua --> python

No commit message

order 7->9 seem to be wrong : the way shape functions are
computed is maybe numerically unstable
 

removed so benchmarks has to update)
Add function on fullMatrix to assemble matrix by perturbations

Integration on bulk and virtual interface terms are performed by a
matrix vector product. Add 2 functions in fullMatrix.h and cpp to
realize this operation

the size is correct (so that they can be used on proxies)

now