d3arr4.cpp

Go to the documentation of this file.

#include "fvar.hpp"

double square(const double value)
{
  return value * value;
}
dvector square(const dvector& values)
{
  dvector results;
  results.allocate(values);
  int min = values.indexmin();
  int max = values.indexmax();

  double* pvaluesi = values.get_v() + min;
  double* presultsi = results.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = square(*pvaluesi);

    ++pvaluesi;
    ++presultsi;
  }
  return results;
}
ivector square(const ivector& values)
{
  ivector results;
  results.allocate(values);
  for (int i = values.indexmin(); i <= values.indexmax(); ++i)
  {
    results(i) = values(i) * values(i);
  }
  return results;
}
dmatrix square(const dmatrix& a)
{
  dmatrix results;
  results.allocate(a);
  for (int i = results.rowmin(); i <= results.rowmax(); ++i)
  {
    results(i) = square(a(i));
  }
  return results;
}
d3_array square(const d3_array& a)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = square(a(i));
  }
  return results;
}