main.c File Reference

contains the main simulation loop More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <mpi.h>
#include <omp.h>
#include "global.h"

Include dependency graph for main.c:

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

contains the main simulation loop

Definition in file main.c.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

This function intializes MPI, calls Timothy initialization and allocation functions. It also contains the main simulation loop where all important simulation steps are called.

Definition at line 40 of file main.c.

References cellsCleanup(), commCleanup(), computePotential(), computePotentialGradient(), createExportList(), decompositionExecute(), decompositionFinalize(), fieldsSolve(), initCellsToGridExchange(), interpolateFieldsToCells(), ioWriteFields(), ioWriteStepPovRay(), ioWriteStepVTK(), MPIrank, MPIsize, nsteps, OMPthreads, povout, printStepNum(), randomStreamFree(), rstOutStep, saveRstFile(), secondsPerStep, simStart, simTime, simulationInit(), statisticsPrint(), statOutStep, step, treeBuild(), treeFree(), updateCellPositions(), updateCellStates(), vnfout, vtkout, vtkOutStep, and waitCellsToGridExchange().

{

  MPI_Init(&argc, &argv);
  MPI_Comm_size(MPI_COMM_WORLD, &MPIsize);
  MPI_Comm_rank(MPI_COMM_WORLD, &MPIrank);

  OMPthreads = omp_get_max_threads();

  simulationInit(argc, argv);

  for (step = 0; step < nsteps; step++) {

    if (!(step % statOutStep))
      printStepNum();

    decompositionExecute();
    treeBuild();
    createExportList();
    initCellsToGridExchange();
    computePotential();
    computePotentialGradient();

    waitCellsToGridExchange();
    fieldsSolve();
    interpolateFieldsToCells();

    if (!(step % statOutStep))
      statisticsPrint();

    if (simStart)
      simTime += secondsPerStep / 3600.0;   /* biological process time in hours */

    if (!(step % vtkOutStep)) {
      if (vtkout)
    ioWriteStepVTK(step);
      if (povout)
    ioWriteStepPovRay(step, 0);
      if (vnfout)
    ioWriteFields(step);
    }

    updateCellPositions();
    updateCellStates();
    commCleanup();
    treeFree();

    if (!(step % rstOutStep))
      saveRstFile();
  }

  MPI_Barrier(MPI_COMM_WORLD);

  decompositionFinalize();
  randomStreamFree();
  cellsCleanup();

  if (MPIrank == 0)
    printf("\nEnd of simulation run.\n");

  MPI_Finalize();

  return 0;
}

Here is the call graph for this function: