GroupOfElement in FunctionSpace can be modified -- Reordering

FunctionSpace/FunctionSpace.cpp

@@ -42,7 +42,7 @@ FunctionSpace::~FunctionSpace(void){
     delete eToGod;
 }
 
-void FunctionSpace::build(const GroupOfElement& goe,
+void FunctionSpace::build(GroupOfElement& goe,
                           const Basis& basis){
 
   // Save GroupOfElement & Mesh //

FunctionSpace/FunctionSpace.h

@@ -40,8 +40,8 @@
 class FunctionSpace{
  protected:
   // Geometry //
-  const Mesh*           mesh;
-  const GroupOfElement* goe;
+  const Mesh*     mesh;
+  GroupOfElement* goe;
 
   // Basis //
   std::vector<const Basis*>* basis;
@@ -68,8 +68,8 @@ class FunctionSpace{
   const Basis&                     getBasis(unsigned int i) const;
   unsigned int                     getNBasis(void) const;
 
-  const GroupOfElement& getSupport(void) const;
-  bool                  isScalar(void) const;
+  GroupOfElement& getSupport(void) const;
+  bool            isScalar(void) const;
 
   std::vector<Dof> getKeys(const MElement& element) const;
   std::vector<Dof> getKeys(const MVertex& vertex) const;
@@ -88,7 +88,7 @@ class FunctionSpace{
   // Init
   FunctionSpace(void);
 
-  void build(const GroupOfElement& goe,
+  void build(GroupOfElement& goe,
 	     const Basis& basis);
 
   // Dof
@@ -187,7 +187,7 @@ inline unsigned int FunctionSpace::getNBasis(void) const{
   return nBasis;
 }
 
-inline const GroupOfElement& FunctionSpace::getSupport(void) const{
+inline GroupOfElement& FunctionSpace::getSupport(void) const{
   return *goe;
 }
 

FunctionSpace/FunctionSpaceScalar.cpp

 #include "FunctionSpaceScalar.h"
 
-FunctionSpaceScalar::FunctionSpaceScalar(const GroupOfElement& goe,
+FunctionSpaceScalar::FunctionSpaceScalar(GroupOfElement& goe,
                                          const Basis& basis){
   scalar = true;
   build(goe, basis);

FunctionSpace/FunctionSpaceScalar.h

@@ -18,7 +18,7 @@
 
 class FunctionSpaceScalar : public FunctionSpace{
  public:
-  FunctionSpaceScalar(const GroupOfElement& goe, const Basis& basis);
+  FunctionSpaceScalar(GroupOfElement& goe, const Basis& basis);
   virtual ~FunctionSpaceScalar(void);
 
   double

FunctionSpace/FunctionSpaceVector.cpp

 #include "Mapper.h"
 #include "FunctionSpaceVector.h"
 
-FunctionSpaceVector::FunctionSpaceVector(const GroupOfElement& goe,
+FunctionSpaceVector::FunctionSpaceVector(GroupOfElement& goe,
                                          const Basis& basis){
   scalar = false;
   build(goe, basis);

FunctionSpace/FunctionSpaceVector.h

@@ -19,7 +19,7 @@
 
 class FunctionSpaceVector : public FunctionSpace{
  public:
-  FunctionSpaceVector(const GroupOfElement& goe, const Basis& basis);
+  FunctionSpaceVector(GroupOfElement& goe, const Basis& basis);
   virtual ~FunctionSpaceVector(void);
 
   fullVector<double>