1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
/*
This file is part of p4est.
p4est is a C library to manage a collection (a forest) of multiple
connected adaptive quadtrees or octrees in parallel.
Copyright (C) 2010 The University of Texas System
Additional copyright (C) 2011 individual authors
Written by Carsten Burstedde, Lucas C. Wilcox, and Tobin Isaac
p4est is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
p4est is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with p4est; if not, write to the Free Software Foundation, Inc.,
51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
*/
/*
* This is an internal header file.
* It is included from a c file that may be preceded by p4est_to_p8est.h.
*/
#ifndef P4EST_PARTICLES_GLOBAL_H
#define P4EST_PARTICLES_GLOBAL_H
#include <sc_statistics.h>
#ifndef P4_TO_P8
#include <p4est.h>
#else
#include <p8est.h>
#endif /* P4_TO_P8 */
SC_EXTERN_C_BEGIN;
typedef enum part_stats
{
PART_STATS_FIRST = 0,
PART_STATS_NOTIFY = PART_STATS_FIRST,
PART_STATS_BINARY,
PART_STATS_NARY,
PART_STATS_COMM,
PART_STATS_WAIT1,
PART_STATS_WAIT2,
PART_STATS_PEERS,
PART_STATS_SEARCHP,
PART_STATS_SEARCHL,
PART_STATS_TRANSF,
PART_STATS_TRANSC,
PART_STATS_BUILD,
PART_STATS_PERTREEF,
PART_STATS_PERTREEB,
PART_STATS_RK,
PART_STATS_LAST
}
part_stats_t;
typedef double (*part_init_density_t) (double x, double y, double z,
void *data);
typedef struct part_global
{
sc_MPI_Comm mpicomm;
int mpisize, mpirank;
/* command line option variables */
int minlevel;
int maxlevel;
int bricklev;
int order;
int olap_notify;
#if 0
int checkp;
#endif
int printn;
int scaling;
int collapse;
double num_particles;
double elem_particles;
double deltat;
double finaltime;
const char *prefix;
/* variables derived from command line */
int ntop;
int nint;
int nbot;
int vtk;
int mpiwrap;
int build_step;
int build_part;
int build_wrap;
/* runtime variables */
int bricklength;
int stage;
int verylast;
double global_density;
double t, lxyz[3], hxyz[3], dxyz[3];
p4est_locidx_t add_count; /**< Number of particles for build_add */
sc_array_t *padata; /**< pa_data_t Numerical data of local particles */
sc_array_t *pfound; /**< int Target rank for local particles */
sc_array_t *iremain; /**< locidx_t Index into padata of stay-local particles */
p4est_locidx_t ireindex, ire2; /**< Running index into iremain */
p4est_locidx_t qremain; /**< Number of particles remaining in quadrant */
sc_array_t *ireceive; /**< Index into particle receive buffer */
p4est_locidx_t irvindex, irv2; /**< Running index into ireceive */
p4est_locidx_t qreceive; /**< Number of particles received for quadrant */
sc_array_t *recevs; /**< comm_prank_t with one entry per receiver, sorted */
sc_array_t *recv_req; /**< sc_MPI_Request for receiving */
sc_array_t *send_req; /**< sc_MPI_Request for sending */
sc_array_t *prebuf; /**< pa_data_t All received particles */
sc_array_t *cfound; /**< char Flag for received particles */
sc_hash_t *psend; /**< comm_psend_t with one entry per receiver */
sc_mempool_t *psmem; /**< comm_psend_t to use as hash table entries */
p4est_locidx_t prevlp; /**< Count local particles in partition callback */
p4est_locidx_t qcount; /**< Count local quadrants in partition callback */
sc_array_t *src_fixed; /**< int Particle counts per quadrant */
sc_array_t *dest_fixed; /**< int Particle counts per quadrant */
part_init_density_t pidense;
void *piddata;
p4est_gloidx_t gpnum, gplost;
sc_array_t *ilh[2], *jlh[2], *klh[2];
/* parallel statistics */
sc_statinfo_t si[PART_STATS_LAST];
/* the mesh */
p4est_connectivity_t *conn;
p4est_t *p4est;
}
part_global_t;
SC_EXTERN_C_END;
#endif /* !P4EST_PARTICLES_GLOBAL_H */
|
