cmpdif8.cpp

Go to the documentation of this file.

/*
 * $Id$
 *
 * Author: David Fournier
 * Copyright (c) 2008-2012 Regents of the University of California
 */
#include <fvar.hpp>

#ifdef __TURBOC__
  #pragma hdrstop
  #include <iostream.h>
#endif

#ifdef __ZTC__
  #include <iostream.hpp>
#endif

#include <string.h>

void save_double_derivative(const double x, const prevariable_position& _pos)
{
  prevariable_position& pos= (prevariable_position&) _pos;
  // adds the derivative value in a double to a prevariables' guts
  pos.xval()+=x;
}

prevariable_position restore_prevariable_position()
{
  return gradient_structure::fp->restore_prevariable_position();
}
prevariable_position DF_FILE::restore_prevariable_position()
{
  constexpr size_t wsize = sizeof(double_and_int*);

  double_and_int* tmp;
  fread(&tmp, wsize);
  return prevariable_position(tmp);
}

void prevariable::save_prevariable_position() const
{
  gradient_structure::fp->save_prevariable_position(*this);
}
void DF_FILE::save_prevariable_position(const prevariable& v)
{
  constexpr size_t wsize = sizeof(double_and_int*);

  double_and_int* tmp = v.get_v();
  fwrite(&tmp, wsize);
}

void prevariable::save_prevariable_value() const
{
  gradient_structure::fp->save_prevariable_value(*this);
}
void DF_FILE::save_prevariable_value(const prevariable& v)
{
  //double_and_int * tmp=get_v();
  //const unsigned wsize=sizeof(double_and_int*);
  //gradient_structure::get_fp()->fwrite(&tmp,wsize);
  double x=value(v);
  //const unsigned dsize=sizeof(double);
  fwrite(x);
}

void save_double_value(const double x)
{
  gradient_structure::fp->save_double_value(x);
}
void DF_FILE::save_double_value(const double x)
{
  //const unsigned wsize=sizeof(double);
  fwrite(x);
}

void save_int_value(int x)
{
  gradient_structure::fp->save_int_value(x);
}
void DF_FILE::save_int_value(int x)
{
  //const unsigned wsize=sizeof(double);
  fwrite(x);
}

void save_pointer_value(void *ptr)
{
  gradient_structure::fp->save_pointer_value(ptr);
}
void DF_FILE::save_pointer_value(void *ptr)
{
  fwrite(ptr);
}

double restore_prevariable_derivative(const prevariable_position& _pos)
{
  prevariable_position& pos= (prevariable_position&) _pos;
  double tmpout=pos.xval();
  pos.xval()=0.0;
  return tmpout;
}

double restore_prevariable_derivative()
{
  return gradient_structure::fp->restore_prevariable_derivative();
}
double DF_FILE::restore_prevariable_derivative()
{
  // Back up the stream and read the number of bytes written in the
  // ``write function'' corresponding to this ``read function''
  constexpr size_t wsize = sizeof(double_and_int*);

  double_and_int* tmp;
  fread(&tmp, wsize);
  double tmpout=tmp->x;
  tmp->x=0.0;
  return tmpout;
}

double restore_prevariable_value()
{
  return gradient_structure::fp->restore_prevariable_value();
}
double DF_FILE::restore_prevariable_value()
{
  // Back up the stream and read the number of bytes written in the
  // ``write function'' corresponding to this ``read function''
  constexpr size_t wsize = sizeof(double);
  double tmpout = 0;
  fread(&tmpout, wsize);
  return tmpout;
}

double restore_double_value()
{
  return gradient_structure::fp->restore_double_value();
}
double DF_FILE::restore_double_value()
{
  double tmpout = 0;
  fread(tmpout);
  return tmpout;
}

int restore_int_value(void)
{
  return gradient_structure::fp->restore_int_value();
}
int DF_FILE::restore_int_value()
{
  int tmpout = 0;
  fread(tmpout);
  return tmpout;
}

void* restore_pointer_value()
{
  return gradient_structure::fp->restore_pointer_value();
}
void* DF_FILE::restore_pointer_value()
{
  void* tmpout = NULL;
  fread(tmpout);
  return tmpout;
}