cmpdif11.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>

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

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

#include <string.h>

void banded_symmetric_dvar_matrix::save_dvar_matrix_value() const
{
  gradient_structure::fp->save_dvar_matrix_value(d);
}

void banded_lower_triangular_dvar_matrix::save_dvar_matrix_value() const
{
  gradient_structure::fp->save_dvar_matrix_value(d);
}

void banded_symmetric_dmatrix::save_dmatrix_value() const
{
  gradient_structure::fp->save_dmatrix_value(d);
}

banded_symmetric_dmatrix restore_banded_symmetric_dvar_matrix_value(
  const dvar_matrix_position& mpos)
{
  DF_FILE* fp = gradient_structure::fp;

  // restores the size, address, and value information for a dvar_matrix
  banded_symmetric_dmatrix out((const dvar_matrix_position&)mpos);
  //int ierr;
  int min=out.rowmin();
  int max=out.rowmax();
  for (int i=max;i>=min;i--)
  {
    dvar_vector_position vpos=fp->restore_dvar_vector_position();
    out.d(i)=restore_dvar_vector_value(vpos);
  }
  return out;
}

banded_lower_triangular_dmatrix
restore_banded_lower_triangular_dvar_matrix_value(
  const dvar_matrix_position& mpos)
{
  DF_FILE* fp = gradient_structure::fp;

  // restores the size, address, and value information for a dvar_matrix
  banded_lower_triangular_dmatrix out((const dvar_matrix_position&)mpos);
  //int ierr;
  int min=out.rowmin();
  int max=out.rowmax();
  for (int i=max;i>=min;i--)
  {
    dvar_vector_position vpos=fp->restore_dvar_vector_position();
    out.d(i)=restore_dvar_vector_value(vpos);
  }
  return out;
}

void banded_symmetric_dvar_matrix::save_dvar_matrix_position() const
{
  gradient_structure::fp->save_dvar_matrix_position(d);
}

void banded_lower_triangular_dvar_matrix::save_dvar_matrix_position() const
{
  gradient_structure::fp->save_dvar_matrix_position(d);
}

void banded_symmetric_dmatrix::save_dmatrix_derivatives(
  const dvar_matrix_position& pos) const
{
  // puts the derivative values in a dvector into a dvar_vector's guts
  int min=rowmin();
  int max=rowmax();
  if (min!=pos.row_min||max!=pos.row_max)
  {
    cerr << "Incompatible array sizes in " <<
     "void dmatrix::save_dmatrix__derivatives(const dvar_matrix_position& pos)"
     << endl;
  }
  for (int i=min;i<=max;i++)
  {
    const dvector& xx=(*this)(i);
    dvector& x=(dvector&) xx;
    dvar_matrix_position& pos1=(dvar_matrix_position&)pos;
    //((const dvector&)x).save_dvector_derivatives(pos1(i));
    x.save_dvector_derivatives(pos1(i));
  }
}

void banded_lower_triangular_dmatrix::save_dmatrix_derivatives(
  const dvar_matrix_position& pos) const
{
  // puts the derivative values in a dvector into a dvar_vector's guts
  int min=rowmin();
  int max=rowmax();
  if (min!=pos.row_min||max!=pos.row_max)
  {
    cerr << "Incompatible array sizes in " <<
    "void dmatrix::save_dmatrix__derivatives(const dvar_matrix_position& pos)"
    << endl;
  }
  for (int i=min;i<=max;i++)
  {
    dvector& x=(dvector&)d(i);
    dvar_matrix_position& pos1= (dvar_matrix_position&) pos;
    x.save_dvector_derivatives(pos1(i));
  }
}

banded_lower_triangular_dmatrix
  restore_banded_lower_triangular_dvar_matrix_derivatives
  (const dvar_matrix_position& _pos)
{
  dvar_matrix_position& pos= (dvar_matrix_position&) _pos;
  banded_lower_triangular_dmatrix tmp(pos);
  for (int i=pos.row_max;i>=pos.row_min;i--)
  {
    tmp(i)=restore_dvar_vector_derivatives(pos(i));
  }
  return tmp;
}

banded_symmetric_dmatrix restore_bounded_symmetric_dvar_matrix_derivatives
  (const dvar_matrix_position& _pos)
{
  dvar_matrix_position& pos= (dvar_matrix_position&) _pos;
  banded_symmetric_dmatrix tmp(pos);
  for (int i=pos.row_max;i>=pos.row_min;i--)
  {
    tmp(i)=restore_dvar_vector_derivatives(pos(i));
  }
  return tmp;
}