fvar_a48.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 DF_cdble_dv_diveq(void);

dvar_vector& dvar_vector::operator/=(const double x)
  {
    RETURN_ARRAYS_INCREMENT();
    save_identifier_string("wctf");
    double xinv=1./x;
    for (int i=indexmin(); i<=indexmax(); i++)
    {
      elem_value(i)*=xinv;
    }
    save_dvar_vector_position();
    save_double_value(x);
    save_identifier_string("cmtu");
    RETURN_ARRAYS_DECREMENT();
    gradient_structure::GRAD_STACK1->set_gradient_stack(DF_cdble_dv_diveq);
    return(*this);
  }

 void DF_cdble_dv_diveq(void)
 {
    verify_identifier_string("cmtu");
    double x=restore_double_value();
    dvar_vector_position this_pos=restore_dvar_vector_position();
    dvector dfthis=restore_dvar_vector_derivatives(this_pos);
    verify_identifier_string("wctf");
    double xinv=1./x;
    for (int i=dfthis.indexmax(); i>=dfthis.indexmin(); i--)
    {
      // elem_value(i)=elem_value(i)/x;
      dfthis(i)*=xinv;
    }
    dfthis.save_dvector_derivatives(this_pos);
 }


  void DF_vdble_dv_diveq(void);

dvar_vector& dvar_vector::operator/=(const prevariable& x)
  {
    RETURN_ARRAYS_INCREMENT();
    double xinv=1./value(x);
    for (int i=indexmin(); i<=indexmax(); i++)
    {
      elem_value(i)=elem_value(i)*xinv;
    }
    save_identifier_string("wctg");
    save_dvar_vector_value();
    save_dvar_vector_position();
    x.save_prevariable_value();
    x.save_prevariable_position();
    save_identifier_string("cmtu");
    RETURN_ARRAYS_DECREMENT();
    gradient_structure::GRAD_STACK1->set_gradient_stack(DF_vdble_dv_diveq);
    return(*this);
  }

 void DF_vdble_dv_diveq(void)
 {
    verify_identifier_string("cmtu");
    prevariable_position x_pos=restore_prevariable_position();
    double dfx=restore_prevariable_derivative(x_pos);
    double x=restore_prevariable_value();
    dvar_vector_position this_pos=restore_dvar_vector_position();
    dvector tmp=restore_dvar_vector_value(this_pos);
    dvector dfthis=restore_dvar_vector_derivatives(this_pos);
    verify_identifier_string("wctg");
    double tmp1=0.;
    double xinv=1./x;
    int i;
    for (i=dfthis.indexmax(); i>=dfthis.indexmin(); i--)
    {
      // elem_value(i)=elem_value(i)/x;
      tmp1+=dfthis(i)*tmp(i);
    }
    tmp1*=-xinv;

    for (i=dfthis.indexmax(); i>=dfthis.indexmin(); i--)
    {
      // elem_value(i)=elem_value(i)/x;
      dfthis(i)*=xinv;
    }
    dfthis.save_dvector_derivatives(this_pos);
    save_double_derivative(dfx+tmp1,x_pos);
 }