\(\newcommand{\B}[1]{ {\bf #1} }\) \(\newcommand{\R}[1]{ {\rm #1} }\)

a_double_binary_xam.py

View page source

Python: a_double Binary Operators With AD Result: Example and Test

def a_double_binary_xam() :
   #
   import numpy
   from cppad_py import a_double
   from cppad_py import pow_int
   ok = True
   a2 = a_double(2.0)
   a3 = a_double(3.0)
   # ---------------------------------------------------------------------
   # a_double op a_double
   a5       = a2 + a3
   a6       = a2 * a3
   a1_minus = a2 - a3
   a23      = a2 / a3
   #
   ok = ok and a5 == 5.0
   ok = ok and a6 == 6.0
   ok = ok and a1_minus == -1.0
   ok = ok and a23.near_equal( a_double(2.0 / 3.0) )
   # ---------------------------------------------------------------------
   # a_double op double
   a5       = a2 + 3.0
   a6       = a2 * 3.0
   a1_minus = a2 - 3.0
   a23      = a2 / 3.0
   #
   ok = ok and a5 == 5.0
   ok = ok and a6 == 6.0
   ok = ok and a1_minus == -1.0
   ok = ok and a23.near_equal( a_double(2.0 / 3.0) )
   # ---------------------------------------------------------------------
   # double op a_double
   a5       = 3.0 + a2
   a6       = 3.0 * a2
   a1       = 3.0 - a2
   a32      = 3.0 / a2
   #
   ok = ok and a5 == 5.0
   ok = ok and a6 == 6.0
   ok = ok and a1 == 1.0
   ok = ok and a32.near_equal( a_double(3.0 / 2.0) )
   # ---------------------------------------------------------------------
   # pow
   a8 = a2  ** a3
   a9 = a3  ** 2.0
   a4 = 2.0 ** a2
   ok = ok and a8.near_equal( a_double(8.0) )
   ok = ok and a9.near_equal( a_double(9.0) )
   ok = ok and a4.near_equal( a_double(4.0) )
   # ---------------------------------------------------------------------
   # pow_int
   a8 = pow_int(a2, 3)
   a9 = pow_int(a3, 2)
   ok = ok and a8.near_equal( a_double(8.0) )
   ok = ok and a9.near_equal( a_double(9.0) )
   #
   return( ok )
#