/*
  This is a demo for a NG Appication programme
 */


#include <iostream>
using namespace std;

// for tcltk ...
#include "../libsrc/include/incvis.hpp"

// netgen interface
#include <nginterface.h>




int DA_ChooseMe (ClientData clientData,
		 Tcl_Interp * interp,
		 int argc, tcl_const char *argv[])
{
  cout << "Hi" << endl;
  return TCL_OK;
}

int DA_PrintMesh (ClientData clientData,
		  Tcl_Interp * interp,
		  int argc, tcl_const char *argv[])
{
  int i;

  int np = Ng_GetNP();
  int ne = Ng_GetNE();

  cout << "Points: " << np << endl;
  cout << "Tets:   " << ne << endl;

  double point[3];
  int tet[4];
  
  for (i = 1; i <= np; i++)
    {
      Ng_GetPoint (i, point);
      cout << "Point " << i << ": ";
      cout << point[0] << " " << point[1] << " " 
	   << point[2] << endl;
    }

  for (i = 1; i <= ne; i++)
    {
      Ng_GetElement (i, tet);
      cout << "Tet " << i << ": ";
      cout << tet[0] << " " << tet[1] << " " 
	   << tet[2] << " " << tet[3] << endl;
    }

  return TCL_OK;
}


int DA_SetSolution (ClientData clientData,
		    Tcl_Interp * interp,
		    int argc, tcl_const char *argv[])
{
  int i;

  int np = Ng_GetNP();

  double point[3];
  
  double * sol = new double[np];

  for (i = 1; i <= np; i++)
    {
      Ng_GetPoint (i, point);
      sol[i-1] = point[0];
    }

  
  Ng_SolutionData soldata;
  Ng_InitSolutionData (&soldata);

  soldata.data = sol;
  soldata.name = "My Solution";

  Ng_SetSolutionData (&soldata);

  return TCL_OK;
}



// initialize Tcl commands
int DemoApp_Init (Tcl_Interp * interp)
{

  Tcl_CreateCommand (interp, "DA_ChooseMe", DA_ChooseMe,
		     (ClientData)NULL,
		     (Tcl_CmdDeleteProc*) NULL);

  Tcl_CreateCommand (interp, "DA_PrintMesh", DA_PrintMesh,
		     (ClientData)NULL,
		     (Tcl_CmdDeleteProc*) NULL);

  Tcl_CreateCommand (interp, "DA_SetSolution", DA_SetSolution,
		     (ClientData)NULL,
		     (Tcl_CmdDeleteProc*) NULL);

  return TCL_OK;
}