ranfill1.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 <math.h>

void dvar_vector::fill_seqadd(const double base, const double offset)
  {
    gradient_structure* gs = gradient_structure::_instance;
    gs->RETURN_ARRAYS_INCREMENT();

    double temp=0;
    for (int i=indexmin(); i<=indexmax(); i++)
    {
       elem(i) =base+temp;
       temp    =temp+offset;
    }
    gs->RETURN_ARRAYS_DECREMENT();
  }
void imatrix::colfill_seqadd(int j, int base, int offset)
  {
    int temp=0;
    for (int i=rowmin(); i<=rowmax(); i++)
    {
      (*this)[i][j] = base+temp;
      temp          = temp+offset;
    }
  }
void dvar_matrix::colfill_seqadd(const int& j, const double base,
  const double offset)
  {
    double temp=0;
    gradient_structure* gs = gradient_structure::_instance;
    gs->RETURN_ARRAYS_INCREMENT();
    for (int i=rowmin(); i<=rowmax(); i++)
    {
      elem(i,j) =base+temp;
      temp      =temp+offset;
    }
    gs->RETURN_ARRAYS_DECREMENT();
  }

void dvar_matrix::rowfill_seqadd(const int& i, const double base,
  const double offset)
  {
    double temp=0;
    gradient_structure* gs = gradient_structure::_instance;
    gs->RETURN_ARRAYS_INCREMENT();
    for (int j=colmin(); j<=colmax(); j++)
    {
      elem(i,j) =base+temp;
      temp      =temp+offset;
    }
    gs->RETURN_ARRAYS_DECREMENT();
  }

void dvar_matrix::colfill(int j, const dvar_vector& v)
  {
    gradient_structure* gs = gradient_structure::_instance;
    gs->RETURN_ARRAYS_INCREMENT();
    for (int i=rowmin(); i<=rowmax(); i++)
    {
      (*this)[i][j]=v[i];
    }
    gs->RETURN_ARRAYS_DECREMENT();
  }

void dvar_matrix::rowfill(int i, const dvar_vector& v)
  {
    gradient_structure* gs = gradient_structure::_instance;
    gs->RETURN_ARRAYS_INCREMENT();
    // for (int j=colmin(); j<=colmax(); j++)
    // {
    //   (*this)[i][j]=v[j];
    // }
    (*this)[i]=v;
    gs->RETURN_ARRAYS_DECREMENT();
  }