commit 88062b7feddcf3e8354c6857753523c2a80d864f
author Srinivas Vasudevan <srvasude@gmail.com> Thu Dec 12 01:56:54 2019 +0000
committer Rasmus Munk Larsen <rmlarsen@google.com> Thu Dec 12 01:56:54 2019 +0000
tree 39a2b3928922d35946a324bd4828040ce6237709
parent 381f8f3139533739c4a428843e2495b442989fcc

Fix implementation of complex expm1. Add tests that fail with previous implementation, but pass with the current one.

test/array_cwise.cpp

diff --git a/test/array_cwise.cpp b/test/array_cwise.cpp
index ae0f9f9..48ebcc8 100644
--- a/test/array_cwise.cpp
+++ b/test/array_cwise.cpp
@@ -429,6 +429,11 @@
   VERIFY_IS_APPROX(m1.exp(), exp(m1));
   VERIFY_IS_APPROX(m1.exp() / m2.exp(),(m1-m2).exp());
 
+  VERIFY_IS_APPROX(m1.expm1(), expm1(m1));
+  VERIFY_IS_APPROX(expm1(m1), exp(m1) - 1.);
+  // Check for larger magnitude complex numbers that expm1 matches exp - 1.
+  VERIFY_IS_APPROX(expm1(10. * m1), exp(10. * m1) - 1.);
+
   VERIFY_IS_APPROX(sinh(m1), 0.5*(exp(m1)-exp(-m1)));
   VERIFY_IS_APPROX(cosh(m1), 0.5*(exp(m1)+exp(-m1)));
   VERIFY_IS_APPROX(tanh(m1), (0.5*(exp(m1)-exp(-m1)))/(0.5*(exp(m1)+exp(-m1))));