Add missing const

(cherry picked from commit aa9bcb1d)

solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.cpp

@@ -739,8 +739,8 @@ namespace Nektar
      * Explicit terms of the mapping
      */
     void   VCSMapping::ApplyIncNSMappingForcing(
-        const Array<OneD, Array<OneD, NekDouble> >    &inarray,
-        Array<OneD, Array<OneD, NekDouble> >          &outarray)
+        const Array<OneD, const Array<OneD, NekDouble> > &inarray,
+        Array<OneD, Array<OneD, NekDouble> >             &outarray)
     {
         int physTot = m_fields[0]->GetTotPoints();
         Array<OneD, Array<OneD, NekDouble> >       vel(m_nConvectiveFields);
@@ -822,7 +822,7 @@ namespace Nektar
     }
     
         void VCSMapping::MappingAdvectionCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray)
         {
             int physTot = m_fields[0]->GetTotPoints();
@@ -847,8 +847,8 @@ namespace Nektar
         }
         
         void VCSMapping::MappingAccelerationCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &vel,
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &vel,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray)
         {
             int physTot = m_fields[0]->GetTotPoints();
@@ -953,7 +953,7 @@ namespace Nektar
         }
 
         void VCSMapping::MappingViscousCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray)
         {
             // L(U) - 1.0*d^2(u^i)/dx^jdx^j

solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.h

@@ -64,7 +64,7 @@ namespace Nektar
         
         // 
         void ApplyIncNSMappingForcing (
-                const Array<OneD, Array<OneD, NekDouble> >        &inarray,
+                const Array<OneD, const Array<OneD, NekDouble> >  &inarray,
                 Array<OneD, Array<OneD, NekDouble> >              &outarray);
 
         virtual ~VCSMapping();
@@ -124,7 +124,7 @@ namespace Nektar
         // Correction needed for convective terms = N(u) - ( -(u \nabla) u)
         //     velPhys is the velocity field (transformed for physical space)
         void MappingAdvectionCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray);
 
         // Correction needed for time-derivative terms 
@@ -132,8 +132,8 @@ namespace Nektar
         //     vel     is the velocity field (can be in wavespace)
         //     velPhys is the velocity field (transformed for physical space)
         void MappingAccelerationCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &vel,
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &vel,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray);
 
         // Correction needed for pressure terms   
@@ -147,7 +147,7 @@ namespace Nektar
         //     vel     is the velocity field (can be in wavespace)
         //     velPhys is the velocity field (transformed for physical space)
         void MappingViscousCorrection(
-            const Array<OneD, Array<OneD, NekDouble> >        &velPhys,
+            const Array<OneD, const Array<OneD, NekDouble> >  &velPhys,
             Array<OneD, Array<OneD, NekDouble> >              &outarray);
         
     };