\(\newcommand{\W}[1]{ \; #1 \; }\) \(\newcommand{\R}[1]{ {\rm #1} }\) \(\newcommand{\B}[1]{ {\bf #1} }\) \(\newcommand{\D}[2]{ \frac{\partial #1}{\partial #2} }\) \(\newcommand{\DD}[3]{ \frac{\partial^2 #1}{\partial #2 \partial #3} }\) \(\newcommand{\Dpow}[2]{ \frac{\partial^{#1}}{\partial {#2}^{#1}} }\) \(\newcommand{\dpow}[2]{ \frac{ {\rm d}^{#1}}{{\rm d}\, {#2}^{#1}} }\)

near_equal_ext.cpp

View page source

Compare AD with Base Objects: Example and Test

# include <cppad/cppad.hpp>
# include <complex>

bool near_equal_ext(void)
{   bool ok = true;
    using CppAD::AD;
    using CppAD::NearEqual;

    // double
    double x    = 1.00000;
    double y    = 1.00001;
    double a    =  .00005;
    double r    =  .00005;
    double zero = 0.;

    // AD<double>
    AD<double> ax(x);
    AD<double> ay(y);

    ok &= NearEqual(ax, ay, zero, a);
    ok &= NearEqual(ax, y,  r, zero);
    ok &= NearEqual(x, ay,  r,    a);

    // std::complex<double>
    AD<double> cx(x);
    AD<double> cy(y);

    // AD< std::complex<double> >
    AD<double> acx(x);
    AD<double> acy(y);

    ok &= NearEqual(acx, acy, zero, a);
    ok &= NearEqual(acx,  cy, r, zero);
    ok &= NearEqual(acx,   y, r,    a);
    ok &= NearEqual( cx, acy, r,    a);
    ok &= NearEqual(  x, acy, r,    a);

    return ok;
}