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
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
|
/* Copyright 2019 Weiqun Zhang
*
* This file is part of WarpX.
*
* License: BSD-3-Clause-LBNL
*/
#include <algorithm>
#include "WarpXParser.H"
WarpXParser::WarpXParser (std::string const& func_body)
{
define(func_body);
}
void
WarpXParser::define (std::string const& func_body)
{
clear();
m_expression = func_body;
m_expression.erase(std::remove(m_expression.begin(),m_expression.end(),'\n'),
m_expression.end());
std::string f = m_expression + "\n";
#ifdef _OPENMP
int nthreads = omp_get_max_threads();
m_variables.resize(nthreads);
m_parser.resize(nthreads);
m_parser[0] = wp_c_parser_new(f.c_str());
#pragma omp parallel
{
int tid = omp_get_thread_num();
if (tid > 0) {
m_parser[tid] = wp_parser_dup(m_parser[0]);
}
}
#else
m_parser = wp_c_parser_new(f.c_str());
#endif
}
WarpXParser::~WarpXParser ()
{
clear();
}
void
WarpXParser::clear ()
{
m_expression.clear();
#ifdef _OPENMP
if (!m_parser.empty())
{
#pragma omp parallel
{
wp_parser_delete(m_parser[omp_get_thread_num()]);
}
}
m_parser.clear();
m_variables.clear();
#else
if (m_parser) wp_parser_delete(m_parser);
m_parser = nullptr;
#endif
}
void
WarpXParser::registerVariable (std::string const& name, amrex::Real& var)
{
// We assume this is called inside OMP parallel region
#ifdef _OPENMP
wp_parser_regvar(m_parser[omp_get_thread_num()], name.c_str(), &var);
#else
wp_parser_regvar(m_parser, name.c_str(), &var);
#endif
}
void
WarpXParser::registerVariables (std::vector<std::string> const& names)
{
#ifdef _OPENMP
// This must be called outside OpenMP parallel region.
#pragma omp parallel
{
const int tid = omp_get_thread_num();
struct wp_parser* p = m_parser[tid];
auto& v = m_variables[tid];
for (int j = 0; j < names.size(); ++j) {
wp_parser_regvar(p, names[j].c_str(), &(v[j]));
}
}
#else
for (int j = 0; j < names.size(); ++j) {
wp_parser_regvar(m_parser, names[j].c_str(), &(m_variables[j]));
}
#endif
}
void
WarpXParser::setConstant (std::string const& name, amrex::Real c)
{
#ifdef _OPENMP
bool in_parallel = omp_in_parallel();
// We don't know if this is inside OMP parallel region or not
#pragma omp parallel if (!in_parallel)
{
wp_parser_setconst(m_parser[omp_get_thread_num()], name.c_str(), c);
}
#else
wp_parser_setconst(m_parser, name.c_str(), c);
#endif
}
void
WarpXParser::print () const
{
#ifdef _OPENMP
#pragma omp critical(warpx_parser_pint)
wp_ast_print(m_parser[omp_get_thread_num()]->ast);
#else
wp_ast_print(m_parser->ast);
#endif
}
std::string const&
WarpXParser::expr () const
{
return m_expression;
}
std::set<std::string>
WarpXParser::symbols () const
{
std::set<std::string> results;
#ifdef _OPENMP
wp_ast_get_symbols(m_parser[omp_get_thread_num()]->ast, results);
#else
wp_ast_get_symbols(m_parser->ast, results);
#endif
return results;
}
|
