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

dvar3_array elem_prod(const dvar3_array& m1, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  dvar3_array result;
  result.allocate(m1);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = elem_prod(m1(i), m2(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}
dvar3_array elem_div(const dvar3_array& m1, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar3_array result;
  result.allocate(m1);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = elem_div(m1(i), m2(i));
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}
dvar3_array operator+(const dvar3_array& m1, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar3_array result;
  result.allocate(m1);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = m1(i) + m2(i);
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}
dvar3_array operator-(const dvar3_array& m1, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar3_array result;
  result.allocate(m1);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = m1(i) - m2(i);
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}
dvar3_array operator+(const dvariable& d, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();
  dvar3_array result;
  result.allocate(m2);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = d + m2(i);
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}
dvar3_array operator-(const dvariable& d, const d3_array& m2)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  dvar3_array result;
  result.allocate(m2);
  for (int i = result.slicemin(); i <= result.slicemax(); ++i)
  {
    result(i) = d - m2(i);
  }
  gs->RETURN_ARRAYS_DECREMENT();

  return result;
}

dvar3_array operator*(const dvariable& d, const d3_array& m2)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();

     dvar3_array tmp;
     tmp.allocate(m2);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=d*m2(i);
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator*(const dvariable& d, const dvar3_array& m2)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m2);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=d*m2(i);
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator*(double d, const dvar3_array& m2)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m2);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=d*m2(i);
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator-(const dvar3_array& m1, const dvariable& d)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m1);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=m1(i)-d;
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator+(const dvar3_array& m1, const dvariable& d)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m1);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=m1(i)+d;
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator+(double d, const dvar3_array& m1)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m1);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=d+m1(i);
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator*(const dvar3_array& m1, const dvariable& d)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m1);
     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=m1(i)*d;
     }
     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }

dvar3_array operator+(const dvariable& d1, const dvar3_array& m1)
   {
     gradient_structure* gs = gradient_structure::_instance;
     gs->RETURN_ARRAYS_INCREMENT();
     dvar3_array tmp;
     tmp.allocate(m1);

     for (int i=tmp.slicemin();i<=tmp.slicemax();i++)
     {
       tmp(i)=m1(i)+d1;
     }

     gs->RETURN_ARRAYS_DECREMENT();
     return tmp;
   }