Update waves

Created by: acrovato

Context

This pull request brings the last developments from the flow module as well as some upgrade to waves.

Major changes in flow

• flow is now able to compute 3D flows
• The Jacobian (tangent) matrix for the Netwton method is now fully analytical
• The adjacent elements identification is now based on ```std::unordered_set````

Major changes in waves

• tbox and fwk are now included in a namespace
• waves is now fully included in a namespace
• All elements now have a private memory (Mem) which can be publicly accessed
• All elements now have a cache including Gauss points (Gauss) which can be publicly accessed
• MshDeformation now uses a linear elasticity solver
• MshCrack class implemented (creates duplicate internal boundaries in meshes)

Test

The ctest battery has been run on gaston, 100% tests passed under 110s.

Future work

• Use automated indentation rules for the whole code
• Move methods using VTK outside of MshData
• Clean tbox (mainly MshData, Tri3, Tera4 and Hex8)

mrstlnos/src/ElementMatrices.hpp

@@ -35,7 +35,7 @@ template<typename scalar, int element_type, int element_size> void ElementMatric
 
     Cache * cache;
 
-
+    /* @todo could probably be refactored into cache = element->getVCache(); */
     if (element_type==static_cast<int>(ELTYPE::HEX8))
     {
         CacheHex8 &cachehex8= this->trilinosHex8GetCache();
@@ -47,8 +47,8 @@ template<typename scalar, int element_type, int element_size> void ElementMatric
         cache = &cachetetra4;
     }
 
-    std::vector<Pt> gauss_p = cache->getVGauss().p; //(static_cast<CacheHex8*>(cache)->gauss).getP();
-    std::vector<double> gauss_w = cache->getVGauss().w; //(static_cast<CacheHex8*>(cache)->gauss).getW();
+    std::vector<Pt> gauss_p = cache->getVGauss().p;
+    std::vector<double> gauss_w = cache->getVGauss().w;
 
     tVector<scalar,element_size*4> tValues;