enzo-3.0 / src / enzo / grid / particles / Grid_DebugActiveParticles.C

The active_particles branch has multiple heads

/***********************************************************************
/
/  GRID CLASS (DEBUG ACTIVE PARTICLE DATA)
/
/  written by: John Wise
/  date:       Februrary, 2012
/  modified1:  
/
/  PURPOSE:
/
************************************************************************/
#include "preincludes.h"

#include "ErrorExceptions.h"
#include "macros_and_parameters.h"
#include "typedefs.h"
#include "global_data.h"
#include "Fluxes.h"
#include "GridList.h"
#include "ExternalBoundary.h"
#include "Grid.h"
#include "Hierarchy.h"
#include "TopGridData.h"
#include "fortran.def"
#include "CosmologyParameters.h"

#include "ActiveParticle.h"

int grid::DebugActiveParticles(int level)
{

  if (ProcessorNumber != MyProcessorNumber || NumberOfActiveParticles == 0)
    return SUCCESS;

  int i, inside;
  FLOAT *pos;

  /* Check if the active particles are within the grid */

  for (i = 0; i < NumberOfActiveParticles; i++) {
    pos = this->ActiveParticles[i]->ReturnPosition();
    inside = this->PointInGrid(pos);
    if (inside == FALSE) {
      printf("Active particle outside grid!  level %d, grid %d\n", level, this->ID);
      printf("\t pos       = %"GOUTSYM " %"GOUTSYM " %"GOUTSYM "\n", pos[0], pos[1], pos[2]);
      printf("\t left edge = %"GOUTSYM " %"GOUTSYM " %"GOUTSYM "\n",
	     GridLeftEdge[0], GridLeftEdge[1], GridLeftEdge[2]);
      printf("\t right edge = %"GOUTSYM " %"GOUTSYM " %"GOUTSYM "\n",
	     GridRightEdge[0], GridRightEdge[1], GridRightEdge[2]);
      ENZO_FAIL("");
    }

  }

  return SUCCESS;

}
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.