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)
{
  //check that index i is in range
  assert((index_min <= i && i <= index_max) || is_valid_row(i));

  return *(m + i);
}
#endif

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

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

  dvector& dvi = elem(i);

  //check that index j is in range
  assert(dvi.is_valid_index(j));

  return *(dvi.v + j);
}
const double& dmatrix::operator()(int i, int j) const
{
  //check that index i is in range
  assert((index_min <= i && i <= index_max) || is_valid_row(i));

  const dvector& dvi = elem(i);

  //check that index j is in range
  assert(dvi.is_valid_index(j));

  return *(dvi.v + j);
}
#endif