159599516SKenneth E. Jansen #include <stdio.h>
259599516SKenneth E. Jansen #include <stdlib.h>
359599516SKenneth E. Jansen #include <string.h>
459599516SKenneth E. Jansen
5*74490262SBen Matthews #include <unistd.h>
659599516SKenneth E. Jansen #include <sys/types.h>
759599516SKenneth E. Jansen #include <sys/stat.h>
859599516SKenneth E. Jansen
959599516SKenneth E. Jansen #include <iostream>
103de22adcSCameron Smith #include <sstream>
1159599516SKenneth E. Jansen using namespace std;
1259599516SKenneth E. Jansen
13a5c8caf2SBenjamin Matthews #define OMPI_SKIP_MPICXX 1
14a5c8caf2SBenjamin Matthews #include <mpi.h>
15a5c8caf2SBenjamin Matthews
1659599516SKenneth E. Jansen #include <FCMangle.h>
1759599516SKenneth E. Jansen #define input FortranCInterface_GLOBAL_(input,INPUT)
1859599516SKenneth E. Jansen #define timer FortranCInterface_GLOBAL_(timer,TIMER)
1959599516SKenneth E. Jansen
20a5c8caf2SBenjamin Matthews #include "commonM2NFixBnd_c.h"
2159599516SKenneth E. Jansen
2259599516SKenneth E. Jansen extern "C" char phasta_iotype[80];
2359599516SKenneth E. Jansen char phasta_iotype[80];
2459599516SKenneth E. Jansen
2559599516SKenneth E. Jansen extern "C" void input();
2659599516SKenneth E. Jansen //MR CHANGE
2759599516SKenneth E. Jansen extern void setIOparam(); // For SyncIO
2859599516SKenneth E. Jansen //MR CHANGE END
2959599516SKenneth E. Jansen
3059599516SKenneth E. Jansen int myrank; /* made file global for ease in debugging */
3159599516SKenneth E. Jansen
3259599516SKenneth E. Jansen void
catchDebugger()3359599516SKenneth E. Jansen catchDebugger() {
3459599516SKenneth E. Jansen while (1) {
3559599516SKenneth E. Jansen int debuggerPresent=0;
3659599516SKenneth E. Jansen int fakeSTOP = 1; // please stop HERE and assign as next line
3759599516SKenneth E. Jansen // assign or set debuggerPresent=1
3859599516SKenneth E. Jansen if(debuggerPresent) {
3959599516SKenneth E. Jansen break;
4059599516SKenneth E. Jansen }
4159599516SKenneth E. Jansen }
4259599516SKenneth E. Jansen }
4359599516SKenneth E. Jansen
4459599516SKenneth E. Jansen // some useful debugging functions
4559599516SKenneth E. Jansen
4659599516SKenneth E. Jansen void
pdarray(void * darray,int start,int end)4759599516SKenneth E. Jansen pdarray( void* darray , int start, int end ) {
4859599516SKenneth E. Jansen for( int i=start; i < end; i++ ){
4959599516SKenneth E. Jansen cout << ((double*)darray)[i] << endl;
5059599516SKenneth E. Jansen }
5159599516SKenneth E. Jansen }
5259599516SKenneth E. Jansen
5359599516SKenneth E. Jansen void
piarray(void * iarray,int start,int end)5459599516SKenneth E. Jansen piarray( void* iarray , int start, int end ) {
5559599516SKenneth E. Jansen for( int i=start; i < end; i++ ){
5659599516SKenneth E. Jansen cout << ((int*)iarray)[i] << endl;
5759599516SKenneth E. Jansen }
5859599516SKenneth E. Jansen }
5959599516SKenneth E. Jansen
6059599516SKenneth E. Jansen extern "C" int
phasta(int argc,char * argv[])6159599516SKenneth E. Jansen phasta( int argc,
6259599516SKenneth E. Jansen char *argv[] ) {
6359599516SKenneth E. Jansen
6459599516SKenneth E. Jansen int size,ierr;
6559599516SKenneth E. Jansen char inpfilename[100];
6659599516SKenneth E. Jansen char* pauseDebugger = getenv("catchDebugger");
6759599516SKenneth E. Jansen //cout << "pauseDebugger" << pauseDebugger << endl;
6859599516SKenneth E. Jansen
6959599516SKenneth E. Jansen MPI_Init(&argc,&argv);
7059599516SKenneth E. Jansen MPI_Comm_size (MPI_COMM_WORLD, &size);
7159599516SKenneth E. Jansen MPI_Comm_rank (MPI_COMM_WORLD, &myrank);
7259599516SKenneth E. Jansen
7359599516SKenneth E. Jansen workfc.numpe = size;
7459599516SKenneth E. Jansen workfc.myrank = myrank;
7559599516SKenneth E. Jansen
7659599516SKenneth E. Jansen #if (defined WIN32)
7759599516SKenneth E. Jansen if(argc > 2 ){
7859599516SKenneth E. Jansen catchDebugger();
7959599516SKenneth E. Jansen }
8059599516SKenneth E. Jansen #endif
8159599516SKenneth E. Jansen #if (1) // ALWAYS ( defined LAUNCH_GDB ) && !( defined WIN32 )
8259599516SKenneth E. Jansen
8359599516SKenneth E. Jansen if ( pauseDebugger ) {
8459599516SKenneth E. Jansen
8559599516SKenneth E. Jansen int parent_pid = getpid();
8659599516SKenneth E. Jansen int gdb_child = fork();
8759599516SKenneth E. Jansen cout << "gdb_child" << gdb_child << endl;
8859599516SKenneth E. Jansen
8959599516SKenneth E. Jansen if( gdb_child == 0 ) {
9059599516SKenneth E. Jansen
9159599516SKenneth E. Jansen cout << "Debugger Process initiating" << endl;
923de22adcSCameron Smith std::stringstream exec_string;
9359599516SKenneth E. Jansen
9459599516SKenneth E. Jansen #if ( defined decalp )
9559599516SKenneth E. Jansen exec_string <<"xterm -e idb "
9659599516SKenneth E. Jansen << " -pid "<< parent_pid <<" "<< argv[0] << endl;
9759599516SKenneth E. Jansen #endif
9859599516SKenneth E. Jansen #if ( defined LINUX )
9959599516SKenneth E. Jansen exec_string <<"xterm -e gdb"
10059599516SKenneth E. Jansen << " -pid "<< parent_pid <<" "<< argv[0] << endl;
10159599516SKenneth E. Jansen #endif
10259599516SKenneth E. Jansen #if ( defined SUN4 )
10359599516SKenneth E. Jansen exec_string <<"xterm -e dbx "
10459599516SKenneth E. Jansen << " - "<< parent_pid <<" "<< argv[0] << endl;
10559599516SKenneth E. Jansen #endif
10659599516SKenneth E. Jansen #if ( defined IRIX )
10759599516SKenneth E. Jansen exec_string <<"xterm -e dbx "
10859599516SKenneth E. Jansen << " -p "<< parent_pid <<" "<< argv[0] << endl;
10959599516SKenneth E. Jansen #endif
1103de22adcSCameron Smith std::string s = exec_string.str();
1113de22adcSCameron Smith system( s.c_str() );
11259599516SKenneth E. Jansen exit(0);
11359599516SKenneth E. Jansen }
11459599516SKenneth E. Jansen catchDebugger();
11559599516SKenneth E. Jansen }
11659599516SKenneth E. Jansen
11759599516SKenneth E. Jansen #endif
11859599516SKenneth E. Jansen
11959599516SKenneth E. Jansen //MR CHANGE
12059599516SKenneth E. Jansen setIOparam();
12159599516SKenneth E. Jansen //MR CHANGE END
12259599516SKenneth E. Jansen
12359599516SKenneth E. Jansen input();
12459599516SKenneth E. Jansen //MR CHANGE
12559599516SKenneth E. Jansen
12659599516SKenneth E. Jansen MPI_Barrier(MPI_COMM_WORLD);
12759599516SKenneth E. Jansen if ( myrank == 0 ) {
12859599516SKenneth E. Jansen printf("phasta.cc - last call before finalize!\n");
12959599516SKenneth E. Jansen }
13059599516SKenneth E. Jansen //MR CHANGE
13159599516SKenneth E. Jansen
13259599516SKenneth E. Jansen MPI_Finalize();
13359599516SKenneth E. Jansen return 0;
13459599516SKenneth E. Jansen }
135