dvect16.cpp

Go to the documentation of this file.

#include "fvar.hpp"

dvector& dvector::operator+=(const dvector& values)
{
  int min = indexmin();
  int max = indexmax();
#ifndef OPT_LIB
  if (min != values.indexmin() || max != values.indexmax())
  {
    cerr << " Incompatible array bounds in "
         << "dvector& operator+=(const dvector&)\n";
    ad_exit(1);
  }
#endif

  double* pvi = v + min;
  const double* pvaluesi = values.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *pvi += *pvaluesi;
    ++pvi;
    ++pvaluesi;
  }
  return *this;
}
dvector& dvector::operator-=(const dvector& values)
{
  int min = indexmin();
  int max = indexmax();
#ifndef OPT_LIB
  if (min != values.indexmin() || max != values.indexmax())
  {
    cerr << " Incompatible array bounds in "
         << "dvector& operator-=(const dvector&)\n";
    ad_exit(1);
  }
#endif
  double* pvi = v + min;
  const double* pvaluesi = values.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *pvi -= *pvaluesi;
    ++pvi;
    ++pvaluesi;
  }
  return *this;
}
dvector& dvector::operator+=(const double value)
{
  int min = indexmin();
  int max = indexmax();
  double* pvi = v + min;
  for (int i = min; i <= max; ++i)
  {
    *pvi += value;
    ++pvi;
  }
  return *this;
}
dvector& dvector::operator-=(const double value)
{
  int min = indexmin();
  int max = indexmax();
  double* pvi = v + min;
  for (int i = min; i <= max; ++i)
  {
    *pvi -= value;
    ++pvi;
  }
  return *this;
}