diff --git a/src/Function1D.cpp b/src/Function1D.cpp
index a886db4..8a2368c 100644
--- a/src/Function1D.cpp
+++ b/src/Function1D.cpp
@@ -629,4 +629,7 @@ Matrix Aurora::reshape(const Matrix& aMatrix, int aRows, int aColumns, int aSlic
     }
     return Matrix::copyFromRawData(aMatrix.getData(),aRows,aColumns,aSlices);
 }
-
+void Aurora::nantoval(Matrix& aMatrix, double val2) {
+    Eigen::Map<Eigen::VectorXd> srcV(aMatrix.getData(),aMatrix.getDataSize());
+    srcV = srcV.array().isNaN().select(val2,srcV);
+}