fvar_m30.cpp

Go to the documentation of this file.

/*
 * $Id$
 *
 * Author: David Fournier
 * Copyright (c) 2008-2012 Regents of the University of California
 */
// file fvar.cpp
// constructors, destructors and misc functions involving class dvariable

#include "fvar.hpp"

#ifdef __TURBOC__
  #pragma hdrstop
  #include <iostream.h>
#endif

#ifdef __ZTC__
  #include <iostream.hpp>
#endif


#include <stdio.h>
#ifndef __SUN__
#endif
#include <math.h>

dvar_matrix operator/(const dvar_matrix& m, const double e)
{
  dvar_matrix tmp;
  int min = m.indexmin();
  int max = m.indexmax();
  if (min <= max)
  {
    tmp.allocate(min, max);
    dvar_vector* ptmp = &tmp(min);
    const dvar_vector* pm = &m(min);
    for (int i = min; i <= max; ++i)
    {
      *ptmp = *pm / e;
      ++pm;
      ++ptmp;
    }
  }
  return tmp;
}

dvar_matrix operator/(const dvar_matrix& m, const prevariable& e)
{
  dvar_matrix tmp;
  int min = m.indexmin();
  int max = m.indexmax();
  if (min <= max)
  {
    tmp.allocate(min, max);
    dvar_vector* ptmp = &tmp(min);
    const dvar_vector* pm = &m(min);
    for (int i = min; i <= max; ++i)
    {
      *ptmp = *pm / e;
      ++pm;
      ++ptmp;
    }
  }
  return tmp;
}

dvar_matrix operator/(const dmatrix& m, const prevariable& e)
{
  dvar_matrix tmp;
  int min = m.indexmin();
  int max = m.indexmax();
  tmp.allocate(min, max);
  if (min <= max)
  {
    dvar_vector* ptmp = &tmp(min);
    const dvector* pm = &m(min);
    for (int i = min; i <= max; ++i)
    {
      *ptmp = *pm / e;
      ++pm;
      ++ptmp;
    }
  }
  return tmp;
}

dvar_matrix symmetrize(const dvar_matrix& m)
{
  if (m.rowmin() != m.colmin() || m.rowmax() != m.colmax() )
  {
    cerr << " Non square matrix passed to dmatrix symmetrize\n";
    ad_exit(1);
  }
  int rmin=m.rowmin();
  int rmax=m.rowmax();

  dvar_matrix s(rmin,rmax,rmin,rmax);


  for (int i=rmin;i<=rmax;i++)
  {
    s(i,i)=m(i,i);
    for (int j=rmin;j<i;j++)
    {
      s(i,j)=(m(i,j)+m(j,i))/2.;
      s(j,i)=s(i,j);
    }
  }

  return s;
}