-
Notifications
You must be signed in to change notification settings - Fork 2
/
ChomboDriver.cpp
66 lines (55 loc) · 1.52 KB
/
ChomboDriver.cpp
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
/*
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
* Copyright (c) 2019, Lawrence Livermore National Security, LLC.
*
* Produced at the Lawrence Livermore National Laboratory
*
* LLNL-CODE-746361
*
* All rights reserved. See COPYRIGHT for details.
*
* This file is part of the GEOSX Simulation Framework.
*
* GEOSX is a free software; you can redistribute it and/or modify it under
* the terms of the GNU Lesser General Public License (as published by the
* Free Software Foundation) version 2.1 dated February 1999.
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*/
#include "coupler.H"
int stop( int, void * ){ abort(); return 0;}
int main( int argc, char * argv[] )
{
MPI_Init( &argc, &argv );
int rank, N;
MPI_Comm_rank( MPI_COMM_WORLD, &rank );
MPI_Comm_size( MPI_COMM_WORLD, &N );
H5Eset_auto ( H5E_DEFAULT, stop, NULL );
#if 0
int VCOUNT = 0;
int QCOUNT = 0;
int * VID = nullptr;
P * x = nullptr;
V * v = nullptr;
double * pressure = nullptr;
int * QID = nullptr;
int * quads = nullptr;
double dt;
if( rank == 0 )
{
readBoundaryFile( MPI_COMM_NULL, "GEOSboundary.hdf5", dt, 0, VCOUNT, VID, x,
v, 0, QCOUNT, quads, QID, pressure );
}
for( int i = 0; i < QCOUNT; ++i )
{
pressure[i] += 0.055;
}
dt = 1.155;
if( rank==0 )
{
writeBoundaryFile( MPI_COMM_NULL, "GEOSboundary.hdf5", dt, 0, VCOUNT, x, v,
0, QCOUNT, pressure );
}
#endif
MPI_Finalize();
return 0;
}