dmat6.cpp

Go to the documentation of this file.

#include "fvar.hpp"

dvector column(const dmatrix& matrix, int j)
{
  return extract_column(matrix, j);
}
dvector row(const dmatrix& matrix, int i)
{
  return extract_row(matrix, i);
}
dvector extract_column(const dmatrix& matrix, int j)
{
#ifndef OPT_LIB
  if (j < matrix.colmin() || j > matrix.colmax())
  {
    cerr << "Error: Invalid matrix column index specified in "
         << "dvector extract_column(const dmatrix&, int).\n",
    ad_exit(1);
  }
#endif
  int min = matrix.rowmin();
  int max = matrix.rowmax();
  dvector column(min, max);

  double* pcolumni = column.get_v() + min;
  const dvector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    *pcolumni = *(pmatrixi->get_v() + j);

    ++pmatrixi;
    ++pcolumni;
  }
  return column;
}
dvector extract_row(const dmatrix& matrix, int i)
{
  if (i < matrix.rowmin() || i > matrix.rowmax())
  {
    cerr << "Error: Invalid matrix row index specified in "
         << "dvector extract_row(const dmatrix&, int).\n",
    ad_exit(1);
  }

  dvector row(matrix(i).indexmin(), matrix(i).indexmax());
  row = matrix.elem(i);

  return row;
}