cmpdif9.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
#endif

void imatrix::save_imatrix_value()
{
  gradient_structure::fp->save_imatrix_value(*this);
}
void DF_FILE::save_imatrix_value(const imatrix& m)
{
  for (int i=m.rowmin();i<=m.rowmax();++i)
  {
    save_ivector_value(m(i));
    save_ivector_position(m(i));
  }
}
void imatrix::save_imatrix_position()
{
  gradient_structure::fp->save_imatrix_position(*this);
}
void DF_FILE::save_imatrix_position(const imatrix& m)
{
  imatrix_position tmp(m);
  constexpr size_t wsize=sizeof(int);
  constexpr size_t wsize1=sizeof(void*);

  int min=m.rowmin();
  int max=m.rowmax();
  for (int i=min;i<=max;++i)
  {
    fwrite(&(tmp.lb(i)),wsize);
    fwrite(&(tmp.ub(i)),wsize);
    fwrite(&(tmp.ptr(i)),wsize1);
  }
  fwrite(&(tmp.row_min),wsize);
  fwrite(&(tmp.row_max),wsize);
}
imatrix_position restore_imatrix_position()
{
  return gradient_structure::fp->restore_imatrix_position();
}
imatrix_position DF_FILE::restore_imatrix_position()
{
  constexpr size_t wsize=sizeof(int);
  constexpr size_t wsize1=sizeof(void*);

  int min;
  int max;
  fread(&max,wsize);
  fread(&min,wsize);
  imatrix_position tmp(min,max);
  for (int i=max;i>=min;i--)
  {
    fread(&(tmp.ptr(i)),wsize1);
    fread(&(tmp.ub(i)),wsize);
    fread(&(tmp.lb(i)),wsize);
  }
  return tmp;
}
imatrix restore_imatrix_value(const imatrix_position& mpos)
{
  return gradient_structure::fp->restore_imatrix_value(mpos);
}
imatrix DF_FILE::restore_imatrix_value(const imatrix_position& mpos)
{
  // restores the size, address, and value information for a dvar_matrix
  //  the size, address, and value information for a dvar_matrix
  imatrix out(mpos);
  int min=out.rowmin();
  int max=out.rowmax();
  for (int i=max;i>=min;i--)
  {
    ivector_position vpos=restore_ivector_position();
    out(i)=restore_ivector_value(vpos);
  }
  return out;
}