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