gravtree_forcetest.c File Reference

routines for direct summation forces More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <float.h>
#include <mpi.h>
#include "allvars.h"
#include "proto.h"

Go to the source code of this file.

Functions
void	gravity_forcetest (void)

---

Detailed Description

routines for direct summation forces

The code in this file allows to compute checks of the force accuracy by an independent direct summation computation. To this end, one can instruct GADGET in the Makefile to coompute direct summation forces for a certain random subfraction of particles.

Definition in file gravtree_forcetest.c.

---

Function Documentation

void gravity_forcetest (  void   )

This routine does the test of the gravitational tree force by computing the force for a random subset of particles with direct summation. Definition at line 28 of file gravtree_forcetest.c. References gravdata_in::Acc, All, global_data_all_processes::BunchSizeForce, global_data_all_processes::ComovingIntegrationOn, endrun(), Exportflag, FdForceTest, FdTimings, force_treeevaluate_direct(), global_data_all_processes::G, get_random_number(), grav_tree_compare_key(), particle_data::GravAccel, particle_data::GravAccelDirect, GravDataGet, GravDataIn, GravDataIndexTable, GravDataOut, GravDataResult, particle_data::GravPM, global_data_all_processes::Hubble, gravdata_index::Index, NTask, NumForceUpdate, NumPart, gravdata_in::OldAcc, particle_data::OldAcc, global_data_all_processes::Omega0, global_data_all_processes::OmegaLambda, global_data_all_processes::OutputDir, P, global_data_all_processes::PM_Ti_endstep, gravdata_in::Pos, particle_data::Pos, second(), set_softenings(), gravdata_index::SortIndex, TAG_DIRECT_A, TAG_DIRECT_B, gravdata_index::Task, ThisTask, global_data_all_processes::Ti_Current, particle_data::Ti_endstep, global_data_all_processes::Time, timediff(), TimeOfLastTreeConstruction, global_data_all_processes::TotNumPart, particle_data::Type, gravdata_in::u, and gravdata_in::w. Referenced by compute_accelerations().

---