#include <cmath>
#include "Pooma/NewFields.h"

int main(int argc, char *argv[])
{
  Pooma::initialize(argc,argv);

  Array<2,Vector<2> > c(5,4);

  // How do I write a data-parallel statement equivalent to?
  //  for (unsigned i = all values in c's domain's first component)
  //    for (unsigned j = all values in c's domain's second component)
  //      c(i,j) = Vector<2>(i * cos (j), i * sin (j));

  Iota<2>::Iota_t ij(c.domain());
  c.comp(0) = ij.comp(0) * cos(ij.comp(1));
  c.comp(1) = ij.comp(0) * sin(ij.comp(1));
  std::cout << c << std::endl;

  Pooma::finalize();
  return 0;
}