From: Wolfgang Bangerth <bangerth@colostate.edu>
Date: Thu, 19 Jan 2017 01:12:05 +0000 (-0700)
Subject: Implement numbers::signaling_nan() also for Point<dim>.
X-Git-Tag: v8.5.0-rc1~222^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=9444398a57edb0ac5621423f4286eff2fda27f24;p=dealii.git

Implement numbers::signaling_nan() also for Point<dim>.
---

diff --git a/include/deal.II/base/signaling_nan.h b/include/deal.II/base/signaling_nan.h
index 6c9a42354f..5011dd5275 100644
--- a/include/deal.II/base/signaling_nan.h
+++ b/include/deal.II/base/signaling_nan.h
@@ -20,6 +20,7 @@
 #include <deal.II/base/tensor.h>
 #include <deal.II/base/symmetric_tensor.h>
 #include <deal.II/base/derivative_form.h>
+#include <deal.II/base/point.h>
 
 #include <limits>
 
@@ -122,6 +123,27 @@ namespace numbers
 
 
 
+      /**
+       * A specialization of the general NaNInitializer class that provides a
+       * function that returns a Tensor<rank,dim> value whose components are
+       * invalid signaling NaN values.
+       */
+      template <int dim, typename T>
+      struct NaNInitializer<Point<dim,T> >
+      {
+        static Point<dim,T> invalid_element ()
+        {
+          Point<dim,T> nan_point;
+
+          for (unsigned int i=0; i<dim; ++i)
+            nan_point[i] = NaNInitializer<T>::invalid_element();
+
+          return nan_point;
+        }
+      };
+
+
+
       /**
        * A specialization of the general NaNInitializer class that provides a
        * function that returns a SymmetricTensor<rank,dim> value whose