one_over_1_minus_x3.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.);
 
  for(int ii = 0; ii < 1000000000; ii++)
    {
      func3(y, a1, a2, a3);
    }
  std::cout << y(0) << " " << y(1) << " " << y(2) << std::endl;
}