tradoptest_47.cpp

Go to the documentation of this file.

 
 
/* Try to test all combinations of types and operations */
 
 
 
#define ADT_RAD Sacado::Rad::
 
 
 
#include "Sacado_trad.hpp"
#include <cstdio>
 
using std::printf;
 
 
 
typedef ADT_RAD IndepADvar<double> AI;
 
typedef ADT_RAD ADvar<double> A;
 
typedef ADT_RAD ConstADvar<double> C;
 
typedef ADT_RAD ADvari<double> Ai;
 
typedef const ADT_RAD IndepADvar<double> cAI;
 
typedef const ADT_RAD ADvar<double> cA;
 
typedef const ADT_RAD ConstADvar<double> cC;
 
typedef const ADT_RAD ADvari<double> cAi;
 
static int rc;
 
 
 
/* This is to be run through an awk program that changes lines */
 
/* with "BINTEST" or "UNOPTEST" at the beginning of the line into */
 
/* a the desired C++ (which we can then inspect). */
 
 
 
 void
 
botch(const char *what, double wanted, double got)
 
{
 
  printf("%s: expected %g, got %g, diff = %.2g\n", what, wanted, got, wanted-got);
 
  rc = 1;
 
  }
 
 
 
 const double tol = 5e-16;
 
 
 
 int
 
differ(double a, double b)
 
{
 
  double d = a - b;
 
  if (d < 0.)
 
    d = -d;
 
  if (a < 0.)
 
    a = -a;
 
  if (b < 0.)
 
    b = -b;
 
  if (a < b)
 
    a = b;
 
  if (a > 0.)
 
    d /= a;
 
  return d > tol;
 
  }
 
 
 
#ifndef RAD_EQ_ALIAS
 
#define Plus_dx 1.
 
#else
 
#ifdef RAD_AUTO_AD_Const
 
#define Plus_dx 1.
 
#else
 
#define Plus_dx 0.
 
#endif
 
#endif
 
 
 
 int
 
main(void)
 
{
 
  AI xAI, yAI;
 
  A fA, xA, yA;
 
  C xC, yC;
 
  double dx, f, xd;
 
 
 
  rc = 0;
 
 
  /**** Test of atanh ****/
 
  xd = .6; f = atanh(.6); dx = 1.5625;
  xAI = xd;
  fA = atanh(xAI);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xAI)", f, fA.val());
  else if (differ(xAI.adj(), dx)) botch("d atanh(xAI)/dx", dx, xAI.adj());
  {
  A::aval_reset();
  cAI xcAI(xd);
  fA = atanh(xcAI);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xcAI)", f, fA.val());
  else if (differ(xcAI.adj(), dx)) botch("d atanh(xcAI)/dx", dx, xcAI.adj());
  }
  xA = xd;
  fA = atanh(xA);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xA)", f, fA.val());
  else if (differ(xA.adj(), dx)) botch("d atanh(xA)/dx", dx, xA.adj());
  {
  A::aval_reset();
  cA xcA(xd);
  fA = atanh(xcA);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xcA)", f, fA.val());
  else if (differ(xcA.adj(), dx)) botch("d atanh(xcA)/dx", dx, xcA.adj());
  }
  xC = xd;
  fA = atanh(xC);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xC)", f, fA.val());
  else if (differ(xC.adj(), dx)) botch("d atanh(xC)/dx", dx, xC.adj());
  {
  A::aval_reset();
  cC xcC(xd);
  fA = atanh(xcC);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xcC)", f, fA.val());
  else if (differ(xcC.adj(), dx)) botch("d atanh(xcC)/dx", dx, xcC.adj());
  }
  {
  cAi xcAi(xd);
  fA = atanh(xcAi);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xcAi)", f, fA.val());
  else if (differ(xcAi.aval, dx)) botch("d atanh(xcAi)/dx", dx, xcAi.aval);
  }
  {
  A::aval_reset();
  cAi xcAi(xd);
  fA = atanh(xcAi);
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = atanh(xcAi)", f, fA.val());
  else if (differ(xcAi.aval, dx)) botch("d atanh(xcAi)/dx", dx, xcAi.aval);
  }
 
 
  if (!rc) // chatter for cppunit test, which cannot tolerate silence
 
    printf("OK\n");
 
  return rc;
 
  }