diff --git a/src/Matrix.cpp b/src/Matrix.cpp
index 7de7d2e..bfbfbc2 100644
--- a/src/Matrix.cpp
+++ b/src/Matrix.cpp
@@ -252,6 +252,18 @@ namespace Aurora {
         return !mData || mInfo.empty();
     }
 
+    bool Matrix::isNan() const
+    {
+        for(size_t i=0; i<getDataSize(); ++i)
+        {
+            if(mData.get()[i] != NAN)
+            {
+                return false;
+            }
+        }
+        return true;
+    }
+
     bool Matrix::isScalar() const {
         return (getDimSize(0) == 1 &&
                 getDimSize(1) == 1 &&
diff --git a/src/Matrix.h b/src/Matrix.h
index dded537..2f678b9 100644
--- a/src/Matrix.h
+++ b/src/Matrix.h
@@ -178,6 +178,8 @@ namespace Aurora {
          */
         bool isNull() const;
 
+        bool isNan() const;
+
         bool isVector() const;
 
         /**