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
|
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
|
+
-
+
|
MTT::Hybrd_Solver::static_ptr->_yz = MTT::Hybrd_Solver::static_ptr->eval(tryUi);
return MTT::Hybrd_Solver::static_ptr->_yz;
}
void
MTT::Hybrd_Solver::Solve (void)
{
static std::fstream ferr ("MTT.Hybrd_messages", ios::out | ios::trunc | ios::app);
int info;
static int input_errors;
static int user_errors;
static int convergences;
static int progress_errors;
static int limit_errors;
static int unknown_errors;
NLFunc fcn(&Hybrd_Solver::f_hybrd);
NLEqn eqn(Solver::_ui,fcn);
eqn.set_tolerance(0.000001);
eqn.set_tolerance(1.0e-3);
Solver::_ui = eqn.solve(info);
switch (info)
{
case 1:
convergences++;
break;
|
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
|
-
+
-
-
+
-
-
-
-
+
+
+
+
+
|
break;
default:
unknown_errors++;
break;
}
if (1 != info)
{
std::cerr << "input (" << input_errors << ") "
std::cerr
<< " user (" << user_errors << ") "
<< " converge (" << convergences << ") "
<< " progress (" << progress_errors << ") "
<< " limit (" << limit_errors << ")"
<< " unknown (" << unknown_errors << ")"
<< " (max error = " << std::abs(eval(_ui).max()) << ")" << std::endl;
|