d3arr2a.cpp

Go to the documentation of this file.

#include "fvar.hpp"

d3_array log(const d3_array& arr3)
{
  d3_array results;
  results.allocate(arr3);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = log(arr3(i));
  }
  return results;
}
d3_array exp(const d3_array& arr3)
{
  d3_array results;
  results.allocate(arr3);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = exp(arr3(i));
  }
  return results;
}
d3_array sin(const d3_array& arr3)
{
  d3_array results;
  results.allocate(arr3);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = sin(arr3(i));
  }
  return results;
}
d3_array cos(const d3_array& arr3)
{
  d3_array results;
  results.allocate(arr3);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = cos(arr3(i));
  }
  return results;
}
d3_array tan(const d3_array& arr3)
{
  d3_array results;
  results.allocate(arr3);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = tan(arr3(i));
  }
  return results;
}
d3_array elem_prod(const d3_array& a, const d3_array& b)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = elem_prod(a(i), b(i));
  }
  return results;
}
d3_array elem_div(const d3_array& a, const d3_array& b)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = elem_div(a(i), b(i));
  }
  return results;
}
d3_array operator+(const d3_array& a, const d3_array& b)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = a(i) + b(i);
  }
  return results;
}
d3_array operator-(const d3_array& a, const d3_array& b)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = a(i) - b(i);
  }
  return results;
}
double norm2(const d3_array& a)
{
  double sum = 0.0;
  for (int i = a.slicemin(); i <= a.slicemax(); ++i)
  {
    sum += norm2(a(i));
  }
  return sum;
}
/*
Return sum of squared elements in a.

\param a d3_array
*/
double sumsq(const d3_array& a)
{
  return norm2(a);
}
double norm(const d3_array& a)
{
  return sqrt(norm2(a));
}
d3_array operator+(double value, const d3_array& a)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = value + a(i);
  }
  return results;
}
d3_array operator-(double value, const d3_array& a)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = value - a(i);
  }
  return results;
}
d3_array operator+(const d3_array& a, double value)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = a(i) + value;
  }
  return results;
}
d3_array operator-(const d3_array& a, double value)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = a(i) - value;
  }
  return results;
}
d3_array operator*(const d3_array& a, double value)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = a(i) * value;
  }
  return results;
}
d3_array operator*(double value, const d3_array& a)
{
  d3_array results;
  results.allocate(a);
  for (int i = results.slicemin(); i <= results.slicemax(); ++i)
  {
    results(i) = value * a(i);
  }
  return results;
}