75 | | <a name="l00017"></a><a class="code" href="classbdm_1_1mexEpdf.html">00017</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mexEpdf.html">mexEpdf</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>{ |
76 | | <a name="l00018"></a>00018 <span class="keyword">protected</span>: |
77 | | <a name="l00019"></a>00019 <span class="keywordtype">string</span> name; |
78 | | <a name="l00020"></a>00020 mxArray *data; |
79 | | <a name="l00021"></a>00021 <span class="keyword">public</span>: |
80 | | <a name="l00022"></a>00022 <a class="code" href="classbdm_1_1mexEpdf.html">mexEpdf</a>() {}; |
81 | | <a name="l00023"></a><a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298">00023</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298" title="Load from structure with elements:.">from_setting</a>(<span class="keyword">const</span> Setting &S) { |
82 | | <a name="l00024"></a>00024 name = (<span class="keyword">const</span> <span class="keywordtype">char</span> *) S[<span class="stringliteral">"name"</span>]; |
83 | | <a name="l00025"></a>00025 <a class="code" href="classUImxConfig.html" title="Reimplementation of libconfig&#39;s Config class for Matlab mxArray structures.">UImxConfig</a> conf(S); |
84 | | <a name="l00026"></a>00026 data = mxDuplicateArray(conf.<a class="code" href="classUImxConfig.html#accf378aab12b7709d644441fd2e30d00" title="Matlab structure where the info is stored.">mxconfig</a>); |
85 | | <a name="l00027"></a>00027 <span class="comment">//mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());</span> |
86 | | <a name="l00028"></a>00028 <span class="comment">//TODO (future...):</span> |
87 | | <a name="l00029"></a>00029 <span class="comment">//mxArray * init_data = setting2mxarray S["init_data"];</span> |
88 | | <a name="l00030"></a>00030 <span class="comment">//mexCallMATLAB(1, &data, 1, &init_data, name+"_from_setting");</span> |
89 | | <a name="l00031"></a>00031 <span class="comment">//delete init_data;</span> |
90 | | <a name="l00032"></a>00032 } |
91 | | <a name="l00033"></a><a class="code" href="classbdm_1_1mexEpdf.html#a028dd3432d55e1fcc49f786a27f5afe8">00033</a> vec <a class="code" href="classbdm_1_1mexEpdf.html#a028dd3432d55e1fcc49f786a27f5afe8" title="return expected value">mean</a>()<span class="keyword"> const </span>{ |
92 | | <a name="l00034"></a>00034 mxArray *tmp; |
93 | | <a name="l00035"></a>00035 <span class="keywordtype">string</span> fname = name+<span class="stringliteral">"_mean"</span>; |
94 | | <a name="l00036"></a>00036 mexCallMATLAB(1, &tmp, 1, (mxArray **) &data, fname.c_str()); |
95 | | <a name="l00037"></a>00037 <span class="keywordflow">return</span> mxArray2vec(tmp); |
96 | | <a name="l00038"></a>00038 } |
97 | | <a name="l00039"></a>00039 }; |
98 | | <a name="l00040"></a>00040 <a class="code" href="user__info_8h.html#a4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(mexEpdf); |
99 | | <a name="l00041"></a>00041 |
100 | | <a name="l00042"></a><a class="code" href="classbdm_1_1mexBM.html">00042</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mexBM.html">mexBM</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of a system, i.e. all uncertainty is modeled by probabilities.">BM</a>{ |
101 | | <a name="l00043"></a>00043 <span class="keyword">protected</span> : |
102 | | <a name="l00044"></a>00044 <span class="keywordtype">string</span> name; |
103 | | <a name="l00045"></a>00045 <a class="code" href="classbdm_1_1mexEpdf.html">mexEpdf</a> est; |
104 | | <a name="l00046"></a>00046 mxArray *data; |
105 | | <a name="l00047"></a>00047 <span class="keyword">public</span>: |
106 | | <a name="l00048"></a>00048 <a class="code" href="classbdm_1_1mexBM.html">mexBM</a>() {} |
107 | | <a name="l00049"></a>00049 |
108 | | <a name="l00050"></a>00050 mxArray *get_data() { |
109 | | <a name="l00051"></a>00051 <span class="comment">//mexCallMATLAB(0, NULL, 1, &data, "dump");</span> |
110 | | <a name="l00052"></a>00052 <span class="keywordflow">return</span> mxDuplicateArray(data); |
111 | | <a name="l00053"></a>00053 } |
112 | | <a name="l00054"></a>00054 |
113 | | <a name="l00055"></a><a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd">00055</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd" title="This method arrange instance properties according the data stored in the Setting...">from_setting</a>(<span class="keyword">const</span> Setting &S) { |
114 | | <a name="l00056"></a>00056 <a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd" title="This method arrange instance properties according the data stored in the Setting...">BM::from_setting</a>(S); |
115 | | <a name="l00057"></a>00057 name = (<span class="keyword">const</span> <span class="keywordtype">char</span> *) S[<span class="stringliteral">"name"</span>]; |
116 | | <a name="l00058"></a>00058 <a class="code" href="classUImxConfig.html" title="Reimplementation of libconfig&#39;s Config class for Matlab mxArray structures.">UImxConfig</a> conf(S); |
117 | | <a name="l00059"></a>00059 data = mxDuplicateArray(conf.<a class="code" href="classUImxConfig.html#accf378aab12b7709d644441fd2e30d00" title="Matlab structure where the info is stored.">mxconfig</a>); |
118 | | <a name="l00060"></a>00060 <span class="comment">//string fname = name+"_new";</span> |
119 | | <a name="l00061"></a>00061 <span class="comment">//mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());</span> |
120 | | <a name="l00062"></a>00062 <span class="comment">//the following works as long as the posterior is the</span> |
121 | | <a name="l00063"></a>00063 <span class="comment">//only member object there could be a structure of </span> |
122 | | <a name="l00064"></a>00064 <span class="comment">//member objects in the setting and a for cycle over</span> |
123 | | <a name="l00065"></a>00065 <span class="comment">//all of them right here in the code</span> |
124 | | <a name="l00066"></a>00066 Setting &posterior = S[<span class="stringliteral">"posterior"</span>]; |
125 | | <a name="l00067"></a>00067 est.<a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298" title="Load from structure with elements:.">from_setting</a>(posterior); |
126 | | <a name="l00068"></a>00068 } |
127 | | <a name="l00069"></a><a class="code" href="classbdm_1_1mexBM.html#aa99452b52f731f561b01029af18c68f7">00069</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexBM.html#aa99452b52f731f561b01029af18c68f7" title="Incremental Bayes rule.">bayes</a>(<span class="keyword">const</span> vec &dt) { |
128 | | <a name="l00070"></a>00070 <span class="comment">//void bayes() {</span> |
129 | | <a name="l00071"></a>00071 mxArray *tmp, *old; |
130 | | <a name="l00072"></a>00072 mxArray *in[2]; |
131 | | <a name="l00073"></a>00073 in[0] = data; |
132 | | <a name="l00074"></a>00074 in[1] = mxCreateDoubleMatrix(dt.size(), 1, mxREAL); |
133 | | <a name="l00075"></a>00075 vec2mxArray(dt, in[1]); |
134 | | <a name="l00076"></a>00076 mexCallMATLAB(1, &tmp, 2, in, (name+<span class="stringliteral">"_bayes"</span>).c_str()); |
135 | | <a name="l00077"></a>00077 old = data; |
136 | | <a name="l00078"></a>00078 data = mxDuplicateArray(tmp); |
137 | | <a name="l00079"></a>00079 <span class="keywordflow">if</span> (old) mxDestroyArray(old); |
138 | | <a name="l00080"></a>00080 <span class="keywordflow">if</span> (tmp) mxDestroyArray(tmp); |
139 | | <a name="l00081"></a>00081 <span class="comment">//mexCallMATLAB(0, NULL, 1, &data, "dump");</span> |
140 | | <a name="l00082"></a>00082 } |
141 | | <a name="l00083"></a>00083 <span class="keyword">const</span> <a class="code" href="classbdm_1_1mexEpdf.html">mexEpdf</a>& posterior()<span class="keyword"> const </span>{ |
142 | | <a name="l00084"></a>00084 <span class="keywordflow">return</span> est; |
143 | | <a name="l00085"></a>00085 } <span class="comment">//tohle by melo zustat!!</span> |
144 | | <a name="l00086"></a>00086 <span class="keyword">const</span> mexEpdf* _e()<span class="keyword"> const </span>{ |
145 | | <a name="l00087"></a>00087 <span class="keywordflow">return</span> &est; |
146 | | <a name="l00088"></a>00088 } <span class="comment">//tohle by melo zustat!!</span> |
147 | | <a name="l00089"></a>00089 }; |
148 | | <a name="l00090"></a>00090 <a class="code" href="user__info_8h.html#a4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(mexBM); |
149 | | <a name="l00091"></a>00091 |
150 | | <a name="l00092"></a>00092 } |
| 75 | <a name="l00018"></a><a class="code" href="classbdm_1_1mexEpdf.html">00018</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mexEpdf.html" title="epdf with functions implemented in matlab">mexEpdf</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>{ |
| 76 | <a name="l00019"></a>00019 <span class="keyword">protected</span>: |
| 77 | <a name="l00021"></a><a class="code" href="classbdm_1_1mexEpdf.html#a287414259db205e5bc82c0023c09e633">00021</a> <span class="keywordtype">string</span> <a class="code" href="classbdm_1_1mexEpdf.html#a287414259db205e5bc82c0023c09e633" title="prefix of matlab functions">name</a>; |
| 78 | <a name="l00023"></a><a class="code" href="classbdm_1_1mexEpdf.html#a5190ab6b434834d4efdecaa7c264df67">00023</a> mxArray *<a class="code" href="classbdm_1_1mexEpdf.html#a5190ab6b434834d4efdecaa7c264df67" title="pointer to storage structure">data</a>; |
| 79 | <a name="l00024"></a>00024 <span class="keyword">public</span>: |
| 80 | <a name="l00025"></a>00025 <a class="code" href="classbdm_1_1mexEpdf.html" title="epdf with functions implemented in matlab">mexEpdf</a>() {}; |
| 81 | <a name="l00026"></a><a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298">00026</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298">from_setting</a>(<span class="keyword">const</span> Setting &S) { |
| 82 | <a name="l00027"></a>00027 <a class="code" href="classbdm_1_1mexEpdf.html#a287414259db205e5bc82c0023c09e633" title="prefix of matlab functions">name</a> = (<span class="keyword">const</span> <span class="keywordtype">char</span> *) S[<span class="stringliteral">"name"</span>]; |
| 83 | <a name="l00028"></a>00028 <a class="code" href="classUImxConfig.html" title="Reimplementation of libconfig&#39;s Config class for Matlab mxArray structures.">UImxConfig</a> conf(S); |
| 84 | <a name="l00029"></a>00029 <a class="code" href="classbdm_1_1mexEpdf.html#a5190ab6b434834d4efdecaa7c264df67" title="pointer to storage structure">data</a> = mxDuplicateArray(conf.<a class="code" href="classUImxConfig.html#accf378aab12b7709d644441fd2e30d00" title="Matlab structure where the info is stored.">mxconfig</a>); |
| 85 | <a name="l00030"></a>00030 <span class="comment">//mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());</span> |
| 86 | <a name="l00031"></a>00031 <span class="comment">//TODO (future...):</span> |
| 87 | <a name="l00032"></a>00032 <span class="comment">//mxArray * init_data = setting2mxarray S["init_data"];</span> |
| 88 | <a name="l00033"></a>00033 <span class="comment">//mexCallMATLAB(1, &data, 1, &init_data, name+"_from_setting");</span> |
| 89 | <a name="l00034"></a>00034 <span class="comment">//delete init_data;</span> |
| 90 | <a name="l00035"></a>00035 } |
| 91 | <a name="l00036"></a><a class="code" href="classbdm_1_1mexEpdf.html#a028dd3432d55e1fcc49f786a27f5afe8">00036</a> vec <a class="code" href="classbdm_1_1mexEpdf.html#a028dd3432d55e1fcc49f786a27f5afe8" title="return expected value">mean</a>()<span class="keyword"> const </span>{ |
| 92 | <a name="l00037"></a>00037 mxArray *tmp; |
| 93 | <a name="l00038"></a>00038 <span class="keywordtype">string</span> fname = <a class="code" href="classbdm_1_1mexEpdf.html#a287414259db205e5bc82c0023c09e633" title="prefix of matlab functions">name</a>+<span class="stringliteral">"_mean"</span>; |
| 94 | <a name="l00039"></a>00039 mexCallMATLAB(1, &tmp, 1, (mxArray **) &<a class="code" href="classbdm_1_1mexEpdf.html#a5190ab6b434834d4efdecaa7c264df67" title="pointer to storage structure">data</a>, fname.c_str()); |
| 95 | <a name="l00040"></a>00040 <span class="keywordflow">return</span> mxArray2vec(tmp); |
| 96 | <a name="l00041"></a>00041 } |
| 97 | <a name="l00042"></a>00042 }; |
| 98 | <a name="l00043"></a>00043 <a class="code" href="user__info_8h.html#a4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(mexEpdf); |
| 99 | <a name="l00044"></a>00044 |
| 100 | <a name="l00046"></a><a class="code" href="classbdm_1_1mexBM.html">00046</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mexBM.html" title="BM with functions implemented in matlab.">mexBM</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of a system, i.e. all uncertainty is modeled by probabilities.">BM</a>{ |
| 101 | <a name="l00047"></a>00047 <span class="keyword">protected</span> : |
| 102 | <a name="l00049"></a><a class="code" href="classbdm_1_1mexBM.html#a61122bf54322e51a965ce3eacb9ee3ac">00049</a> <span class="keywordtype">string</span> <a class="code" href="classbdm_1_1mexBM.html#a61122bf54322e51a965ce3eacb9ee3ac" title="prefix of matlab functions">name</a>; |
| 103 | <a name="l00051"></a><a class="code" href="classbdm_1_1mexBM.html#adfaf13987bb1263bc32426da8bc94496">00051</a> <a class="code" href="classbdm_1_1mexEpdf.html" title="epdf with functions implemented in matlab">mexEpdf</a> <a class="code" href="classbdm_1_1mexBM.html#adfaf13987bb1263bc32426da8bc94496" title="internal estimator">est</a>; |
| 104 | <a name="l00053"></a><a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f">00053</a> mxArray *<a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a>; |
| 105 | <a name="l00054"></a>00054 <span class="keyword">public</span>: |
| 106 | <a name="l00055"></a>00055 <a class="code" href="classbdm_1_1mexBM.html" title="BM with functions implemented in matlab.">mexBM</a>() {} |
| 107 | <a name="l00056"></a>00056 |
| 108 | <a name="l00058"></a><a class="code" href="classbdm_1_1mexBM.html#a4a35ffc1eedee492e54c031623fd72a8">00058</a> mxArray *<a class="code" href="classbdm_1_1mexBM.html#a4a35ffc1eedee492e54c031623fd72a8" title="duplicate internal data pointer?">get_data</a>() { |
| 109 | <a name="l00059"></a>00059 <span class="comment">//mexCallMATLAB(0, NULL, 1, &data, "dump");</span> |
| 110 | <a name="l00060"></a>00060 <span class="keywordflow">return</span> mxDuplicateArray(<a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a>); |
| 111 | <a name="l00061"></a>00061 } |
| 112 | <a name="l00062"></a>00062 |
| 113 | <a name="l00063"></a><a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd">00063</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd" title="This method arrange instance properties according the data stored in the Setting...">from_setting</a>(<span class="keyword">const</span> Setting &S) { |
| 114 | <a name="l00064"></a>00064 <a class="code" href="classbdm_1_1mexBM.html#a1dc5937132e070894ca6002b05512afd" title="This method arrange instance properties according the data stored in the Setting...">BM::from_setting</a>(S); |
| 115 | <a name="l00065"></a>00065 <a class="code" href="classbdm_1_1mexBM.html#a61122bf54322e51a965ce3eacb9ee3ac" title="prefix of matlab functions">name</a> = (<span class="keyword">const</span> <span class="keywordtype">char</span> *) S[<span class="stringliteral">"name"</span>]; |
| 116 | <a name="l00066"></a>00066 <a class="code" href="classUImxConfig.html" title="Reimplementation of libconfig&#39;s Config class for Matlab mxArray structures.">UImxConfig</a> conf(S); |
| 117 | <a name="l00067"></a>00067 <a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a> = mxDuplicateArray(conf.<a class="code" href="classUImxConfig.html#accf378aab12b7709d644441fd2e30d00" title="Matlab structure where the info is stored.">mxconfig</a>); |
| 118 | <a name="l00068"></a>00068 <span class="comment">//string fname = name+"_new";</span> |
| 119 | <a name="l00069"></a>00069 <span class="comment">//mexCallMATLAB(1, &data, 0, 0, (name+"_new").c_str());</span> |
| 120 | <a name="l00070"></a>00070 <span class="comment">//the following works as long as the posterior is the</span> |
| 121 | <a name="l00071"></a>00071 <span class="comment">//only member object there could be a structure of </span> |
| 122 | <a name="l00072"></a>00072 <span class="comment">//member objects in the setting and a for cycle over</span> |
| 123 | <a name="l00073"></a>00073 <span class="comment">//all of them right here in the code</span> |
| 124 | <a name="l00074"></a>00074 Setting &<a class="code" href="classbdm_1_1mexBM.html#a76b0e1e8f86a90a829e6821d013d7683" title="return correctly typed posterior (covariant return)">posterior</a> = S[<span class="stringliteral">"posterior"</span>]; |
| 125 | <a name="l00075"></a>00075 <a class="code" href="classbdm_1_1mexBM.html#adfaf13987bb1263bc32426da8bc94496" title="internal estimator">est</a>.<a class="code" href="classbdm_1_1mexEpdf.html#af371482ea8d21c8ac01f011b3d942298">from_setting</a>(posterior); |
| 126 | <a name="l00076"></a>00076 } |
| 127 | <a name="l00077"></a><a class="code" href="classbdm_1_1mexBM.html#aa99452b52f731f561b01029af18c68f7">00077</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mexBM.html#aa99452b52f731f561b01029af18c68f7" title="Incremental Bayes rule.">bayes</a>(<span class="keyword">const</span> vec &dt) { |
| 128 | <a name="l00078"></a>00078 <span class="comment">//void bayes() {</span> |
| 129 | <a name="l00079"></a>00079 mxArray *tmp, *old; |
| 130 | <a name="l00080"></a>00080 mxArray *in[2]; |
| 131 | <a name="l00081"></a>00081 in[0] = <a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a>; |
| 132 | <a name="l00082"></a>00082 in[1] = mxCreateDoubleMatrix(dt.size(), 1, mxREAL); |
| 133 | <a name="l00083"></a>00083 vec2mxArray(dt, in[1]); |
| 134 | <a name="l00084"></a>00084 mexCallMATLAB(1, &tmp, 2, in, (<a class="code" href="classbdm_1_1mexBM.html#a61122bf54322e51a965ce3eacb9ee3ac" title="prefix of matlab functions">name</a>+<span class="stringliteral">"_bayes"</span>).c_str()); |
| 135 | <a name="l00085"></a>00085 old = <a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a>; |
| 136 | <a name="l00086"></a>00086 <a class="code" href="classbdm_1_1mexBM.html#af3f5f0313085c2f8db8249414644a62f" title="mxArray with attributes of this object">data</a> = mxDuplicateArray(tmp); |
| 137 | <a name="l00087"></a>00087 <span class="keywordflow">if</span> (old) mxDestroyArray(old); |
| 138 | <a name="l00088"></a>00088 <span class="keywordflow">if</span> (tmp) mxDestroyArray(tmp); |
| 139 | <a name="l00089"></a>00089 <span class="comment">//mexCallMATLAB(0, NULL, 1, &data, "dump");</span> |
| 140 | <a name="l00090"></a>00090 } |
| 141 | <a name="l00092"></a><a class="code" href="classbdm_1_1mexBM.html#a76b0e1e8f86a90a829e6821d013d7683">00092</a> <span class="keyword">const</span> <a class="code" href="classbdm_1_1mexEpdf.html" title="epdf with functions implemented in matlab">mexEpdf</a>& <a class="code" href="classbdm_1_1mexBM.html#a76b0e1e8f86a90a829e6821d013d7683" title="return correctly typed posterior (covariant return)">posterior</a>()<span class="keyword"> const </span>{ |
| 142 | <a name="l00093"></a>00093 <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mexBM.html#adfaf13987bb1263bc32426da8bc94496" title="internal estimator">est</a>; |
| 143 | <a name="l00094"></a>00094 } <span class="comment">//tohle by melo zustat!!</span> |
| 144 | <a name="l00095"></a>00095 }; |
| 145 | <a name="l00096"></a>00096 <a class="code" href="user__info_8h.html#a4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(mexBM); |
| 146 | <a name="l00097"></a>00097 |
| 147 | <a name="l00098"></a>00098 } |