#include <vtkm/Matrix.h>

#include <vtkm/testing/Testing.h>

namespace
{

VTKM_CONT
void BuildMatrix()
{
  std::cout << "Building matrix containing " << std::endl
            << "|  0  1  2 |" << std::endl
            << "| 10 11 12 |" << std::endl;

  ////
  //// BEGIN-EXAMPLE BuildMatrix.cxx
  ////
  vtkm::Matrix<vtkm::Float32, 2, 3> matrix;

  // Using parenthesis notation.
  matrix(0, 0) = 0.0f;
  matrix(0, 1) = 1.0f;
  matrix(0, 2) = 2.0f;

  // Using bracket notation.
  matrix[1][0] = 10.0f;
  matrix[1][1] = 11.0f;
  matrix[1][2] = 12.0f;
  ////
  //// END-EXAMPLE BuildMatrix.cxx
  ////

  vtkm::Vec2f_32 termVec(1.0f, 0.1f);
  vtkm::Vec3f_32 multVec = vtkm::MatrixMultiply(termVec, matrix);
  //  std::cout << multVec << std::endl;
  VTKM_TEST_ASSERT(test_equal(multVec, vtkm::make_Vec(1.0, 2.1, 3.2)),
                   "Unexpected product.");
}

void Run()
{
  BuildMatrix();
}

} // anonymous namespace

int Matrix(int argc, char* argv[])
{
  return vtkm::testing::Testing::Run(Run, argc, argv);
}
