f5arr1.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"

dvariable sum(const dvar5_array& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  dvariable tmp=0.;

  int min = m.indexmin();
  int max = m.indexmax();
  const dvar4_array* pmi = &m(min);
  for (int i = min; i <= max; ++i)
  {
    tmp+=sum(*pmi);

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


dvar5_array sqrt(const dvar5_array& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar5_array tmp;
  tmp.allocate(m);
  for (int i=tmp.indexmin();i<=tmp.indexmax();i++)
  {
    tmp(i)=sqrt(m(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvar5_array exp(const dvar5_array& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar5_array tmp;
  tmp.allocate(m);
  for (int i=tmp.indexmin();i<=tmp.indexmax();i++)
  {
    tmp(i)=exp(m(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvar5_array mfexp(const dvar5_array& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar5_array tmp;
  tmp.allocate(m);
  for (int i=tmp.indexmin();i<=tmp.indexmax();i++)
  {
    tmp(i)=mfexp(m(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvar5_array log(const dvar5_array& m)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar5_array tmp;
  tmp.allocate(m);
  for (int i=tmp.indexmin();i<=tmp.indexmax();i++)
  {
    tmp(i)=log(m(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}