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

#ifndef OPT_LIB

prevariable dvar_matrix::operator () (int i, int j)
{
  if (i < rowmin())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is lower than " << rowmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is higher than " << rowmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  dvar_vector* pmi = m + i;
  if (j < pmi->indexmin())
  {
    cerr << "matrix column bound exceeded -- column index " << i
         << " is lower than " << elem(i).indexmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (j > pmi->indexmax())
  {
    cerr << "matrix column bound exceeded -- column index " << i
         << " is higher than " << elem(i).indexmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  return pmi->va + j;
}
dvar_vector& dvar_matrix::operator[] (int i)
{
  if (i < rowmin())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is lower than " << rowmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is higher than " << rowmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  return m[i];
}
const prevariable dvar_matrix::operator()(int i, int j) const
{
  if (i < rowmin())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is lower than " << rowmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is higher than " << rowmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  dvar_vector* pmi = m + i;
  if (j < pmi->indexmin())
  {
    cerr << "matrix column bound exceeded -- column index " << i
         << " is lower than " << elem(i).indexmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (j > pmi->indexmax())
  {
    cerr << "matrix column bound exceeded -- column index " << i
         << " is higher than " << elem(i).indexmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  return pmi->va + j;
}
const dvar_vector& dvar_matrix::operator[](int i) const
{
  if (i < rowmin())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is lower than " << rowmin() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix row bound exceeded -- row index " << i
         << " is higher than " << rowmax() << " in "
         << __FILE__ << ':' << __LINE__ << ".\n";
    ad_exit(1);
  }
  return m[i];
}
#endif