fvar_opr.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 <stdio.h>
#include <math.h>


void gradfree(dlink *);

void prevariable::operator+=(const double uu)
 {
   (*v).x+=uu;
 }

void prevariable::operator-=(const double uu)
 {
   (*v).x-=uu;
 }

prevariable& operator+(const prevariable& v1, const prevariable& v2)
{
  if (++gradient_structure::RETURN_PTR > gradient_structure::MAX_RETURN)
    gradient_structure::RETURN_PTR = gradient_structure::MIN_RETURN;

 gradient_structure::RETURN_PTR->v->x =(*v1.v).x+ (*v2.v).x;
 gradient_structure::GRAD_STACK1->set_gradient_stack4(default_evaluation4,
    &(gradient_structure::RETURN_PTR->v->x),&((*v1.v).x), &((*v2.v).x));
  return(*gradient_structure::RETURN_PTR);
}

prevariable& operator*(const prevariable& v1, const prevariable& v2)
{
  if (++gradient_structure::RETURN_PTR > gradient_structure::MAX_RETURN)
    gradient_structure::RETURN_PTR = gradient_structure::MIN_RETURN;
  gradient_structure::RETURN_PTR->v->x= (*v1.v).x * (*v2.v).x;
  gradient_structure::GRAD_STACK1->set_gradient_stack(default_evaluation3,
    &(gradient_structure::RETURN_PTR->v->x),
    &((*v1.v).x),(*v2.v).x,&((*v2.v).x),(*v1.v).x );
  return(*gradient_structure::RETURN_PTR);
 }

int prevariable::operator==(const prevariable& v1) const
    {
      return(v->x==value(v1));
    }

int prevariable::operator <=(const prevariable& v1) const
    {
      return(v->x<=value(v1));
    }

int prevariable::operator<(const prevariable& v1) const
    {
      return(v->x<value(v1));
    }

int prevariable::operator>(const prevariable& v1) const
    {
      return(v->x>value(v1));
    }

int prevariable::operator>=(const prevariable& v1) const
    {
      return(v->x >=value(v1));
    }

int prevariable::operator!=(const prevariable& v1) const
    {
      return(v->x!=value(v1));
    }