#include <iostream>
#include <iomanip>
#include "2d_FE.h"
#include "readmesh.cc"

int main()
{
   mesh_initialized = read_mesh("../triangle/data.1");

   if (mesh_initialized)
   {
      for (int i=0; i<num_nodes; i++)
      {
         double x = nodes[i].p.x;
         double y = nodes[i].p.y;
         std::cout << std::setw(10) << x << " " << std::setw(10) << y << " "
                   << std::setw(12) << x*x+y*y << std::endl;
      }
   }
   else
   {
      std::cout << "Error initializing the mesh. Aborting!" << std::endl;
   }
}