-
Notifications
You must be signed in to change notification settings - Fork 0
/
HistorySingleParameter.cpp
60 lines (51 loc) · 1.88 KB
/
HistorySingleParameter.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
#include "HistorySingleParameter.h"
std::ostream& operator<<(std::ostream& os, const HistorySingleParameter& h) {
os << "HistorySingleParameter[" << (h.size()-1) << "..0]" << std::endl;
h.print(os);
return(os);
}
/////////////////////////
///////////////////////// HistorySingleParameter methods
/////////////////////////
void HistorySingleParameter::print(std::ostream& os) const {
for (int i = (size()-1); i >= 0; --i) {
os << " [" << i << "] ";
(((HistorySingleParameter*)this)->history)[i].print(os);
}
}
void HistorySingleParameter::simple_print(std::ostream& os) const {
os << "parameter" << "\t" << "logL" << std::endl;
for (int i = (size()-1); i >= 0; --i) {
os << i << "\t";
(((HistorySingleParameter*)this)->history)[i].simple_print(os);
os << std::endl;
}
}
/////////////////////////
///////////////////////// HistorySingleParameter::HistoryEntry methods
/////////////////////////
inline void HistorySingleParameter::HistoryEntry::print(std::ostream& os) const {
os << " delta=" << delta;
os << " parameter[-delta,-0.5delta,0,+0.5delta,+delta]=["
<< parameter_m_delta << ", "
// << parameter_m_05delta << ", "
<< parameter << ", "
// << parameter_p_05delta << ", "
<< parameter_p_delta << "]";
os << " epsilon=" << epsilon;
os << " logL[-delta,-0.5delta,0,+0.5delta,+delta]=["
<< logL_m_delta << ", "
// << logL_m_05delta << ", "
<< logL << ", "
// << logL_p_05delta << ", "
<< logL_p_delta << "]";
os << " f_lo=" << f_lo;
os << " f_hi=" << f_hi;
os << " f_parameter=" << f_parameter;
os << " fp_parameter=" << fp_parameter;
os << std::endl;
}
inline void HistorySingleParameter::HistoryEntry::simple_print(
std::ostream& os) const {
os << parameter << "\t" << logL;
}