fvar_a31.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>

dvar_vector colsum(const dvar_matrix& v)
{
  int cmin=v.colmin();
  int cmax=v.colmax();
  int rmin=v.rowmin();
  int rmax=v.rowmax();

  dvar_vector tmp(cmin, cmax);
  tmp.initialize();
  for (int j = cmin; j <= cmax; ++j)
  {
    for (int i = rmin; i <= rmax; ++i)
    {
      tmp(j) += v(i, j);
    }
  }
  return tmp;
}

dvar_vector rowsum(const dvar_matrix& v)
{
  //int cmin=v.colmin();
  //int cmax=v.colmax();
  int rmin=v.rowmin();
  int rmax=v.rowmax();

  dvar_vector tmp(rmin, rmax);
  const dvar_vector* pvi = &v(rmin);
  for (int i = rmin; i <= rmax; ++i)
  {
    tmp(i) = sum(*pvi);
    ++pvi;
  }
  return tmp;
}