fvar_m38.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_matrix operator*(const prevariable& x, const dmatrix& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  int rmin=m.rowmin();
  int rmax=m.rowmax();

  ivector cmin(rmin,rmax);
  ivector cmax(rmin,rmax);

  const dvector* pmi = &m(rmin);
  int* pcmini = cmin.get_v() + rmin;
  int* pcmaxi = cmax.get_v() + rmin;
  for (int i=rmin;i<=rmax;++i)
  {
    *pcmini = pmi->indexmin();
    *pcmaxi = pmi->indexmax();

    ++pcmini;
    ++pcmaxi;
    ++pmi;
  }

  dvar_matrix tmp(rmin,rmax,cmin,cmax);

  pmi = &m(rmin);
  dvar_vector *ptmpi = &tmp(rmin);
  for (int i=rmin;i<=rmax;++i)
  {
    *ptmpi = x * *pmi;

    ++pmi;
    ++ptmpi;
  }

  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvar_matrix operator*(const dvar_matrix& m, const double x)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  int rmin=m.rowmin();
  int rmax=m.rowmax();

  ivector cmin(rmin,rmax);
  ivector cmax(rmin,rmax);

  const dvar_vector* pmi = &m(rmin);
  int* pcmini = cmin.get_v() + rmin;
  int* pcmaxi = cmax.get_v() + rmin;
  for (int i=rmin;i<=rmax;++i)
  {
    *pcmini = pmi->indexmin();
    *pcmaxi = pmi->indexmax();

    ++pcmini;
    ++pcmaxi;
    ++pmi;
  }
  dvar_matrix tmp(rmin,rmax,cmin,cmax);

  pmi = &m(rmin);
  dvar_vector *ptmpi = &tmp(rmin);
  for (int i=rmin;i<=rmax;++i)
  {
    *ptmpi = x * *pmi;

    ++pmi;
    ++ptmpi;
  }

  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvar_matrix operator*(const double x, const dvar_matrix& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  int rmin=m.rowmin();
  int rmax=m.rowmax();

  ivector cmin(rmin,rmax);
  ivector cmax(rmin,rmax);

  const dvar_vector* pmi = &m(rmin);
  int* pcmini = cmin.get_v() + rmin;
  int* pcmaxi = cmax.get_v() + rmin;
  for (int i=rmin;i<=rmax;++i)
  {
    *pcmini = pmi->indexmin();
    *pcmaxi = pmi->indexmax();

    ++pcmini;
    ++pcmaxi;
    ++pmi;
  }

  dvar_matrix tmp(rmin,rmax,cmin,cmax);

  pmi = &m(rmin);
  dvar_vector *ptmpi = &tmp(rmin);
  for (int i=rmin;i<=rmax;++i)
  {
    *ptmpi = x * *pmi;

    ++pmi;
    ++ptmpi;
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}