From: Jean-Paul Pelteret <jppelteret@gmail.com>
Date: Sat, 21 Jul 2018 14:41:04 +0000 (+0200)
Subject: Add AD number static asserts for FullMatrix and Vector classes.
X-Git-Tag: v9.1.0-rc1~899^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=99edea03f4e75fd72faf5a4d1b1286af9e3bc2e0;p=dealii.git

Add AD number static asserts for FullMatrix and Vector classes.
---

diff --git a/include/deal.II/lac/full_matrix.h b/include/deal.II/lac/full_matrix.h
index fc2cb0f76a..77600face8 100644
--- a/include/deal.II/lac/full_matrix.h
+++ b/include/deal.II/lac/full_matrix.h
@@ -23,6 +23,8 @@
 #include <deal.II/base/table.h>
 #include <deal.II/base/tensor.h>
 
+#include <deal.II/differentiation/ad/ad_number_traits.h>
+
 #include <deal.II/lac/exceptions.h>
 #include <deal.II/lac/identity_matrix.h>
 
@@ -68,6 +70,12 @@ template <typename number>
 class FullMatrix : public Table<2, number>
 {
 public:
+  // The assertion in full_matrix.templates.h for whether or not a number is
+  // finite is not compatible for AD number types.
+  static_assert(
+    !Differentiation::AD::is_ad_number<number>::value,
+    "The FullMatrix class does not support auto-differentiable numbers.");
+
   /**
    * A type of used to index into this container.
    */
diff --git a/include/deal.II/lac/vector.h b/include/deal.II/lac/vector.h
index 13c67eeaaa..a6dfa1cdc9 100644
--- a/include/deal.II/lac/vector.h
+++ b/include/deal.II/lac/vector.h
@@ -24,6 +24,8 @@
 #include <deal.II/base/logstream.h>
 #include <deal.II/base/subscriptor.h>
 
+#include <deal.II/differentiation/ad/ad_number_traits.h>
+
 #include <deal.II/lac/vector_operation.h>
 #include <deal.II/lac/vector_type_traits.h>
 
@@ -108,6 +110,12 @@ template <typename Number>
 class Vector : public Subscriptor
 {
 public:
+  // The assertion in vector.templates.h for whether or not a number is
+  // finite is not compatible for AD number types.
+  static_assert(
+    !Differentiation::AD::is_ad_number<Number>::value,
+    "The Vector class does not support auto-differentiable numbers.");
+
   /**
    * Declare standard types used in all containers. These types parallel those
    * in the <tt>C++</tt> standard libraries <tt>vector<...></tt> class.