HARM
harm and utilities
 All Data Structures Files Functions Variables Typedefs Macros Pages
Functions
initbase.gridsectioning.c File Reference

General initialization of code related to grid sectioning. More...

#include "decs.h"
Include dependency graph for initbase.gridsectioning.c:

Go to the source code of this file.

Functions

static int check_limitsinbox (FTYPE rlo, FTYPE rhi, int iloabs, int ihiabs, int doreport)
 check whether chosen radii and computed indices suggest should stop calculation More...
 
int findindexfromradius (FTYPE xr, int *xcpupos1, int *xi)
 
int compute_numcompzones (int(*sectiondef)[NDIM], long long int *localnumcompzones)
 

Detailed Description

General initialization of code related to grid sectioning.

USAGEMARK: on using grid sectioning: 1) findandsetactivesection() calls two problem dependent functions. They specify limits of integration range in terms of absolute grid indicies and the update period. 2) Big jumps issue when exposing new cells. If jump is more than MAXBND, then cells are defined by (e.g.) t=0 data instead of being defined by bounds.c because jumped more than MAXBND and bounds only does MAXBND cells: Solution: User should ensure that jumps are <MAXBND or bound region outside standard boundary region in a reasonable way. 3) Jon version of OUTFLOW – field divb problem when adding grid cells since outflowed field doesn't satisfy divb=0. Assume not adding and if do then bounds were written specially so added grid cells will have divb=0.

Definition in file initbase.gridsectioning.c.

Function Documentation

int check_limitsinbox ( FTYPE  rlo,
FTYPE  rhi,
int  iloabs,
int  ihiabs,
int  doreport 
)
static

check whether chosen radii and computed indices suggest should stop calculation

Definition at line 962 of file initbase.gridsectioning.c.

int compute_numcompzones ( int(*)  sectiondef[NDIM],
long long int *  localnumcompzones 
)

Definition at line 499 of file initbase.gridsectioning.c.

int findindexfromradius ( FTYPE  xr,
int *  xcpupos1,
int *  xi 
)
Returns
0 on success Somewhat user-dependent since in general might want to find an arbitrary grid index from an arbitrary physical position Notes: 1) BCtype could change in time, so don't assume fixed – ok so far

Definition at line 260 of file initbase.gridsectioning.c.

Here is the call graph for this function: