dmat8.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"
#include <cassert>

bool dmatrix::is_valid_row(const int i) const
{
  const bool valid = index_min <= i && i <= index_max;
  if (!valid)
  {
    cerr << "Error: Used invalid i = " << i << " for dmatrix rows bounded by ["
         << index_min << ", " << index_max << "].\n";
  }
  return valid;
}
#if !defined(OPT_LIB)

dvector& dmatrix::operator()(int i)
{
#ifndef OPT_LIB
  //check that index i is in range
  assert(index_min <= i && i <= index_max);
#endif

  return *(m + i);
}
#endif

#if !defined(OPT_LIB) || defined(__INTEL_COMPILER)

double& dmatrix::operator()(int i, int j)
{
#ifndef OPT_LIB
  //check that index i is in range
  assert(index_min <= i && i <= index_max);
#endif

  dvector* pmi =  m + i;

#ifndef OPT_LIB
  //check that index j is in range
  assert(pmi->is_valid_index(j));
#endif

  return *(pmi->get_v() + j);
}
const double& dmatrix::operator()(int i, int j) const
{
#ifndef OPT_LIB
  //check that index i is in range
  assert(index_min <= i && i <= index_max);
#endif

  dvector* pmi =  m + i;

#ifndef OPT_LIB
  //check that index j is in range
  assert(pmi->is_valid_index(j));
#endif

  return *(pmi->get_v() + j);
}
#endif