predict_sub.c
b87f34a0ce7emaster

predict_sub.c
View Options

			#include <stdio.h>
			#include <stdlib.h>
			#include <string.h>
			#include <math.h>
			#include <mpi.h>
			#include <gsl/gsl_math.h>

			#include "allvars.h"
			#include "proto.h"


			#ifdef PY_INTERFACE

			/*! This function makes sure that all particle coordinates (Pos) are
			* periodically mapped onto the interval [0, BoxSize]. After this function
			* has been called, a new domain decomposition should be done, which will
			* also force a new tree construction.
			*/
			#ifdef PERIODIC
			void do_box_wrappingQ(void)
			{
			int i, j;
			double boxsize[3];

			for(j = 0; j < 3; j++)
			boxsize[j] = All.BoxSize;

			#ifdef LONG_X
			boxsize[0] *= LONG_X;
			#endif
			#ifdef LONG_Y
			boxsize[1] *= LONG_Y;
			#endif
			#ifdef LONG_Z
			boxsize[2] *= LONG_Z;
			#endif

			for(i = 0; i < NumPartQ; i++)
			for(j = 0; j < 3; j++)
			{
			while(Q[i].Pos[j] < 0)
			Q[i].Pos[j] += boxsize[j];

			while(Q[i].Pos[j] >= boxsize[j])
			Q[i].Pos[j] -= boxsize[j];
			}
			}
			#endif


			#endif /PY_INTERFACE/