netgen/ng/ngappinit.cpp

421 lines
9.3 KiB
C++
Raw Normal View History

2009-01-13 04:40:13 +05:00
/*
The main function of netgen.
This file is a modification of tkAppInit.c from the Tcl/Tk package
*/
#include <mystdlib.h>
#include "incvis.hpp"
#include <meshing.hpp>
2009-01-19 23:19:17 +05:00
2011-06-26 13:35:08 +06:00
namespace netgen
{
int id = 0, ntasks = 1;
}
2009-01-13 04:40:13 +05:00
#ifdef PARALLEL
#include <mpi.h>
2011-07-04 18:29:18 +06:00
extern void ParallelRun();
namespace netgen
{
MPI_Comm mesh_comm;
}
2009-01-13 04:40:13 +05:00
#endif
2009-01-13 04:40:13 +05:00
namespace netgen
{
#include "writeuser.hpp"
extern string ngdir;
2011-01-11 01:18:01 +05:00
Flags parameters;
2009-01-13 04:40:13 +05:00
}
using netgen::parameters;
using netgen::ngdir;
using netgen::verbose;
2009-01-25 17:35:25 +05:00
using netgen::Array;
2009-01-13 04:40:13 +05:00
using netgen::RegisterUserFormats;
/*
* The following variable is a special hack that is needed in order for
* Sun shared libraries to be used for Tcl.
*/
2009-09-09 02:37:33 +06:00
2009-01-13 04:40:13 +05:00
// extern "C" int matherr();
2009-09-09 02:37:33 +06:00
// int *tclDummyMathPtr = (int *) matherr;
2009-01-13 04:40:13 +05:00
extern "C" int Ng_ServerSocketManagerInit (int port);
extern "C" int Ng_ServerSocketManagerRun (void);
2011-07-15 03:36:19 +06:00
2009-01-13 04:40:13 +05:00
bool nodisplay = false;
bool shellmode = false;
/*
*
* The Netgen main function
*
*/
int main(int argc, char ** argv)
{
2012-02-12 08:25:26 +06:00
2009-01-13 04:40:13 +05:00
#ifdef PARALLEL
2012-06-16 18:09:38 +06:00
int mpi_required = MPI_THREAD_MULTIPLE;
#ifdef VTRACE
2012-06-16 18:59:03 +06:00
mpi_required = MPI_THREAD_SINGLE;
2012-06-16 18:09:38 +06:00
#endif
int mpi_provided;
2012-06-13 15:07:11 +06:00
MPI_Init_thread(&argc, &argv, mpi_required, &mpi_provided);
2011-07-15 03:36:19 +06:00
2009-01-13 04:40:13 +05:00
MPI_Comm_size(MPI_COMM_WORLD, &netgen::ntasks);
MPI_Comm_rank(MPI_COMM_WORLD, &netgen::id);
2011-07-15 03:36:19 +06:00
MPI_Comm_dup ( MPI_COMM_WORLD, &netgen::mesh_comm);
2009-01-13 04:40:13 +05:00
#endif
2011-06-26 13:35:08 +06:00
2009-01-13 04:40:13 +05:00
if ( netgen::id == 0 )
{
2009-01-16 01:41:14 +05:00
cout << "NETGEN-" << PACKAGE_VERSION << endl;
2009-01-13 04:40:13 +05:00
2011-01-11 01:18:01 +05:00
cout << "Developed by Joachim Schoeberl at" << endl
2011-02-19 03:50:58 +05:00
<< "2010-xxxx Vienna University of Technology" << endl
2011-01-11 01:18:01 +05:00
<< "2006-2010 RWTH Aachen University" << endl
<< "1996-2006 Johannes Kepler University Linz" << endl;
2009-01-13 04:40:13 +05:00
#ifdef OCCGEOMETRY
cout << "Including OpenCascade geometry kernel" << endl;
#endif
#ifdef ACIS
cout << "Including ACIS geometry kernel" << endl;
#endif
#ifdef LINUX
//feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
//feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW | FE_UNDERFLOW);
//cout << "Handle Floating Point Exceptions: " << fegetexcept() << endl;
#endif
#ifdef DEBUG
cout << "You are running the debug version !" << endl;
#endif
#ifdef PARALLEL
2012-06-13 15:07:11 +06:00
if (netgen::ntasks == 1)
{
cout << "Run parallel Netgen with 'mpirun -np xy netgen'" << endl;
}
else
{
cout << "Running MPI - parallel using "
<< netgen::ntasks << " processor"
<< ((netgen::ntasks > 1) ? "s " : " ") << endl;
cout << "MPI-version = " << MPI_VERSION << '.' << MPI_SUBVERSION << endl;
if (mpi_provided == MPI_THREAD_MULTIPLE)
cout << "multithreaded MPI is supported" << endl;
}
2009-01-13 04:40:13 +05:00
#endif
}
2011-07-04 18:29:18 +06:00
2009-01-13 04:40:13 +05:00
2012-07-06 13:49:05 +06:00
netgen::h_argc = argc;
netgen::h_argv = argv;
2009-01-13 04:40:13 +05:00
// command line arguments:
for (int i = 1; i < argc; i++)
{
if (argv[i][0] == '-')
parameters.SetCommandLineFlag (argv[i]);
else
parameters.SetFlag ("geofile", argv[i]);
}
if (getenv ("NETGENDIR") && strlen (getenv ("NETGENDIR")))
ngdir = getenv ("NETGENDIR");
else
ngdir = ".";
verbose = parameters.GetDefineFlag ("V");
if (verbose)
cout << "NETGENDIR = " << ngdir << endl;
2012-02-12 08:25:26 +06:00
2009-01-13 04:40:13 +05:00
if ( netgen::id == 0 )
{
if (parameters.StringFlagDefined ("testout"))
netgen::testout = new ofstream (parameters.GetStringFlag ("testout", "test.out"));
2009-01-13 04:40:13 +05:00
#ifdef SOCKETS
Ng_ServerSocketManagerInit(static_cast<int>(parameters.GetNumFlag("serversocket",-1)));
if(parameters.GetNumFlag("serversocket",-1) > 0 && !parameters.GetDefineFlag("display"))
nodisplay = true;
#endif
if(parameters.GetDefineFlag("batchmode"))
nodisplay = true;
2011-06-16 22:48:54 +06:00
2009-01-13 04:40:13 +05:00
if(parameters.GetDefineFlag("shellmode"))
{
nodisplay = true;
shellmode = true;
}
Tcl_FindExecutable(NULL);
// initialize application
Tcl_Interp * myinterp = Tcl_CreateInterp ();
if (Tcl_AppInit (myinterp) == TCL_ERROR)
{
cerr << "Exit Netgen due to initialization problem" << endl;
exit (1);
}
// parse tcl-script
int errcode;
2009-01-18 00:37:57 +05:00
bool internaltcl = false;
2009-01-13 04:40:13 +05:00
if (shellmode)
internaltcl = false;
if (verbose)
{
cout << "Tcl header version = " << TCL_PATCH_LEVEL << endl;
Tcl_Eval (myinterp, "puts \"Tcl runtime version = [info patchlevel] \";");
}
if (parameters.GetDefineFlag ("internaltcl"))
internaltcl=true;
if (parameters.GetDefineFlag ("externaltcl"))
internaltcl=false;
2009-01-18 00:37:57 +05:00
2009-01-13 04:40:13 +05:00
if (internaltcl)
{
if (verbose)
cout << "using internal Tcl-script" << endl;
// connect to one string
extern const char * ngscript[];
const char ** hcp = ngscript;
int len = 0;
while (*hcp)
len += strlen (*hcp++);
char * tr1 = new char[len+1];
*tr1 = 0;
hcp = ngscript;
char * tt1 = tr1;
while (*hcp)
{
strcat (tt1, *hcp);
tt1 += strlen (*hcp++);
}
errcode = Tcl_Eval (myinterp, tr1);
delete [] tr1;
}
else
{
string startfile = ngdir + "/ng.tcl";
if (verbose)
cout << "Load Tcl-script from " << startfile << endl;
errcode = Tcl_EvalFile (myinterp, (char*)startfile.c_str());
}
if (errcode)
{
cout << "Error in Tcl-Script:" << endl;
// cout << "result = " << myinterp->result << endl;
cout << "result = " << Tcl_GetStringResult (myinterp) << endl;
// cout << "in line " << myinterp->errorLine << endl;
2009-01-13 04:40:13 +05:00
cout << "\nMake sure to set environment variable NETGENDIR to directory containing ng.tcl" << endl;
2009-01-13 04:40:13 +05:00
exit (1);
}
// lookup user file formats and insert into format list:
2009-01-25 17:35:25 +05:00
Array<const char*> userformats;
2009-09-22 13:12:00 +06:00
Array<const char*> extensions;
RegisterUserFormats (userformats, extensions);
2009-01-13 04:40:13 +05:00
ostringstream fstr;
2009-09-22 13:12:00 +06:00
tcl_const char * exportft = Tcl_GetVar (myinterp, "exportfiletype", 0);
2009-01-13 04:40:13 +05:00
for (int i = 1; i <= userformats.Size(); i++)
2009-09-22 13:12:00 +06:00
{
fstr << ".ngmenu.file.filetype add radio -label \""
<< userformats.Get(i) << "\" -variable exportfiletype\n";
fstr << "lappend meshexportformats { {" << userformats.Get(i) << "} {" << extensions.Get(i) << "} }\n";
}
2009-01-13 04:40:13 +05:00
Tcl_Eval (myinterp, (char*)fstr.str().c_str());
2009-09-22 13:12:00 +06:00
Tcl_SetVar (myinterp, "exportfiletype", exportft, 0);
2009-01-13 04:40:13 +05:00
#ifdef SOCKETS
Ng_ServerSocketManagerRun();
#endif
// start event-loop
Tk_MainLoop();
Tcl_DeleteInterp (myinterp);
Tcl_Exit(0);
}
2011-07-04 18:29:18 +06:00
2009-01-13 04:40:13 +05:00
#ifdef PARALLEL
else
{
ParallelRun();
MPI_Finalize();
}
2011-07-04 18:29:18 +06:00
2009-01-13 04:40:13 +05:00
#endif
return 0;
}
2009-01-18 05:25:35 +05:00
/*
2009-01-13 04:40:13 +05:00
extern "C" int Tix_Init (Tcl_Interp * interp);
extern "C" int Itcl_Init (Tcl_Interp * interp);
extern "C" int Itk_Init (Tcl_Interp * interp);
2009-01-18 05:25:35 +05:00
*/
2009-01-13 04:40:13 +05:00
extern "C" int Ng_Init (Tcl_Interp * interp);
extern "C" int Ng_Vis_Init (Tcl_Interp * interp);
// extern Tcl_PackageInitProc * Tk_SafeInit;
/*
*
* Initialize packages
*
*/
2009-04-05 17:30:37 +06:00
// extern "C" int NGSolve_Init (Tcl_Interp * interp);
2009-01-13 04:40:13 +05:00
int Tcl_AppInit(Tcl_Interp * interp)
{
if (Tcl_Init(interp) == TCL_ERROR) {
cerr << "Problem in Tcl_Init: " << endl;
cout << "result = " << Tcl_GetStringResult (interp) << endl;
// cerr << interp->result << endl;
2009-01-13 04:40:13 +05:00
// return TCL_ERROR;
}
if (!nodisplay && Tk_Init(interp) == TCL_ERROR) {
cerr << "Problem in Tk_Init: " << endl;
cout << "result = " << Tcl_GetStringResult (interp) << endl;
// cerr << interp->result << endl;
2009-01-13 04:40:13 +05:00
// return TCL_ERROR;
}
if (Ng_Init(interp) == TCL_ERROR) {
cerr << "Problem in Ng_Init: " << endl;
cout << "result = " << Tcl_GetStringResult (interp) << endl;
// cerr << interp->result << endl;
2009-01-13 04:40:13 +05:00
// return TCL_ERROR;
}
2011-06-16 23:55:08 +06:00
2009-01-13 04:40:13 +05:00
if (!nodisplay && Ng_Vis_Init(interp) == TCL_ERROR) {
cerr << "Problem in Ng_Vis_Init: " << endl;
cout << "result = " << Tcl_GetStringResult (interp) << endl;
// cerr << interp->result << endl;
2009-01-13 04:40:13 +05:00
// return TCL_ERROR;
}
2009-04-05 17:30:37 +06:00
2011-07-04 18:29:18 +06:00
2009-01-13 04:40:13 +05:00
#ifdef TRAFO
// extern int Trafo_Init (Tcl_Interp * interp);
// if (Trafo_Init(interp) == TCL_ERROR)
// {
// cerr << "Problem in Trafo_Init: " << endl;
// cerr << interp->result << endl;
// return TCL_ERROR;
// }
#endif
#ifdef EBGELAST
extern int EBGElast_Init (Tcl_Interp * interp);
if(EBGElast_Init(interp) == TCL_ERROR)
{
cerr << "Problem in EBGElast_Init: " << endl;
cerr << interp->result << endl;
return TCL_ERROR;
}
#endif
#ifdef SMALLTRAFO
extern int SmallModels_Init (Tcl_Interp * interp);
if(SmallModels_Init(interp) == TCL_ERROR)
{
cerr << "Problem in SmallModel_Init: " << endl;
cerr << interp->result << endl;
return TCL_ERROR;
}
#endif
#ifdef SOCKETS
extern int Ng_Socket_Init (Tcl_Interp * interp);
if ( Ng_Socket_Init(interp) == TCL_ERROR)
{
cerr << "Problem in Ng_Socket_Init: " << endl;
cerr << interp->result << endl;
return TCL_ERROR;
}
#endif
#ifdef ZUGSTANGE
extern int Zugstange_Init (Tcl_Interp * interp);
if (Zugstange_Init(interp) == TCL_ERROR)
{
cerr << "Problem in Zugstange_Init: " << endl;
cerr << interp->result << endl;
return TCL_ERROR;
}
#endif
Tcl_StaticPackage(interp, "Tk", Tk_Init, 0);
return TCL_OK;
}