imat7.cpp

Go to the documentation of this file.

#include "fvar.hpp"
#include "admb_messages.h"

int sum(const imatrix& matrix)
{
  int min = matrix.rowmin();
  int max = matrix.rowmax();

  int total = 0;
  const ivector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    total += sum(*pmatrixi);
    ++pmatrixi;
  }
  return total;
}
int colsum(const imatrix& matrix, int column)
{
#ifndef OPT_LIB
  if (column < matrix.colmin() || column > matrix.colmax())
  {
    ADMB_ARRAY_BOUNDS_ERROR("Specified column is out of bounds",
    "in colsum(const imatrix&,int column)",
    matrix.colmin(), matrix.colmax(), column);
  }
#endif
  int sum = 0;
  int min = matrix.rowmin();
  int max = matrix.rowmax();
  const ivector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    sum += *(pmatrixi->get_v() + column);
    ++pmatrixi;
  }
  return sum;
}
ivector column(const imatrix& matrix, int column)
{
#ifndef OPT_LIB
  if (column < matrix.colmin() || column > matrix.colmax())
  {
    ADMB_ARRAY_BOUNDS_ERROR("Specified column is out of bounds",
    "in column(const imatrix&, int column)",
    matrix.colmin(), matrix.colmax(), column);
  }
#endif
  int min = matrix.rowmin();
  int max = matrix.rowmax();
  ivector vector(min, max);
  int* pvectori = &vector(min);
  const ivector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    *pvectori = *(pmatrixi->get_v() + column);
    ++pmatrixi;
    ++pvectori;
  }
  return vector;
}