fvar_ma2.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 dvar_matrix& m1, const dmatrix& m2)
 {
   gradient_structure* gs = gradient_structure::_instance;
   gs->RETURN_ARRAYS_INCREMENT();

   if (m1.colmin() != m2.colmin() || m1.colmax() != m2.colmax())
   {
     cerr << " Incompatible array bounds in "
     "dvar_matrix operator+(const dvar_matrix& x, const dmatrix& m)\n";
     ad_exit(21);
   }

   kkludge_object kk;
   dvar_matrix tmp(m1.rowmin(),m1.rowmax(),kk);

   for (int i=m1.rowmin(); i<=m1.rowmax(); i++)
   {
     tmp(i)=m1(i)+m2(i);
   }
   gs->RETURN_ARRAYS_DECREMENT();
   return(tmp);
 }

 dvar_matrix operator+(const dvar_matrix& m1, const dvar_matrix& m2)
 {
   gradient_structure* gs = gradient_structure::_instance;
   gs->RETURN_ARRAYS_INCREMENT();

   if (m1.colmin() != m2.colmin() || m1.colmax() != m2.colmax())
   {
     cerr << " Incompatible array bounds in "
     "dvar_matrix operator+(const dvar_matrix& x, const dvar_matrix& m)\n";
     ad_exit(21);
   }
   kkludge_object kk;
   dvar_matrix tmp(m1.rowmin(),m1.rowmax(),kk);
   //dvar_matrix tmp;
   //tmp.allocate(m1);

   for (int i=m1.rowmin(); i<=m1.rowmax(); i++)
   {
     tmp(i)=m1(i)+m2(i);
   }
   gs->RETURN_ARRAYS_DECREMENT();
   return(tmp);
 }

dvar_matrix operator+(const dmatrix& m1, const dvar_matrix& m2)
 {
   gradient_structure* gs = gradient_structure::_instance;
   gs->RETURN_ARRAYS_INCREMENT();

   if (m1.colmin() != m2.colmin() || m1.colmax() != m2.colmax())
   {
     cerr << " Incompatible array bounds in "
     "dvar_matrix operator+(const dmatrix& x, const dvar_matrix& m)\n";
     ad_exit(21);
   }
   kkludge_object kk;
   dvar_matrix tmp(m1.rowmin(),m1.rowmax(),kk);

   for (int i=m1.rowmin(); i<=m1.rowmax(); i++)
   {
     tmp(i)=m1(i)+m2(i);
   }
   gs->RETURN_ARRAYS_DECREMENT();
   return(tmp);
 }