1
1
#include " example1.hpp"
2
2
3
- #include < cmath >
4
- #include < time.h >
3
+ #include < ctime >
4
+ #include < iostream >
5
5
6
6
#include < Model/PhasorDynamics/Branch/Branch.hpp>
7
7
#include < Model/PhasorDynamics/Bus/Bus.hpp>
12
12
#include < Solver/Dynamic/Ida.hpp>
13
13
#include < Utilities/Testing.hpp>
14
14
15
- #define _USE_MATH_DEFINES
16
- #define _CRT_SECURE_NO_WARNINGS
17
-
18
15
int main ()
19
16
{
20
17
using namespace GridKit ::PhasorDynamics;
21
18
using namespace AnalysisManager ::Sundials;
22
19
23
- printf ( " Example 1 version 2\n " ) ;
20
+ std::cout << " Example 1 version 2\n " ;
24
21
25
22
/* Create model parts */
26
23
SystemModel<double , size_t > sys;
@@ -104,7 +101,6 @@ int main()
104
101
/ (1.0 + std::abs (reference_solution[i / 48 ][2 ]));
105
102
if (err > error_V)
106
103
error_V = err;
107
- // std::cout << "t = " << ti << ": Vr = " << Vr << ", Vi = " << Vi << ", dw = " << dw;
108
104
std::cout << " GridKit: t = " << ti
109
105
<< " , |V| = " << std::sqrt (Vr * Vr + Vi * Vi)
110
106
<< " , w = " << (1.0 + dw) << " \n " ;
@@ -138,11 +134,8 @@ int main()
138
134
}
139
135
++j;
140
136
++i;
141
- // if (i > 500)
142
- // break;
143
137
}
144
138
145
- // std::cout << buffer.str();
146
139
int status = 0 ;
147
140
std::cout << " Max error in |V| = " << error_V << " \n " ;
148
141
if (error_V > 2e-4 )
@@ -153,5 +146,5 @@ int main()
153
146
154
147
std::cout << " \n\n Complete in " << (stop - start) / CLOCKS_PER_SEC << " seconds\n " ;
155
148
156
- return 0 ;
149
+ return status ;
157
150
}
0 commit comments