Prev Next near_equal_ext.cpp Headings

@(@\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}} }@)@ This is cppad-20221105 documentation. Here is a link to its current documentation .
Compare AD with Base Objects: Example and Test

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

bool NearEqualExt(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;
}
Input File: example/general/near_equal_ext.cpp