dmat4.cpp

Go to the documentation of this file.

#if (__cplusplus > 199711L)
  #include <iterator>
  #include <algorithm>
#endif

#include "fvar.hpp"

struct dvec_ptr_ptr
{
  void** m;
};

dmatrix& dmatrix::operator=(const dmatrix& other)
{
  if (allocated(*this))
  {
    if (rowmin() != other.rowmin() || rowmax() != other.rowmax()
        || colmin() != other.colmin() || colmax() != other.colmax())
    {
      cerr << "Incompatible array bounds in "
           << "dmatrix& dmatrix::operator=(const dmatrix&)\n";
      ad_exit(1);
    }
    // check for condition that both matrices point to the same object
    if (m != other.m)
    {
#if (__cplusplus <= 199711L)
      for (int i = rowmin(); i <= rowmax(); ++i)
      {
        *(m+i) = other[i];
      }
#else
      dvector* iter = other.begin();
      std::for_each(begin(), end(), [&iter](dvector& v) {
        v = *iter;
        ++iter;
      });
#endif
    }
  }
  else
  {
    index_min = other.index_min;
    index_max = other.index_max;
    shape = other.shape;
    if (shape)
    {
      (shape->ncopies)++;
    }
    m = other.m;
  }
  return *this;
}
dmatrix& dmatrix::operator+=(const dmatrix& other)
{
  if (rowmin() != other.rowmin() || rowmax() != other.rowmax() )
  {
    cerr << "Incompatible array bounds in "
         << "dmatrix& dmatrix::operator+=(const dmatrix&)\n";
    ad_exit(1);
  }
#if (__cplusplus <= 199711L)
  for (int i = rowmin(); i <= rowmax(); ++i)
  {
    elem(i) += other.elem(i);
  }
#else
  dvector* iter = other.begin();
  std::for_each(begin(), end(), [&iter](dvector& v) {
    v += *iter;
    ++iter;
  });
#endif
  return *this;
}
dmatrix& dmatrix::operator-=(const dmatrix& other)
{
  if (rowmin() != other.rowmin() || rowmax() != other.rowmax() )
  {
    cerr << "Incompatible array bounds in "
         << "dmatrix& dmatrix::operator-=(const dmatrix&)\n";
    ad_exit(1);
  }
#if (__cplusplus <= 199711L)
  for (int i = rowmin(); i <= rowmax(); ++i)
  {
    elem(i) -= other.elem(i);
  }
#else
  dvector* iter = other.begin();
  std::for_each(begin(), end(), [&iter](dvector& v) {
    v -= *iter;
    ++iter;
  });
#endif
  return *this;
}