From 9550af5292cfba5487de02dde583876ae31c0890 Mon Sep 17 00:00:00 2001
From: Douglas Serson <d.serson@gmail.com>
Date: Thu, 8 Sep 2016 12:14:14 +0100
Subject: [PATCH] Add missing const

(cherry picked from commit aa9bcb1d97474377b966a57e769afa714cd3baa8)
---
 .../EquationSystems/VCSMapping.cpp                   | 12 ++++++------
 .../EquationSystems/VCSMapping.h                     | 10 +++++-----
 2 files changed, 11 insertions(+), 11 deletions(-)

diff --git a/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.cpp b/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.cpp
index 80d9c94aa8..bb0be7c16f 100644
--- a/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.cpp
+++ b/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
diff --git a/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.h b/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.h
index e5cc870b13..319e15f3d9 100644
--- a/solvers/IncNavierStokesSolver/EquationSystems/VCSMapping.h
+++ b/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);
         
     };
-- 
GitLab