i4arr.cpp

Go to the documentation of this file.

#include <fvar.hpp>
#include "admb_messages.h"

void i4_array::allocate()
{
  t = nullptr;
  shape = nullptr;
}
void i4_array::initialize()
{
  for (int i = hslicemin(); i <= hslicemax(); ++i)
  {
    if (allocated(elem(i)))
    {
      elem(i).initialize();
    }
  }
}
i4_array::i4_array()
{
  allocate();
}
i4_array::i4_array(int hsl, int hsu)
{
  allocate(hsl, hsu);
}
void i4_array::allocate(int hsl, int hsu)
{
  if (hsl > hsu)
  {
    return allocate();
  }
  if ((t = new i3_array[static_cast<unsigned int>(hsu - hsl + 1)]) == 0)
  {
    cerr << " Error: i4_array unable to allocate t in "
         << "i4_array::allocate(int, int).\n";
    ad_exit(1);
  }
  if ((shape=new vector_shapex(hsl, hsu, t)) == 0)
  {
    cerr << " Error: i4_array unable to allocate shape in "
         << "i4_array::allocate(int, int).\n";
    ad_exit(1);
  }
  t -= indexmin();
  for (int i = hsl; i <= hsu; ++i)
  {
    elem(i).allocate();
  }
}
i4_array::i4_array(
  int hsl, int hsu,
  int sl, const ivector& sh,
  int nrl, const imatrix& nrh,
  int ncl, const i3_array& nch)
{
  allocate(hsl,hsu,sl,sh,nrl,nrh,ncl,nch);
}
i4_array::i4_array(
  int hsl, int hsu,
  int sl, int sh,
  int nrl, int nrh,
  int ncl, int nch)
{
  allocate(hsl, hsu, sl, sh, nrl, nrh, ncl, nch);
}
void i4_array::allocate(
  int hsl, int hsu,
  int sl, int sh,
  int nrl, int nrh,
  int ncl, int nch)
{
  if (hsl > hsu)
  {
    return allocate();
  }
  allocate(hsl, hsu);
  for (int i = hslicemin(); i <= hslicemax(); ++i)
  {
    elem(i).allocate(sl, sh, nrl, nrh, ncl, nch);
  }
}
void i4_array::allocate(
  const ad_integer& hsl, const ad_integer& hsu,
  const index_type& sl, const index_type& sh,
  const index_type& nrl, const index_type& nrh,
  const index_type& ncl, const index_type& nch)
{
  unsigned int ss =
    static_cast<unsigned int>(hsu < hsl ? 0 : hsu - hsl + 1);
  if (ss > 0)
  {
    if ((t = new i3_array[ss]) == 0)
    {
      cerr << " Error: i4_array unable to allocate memory in "
           << __FILE__ << ':' << __LINE__ << '\n';
      ad_exit(1);
    }
    if ((shape = new vector_shapex(hsl, hsu, t)) == 0)
    {
       cerr << " Error: i4_array unable to allocate memory in "
            << __FILE__ << ':' << __LINE__ << '\n';
       ad_exit(1);
    }
    t -= indexmin();
    for (int i = hsl; i <= hsu; ++i)
    {
      (*this)(i).allocate(
         ad_integer(sl(i)), ad_integer(sh(i)),
         nrl(i), nrh(i),
         ncl(i), nch(i));
    }
  }
  else
  {
    t = 0;
    shape = 0;
  }
}
i4_array::i4_array(const i4_array& other)
{
  shallow_copy(other);
}
void i4_array::shallow_copy(const i4_array& other)
{
  if (other.shape)
  {
    shape = other.shape;
    ++(shape->ncopies);
    t = other.t;
  }
  else
  {
#ifdef DEBUG
    cerr << "Warning -- Unable to shallow copy an unallocated i4_array.\n";
#endif
    allocate();
  }
}
i4_array::~i4_array()
{
  deallocate();
}
void i4_array::deallocate()
{
  if (shape)
  {
    if (shape->ncopies > 0)
    {
      --(shape->ncopies);
    }
    else
    {
      t = static_cast<i3_array*>(shape->get_truepointer());
      delete [] t;
      delete shape;
    }
    allocate();
  }
#if defined(DIAG)
  else
  {
    cerr << "Warning -- Unable to deallocate an unallocated i4_array.\n";
  }
#endif
}