one_over_1_minus_x4.cpp

Go to the documentation of this file.

#include "one_over_1_minus_x.hpp"
#include <iostream>
 
int main()
{
  Tensor1<double, 3> y(0., 1., 2.);
  Tensor1<double, 3> a1(2., 3., 4.);
  Tensor1<double, 3> a2(5., 6., 7.);
  Tensor1<double, 3> a3(8., 9., 10.);
  Tensor1<double, 3> a4(11., 12., 13.);
 
  for(int ii = 0; ii < 100000000; ii++)
    {
      func4(y, a1, a2, a3, a4);
    }
  std::cout << y(0) << " " << y(1) << " " << y(2) << std::endl;
}