// $Id: sinh.cpp 3785 2016-02-08 12:53:06Z bradbell $
/* --------------------------------------------------------------------------
CppAD: C++ Algorithmic Differentiation: Copyright (C) 2003-15 Bradley M. Bell
CppAD is distributed under multiple licenses. This distribution is under
the terms of the
GNU General Public License Version 3.
A copy of this license is included in the COPYING file of this distribution.
Please visit http://www.coin-or.org/CppAD/ for information on other licenses.
-------------------------------------------------------------------------- */
/*
Old example now just used for validation testing.
*/
# include <cppad/cppad.hpp>
bool Sinh(void)
{ bool ok = true;
using CppAD::sinh;
using CppAD::cosh;
using namespace CppAD;
// independent variable vector
CPPAD_TESTVECTOR(AD<double>) U(1);
U[0] = 1.;
Independent(U);
// dependent variable vector
CPPAD_TESTVECTOR(AD<double>) Z(1);
Z[0] = sinh(U[0]);
// create f: U -> Z and vectors used for derivative calculations
ADFun<double> f(U, Z);
CPPAD_TESTVECTOR(double) v(1);
CPPAD_TESTVECTOR(double) w(1);
// check value
double sin_u = sinh( Value(U[0]) );
double cos_u = cosh( Value(U[0]) );
ok &= NearEqual(sin_u, Value(Z[0]), 1e-10 , 1e-10);
// forward computation of partials w.r.t. u
size_t j;
size_t p = 5;
double jfac = 1.;
v[0] = 1.;
for(j = 1; j < p; j++)
{ w = f.Forward(j, v);
double value;
if( j % 2 == 1 )
value = cos_u;
else value = sin_u;
jfac *= j;
ok &= NearEqual(jfac*w[0], value, 1e-10 , 1e-10); // d^jz/du^j
v[0] = 0.;
}
// reverse computation of partials of Taylor coefficients
CPPAD_TESTVECTOR(double) r(p);
w[0] = 1.;
r = f.Reverse(p, w);
jfac = 1.;
for(j = 0; j < p; j++)
{
double value;
if( j % 2 == 0 )
value = cos_u;
else value = sin_u;
ok &= NearEqual(jfac*r[j], value, 1e-10 , 1e-10); // d^jz/du^j
jfac *= (j + 1);
}
return ok;
}