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 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)
{
  dvariable* RETURN_PTR = gradient_structure::next_RETURN_PTR();

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

  return *RETURN_PTR;
}

prevariable& operator*(const prevariable& v1, const prevariable& v2)
{
  dvariable* RETURN_PTR = gradient_structure::next_RETURN_PTR();

  double_and_int* di1 = v1.v;
  double_and_int* di2 = v2.v;
  RETURN_PTR->v->x = di1->x * di2->x;
  gradient_structure::GRAD_STACK1->set_gradient_stack(default_evaluation3,
    &(RETURN_PTR->v->x),
    &(di1->x), di2->x, &(di2->x), di1->x);

  return *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));
    }