fvar_ma3.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_vector operator*(const dvar_vector& x, const dmatrix& m)
 {
   RETURN_ARRAYS_INCREMENT();

   if (x.indexmin() != m.rowmin() || x.indexmax() != m.rowmax())
   {
     cerr << " Incompatible array bounds in dvar_vector  operator * (const dvar_vector& x,const dmatrix& m)\n";
     ad_exit(21);
   }
   dvar_vector tmp(m.colmin(),m.colmax());

   for (int j=m.colmin(); j<=m.colmax(); j++)
   {
     tmp[j]=0.0;
     for (int i=x.indexmin(); i<=x.indexmax(); i++)
     {
       tmp[j]+=x[i]*m[i][j];
     }
     //cout << "tmp[j]= "<< value(tmp[j])<<"\n";
   }
   RETURN_ARRAYS_DECREMENT();
   return(tmp);
 }

dvar_vector operator*(const dmatrix& m, const dvar_vector& x)
 {
   RETURN_ARRAYS_INCREMENT();

   if (x.indexmin() != m.colmin() || x.indexmax() != m.colmax())
   {
     cerr << " Incompatible array bounds in dvar_vector  operator*(const dvar_vector& x, const dmatrix& m)\n";
     ad_exit(21);
   }

   dvar_vector tmp(m.rowmin(),m.rowmax());

   for (int i=m.rowmin(); i<=m.rowmax(); i++)
   {
     tmp[i]=0.0;
     for (int j=x.indexmin(); j<=x.indexmax(); j++)
     {
       tmp[i]+=m[i][j]*x[j];
     }
   }
   RETURN_ARRAYS_DECREMENT();
   return(tmp);
 }
*/

dvariable norm(const dvar_matrix& m1)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  int min = m1.rowmin();
  int max = m1.rowmax();

  dvariable tmp;
  tmp = 0.0;

  const dvar_vector* pm1i = &m1(min);
  for (int i = min; i <= max; ++i)
  {
    tmp += norm2(*pm1i);

    ++pm1i;
  }

  if (tmp > 0.0)
  {
    tmp = pow(tmp, 0.5);
  }
  else
  {
    tmp = 0.0;
  }
  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}

dvariable norm2(const dvar_matrix& m1)
{
  gradient_structure* gs = gradient_structure::_instance;
  gs->RETURN_ARRAYS_INCREMENT();

  int min = m1.rowmin();
  int max = m1.rowmax();

  dvariable tmp;
  tmp = 0.0;

  const dvar_vector* pm1i = &m1(min);
  for (int i = min; i <= max; ++i)
  {
    tmp += norm2(*pm1i);

    ++pm1i;
  }

  gs->RETURN_ARRAYS_DECREMENT();
  return tmp;
}
dvariable sumsq(const dvar_matrix& m1) {return(norm2(m1));}

   /*
    dvar_matrix trans(const dvar_matrix& m1)
    {
      RETURN_ARRAYS_INCREMENT();

      dvar_matrix tmp(m1.colmin(),m1.colmax(),m1.rowmin(),m1.rowmax());

      for (int i=m1.rowmin(); i<=m1.rowmax(); i++)
      {
        for (int j=m1.colmin(); j<=m1.colmax(); j++)
        {
          tmp.elem(j,i)=m1.elem(i,j);
        }
      }
      RETURN_ARRAYS_DECREMENT();
      return (tmp);
    }
   */