dvect12.cpp

Go to the documentation of this file.

#include <fvar.hpp>

double sum(const dvector& values)
{
  double total = 0.0;
  int min = values.indexmin();
  int max = values.indexmax();
  double* pvalues = values.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    total += *pvalues;

    ++pvalues;
  }
  return total;
}
dvector colsum(const dmatrix& matrix)
{
  int cmin = matrix.colmin();
  int cmax = matrix.colmax();
  int rmin = matrix.rowmin();
  int rmax = matrix.rowmax();

  dvector sums(cmin, cmax);
  sums.initialize();
  double* psumsj = sums.get_v() + cmin;
  for (int j=cmin; j<=cmax; ++j)
  {
    const dvector* pmatrixi = &matrix(rmin);
    for (int i=rmin; i<=rmax; ++i)
    {
      *psumsj += *(pmatrixi->get_v() + j);

      ++pmatrixi;
    }
    ++psumsj;
  }
  return sums;
}
dvector rowsum(const dmatrix& matrix)
{
  int min = matrix.rowmin();
  int max = matrix.rowmax();

  dvector sums(min, max);
  double* psumsi = sums.get_v() + min;
  const dvector* pmatrixi = &matrix(min);
  for (int i = min; i <= max; ++i)
  {
    *psumsi = sum(*pmatrixi);

    ++psumsi;
    ++pmatrixi;
  }
  return sums;
}
double sum(const dmatrix& matrix)
{
  double total = 0.0;
  int min = matrix.rowmin();
  int max = matrix.rowmax();

  const dvector* pmatrixi = &matrix.elem(min);
  for (int i = min; i <= max; ++i)
  {
    total += sum(*pmatrixi);

    ++pmatrixi;
  }
  return total;
}