xgmain.cpp

来自「pic 模拟程序！面向对象」· C++ 代码 · 共 66 行

extern "C++" void Physics(void);
extern "C++" void Maintain_Diagnostics();
extern "C" void XGMainLoop();
extern "C" void Quit(void);
#include "ovector.h"
#include "oopiclist.h"
#include "diagn.h"
#include "oops.h"
extern oopicList<Diagnostics> *theDiagnostics;

#ifdef TIMING
extern "C" {
#include <sys/time.h>
}
double Time0=-1;
long int TimeIndex = 0;
#endif

void XGMainLoop()
{

#ifdef TIMING
  struct timeval TheTime;
  double TimeUsed;

  gettimeofday(&TheTime,0);

  // initialize the first time
  if(Time0==-1)  Time0 = 1.0e-6*(double)TheTime.tv_usec + (double)TheTime.tv_sec;

  TimeUsed =  1.0e-6*(double)TheTime.tv_usec + (double)TheTime.tv_sec - Time0;
  
  fprintf(stderr,"\n%ld %f",TimeIndex,(Scalar)TimeUsed);
  TimeIndex++;

#endif //TIMING

  // We have to catch exceptions here, because XGrafix can't do it
	 
  try {
   Physics();  // fields then particles
	  {
	  oopicListIter <Diagnostics> Diter(*theDiagnostics);
	  for(Diter.restart();!Diter.Done();Diter++)
		 Diter.current()->UpdatePreDiagnostics();
	}
	  {
	  oopicListIter <Diagnostics> Diter(*theDiagnostics);
	  for(Diter.restart();!Diter.Done();Diter++)
		 Diter.current()->UpdatePostDiagnostics();
	}
  }
  catch(Oops& oopsPhys){
    TXSTD::cout << "\nAn error occurred while processing XGMainLoop(). "
		<< "Back trace follows.\n"
		<< oopsPhys.getMessage() << endl;
    Quit();
    exit(1);
  }


}

void taskYield()
{}