imat2.cpp

Go to the documentation of this file.

#include "fvar.hpp"

#ifdef __TURBOC__
  #pragma hdrstop
#endif

#if !defined(OPT_LIB)

int& imatrix::operator()(int i, int j)
{
  if (i < rowmin())
  {
    cerr << "matrix bound exceeded -- row index too low in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << i << "\".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix bound exceeded -- row index too high in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << i << "\".\n";
    ad_exit(1);
  }
  ivector* pmi = m + i;
  if (j < pmi->indexmin())
  {
    cerr << "matrix bound exceeded -- column index too low in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << j << "\".\n";
    ad_exit(1);
  }
  if (j > pmi->indexmax())
  {
    cerr << "matrix bound exceeded -- column index too high in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << j << "\".\n";
    ad_exit(1);
  }
  return *(pmi->v + j);
}
const int& imatrix::operator()(int i, int j) const
{
  if (i < rowmin())
  {
    cerr << "matrix bound exceeded -- row index too low in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << i << "\".\n";
    ad_exit(1);
  }
  if (i > rowmax())
  {
    cerr << "matrix bound exceeded -- row index too high in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << i << "\".\n";
    ad_exit(1);
  }
  ivector* pmi = m + i;
  if (j < pmi->indexmin())
  {
    cerr << "matrix bound exceeded -- column index too low in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << j << "\".\n";
    ad_exit(1);
  }
  if (j > pmi->indexmax())
  {
    cerr << "matrix bound exceeded -- column index too high in "
         << "imatrix::operator()(int, int)\n"
         << "value was \"" << j << "\".\n";
    ad_exit(1);
  }
  return *(pmi->v + j);
}
#endif