dmat27.cpp

Go to the documentation of this file.

#include "fvar.hpp"

double colsum(const dmatrix& m, int col)
{
  int mmin = m.rowmin();
  int mmax = m.rowmax();
#ifndef OPT_LIB
  if (col < mmin || col > mmax)
  {
    cerr << "Row out of bounds in function"
            " colsum(const imatrix& m,int col)" << endl;
    ad_exit(1);
  }
#endif
  double sum = 0;
  const dvector* pmi = &m(mmin);
  for (int i = mmin; i <= mmax; ++i)
  {
    sum += *(pmi->get_v() + col);
    ++pmi;
  }
  return sum;
}