f3arr6.cpp

Go to the documentation of this file.

#include "fvar.hpp"

dvar3_array square(const dvar3_array& vararr3)
{
  dvar3_array result;
  result.allocate(vararr3);

  int min = result.slicemin();
  int max = result.slicemax();
  dvar_matrix* presulti = &result(min);
  const dvar_matrix* pvararr3i = &vararr3(min);
  for (int i = min; i <= max; ++i)
  {
    *presulti = square(*pvararr3i);

    ++presulti;
    ++pvararr3i;
  }
  return result;
}