blob: 96d061f50496259b245ba68c607dbc081326c91d (
plain)
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
|
#include <cassert>
#include <cstdlib>
#include <iomanip>
#include <limits>
#include <cctk.h>
#include <cctk_Parameters.h>
#include <dist.hh>
#include <carpet.hh>
#include <TimerNode.hh>
namespace Carpet {
using namespace std;
int CarpetMultiModelStartup()
{
DECLARE_CCTK_PARAMETERS;
// Increase the default output precision, so that all relevant
// digits are displayed. (The C++ output streams are mostly used
// for debug messages.)
int const precision = numeric_limits<CCTK_REAL>::digits10;
cout << setprecision(precision);
cerr << setprecision(precision);
comm_universe = MPI_COMM_WORLD;
SplitUniverse (comm_universe, model, comm_world, true);
dist::pseudoinit (comm_world);
return 0;
}
int CarpetStartup()
{
CCTK_RegisterBanner ("AMR driver provided by Carpet");
GHExtension = CCTK_RegisterGHExtension("Carpet");
CCTK_RegisterGHExtensionSetupGH (GHExtension, SetupGH);
CCTK_OverloadInitialise (Initialise);
CCTK_OverloadEvolve (Evolve);
CCTK_OverloadShutdown (Shutdown);
CCTK_OverloadOutputGH (OutputGH);
CCTK_OverloadSyncGroupsByDirI (SyncGroupsByDirI);
CCTK_OverloadEnableGroupStorage (EnableGroupStorage);
CCTK_OverloadDisableGroupStorage (DisableGroupStorage);
CCTK_OverloadGroupStorageIncrease (GroupStorageIncrease);
CCTK_OverloadGroupStorageDecrease (GroupStorageDecrease);
CCTK_OverloadEnableGroupComm (EnableGroupComm);
CCTK_OverloadDisableGroupComm (DisableGroupComm);
CCTK_OverloadBarrier (Barrier);
// CCTK_OverloadNamedBarrier (NamedBarrier);
CCTK_OverloadExit ((int (*)(cGH*, int)) Exit);
CCTK_OverloadAbort ((int (*)(cGH*, int)) Abort);
CCTK_OverloadMyProc (MyProc);
CCTK_OverloadnProcs (nProcs);
CCTK_OverloadArrayGroupSizeB (ArrayGroupSizeB);
CCTK_OverloadQueryGroupStorageB (QueryGroupStorageB);
CCTK_OverloadGroupDynamicData (GroupDynamicData);
// This must happen before any Timer objects are created
main_timer_tree.root = new TimerNode(&main_timer_tree,"main");
main_timer_tree.current = 0; // No timer has been started yet
main_timer_tree.root->start();
mode_timer_tree.root = new TimerNode(&mode_timer_tree,"meta mode");
mode_timer_tree.current = 0; // No timer has been started yet
mode_timer_tree.root->start();
return 0;
}
} // namespace Carpet
|
