fvar_m33.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"

void dmatrix::allocate(const dvar_matrix& dm)
{
  int nrl = dm.rowmin();
  int nrh = dm.rowmax();

  index_min = nrl;
  index_max = nrh;

  if ((m = new dvector[rowsize()]) == 0)
  {
    cerr << " Error allocating memory in dmatrix contructor\n";
    ad_exit(21);
  }
  if ((shape = new mat_shapex(m)) == 0)
  {
    cerr << " Error allocating memory in dmatrix contructor\n";
    ad_exit(21);
  }

  m -= nrl;
  dvector* pmi = m + nrl;
  const dvar_vector* pdmi = &dm(nrl);
  for (int i = nrl; i <= nrh; ++i)
  {
    pmi->allocate(*pdmi);

    ++pmi;
    ++pdmi;
  }
}