cranfill.cpp

Go to the documentation of this file.

#include <fvar.hpp>

#ifdef __TURBOC__
  #pragma hdrstop
  #include <iostream.h>
#endif

#ifdef __ZTC__
  #include <iostream.hpp>
#endif

#include <math.h>

double var(const dvector& vec)
{
  double x = norm(vec) / sqrt(static_cast<double>(vec.size()));
  double y = mean(vec);
  return x * x - y * y;
}
double std_dev(const dvector& vec)
{
  return sqrt(var(vec));
}
double mean(const dvector& vec)
{
  double sum = 0;
  for (int i = vec.indexmin(); i <= vec.indexmax(); ++i)
  {
    sum += vec.elem(i);
  }
  return sum / vec.size();
}
void dvector::fill_seqadd(const double base, const double offset)
{
  double sequence = base;
  for (int i = indexmin(); i <= indexmax(); ++i)
  {
    elem(i) = sequence;
    sequence += offset;
  }
}
void ivector::fill_seqadd(int base, int offset)
{
  int sequence = base;
  for (int i = indexmin(); i <= indexmax(); ++i)
  {
    elem(i) = sequence;
    sequence += offset;
  }
}
void lvector::fill_seqadd(const AD_LONG_INT& base, const AD_LONG_INT& offset)
{
  AD_LONG_INT sequence = base;
  for (int i = indexmin(); i <= indexmax(); ++i)
  {
    elem(i) = sequence;
    sequence += offset;
  }
}
void dmatrix::colfill_seqadd(
  const int& column,
  const double base,
  const double offset)
{
  double sequence = base;
  for (int i = rowmin(); i <= rowmax(); ++i)
  {
    elem(i, column) = sequence;
    sequence += offset;
  }
}
void dmatrix::colfill(int column, const dvector& vec)
{
  for (int i = rowmin(); i <= rowmax(); ++i)
  {
    (*this)[i][column] = vec[i];
  }
}
void dmatrix::rowfill(int row, const dvector& vec)
{
  elem(row) = vec;
}
void dmatrix::rowfill_seqadd(
  const int& row,
  const double base,
  const double offset)
{
  elem(row).fill_seqadd(base, offset);
}