fvar_a24.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"

  void dvar_vector::initialize(void)
  {
    if (!(!(*this)))  // only initialize allocated objects
    {
      for (int i=indexmin();i<=indexmax();i++)
      {
        //this->elem(i)=0.0;
        va[i].x=0.0;
      }
      grad_stack* GRAD_STACK1 = gradient_structure::GRAD_STACK1;
      DF_FILE* fp = gradient_structure::fp;
      save_identifier_string("b");
      fp->save_dvar_vector_position(*this);
      save_identifier_string("a");
      GRAD_STACK1->set_gradient_stack(dv_init);
    }
  }


void dv_init(void)
{
  DF_FILE* fp = gradient_structure::fp;

  // int ierr=fsetpos(gradient_structure::get_fp(),&filepos);
  verify_identifier_string("a");
  dvar_vector_position tmp_pos=fp->restore_dvar_vector_position();
  verify_identifier_string("b");
  dvector dftmp(tmp_pos.indexmin(),tmp_pos.indexmax());
  for (int i=dftmp.indexmin();i<=dftmp.indexmax();i++)
  {
    //vtmp.elem(i)=value(v1.elem(i))+value(v2.elem(i));
    dftmp.elem(i)=0.0;
  }
  dftmp.save_dvector_derivatives_na(tmp_pos);
  //ierr=fsetpos(gradient_structure::get_fp(),&filepos);
}


void dvar_vector::initialize(const dvector& ww)
  {
      if (indexmin() != ww.indexmin() ||  indexmax() != ww.indexmax())
      {
        cerr << "Index bounds do not match in "
              "void dvar_vector::initialize(const dvector& ww)"<<endl;
        ad_exit(1);
      }

    for (int i=indexmin();i<=indexmax();i++)
    {
      //this->elem(i)=0.0;
      va[i].x=ww.elem(i);
    }
    grad_stack* GRAD_STACK1 = gradient_structure::GRAD_STACK1;
    DF_FILE* fp = gradient_structure::fp;
    save_identifier_string("b");
    fp->save_dvar_vector_position(*this);
    save_identifier_string("a");
    GRAD_STACK1->set_gradient_stack(dv_init);
  }