From 1a61f09b6a6698a039927c8bf5f3c5be4471b17a Mon Sep 17 00:00:00 2001
From: Alex Lindsay <alexlindsay239@gmail.com>
Date: Mon, 3 Dec 2018 10:36:29 -0700
Subject: [PATCH] Add copy assignment operators for T2 and DualNumber<T2,D2>

---
 src/numerics/include/metaphysicl/dualnumber.h | 22 +++++++++++++++++++
 .../include/metaphysicl/dualnumber_decl.h     |  6 +++++
 2 files changed, 28 insertions(+)

diff --git a/src/numerics/include/metaphysicl/dualnumber.h b/src/numerics/include/metaphysicl/dualnumber.h
index d3345944..b6067bd4 100644
--- a/src/numerics/include/metaphysicl/dualnumber.h
+++ b/src/numerics/include/metaphysicl/dualnumber.h
@@ -72,6 +72,17 @@ inline
 DualNumber<T,D>
 DualNumber<T,D>::operator! () const { return DualNumber<T,D>(!_val, !_deriv); }
 
+template <typename T, typename D>
+template <typename T2, typename D2>
+inline
+DualNumber<T,D> &
+DualNumber<T,D>::operator=(const DualNumber<T2,D2> & dn)
+{
+  _val = dn.value();
+  _deriv = dn.derivatives();
+  return *this;
+}
+
 template <typename T, typename D>
 template <typename T2, typename D2>
 inline
@@ -106,6 +117,17 @@ DualNumber<T,D>::operator=(const DualNumberSurrogate<T2,D2> & dns)
   return *this;
 }
 
+template <typename T, typename D>
+template <typename T2>
+inline
+DualNumber<T,D> &
+DualNumber<T,D>::operator=(const T2 & scalar)
+{
+  _val = scalar;
+  _deriv = 0;
+  return *this;
+}
+
 //
 // Member function definitions
 //
diff --git a/src/numerics/include/metaphysicl/dualnumber_decl.h b/src/numerics/include/metaphysicl/dualnumber_decl.h
index 2da2326e..d2fa49cb 100644
--- a/src/numerics/include/metaphysicl/dualnumber_decl.h
+++ b/src/numerics/include/metaphysicl/dualnumber_decl.h
@@ -74,6 +74,9 @@ class DualNumber : public safe_bool<DualNumber<T,D> >
   DualNumber& operator= (const DualNumber<T, D> & /*src*/) = default;
 #endif
 
+  template <typename T2, typename D2>
+  DualNumber & operator=(const DualNumber<T2,D2> & dn);
+
   template <typename T2, typename D2>
   DualNumber & operator=(const NotADuckDualNumber<T2,D2> & nd_dn);
 
@@ -83,6 +86,9 @@ class DualNumber : public safe_bool<DualNumber<T,D> >
   template <typename T2, typename D2>
   DualNumber & operator= (const DualNumberSurrogate<T2, D2> & dns);
 
+  template <typename T2>
+  DualNumber & operator= (const T2 & scalar);
+
   T& value();
 
   const T& value() const;