tradoptest_62.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 + ****/
 
  xd = 17.; f = 17.; dx = Plus_dx;
  xAI = xd;
  fA = +xAI;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xAI", f, fA.val());
  else if (differ(xAI.adj(), dx)) botch("d +xAI/dx", dx, xAI.adj());
  {
  A::aval_reset();
  cAI xcAI(xd);
  fA = +xcAI;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xcAI", f, fA.val());
  else if (differ(xcAI.adj(), dx)) botch("d +xcAI/dx", dx, xcAI.adj());
  }
  xA = xd;
  fA = +xA;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xA", f, fA.val());
  else if (differ(xA.adj(), dx)) botch("d +xA/dx", dx, xA.adj());
  {
  A::aval_reset();
  cA xcA(xd);
  fA = +xcA;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xcA", f, fA.val());
  else if (differ(xcA.adj(), dx)) botch("d +xcA/dx", dx, xcA.adj());
  }
  xC = xd;
  fA = +xC;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xC", f, fA.val());
  else if (differ(xC.adj(), dx)) botch("d +xC/dx", dx, xC.adj());
  {
  A::aval_reset();
  cC xcC(xd);
  fA = +xcC;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xcC", f, fA.val());
  else if (differ(xcC.adj(), dx)) botch("d +xcC/dx", dx, xcC.adj());
  }
  {
  cAi xcAi(xd);
  fA = +xcAi;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xcAi", f, fA.val());
  else if (differ(xcAi.aval, dx)) botch("d +xcAi/dx", dx, xcAi.aval);
  }
  {
  A::aval_reset();
  cAi xcAi(xd);
  fA = +xcAi;
  A::Gradcomp();
  if (differ(fA.val(), f)) botch("fA = +xcAi", f, fA.val());
  else if (differ(xcAi.aval, dx)) botch("d +xcAi/dx", dx, xcAi.aval);
  }
 
 
  if (!rc) // chatter for cppunit test, which cannot tolerate silence
 
    printf("OK\n");
 
  return rc;
 
  }