fixing av for ip

solvers/CompressibleFlowSolver/EquationSystems/NavierStokesCFE.h

@@ -388,6 +388,14 @@ namespace Nektar
         constexpr unsigned short nVarMax = 5;
         constexpr unsigned short nDimMax = 3;
 
+        // Update viscosity and thermal conductivity
+        // unfortunately this is very difficult to vectorize
+        // with the current implementation
+        Array<OneD, NekDouble> temperature (nPts, 0.0);
+        m_varConv->GetTemperature(inarray, temperature);
+        GetViscosityAndThermalCondFromTemp(temperature, m_mu,
+            m_thermalConductivity);
+
         // vector loop
         using namespace tinysimd;
         using vec_t = simd<NekDouble>;
@@ -429,10 +437,11 @@ namespace Nektar
             }
 
             // get temp
-            vec_t temperature = m_varConv->GetTemperature(inTmp.data());
+            // vec_t temperature = m_varConv->GetTemperature(inTmp.data());
             // get viscosity
             vec_t mu;
-            GetViscosityFromTempKernel(temperature, mu);
+            // GetViscosityFromTempKernel(temperature, mu);
+            mu.load(&(m_mu[p]), is_not_aligned);
 
             for (size_t nderiv = 0; nderiv < nDim; ++nderiv)
             {
@@ -519,10 +528,11 @@ namespace Nektar
             }
 
             // get temp
-            NekDouble temperature = m_varConv->GetTemperature(inTmp.data());
+            // NekDouble temperature = m_varConv->GetTemperature(inTmp.data());
             // get viscosity
             NekDouble mu;
-            GetViscosityFromTempKernel(temperature, mu);
+            // GetViscosityFromTempKernel(temperature, mu);
+            mu = m_mu[p];
 
             for (int nderiv = 0; nderiv < nDim; ++nderiv)
             {