diff --git a/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h b/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
index 6b9b5b42e58b705c5318e7d101c595a0a7d7d9af..6f91e8ec7d2ca2eea710759ffaaacb21e55f156e 100644
--- a/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
+++ b/ProcessLib/HeatTransportBHE/LocalAssemblers/HeatTransportBHELocalAssemblerSoil-impl.h
@@ -10,7 +10,6 @@
 
 #pragma once
 
-#include <valarray>
 #include <vector>
 
 #include "HeatTransportBHEProcessAssemblerInterface.h"
diff --git a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
index f79e1721477e228f3bea8b410c458b22611db324..1ffccc3d791a849803d4c68de6f854b66b98dfb7 100644
--- a/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
+++ b/ProcessLib/LIE/SmallDeformation/LocalAssembler/SmallDeformationLocalAssemblerMatrix-impl.h
@@ -11,7 +11,6 @@
 #pragma once
 
 #include <Eigen/Core>
-#include <valarray>
 #include <vector>
 
 #include "IntegrationPointDataMatrix.h"
diff --git a/Tests/NumLib/TestShapeFunctions.cpp b/Tests/NumLib/TestShapeFunctions.cpp
index e84f7709c38f803cf461f995861b1721b950c00b..1fe286985f85646b5c70cbf3240d35ae9dd1dff3 100644
--- a/Tests/NumLib/TestShapeFunctions.cpp
+++ b/Tests/NumLib/TestShapeFunctions.cpp
@@ -12,11 +12,11 @@
 
 #include <gtest/gtest.h>
 
+#include <Eigen/Core>
 #include <algorithm>
 #include <limits>
 #include <numeric>
 #include <type_traits>
-#include <valarray>
 
 #include "MeshLib/Elements/Elements.h"
 #include "NumLib/Fem/ShapeFunction/ShapeHex20.h"
@@ -156,13 +156,13 @@ TEST(NumLib, FemShapeQuad4)
     static const double eps = std::numeric_limits<double>::epsilon();
     static const unsigned NNodes = 4;
     static const unsigned dim = 2;
-    std::valarray<double> r(dim);
-    std::valarray<double> N(NNodes);
-    std::valarray<double> dN(NNodes * dim);
+    Eigen::VectorXd r(dim);
+    Eigen::VectorXd N(NNodes);
+    Eigen::VectorXd dN(NNodes * dim);
 
     // check N, dN at specific location
     {
-        r = .5;  // r = (0,5, 0.5)
+        r << 0.5, 0.5;
         ShapeQuad4::computeShapeFunction(r, N);
         ShapeQuad4::computeGradShapeFunction(r, dN);
         double exp_N[] = {0.5625, 0.1875, 0.0625, 0.1875};
@@ -172,13 +172,13 @@ TEST(NumLib, FemShapeQuad4)
         ASSERT_ARRAY_NEAR(exp_dN, dN, dN.size(), eps);
     }
 
-    std::valarray<double> exp_N(NNodes);
+    Eigen::VectorXd exp_N(NNodes);
     // check N_i(r_j)= {i==j: 1, i!=j: 0}
     for (unsigned i = 0; i < NNodes; i++)
     {
         r[0] = (i == 0 || i == 3) ? 1 : -1;
         r[1] = i < 2 ? 1 : -1;
-        exp_N = .0;
+        exp_N.setZero(NNodes);
         exp_N[i] = 1.0;
         ShapeQuad4::computeShapeFunction(r, N);
         ASSERT_ARRAY_NEAR(exp_N, N, NNodes, eps);
diff --git a/scripts/test/clang_sanitizer_blacklist.txt b/scripts/test/clang_sanitizer_blacklist.txt
index 4aa88c22164eee7a304732288ea8d6facf15bb08..07a89b304065f6f872c66d3412aedb950c25c36c 100644
--- a/scripts/test/clang_sanitizer_blacklist.txt
+++ b/scripts/test/clang_sanitizer_blacklist.txt
@@ -29,7 +29,3 @@ fun:_ZNKSs12find_last_ofEPKcmm
 fun:_ZNSt23mersenne_twister_engineImLm32ELm624ELm397ELm31ELm2567483615ELm11ELm4294967295ELm7ELm2636928640ELm15ELm4022730752ELm18ELm1812433253EE11_M_gen_randEv
 fun:_ZNSt24uniform_int_distributionIiEclISt23mersenne_twister_engineImLm32ELm624ELm397ELm31ELm2567483615ELm11ELm4294967295ELm7ELm2636928640ELm15ELm4022730752ELm18ELm1812433253EEEEiRT_
 fun:_ZNSt24uniform_int_distributionIiEclISt23mersenne_twister_engineImLm32ELm624ELm397ELm31ELm2567483615ELm11ELm4294967295ELm7ELm2636928640ELm15ELm4022730752ELm18ELm1812433253EEEEiRT_RKNS0_10param_typeE
-
-# valarray fill loop is counting down and reaching integer subtraction overflow
-# STL gcc-5.3: called from operator=(double)
-fun:_ZSt15__valarray_fillIdEvPT_mRKS0_