dvect14.cpp

Go to the documentation of this file.

/*
 * $Id$
 *
 * Author: David Fournier
 * Copyright (c) 2008-2012 Regents of the University of California
 */
// file fvar_fn.cpp
// math.h functions involving prevariables
#include "fvar.hpp"
#include <stdio.h>

double sfabs(const double v1)
 {
   #define cutoff 0.001
   #define cutoff2 0.000001
   double a = 2./cutoff;
   double b = 1./cutoff2;

   double return_value;
   double tmp1 = v1;
   double tmp = ::fabs(tmp1);
   if (tmp > cutoff)
   {
        return_value=tmp;
   }
   else
   {
     double t2 = tmp * tmp;

     if (tmp1 >=0)
     {
       return_value=t2*(a-b*tmp1);
     }
     else
     {
       return_value=t2*(a+b*tmp1);
     }
   }
   return(return_value);
 }

dvector sfabs(const dvector& t1)
{
  int min = t1.indexmin();
  int max = t1.indexmax();
  dvector tmp(min, max);
  double* ptmpi = tmp.get_v() + min;
  double* pt1i = t1.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *ptmpi = sfabs(*pt1i);
    ++pt1i;
    ++ptmpi;
  }
  return tmp;
}