imat8.cpp

Go to the documentation of this file.

/*
Author: David Fournier
Copyright (c) 2008-2012 Regents of the University of California
*/

#include <fvar.hpp>

ivector colsum(const imatrix& matrix)
{
  int jmin = matrix.colmin();
  int jmax = matrix.colmax();
  ivector colsums(jmin, jmax);

  int imin = matrix.rowmin();
  int imax = matrix.rowmax();

  int* pcolumnsj = colsums.get_v() + jmin;
  for (int j = jmin; j <= jmax; ++j)
  {
    *pcolumnsj = 0;

    const ivector* pmatrixi = &matrix(imin);
    for (int i = imin; i <= imax; ++i)
    {
      *pcolumnsj += *(pmatrixi->get_v() + j);

      ++pmatrixi;
    }

    ++pcolumnsj;
  }
  return colsums;
}
ivector rowsum(const imatrix& matrix)
{
  int min = matrix.rowmin();
  int max = matrix.rowmax();
  ivector rowsums(min, max);

  int* prowsumsi = rowsums.get_v() + min;
  const ivector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    *prowsumsi = sum(*pmatrixi);
    ++pmatrixi;
    ++prowsumsi;
  }

  return rowsums;
}
void imatrix::fill_seqadd(int start, int increment)
{
  int current = start;

  int imin = indexmin();
  int imax = indexmax();
  for (int i = imin; i <= imax; ++i)
  {
    if (allocated(elem(i)))
    {
      elem(i).fill_seqadd(current, increment);
      current = elem(i, elem(i).indexmax()) + increment;
    }
  }
}