dvect6.cpp

Go to the documentation of this file.

#include "fvar.hpp"
#include <cmath>

#ifdef __TURBOC__
  #pragma hdrstop
#endif

void shape_check(const dvector& v1, const dvector& v2,
  const char *function_name);

dvector sin(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::sin(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector cos(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::cos(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector tan(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::tan(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector sinh(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::sinh(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector cosh(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::cosh(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector tanh(const dvector& vec)
{
  int min = vec.indexmin();
  int max = vec.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pveci = vec.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::tanh(*pveci);

    ++presultsi;
    ++pveci;
  }
  return results;
}
dvector pow(const dvector& bases, const double exponent)
{
  int min = bases.indexmin();
  int max = bases.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pbasesi = bases.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::pow(*pbasesi, exponent);

    ++presultsi;
    ++pbasesi;
  }
  return results;
}
dvector pow(const dvector& bases, int exponent)
{
  int min = bases.indexmin();
  int max = bases.indexmax();
  dvector results(min, max);
  double* presultsi = results.get_v() + min;
  double* pbasesi = bases.get_v() + min;
  for (int i = min; i <= max; ++i)
  {
    *presultsi = std::pow(*pbasesi, exponent);

    ++presultsi;
    ++pbasesi;
  }
  return results;
}
dvector pow(const dvector& bases, const dvector& exponents)
{
  shape_check(bases, exponents,
    "dvector pow(const dvector& bases ,constdvector& exponents)");

  dvector results(bases.indexmin(), bases.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::pow(bases.elem(i), exponents.elem(i));
  }
  return results;
}
dvector pow(const double base, const dvector& exponents)// ***
{
  dvector results(exponents.indexmin(), exponents.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::pow(base, exponents.elem(i));
  }
  return results;
}
dvector asin(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::asin(vec.elem(i));
  }
  return results;
}
dvector acos(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::acos(vec.elem(i));
  }
  return results;
}
dvector atan(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::atan(vec.elem(i));
  }
  return results;
}
dvector log10(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::log10(vec.elem(i));
  }
  return results;
}
dvector sqrt(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    results.elem(i) = std::sqrt(vec.elem(i));
  }
  return results;
}
dvector sqr(const dvector& vec)
{
  dvector results(vec.indexmin(), vec.indexmax());
  for (int i = results.indexmin(); i <= results.indexmax(); ++i)
  {
    double x = vec.elem(i);
    results.elem(i) = x * x;
  }
  return results;
}