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

#include <string.h>
int save_identifier_string(const char*);
void verify_identifier_string(const char*);

//struct dvar_matrix_position;
//struct dvector_position;
long int reset_gs_stack(void);
void reset_gs_stack(long int);

#ifndef OPT_LIB
  #define CHK_ID_STRING
#endif

dvar_vector_position::dvar_vector_position():
  min(0),
  max(-1),
  va(nullptr)
{
}
dvar_vector_position::dvar_vector_position(const dvar_vector_position& dvp):
  min(dvp.min),
  max(dvp.max),
  va(dvp.va)
{
}
dvar_vector_position::dvar_vector_position(const dvar_vector& v):
  min(v.indexmin()),
  max(v.indexmax()),
  va(v.get_va())
{
}
double& dvar_vector_position::operator()(const int& i)
{
  if (i < min || i > max)
  {
    cerr << "Error -- Index out of bounds in "
     "double& dvar_vector_position::operator()(const int& i)\n";
    ad_exit(1);
  }
  return va[i].x;
}
dvector_position::dvector_position(const dvector& vv)
{
  min=vv.indexmin();
  max=vv.indexmax();
  v=vv.get_v();
}
dvector_position::dvector_position(const dvector_position& dvp)
{
  min=dvp.min;
  max=dvp.max;
  v=dvp.v;
}
dvector_position::dvector_position(void)
{
  min=0;
  max=-1;
  v=0;
}

ivector_position::ivector_position(const ivector& iv)
{
  min=iv.indexmin();
  max=iv.indexmax();
  v=iv.get_v();
}
ivector_position::ivector_position(const ivector_position& dvp)
{
  min=dvp.min;
  max=dvp.max;
  v=dvp.v;
}
ivector_position::ivector_position(void)
{
  min=0;
  max=-1;
  v=0;
}
dvar_vector_position dvar_matrix_position::operator()(int i)
{
#ifndef OPT_LIB
  if (i<row_min||i>row_max)
  {
    cerr << "Index out of range in\n"
        "  dvar_vector_position dvar_matrix_position::operator () (int i)"
      << endl;
  }
#endif
  dvar_vector_position tmp;
  tmp.min = *(lb.get_v() + i);
  tmp.max = *(ub.get_v() + i);
  tmp.va=(double_and_int *) ptr(i);
  return tmp;
}

dvector_position dmatrix_position::operator()(int i)
{
#ifndef OPT_LIB
  if (i<row_min||i>row_max)
  {
    cerr << "Index out of range in\n"
        "  dvector_position dmatrix_position::operator () (int i)"
      << endl;
  }
#endif
  dvector_position tmp;
  tmp.min = *(lb.get_v() + i);
  tmp.max = *(ub.get_v() + i);
  tmp.v=(double *) ptr(i);
  return tmp;
}

dvar_matrix_position::dvar_matrix_position(const dvar_matrix& m,int x):
  lb(m.rowmin(),m.rowmax()),
  ub(m.rowmin(),m.rowmax()),
  ptr(m.rowmin(),m.rowmax())
{
  row_min=m.rowmin();
  row_max=m.rowmax();

  if (row_min <= row_max)
  {
    const dvar_vector* pm = &m(row_min);
    int* plb = lb.get_v() + row_min;
    int* pub = ub.get_v() + row_min;
    for (int i=row_min;i<=row_max;i++)
    {
      if (allocated(*pm))
      {
        *plb = pm->indexmin();
        *pub = pm->indexmax();
        ptr(i) = pm->get_va();
      }
      else
      {
        *plb = 0;
        *pub = -1;
        ptr(i) = 0;
      }
      ++pm;
      ++plb;
      ++pub;
    }
  }
}

dmatrix_position::dmatrix_position(const dmatrix& m)
  : lb(m.rowmin(),m.rowmax()), ub(m.rowmin(),m.rowmax()),
  ptr(m.rowmin(),m.rowmax())
{
  row_min=m.rowmin();
  row_max=m.rowmax();
  const dvector* pmi = &m(row_min);
  int* plbi = lb.get_v() + row_min;
  int* pubi = ub.get_v() + row_min;
  for (int i = row_min; i <= row_max; ++i)
  {
    *plbi = pmi->indexmin();
    *pubi = pmi->indexmax();
    ptr(i) = pmi->get_v();

    ++pmi;
    ++plbi;
    ++pubi;
  }
}

dvar_matrix_position::dvar_matrix_position(int min,int max)
  : lb(min,max), ub(min,max), ptr(min,max)
{
  row_min=min;
  row_max=max;
  int* plbi = lb.get_v() + row_min;
  int* pubi = ub.get_v() + row_min;
  for (int i=row_min;i<=row_max;i++)
  {
    *plbi = 0;
    *pubi = -1;
    ptr(i)=0;

    ++plbi;
    ++pubi;
  }
}

dmatrix_position::dmatrix_position(int min,int max)
  : lb(min,max), ub(min,max), ptr(min,max)
{
  row_min=min;
  row_max=max;
  int* plbi = lb.get_v() + row_min;
  int* pubi = ub.get_v() + row_min;
  for (int i=row_min;i<=row_max;i++)
  {
    *plbi = 0;
    *pubi = -1;
    ptr(i)=0;

    ++plbi;
    ++pubi;
  }
}

dvar_matrix_position::dvar_matrix_position(const dvar_matrix_position& p)
  : lb(p.row_min,p.row_max), ub(p.row_min,p.row_max),
    ptr(p.row_min,p.row_max)
{
  row_min=p.row_min;
  row_max=p.row_max;
  lb=p.lb;
  ub=p.ub;
  ptr=p.ptr;
  // cout << "ptr= " << ptr ;
}
dmatrix_position::dmatrix_position(const dmatrix_position& p)
  : lb(p.row_min,p.row_max), ub(p.row_min,p.row_max),
    ptr(p.row_min,p.row_max)
{
  row_min=p.row_min;
  row_max=p.row_max;
  lb=p.lb;
  ub=p.ub;
  ptr=p.ptr;
}
int save_identifier_string(const char* str)
{
#if defined(CHK_ID_STRING)
  //int wsize=sizeof(char);
  size_t length=strlen(str);
  gradient_structure::fp->fwrite(str,length);
#endif
  return 0;
}