#include <assert.h>
#include <stdlib.h>

#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctki_GHExtensions.h"
#include "cctki_ScheduleBindings.h"
#include "cctki_WarnLevel.h"

#include "carpet.hh"

static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/Carpet/src/Initialise.cc,v 1.1 2001/07/04 12:29:47 schnetter Exp $";



namespace Carpet {
  
  using namespace std;
  
  
  
  int Initialise (tFleshConfig* fc)
  {
    DECLARE_CCTK_PARAMETERS;
    
    // Initialise stuff
    const int convlev = 0;
    cGH* const cgh = CCTK_SetupGH (fc, convlev);
    CCTKi_AddGH (fc, convlev, cgh);
    
    // Delay checkpoint until MPI has been initialised
    Checkpoint ("starting Initialise...");
    
    // Initialise stuff
    cgh->cctk_iteration = 0;
    cgh->cctk_time = cctk_initial_time;
    
    // Enable storage and communtication
    CCTKi_ScheduleGHInit (cgh);
    
    // Initialise stuff
    CCTKi_InitGHExtensions (cgh);
    
    // Check parameters
    Checkpoint ("PARAMCHECK");
    CCTK_ScheduleTraverse ("CCTK_PARAMCHECK", cgh, CallFunction);
    CCTKi_FinaliseParamWarn();
    
    BEGIN_REFLEVEL_LOOP(cgh) {
      
      // Checking
      Poison (cgh, allbutlasttime);
      
      // Set up the grid
      Checkpoint ("%*sScheduling BASEGRID", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_BASEGRID", cgh, CallFunction);
      if (reflevel==0) {
	base_delta_time = cgh->cctk_delta_time;
      } else {
// 	assert (abs(cgh->cctk_delta_time - base_delta_time / reflevelfactor)
// 		< 1e-6 * base_delta_time);
	// This circumvents a bug in CactusBase/Time
	cgh->cctk_delta_time = base_delta_time / reflevelfact;
      }
      
      // Set up the initial data
      Checkpoint ("%*sScheduling INITIAL", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_INITIAL", cgh, CallFunction);
      Checkpoint ("%*sScheduling POSTINITIAL", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_POSTINITIAL", cgh, CallFunction);
      
      // Recover
      Checkpoint ("%*sScheduling RECOVER_VARIABLES", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_RECOVER_VARIABLES", cgh, CallFunction);
      Checkpoint ("%*sScheduling CPINITIAL", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_CPINITIAL", cgh, CallFunction);
      
      // Poststep
      Checkpoint ("%*sScheduling POSTSTEP", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_POSTSTEP", cgh, CallFunction);
      
      // Checking
      PoisonCheck (cgh, allbutlasttime);
      
    } END_REFLEVEL_LOOP(cgh);
    
    BEGIN_REVERSE_REFLEVEL_LOOP(cgh) {
      
      // Restrict
      Restrict (cgh);
      
      // Checking
      CalculateChecksums (cgh, allbutlasttime);
      
      // Analysis
      Checkpoint ("%*sScheduling ANALYSIS", 2*reflevel, "");
      CCTK_ScheduleTraverse ("CCTK_ANALYSIS", cgh, CallFunction);
      
      // Output
      Checkpoint ("%*sOutputGH", 2*reflevel, "");
      CCTK_OutputGH (cgh);
      
      // Checking
      CheckChecksums (cgh, allbutlasttime);
      
    } END_REVERSE_REFLEVEL_LOOP(cgh);
    
    Checkpoint ("done with Initialise.");
    
    return 0;
  }
  
} // namespace Carpet