/
om_pm_interface.cpp
123 lines (101 loc) · 4.01 KB
/
om_pm_interface.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
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
/*
* This file is part of OpenModelica.
*
* Copyright (c) 1998-CurrentYear, Linköping University,
* Department of Computer and Information Science,
* SE-58183 Linköping, Sweden.
*
* All rights reserved.
*
* THIS PROGRAM IS PROVIDED UNDER THE TERMS OF GPL VERSION 3
* AND THIS OSMC PUBLIC LICENSE (OSMC-PL).
* ANY USE, REPRODUCTION OR DISTRIBUTION OF THIS PROGRAM CONSTITUTES RECIPIENT'S
* ACCEPTANCE OF THE OSMC PUBLIC LICENSE.
*
* The OpenModelica software and the Open Source Modelica
* Consortium (OSMC) Public License (OSMC-PL) are obtained
* from Linköping University, either from the above address,
* from the URLs: http://www.ida.liu.se/projects/OpenModelica or
* http://www.openmodelica.org, and in the OpenModelica distribution.
* GNU version 3 is obtained from: http://www.gnu.org/copyleft/gpl.html.
*
* This program is distributed WITHOUT ANY WARRANTY; without
* even the implied warranty of MERCHANTABILITY or FITNESS
* FOR A PARTICULAR PURPOSE, EXCEPT AS EXPRESSLY SET FORTH
* IN THE BY RECIPIENT SELECTED SUBSIDIARY LICENSE CONDITIONS
* OF OSMC-PL.
*
* See the full OSMC Public License conditions for more details.
*
*/
/*
Mahder.Gebremedhin@liu.se 2020-10-12
*/
#include <iostream>
#include "om_pm_interface.hpp"
#include "om_pm_model.hpp"
extern "C" {
using namespace openmodelica::parmodelica;
typedef Equation::FunctionType FunctionType;
PMTimer seq_ode_timer;
void* PM_Model_create(const char* model_name, DATA* data, threadData_t* threadData) {
OMModel* pm_om_model = new OMModel(model_name);
pm_om_model->data = data;
pm_om_model->threadData = threadData;
pm_om_model->max_num_threads = tbb::this_task_arena::max_concurrency();
return pm_om_model;
}
void PM_Model_load_ODE_system(void* v_model, FunctionType* ode_system_funcs) {
OMModel& model = *(static_cast<OMModel*>(v_model));
model.ode_system_funcs = ode_system_funcs;
model.load_ODE_system();
}
void PM_evaluate_ODE_system(void* v_model) {
OMModel& model = *(static_cast<OMModel*>(v_model));
model.ODE_scheduler.execute();
// pm_om_model.ODE_scheduler.execution_timer.start_timer();
// for(int i = 0; i < size; ++i)
// functionODE_systems[i](data, threadData);
// pm_om_model.ODE_scheduler.execution_timer.stop_timer();
}
void seq_ode_timer_start() {
seq_ode_timer.start_timer();
}
void seq_ode_timer_stop() {
seq_ode_timer.stop_timer();
}
void seq_ode_timer_reset() {
seq_ode_timer.reset_timer();
}
void seq_ode_timer_get_elapsed_time2() {
std::cerr << seq_ode_timer.get_elapsed_time();
}
double seq_ode_timer_get_elapsed_time() {
return seq_ode_timer.get_elapsed_time();
}
void dump_times(void* v_model) {
OMModel& model = *(static_cast<OMModel*>(v_model));
#ifdef USE_LEVEL_SCHEDULER
utility::log("") << "Using level scheduler" << std::endl;
#else
#ifdef USE_FLOW_SCHEDULER
utility::log("") << "Using flow scheduler" << std::endl;
#else
#error "please specify scheduler. See makefile"
#endif
#endif
utility::log("") << "Nr.of threads " << model.max_num_threads << std::endl;
utility::log("") << "Nr.of ODE evaluations: " << model.ODE_scheduler.total_evaluations << std::endl;
utility::log("") << "Nr.of profiling ODE Evaluations: " << model.ODE_scheduler.sequential_evaluations << std::endl;
// utility::log("") << "Total ODE evaluation time : " << model.ODE_scheduler.total_parallel_cost << std::endl;
utility::log("") << "Total ODE evaluation time : " << model.ODE_scheduler.execution_timer.get_elapsed_time()
<< std::endl;
utility::log("") << "Avg. ODE evaluation time : "
<< model.ODE_scheduler.execution_timer.get_elapsed_time() /
model.ODE_scheduler.parallel_evaluations
<< std::endl;
utility::log("") << "Total ODE loading time: " << model.load_system_timer.get_elapsed_time() << std::endl;
utility::log("") << "Total ODE Clustering time: " << model.ODE_scheduler.clustering_timer.get_elapsed_time()
<< std::endl;
}
} // extern "C"