i3arr.cpp

Go to the documentation of this file.

#include "fvar.hpp"

i3_array::i3_array()
{
  allocate();
}
i3_array::~i3_array()
{
  deallocate();
}
void i3_array::allocate(
  int sl, int sh,
  int nrl, int nrh,
  int ncl, int nch)
{
  allocate(sl, sh);
  for (int i = sl; i <= sh; ++i)
  {
    t[i].allocate(nrl, nrh, ncl, nch);
  }
}
void i3_array::allocate(int sl, int sh)
{
  if (sl > sh)
  {
    return allocate();
  }
  if ((shape = new three_array_shape(sl, sh)) == 0)
  {
    cerr << " Error: Unable to allocate memory in "
         << " i3_array::allocate(int, int).\n";
    ad_exit(1);
  }
  if ((t = new imatrix[slicesize()]) == 0)
  {
    cerr << " Error: Unable to allocate memory in "
         << " i3_array::allocate(int, int).\n";
    ad_exit(1);
  }
  t -= slicemin();
  for (int i = sl; i <= sh; ++i)
  {
    t[i].allocate();
  }
}
i3_array::i3_array(
  int sl, int sh,
  int nrl, int nrh,
  const ivector& ncl, const ivector& nch)
{
  allocate(sl, sh, nrl, nrh, ncl, nch);
}
i3_array::i3_array(int sl, int sh)
{
  allocate(sl, sh);
}
i3_array::i3_array(
  int sl, int sh,
  int nrl, const ivector& nrh,
  int ncl, int nch)
{
  allocate(sl, sh, nrl, nrh, ncl, nch);
}
void i3_array::allocate(
  int sl, int sh,
  int nrl, int nrh,
  const ivector& ncl, const ivector& nch)
{
  allocate(sl, sh);
  for (int i = sl; i <= sh; ++i)
  {
    elem(i).allocate(nrl, nrh, ncl(i), nch(i));
  }
}
void i3_array::allocate(
  int sl,int sh,
  int nrl,const ivector& nrh,
  int ncl, int nch)
{
  allocate(sl, sh);
  for (int i = sl; i <= sh; ++i)
  {
    elem(i).allocate(nrl, nrh(i), ncl, nch);
  }
}
i3_array::i3_array(
  int sl, int sh,
  int nrl, int nrh,
  int ncl, int nch)
{
  allocate(sl, sh, nrl, nrh, ncl, nch);
}
i3_array::i3_array(
  int sl, int sh,
  int nrl, const ivector& nrh,
  int ncl, const imatrix& nch)
{
  allocate(sl, sh, nrl, nrh, ncl, nch);
}
void i3_array::allocate(
  int sl, int sh,
  int nrl, const ivector& nrh,
  int ncl, const imatrix& nch)
{
  allocate(sl, sh);
  for (int i = slicemin(); i <= slicemax(); ++i)
  {
    t[i].allocate(nrl, nrh(i), ncl, nch(i));
  }
}
i3_array::i3_array(int sl, int sh, const imatrix& m1)
{
  allocate(sl, sh);
  for (int i = slicemin(); i <= slicemax(); ++i)
  {
    t[i].shallow_copy(m1);
  }
}
void i3_array::allocate(const i3_array& other)
{
  allocate(other.slicemin(), other.slicemax());
  for (int i = slicemin(); i <= slicemax(); ++i)
  {
    elem(i).allocate(other.elem(i));
  }
}
void i3_array::allocate()
{
  shape = nullptr;
  t = nullptr;
}
void i3_array::initialize()
{
  for (int i = slicemin(); i <= slicemax(); ++i)
  {
    elem(i).initialize();
  }
}
i3_array::i3_array(const i3_array& other)
{
  shallow_copy(other);
}
void i3_array::shallow_copy(const i3_array& other)
{
  if (other.shape)
  {
    shape = other.shape;
    ++(shape->ncopies);
    t = other.t;
  }
  else
  {
#ifdef DEBUG
    cerr << "Warning -- Unable to shallow copy an unallocated i3_array.\n";
#endif
    allocate();
  }
}
void i3_array::deallocate()
{
  if (shape)
  {
    if (shape->ncopies > 0)
    {
      --(shape->ncopies);
    }
    else
    {
      t += indexmin();
      delete [] t;
      delete shape;
    }
    allocate();
  }
#if defined(DIAG)
  else
  {
    cerr << "Warning -- Unable to deallocate an unallocated i3_array.\n";
  }
#endif
}