11#include < iostream>
2- #include < vector>
32
43#include " gsl_wrapper.h"
54#include " recon.h"
@@ -25,18 +24,18 @@ void example_VDI2048_1(bool increase_var = false, bool keep_free = true) {
2524 system.add_var (" m_D" , 2.092 , 0.272 );
2625
2726 if (keep_free) {
28- system.add_var (" m_FD1" , 0 , -1 );
29- system.add_var (" m_FD2" , 0 , -1 );
30- system.add_var (" m_FD3" , 0 , -1 );
31- system.add_var (" m_HDAnz" , 0 , -1 );
27+ system.add_var (" m_FD1" , 0 , -1 );
28+ system.add_var (" m_FD2" , 0 , -1 );
29+ system.add_var (" m_FD3" , 0 , -1 );
30+ system.add_var (" m_HDAnz" , 0 , -1 );
3231 }
3332
3433 system.add_covariance_coeff (" m_FDKeI" , " m_FDKeII" , 0.2 );
3534 system.add_covariance_coeff (" m_SpI" , " m_SpII" , 0.4 );
3635
3736 if (increase_var) {
38- system.change_var (" m_FDKeI" , 46.241 , 2.500 );
39- system.change_var (" m_FDKeII" , 45.668 , 2.500 );
37+ system.change_var (" m_FDKeI" , 46.241 , 2.500 );
38+ system.change_var (" m_FDKeII" , 45.668 , 2.500 );
4039 }
4140
4241 // get and print the covariance matrix
@@ -54,25 +53,25 @@ void example_VDI2048_1(bool increase_var = false, bool keep_free = true) {
5453 // M_FD2 - M_FD3 = 0
5554 // M_HDAnz - M_HDNK = 0
5655 auto F_with_free = []() -> matrix {
57- double Fc[][15 ] = {
58- {1 , 1 , 0 , 0 , -0.2 , 0 , 0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 , 0 },
59- {0 , 0 , 1 , 1 , -0.6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 },
60- {0 , 0 , 0 , 0 , 0.4 , 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , -1 , 0 },
61- {0 , 0 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 , -1 },
62- {0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , -1 , 0 , 0 },
63- {0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , -1 , 0 },
64- {0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 , 0 , 0 , 1 },
65- };
66- return matrix (7 , 15 , Fc);
56+ double Fc[][15 ] = {
57+ {1 , 1 , 0 , 0 , -0.2 , 0 , 0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 , 0 },
58+ {0 , 0 , 1 , 1 , -0.6 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 },
59+ {0 , 0 , 0 , 0 , 0.4 , 1 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , -1 , 0 },
60+ {0 , 0 , 0 , 0 , 0.0 , 0 , 1 , 1 , 1 , 0 , 0 , 0 , 0 , 0 , -1 },
61+ {0 , 0 , 0 , 0 , 0.0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , -1 , 0 , 0 },
62+ {0 , 0 , 0 , 0 , 0.0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 1 , -1 , 0 },
63+ {0 , 0 , 0 , 0 , 0.0 , 0 , 0 , 0 , 0 , -1 , 0 , 0 , 0 , 0 , 1 },
64+ };
65+ return matrix (7 , 15 , Fc);
6766 };
6867
6968 auto F_without_free = []() -> matrix {
70- double Fc[][11 ] = {
71- {1 , 1 , -1 , -1 , 0.4 , 0 , 0 , 0 , 0 , 0 , 0 },
72- {0 , 0 , 1 , 1 , -1 , -1 , -1 , -1 , -1 , 0 , 0 },
73- {0 , 0 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , -1 , 0 },
74- };
75- return matrix (3 , 11 , Fc);
69+ double Fc[][11 ] = {
70+ {1 , 1 , -1 , -1 , 0.4 , 0 , 0 , 0 , 0 , 0 , 0 },
71+ {0 , 0 , 1 , 1 , -1 , -1 , -1 , -1 , -1 , 0 , 0 },
72+ {0 , 0 , 0 , 0 , 0 , 0 , 1 , 1 , 1 , -1 , 0 },
73+ };
74+ return matrix (3 , 11 , Fc);
7675 };
7776
7877 matrix F = keep_free ? F_with_free () : F_without_free ();
@@ -137,9 +136,8 @@ void example_VDI2048_1(bool increase_var = false, bool keep_free = true) {
137136 PRINT_TITLE (" Updated variables and confidences" );
138137 auto updated_sys = system.updated (x_new, confint_new);
139138 updated_sys.print_vars ();
140-
141139}
142140
143141void example_VDI2048 () {
144142 example_VDI2048_1 (true , true );
145- }
143+ }
0 commit comments