grid.c File Reference

contains functions that build the computational grid More...

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <inttypes.h>
#include "global.h"
#include "fields.h"

Include dependency graph for grid.c:

Go to the source code of this file.

Functions
void	computeGridSize ()
void	allocateGrid ()

Detailed Description

contains functions that build the computational grid

Definition in file grid.c.

Function Documentation

void allocateGrid

(

)

This function allocates grid arrays.

Definition at line 81 of file grid.c.

References grid, gridBuffer, gridEndIdx, gridResolution, gridSize, gridStartIdx, MPIcoords, MPIrank, MPIsize, sdim, stopRun(), int64Vector3d::x, int64Vector3d::y, and int64Vector3d::z.

{
  int i, j, k;

  gridStartIdx =
      (struct int64Vector3d *) malloc(MPIsize *
                      sizeof(struct int64Vector3d));
  gridEndIdx =
      (struct int64Vector3d *) malloc(MPIsize *
                      sizeof(struct int64Vector3d));

  for (i = 0; i < MPIsize; i++) {
    gridStartIdx[i].x = gridSize.x * MPIcoords[i][0];
    gridStartIdx[i].y = gridSize.y * MPIcoords[i][1];
    if (sdim == 3)
      gridStartIdx[i].z = gridSize.z * MPIcoords[i][2];
    else
      gridStartIdx[i].z = 0;
    gridEndIdx[i].x = gridStartIdx[i].x + gridSize.x - 1;
    gridEndIdx[i].y = gridStartIdx[i].y + gridSize.y - 1;
    if (sdim == 3)
      gridEndIdx[i].z = gridStartIdx[i].z + gridSize.z - 1;
    else
      gridEndIdx[i].z = 0;
  }

  if (!
      (gridBuffer =
       (struct doubleVector3d *) calloc(gridSize.x * gridSize.y *
                    gridSize.z,
                    sizeof(struct doubleVector3d))))
    stopRun(106, "gridBuffer", __FILE__, __LINE__);


  for (i = 0; i < gridSize.x; i++)
    for (j = 0; j < gridSize.y; j++)
      for (k = 0; k < gridSize.z; k++) {
    grid(i, j, k).x =
        lowerGridCorner.x + gridResolution * (gridStartIdx[MPIrank].x +
                          i);
    grid(i, j, k).y =
        lowerGridCorner.y + gridResolution * (gridStartIdx[MPIrank].y +
                          j);
    grid(i, j, k).z =
        lowerGridCorner.z + gridResolution * (gridStartIdx[MPIrank].z +
                          k);
      }
}

Here is the call graph for this function:

Here is the caller graph for this function:

void computeGridSize

(

)

This function computes the sizes of the grid.

Definition at line 38 of file grid.c.

References boxVolume, csize, csizeInUnits, gfH, globalGridSize, gridI, gridJ, gridK, gridResolution, gridSize, MPIdim, MPIrank, nx, ny, nz, sdim, upperGridCorner, doubleVector3d::x, int64Vector3d::x, doubleVector3d::y, int64Vector3d::y, doubleVector3d::z, and int64Vector3d::z.

{
  gridResolution = csize * gfH;
  boxVolume = pow((gridResolution / csize) * csizeInUnits * 0.0001, 3);
  lowerGridCorner.x = 0.0;
  lowerGridCorner.y = 0.0;
  lowerGridCorner.z = 0.0;
  upperGridCorner.x = (double) nx - 1;
  upperGridCorner.y = (double) ny - 1;
  upperGridCorner.z = (double) nz - 1;


  globalGridSize.x = upperGridCorner.x - lowerGridCorner.x + 1;
  globalGridSize.y = upperGridCorner.y - lowerGridCorner.y + 1;
  globalGridSize.z = upperGridCorner.z - lowerGridCorner.z + 1;

  gridI = (int64_t) ((globalGridSize.x + 1) / gridResolution);
  gridJ = (int64_t) ((globalGridSize.y + 1) / gridResolution);
  if (sdim == 3)
    gridK = (int64_t) ((globalGridSize.z + 1) / gridResolution);
  else
    gridK = 0;

  gridI = gridI + (MPIdim[0] - gridI % MPIdim[0]);
  gridJ = gridJ + (MPIdim[1] - gridJ % MPIdim[1]);
  if (sdim == 3)
    gridK = gridK + (MPIdim[2] - gridK % MPIdim[2]);

  gridSize.x = gridI / MPIdim[0];
  gridSize.y = gridJ / MPIdim[1];
  if (sdim == 3)
    gridSize.z = gridK / MPIdim[2];
  else
    gridSize.z = 1;

  if (MPIrank == 0)
    printf("Grid size:%" PRId64 "x%" PRId64 "x%" PRId64 "\n", gridI, gridJ,
       gridK);
}

Here is the caller graph for this function: