Context Navigation

bdmbase_8h_source.html @ 472

Revision 472, 100.0 kB (checked in by mido, 15 years ago)
obnovena dokumentace, vcetne pridani mnoha novych doc-souboru do svn

Line
1	<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2	<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
3	<title>mixpp: bdmbase.h Source File</title>
4	<link href="tabs.css" rel="stylesheet" type="text/css">
5	<link href="doxygen.css" rel="stylesheet" type="text/css">
6	</head><body>
7	<!-- Generated by Doxygen 1.5.9 -->
8	<script type="text/javascript">
9	<!--
10	function changeDisplayState (e){
11	var num=this.id.replace(/[^[0-9]/g,'');
12	var button=this.firstChild;
13	var sectionDiv=document.getElementById('dynsection'+num);
14	if (sectionDiv.style.display=='none'\|\|sectionDiv.style.display==''){
15	sectionDiv.style.display='block';
16	button.src='open.gif';
17	}else{
18	sectionDiv.style.display='none';
19	button.src='closed.gif';
20	}
21	}
22	function initDynSections(){
23	var divs=document.getElementsByTagName('div');
24	var sectionCounter=1;
25	for(var i=0;i<divs.length-1;i++){
26	if(divs[i].className=='dynheader'&&divs[i+1].className=='dynsection'){
27	var header=divs[i];
28	var section=divs[i+1];
29	var button=header.firstChild;
30	if (button!='IMG'){
31	divs[i].insertBefore(document.createTextNode(' '),divs[i].firstChild);
32	button=document.createElement('img');
33	divs[i].insertBefore(button,divs[i].firstChild);
34	}
35	header.style.cursor='pointer';
36	header.onclick=changeDisplayState;
37	header.id='dynheader'+sectionCounter;
38	button.src='closed.gif';
39	section.id='dynsection'+sectionCounter;
40	section.style.display='none';
41	section.style.marginLeft='14px';
42	sectionCounter++;
43	}
44	}
45	}
46	window.onload = initDynSections;
47	-->
48	</script>
49	<div class="navigation" id="top">
50	<div class="tabs">
51	<ul>
52	<li><a href="main.html"><span>Main Page</span></a></li>
53	<li><a href="pages.html"><span>Related Pages</span></a></li>
54	<li><a href="annotated.html"><span>Classes</span></a></li>
55	<li class="current"><a href="files.html"><span>Files</span></a></li>
56	</ul>
57	</div>
58	<div class="tabs">
59	<ul>
60	<li><a href="files.html"><span>File List</span></a></li>
61	<li><a href="globals.html"><span>File Members</span></a></li>
62	</ul>
63	</div>
64	<h1>bdmbase.h</h1><a href="bdmbase_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
65	<a name="l00013"></a>00013 <span class="preprocessor">#ifndef BDMBASE_H</span>
66	<a name="l00014"></a>00014 <span class="preprocessor"></span><span class="preprocessor">#define BDMBASE_H</span>
67	<a name="l00015"></a>00015 <span class="preprocessor"></span>
68	<a name="l00016"></a>00016 <span class="preprocessor">#include <map></span>
69	<a name="l00017"></a>00017
70	<a name="l00018"></a>00018 <span class="preprocessor">#include "../itpp_ext.h"</span>
71	<a name="l00019"></a>00019 <span class="preprocessor">#include "../bdmroot.h"</span>
72	<a name="l00020"></a>00020 <span class="preprocessor">#include "../shared_ptr.h"</span>
73	<a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="user__info_8h.html" title="UI (user info) class for loading/saving objects from/to configuration files. It is...">user_info.h</a>"</span>
74	<a name="l00022"></a>00022
75	<a name="l00023"></a>00023 <span class="keyword">using namespace </span>libconfig;
76	<a name="l00024"></a>00024 <span class="keyword">using namespace </span>itpp;
77	<a name="l00025"></a>00025 <span class="keyword">using namespace </span>std;
78	<a name="l00026"></a>00026
79	<a name="l00027"></a>00027 <span class="keyword">namespace </span>bdm
80	<a name="l00028"></a>00028 {
81	<a name="l00029"></a>00029
82	<a name="l00030"></a>00030 <span class="keyword">typedef</span> std::map<string, int> RVmap;
83	<a name="l00031"></a>00031 <span class="keyword">extern</span> ivec RV_SIZES;
84	<a name="l00032"></a>00032 <span class="keyword">extern</span> Array<string> RV_NAMES;
85	<a name="l00033"></a>00033
86	<a name="l00035"></a><a class="code" href="classbdm_1_1str.html">00035</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1str.html" title="Structure of RV, i.e. RVs expanded into a flat list of IDs, used for debugging.">str</a>
87	<a name="l00036"></a>00036 {
88	<a name="l00037"></a>00037 <span class="keyword">public</span>:
89	<a name="l00039"></a><a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0">00039</a> ivec <a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0" title="vector id ids (non-unique!)">ids</a>;
90	<a name="l00041"></a><a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a">00041</a> ivec <a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a" title="vector of times">times</a>;
91	<a name="l00043"></a><a class="code" href="classbdm_1_1str.html#7c6b9dccc6814d7bd1e8e2208a0d3f3f">00043</a> <a class="code" href="classbdm_1_1str.html#7c6b9dccc6814d7bd1e8e2208a0d3f3f" title="Default constructor.">str</a>(ivec ids0, ivec times0) : <a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0" title="vector id ids (non-unique!)">ids</a>(ids0), <a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a" title="vector of times">times</a>(times0) {
92	<a name="l00044"></a>00044 it_assert_debug(times0.length() == ids0.length(), <span class="stringliteral">"Incompatible input"</span>);
93	<a name="l00045"></a>00045 };
94	<a name="l00046"></a>00046 };
95	<a name="l00047"></a>00047
96	<a name="l00086"></a><a class="code" href="classbdm_1_1RV.html">00086</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
97	<a name="l00087"></a>00087 {
98	<a name="l00088"></a>00088 <span class="keyword">protected</span>:
99	<a name="l00090"></a><a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb">00090</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb" title="size of the data vector">dsize</a>;
100	<a name="l00092"></a><a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385">00092</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="number of individual rvs">len</a>;
101	<a name="l00094"></a><a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b">00094</a> ivec <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>;
102	<a name="l00096"></a><a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896">00096</a> ivec <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>;
103	<a name="l00097"></a>00097
104	<a name="l00098"></a>00098 <span class="keyword">private</span>:
105	<a name="l00100"></a>00100 <span class="keywordtype">void</span> init(<span class="keyword">const</span> Array<std::string> &in_names, <span class="keyword">const</span> ivec &in_sizes, <span class="keyword">const</span> ivec &in_times);
106	<a name="l00101"></a>00101 <span class="keywordtype">int</span> init(<span class="keyword">const</span> <span class="keywordtype">string</span> &name, <span class="keywordtype">int</span> size);
107	<a name="l00102"></a>00102 <span class="keyword">public</span>:
108	<a name="l00105"></a>00105
109	<a name="l00107"></a><a class="code" href="classbdm_1_1RV.html#3fa35c2744e687c74f7742990e6aa264">00107</a> <a class="code" href="classbdm_1_1RV.html#3fa35c2744e687c74f7742990e6aa264" title="Full constructor.">RV</a>(<span class="keyword">const</span> Array<std::string> &in_names, <span class="keyword">const</span> ivec &in_sizes, <span class="keyword">const</span> ivec &in_times) { init(in_names, in_sizes, in_times); }
110	<a name="l00108"></a>00108
111	<a name="l00110"></a><a class="code" href="classbdm_1_1RV.html#381f5d963bb54efd4e76815dc93d50ce">00110</a> <a class="code" href="classbdm_1_1RV.html#381f5d963bb54efd4e76815dc93d50ce" title="Constructor with times=0.">RV</a>(<span class="keyword">const</span> Array<std::string> &in_names, <span class="keyword">const</span> ivec &in_sizes) { init(in_names, in_sizes, zeros_i(in_names.length())); }
112	<a name="l00111"></a>00111
113	<a name="l00113"></a><a class="code" href="classbdm_1_1RV.html#cb47d71d82639ec7ddf2635b21aafd00">00113</a> <a class="code" href="classbdm_1_1RV.html#cb47d71d82639ec7ddf2635b21aafd00" title="Constructor with sizes=1, times=0.">RV</a>(<span class="keyword">const</span> Array<std::string> &in_names) { init(in_names, ones_i(in_names.length()), zeros_i(in_names.length())); }
114	<a name="l00114"></a>00114
115	<a name="l00116"></a><a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412">00116</a> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a>() : <a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb" title="size of the data vector">dsize</a>(0), <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="number of individual rvs">len</a>(0), <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>(0), <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>(0) {}
116	<a name="l00118"></a>00118 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a>(<span class="keywordtype">string</span> name, <span class="keywordtype">int</span> sz, <span class="keywordtype">int</span> tm = 0);
117	<a name="l00120"></a>00120
118	<a name="l00123"></a>00123
119	<a name="l00125"></a>00125 <span class="keyword">friend</span> std::ostream &<a class="code" href="classbdm_1_1RV.html#cad6f15535d10437f6690cb4492176a8" title="State output, e.g. for debugging.">operator<< </a>(std::ostream &os, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv);
120	<a name="l00126"></a>00126
121	<a name="l00127"></a>00127 <span class="keywordtype">int</span> _dsize()<span class="keyword"> const </span>{ <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb" title="size of the data vector">dsize</a>; }
122	<a name="l00128"></a>00128
123	<a name="l00130"></a>00130 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#39f9b44752376f07e346579b5ff8cbbd" title="Recount size of the corresponding data vector.">countsize</a>() <span class="keyword">const</span>;
124	<a name="l00131"></a>00131 ivec cumsizes() <span class="keyword">const</span>;
125	<a name="l00132"></a>00132 <span class="keywordtype">int</span> length()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="number of individual rvs">len</a>;}
126	<a name="l00133"></a>00133 <span class="keywordtype">int</span> id(<span class="keywordtype">int</span> at)<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>(at);}
127	<a name="l00134"></a>00134 <span class="keywordtype">int</span> size(<span class="keywordtype">int</span> at)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> RV_SIZES(<a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>(at)); }
128	<a name="l00135"></a>00135 <span class="keywordtype">int</span> time(<span class="keywordtype">int</span> at)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>(at); }
129	<a name="l00136"></a>00136 std::string name(<span class="keywordtype">int</span> at)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> RV_NAMES(<a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>(at)); }
130	<a name="l00137"></a>00137 <span class="keywordtype">void</span> set_time(<span class="keywordtype">int</span> at, <span class="keywordtype">int</span> time0) { <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>(at) = time0; }
131	<a name="l00139"></a>00139
132	<a name="l00140"></a>00140 <span class="comment">//TODO why not inline and later??</span>
133	<a name="l00141"></a>00141
134	<a name="l00144"></a>00144
135	<a name="l00146"></a>00146 ivec <a class="code" href="classbdm_1_1RV.html#1367055b2608ad9d50d01a2da6cec869" title="Find indices of self in another rv,.">findself</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> &rv2) <span class="keyword">const</span>;
136	<a name="l00148"></a>00148 <span class="keywordtype">bool</span> <a class="code" href="classbdm_1_1RV.html#3bb011af88ae98f5c97955bb2c95cfee" title="Compare if rv2 is identical to this RV.">equal</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> &rv2) <span class="keyword">const</span>;
137	<a name="l00150"></a>00150 <span class="keywordtype">bool</span> <a class="code" href="classbdm_1_1RV.html#87841b5ee43997b79789c0c22047e224" title="Add (concat) another variable to the current one,.">add</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> &rv2);
138	<a name="l00152"></a>00152 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> <a class="code" href="classbdm_1_1RV.html#aec44dabdf0a6d90fbae95e1356eda39" title="Subtract another variable from the current one.">subt</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> &rv2) <span class="keyword">const</span>;
139	<a name="l00154"></a>00154 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> <a class="code" href="classbdm_1_1RV.html#c0c04ab99ded62e1e05bcefb151ebeb9" title="Select only variables at indices ind.">subselect</a>(<span class="keyword">const</span> ivec &ind) <span class="keyword">const</span>;
140	<a name="l00155"></a>00155
141	<a name="l00157"></a><a class="code" href="classbdm_1_1RV.html#6f3daf0f857fff53a1c0a445fe160966">00157</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1RV.html#6f3daf0f857fff53a1c0a445fe160966" title="Select only variables at indices ind.">operator()</a>(<span class="keyword">const</span> ivec &ind)<span class="keyword"> const </span>{ <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#c0c04ab99ded62e1e05bcefb151ebeb9" title="Select only variables at indices ind.">subselect</a>(ind); }
142	<a name="l00158"></a>00158
143	<a name="l00160"></a>00160 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1RV.html#6f3daf0f857fff53a1c0a445fe160966" title="Select only variables at indices ind.">operator()</a>(<span class="keywordtype">int</span> di1, <span class="keywordtype">int</span> di2) <span class="keyword">const</span>;
144	<a name="l00161"></a>00161
145	<a name="l00163"></a>00163 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#bd77cb0899f7fba1d0e34396137b7f1f" title="Shift time by delta.">t</a>(<span class="keywordtype">int</span> delta);
146	<a name="l00165"></a>00165
147	<a name="l00168"></a>00168
148	<a name="l00170"></a>00170 <a class="code" href="classbdm_1_1str.html" title="Structure of RV, i.e. RVs expanded into a flat list of IDs, used for debugging.">str</a> <a class="code" href="classbdm_1_1RV.html#8bdd7f4072749c7cc864b78a418ca388" title="generate str from rv, by expanding sizes TODO to_string..">tostr</a>() <span class="keyword">const</span>;
149	<a name="l00173"></a>00173 ivec <a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &crv) <span class="keyword">const</span>;
150	<a name="l00176"></a>00176 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv2, ivec &selfi, ivec &rv2i) <span class="keyword">const</span>;
151	<a name="l00178"></a><a class="code" href="classbdm_1_1RV.html#9a5c5aeeaf39299d07763377b05a8241">00178</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#9a5c5aeeaf39299d07763377b05a8241" title="Minimum time-offset.">mint</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> min(<a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>);}
152	<a name="l00180"></a>00180
153	<a name="l00181"></a>00181 <span class="comment">// TODO aktualizovat dle soucasneho UI</span>
154	<a name="l00196"></a>00196 <span class="comment"></span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#6564bdd36b6845a51446a192d33e1194" title="UI for class RV (description of data vectors).">from_setting</a>(<span class="keyword">const</span> Setting &<span class="keyword">set</span>);
155	<a name="l00197"></a>00197
156	<a name="l00198"></a>00198 <span class="comment">// TODO dodelat void to_setting( Setting &set ) const;</span>
157	<a name="l00199"></a>00199
158	<a name="l00201"></a>00201 <span class="keyword">static</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#ab92b361335a644df654cf7cf7e5ffc0" title="Invalidate all named RVs. Use before initializing any RV instances, with care...">clear_all</a>();
159	<a name="l00202"></a>00202 };
160	<a name="l00203"></a>00203 <a class="code" href="user__info_8h.html#4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(RV);
161	<a name="l00204"></a>00204
162	<a name="l00206"></a>00206 RV concat(<span class="keyword">const</span> RV &rv1, <span class="keyword">const</span> RV &rv2);
163	<a name="l00207"></a>00207
164	<a name="l00209"></a>00209 <span class="keyword">extern</span> RV RV0;
165	<a name="l00210"></a>00210
166	<a name="l00212"></a>00212
167	<a name="l00213"></a><a class="code" href="classbdm_1_1fnc.html">00213</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1fnc.html" title="Class representing function of variable represented by rv.">fnc</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
168	<a name="l00214"></a>00214 {
169	<a name="l00215"></a>00215 <span class="keyword">protected</span>:
170	<a name="l00217"></a><a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62">00217</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62" title="Length of the output vector.">dimy</a>;
171	<a name="l00218"></a>00218 <span class="keyword">public</span>:
172	<a name="l00220"></a><a class="code" href="classbdm_1_1fnc.html#65e9fb11f6cc7884018d8f229bcabbb7">00220</a> <a class="code" href="classbdm_1_1fnc.html#65e9fb11f6cc7884018d8f229bcabbb7" title="default constructor">fnc</a>() {};
173	<a name="l00222"></a><a class="code" href="classbdm_1_1fnc.html#6277b11d7fffc7ef8a2fa3e84ae5bad4">00222</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1fnc.html#6277b11d7fffc7ef8a2fa3e84ae5bad4" title="function evaluates numerical value of at cond ">eval</a>(<span class="keyword">const</span> vec &cond) {
174	<a name="l00223"></a>00223 <span class="keywordflow">return</span> vec(0);
175	<a name="l00224"></a>00224 };
176	<a name="l00225"></a>00225
177	<a name="l00227"></a><a class="code" href="classbdm_1_1fnc.html#0786e40fade2663a70d654c1dda5d73e">00227</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1fnc.html#0786e40fade2663a70d654c1dda5d73e" title="function substitutes given value into an appropriate position">condition</a>(<span class="keyword">const</span> vec &val) {};
178	<a name="l00228"></a>00228
179	<a name="l00230"></a><a class="code" href="classbdm_1_1fnc.html#083832294da9d1e40804158b979c4341">00230</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1fnc.html#083832294da9d1e40804158b979c4341" title="access function">dimension</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62" title="Length of the output vector.">dimy</a>;}
180	<a name="l00231"></a>00231 };
181	<a name="l00232"></a>00232
182	<a name="l00233"></a>00233 <span class="keyword">class </span>mpdf;
183	<a name="l00234"></a>00234
184	<a name="l00236"></a>00236
185	<a name="l00237"></a><a class="code" href="classbdm_1_1epdf.html">00237</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
186	<a name="l00238"></a>00238 {
187	<a name="l00239"></a>00239 <span class="keyword">protected</span>:
188	<a name="l00241"></a><a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce">00241</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>;
189	<a name="l00243"></a><a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8">00243</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>;
190	<a name="l00244"></a>00244
191	<a name="l00245"></a>00245 <span class="keyword">public</span>:
192	<a name="l00257"></a>00257 <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>() : <a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>(0), <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>() {};
193	<a name="l00258"></a>00258 <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a> &e) : <a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>(e.<a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>), <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>(e.<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>) {};
194	<a name="l00259"></a>00259 epdf(<span class="keyword">const</span> RV &rv0):<a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>(rv0._dsize()) {<a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff" title="Name its rv.">set_rv</a>(rv0);};
195	<a name="l00260"></a>00260 <span class="keywordtype">void</span> set_parameters(<span class="keywordtype">int</span> dim0) {<a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a> = dim0;}
196	<a name="l00262"></a>00262
197	<a name="l00265"></a>00265
198	<a name="l00267"></a><a class="code" href="classbdm_1_1epdf.html#75ba57a89681a8efec488f5796c7cc1b">00267</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#75ba57a89681a8efec488f5796c7cc1b" title="Returns a sample, from density .">sample</a>()<span class="keyword"> const </span>{it_error(<span class="stringliteral">"not implemneted"</span>); <span class="keywordflow">return</span> vec(0);};
199	<a name="l00269"></a>00269 <span class="keyword">virtual</span> mat <a class="code" href="classbdm_1_1epdf.html#b4cf45fd83cc7573ede9fe1215256058" title="Returns N samples, from density .">sample_m</a>(<span class="keywordtype">int</span> N) <span class="keyword">const</span>;
200	<a name="l00272"></a><a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249">00272</a> <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249">evallog</a>(<span class="keyword">const</span> vec &val)<span class="keyword"> const </span>{it_error(<span class="stringliteral">"not implemneted"</span>); <span class="keywordflow">return</span> 0.0;};
201	<a name="l00274"></a><a class="code" href="classbdm_1_1epdf.html#34956d4dd3176eeb5937cf48a1546b62">00274</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#34956d4dd3176eeb5937cf48a1546b62" title="Compute log-probability of multiple values argument val.">evallog_m</a>(<span class="keyword">const</span> mat &Val)<span class="keyword"> const </span>{
202	<a name="l00275"></a>00275 vec x(Val.cols());
203	<a name="l00276"></a>00276 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < Val.cols(); i++) {x(i) = <a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249">evallog</a>(Val.get_col(i)) ;}
204	<a name="l00277"></a>00277 <span class="keywordflow">return</span> x;
205	<a name="l00278"></a>00278 }
206	<a name="l00280"></a><a class="code" href="classbdm_1_1epdf.html#512dbf057299e871fbc61ed55acb5888">00280</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#512dbf057299e871fbc61ed55acb5888" title="Compute log-probability of multiple values argument val.">evallog_m</a>(<span class="keyword">const</span> Array<vec> &Avec)<span class="keyword"> const </span>{
207	<a name="l00281"></a>00281 vec x(Avec.size());
208	<a name="l00282"></a>00282 <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < Avec.size(); i++) {x(i) = <a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249">evallog</a>(Avec(i)) ;}
209	<a name="l00283"></a>00283 <span class="keywordflow">return</span> x;
210	<a name="l00284"></a>00284 }
211	<a name="l00286"></a><a class="code" href="classbdm_1_1epdf.html#e584eac5579c1b6384947ecf66166c77">00286</a> <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>* <a class="code" href="classbdm_1_1epdf.html#e584eac5579c1b6384947ecf66166c77" title="Return conditional density on the given RV, the remaining rvs will be in conditioning...">condition</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>)<span class="keyword"> const </span>{it_warning(<span class="stringliteral">"Not implemented"</span>); <span class="keywordflow">return</span> NULL;}
212	<a name="l00287"></a>00287
213	<a name="l00289"></a><a class="code" href="classbdm_1_1epdf.html#3fb2ece54f720b62ad325e61214fa0a1">00289</a> <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classbdm_1_1epdf.html#3fb2ece54f720b62ad325e61214fa0a1" title="Return marginal density on the given RV, the remainig rvs are intergrated out.">marginal</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>)<span class="keyword"> const </span>{it_warning(<span class="stringliteral">"Not implemented"</span>); <span class="keywordflow">return</span> NULL;}
214	<a name="l00290"></a>00290
215	<a name="l00292"></a><a class="code" href="classbdm_1_1epdf.html#0c8d13d788094605b3cdab20a6e22b75">00292</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#0c8d13d788094605b3cdab20a6e22b75" title="return expected value">mean</a>()<span class="keyword"> const </span>{it_error(<span class="stringliteral">"not implemneted"</span>); <span class="keywordflow">return</span> vec(0);};
216	<a name="l00293"></a>00293
217	<a name="l00295"></a><a class="code" href="classbdm_1_1epdf.html#4543d40ccfc5d549b0c831a57a3ed84c">00295</a> <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#4543d40ccfc5d549b0c831a57a3ed84c" title="return expected variance (not covariance!)">variance</a>()<span class="keyword"> const </span>{it_error(<span class="stringliteral">"not implemneted"</span>); <span class="keywordflow">return</span> vec(0);};
218	<a name="l00297"></a><a class="code" href="classbdm_1_1epdf.html#0a7705e80601957f7226d333a0a6d539">00297</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1epdf.html#0a7705e80601957f7226d333a0a6d539" title="Lower and upper bounds of percentage % quantile, returns mean-2*sigma as default...">qbounds</a>(vec &lb, vec &ub, <span class="keywordtype">double</span> percentage = 0.95)<span class="keyword"> const </span>{
219	<a name="l00298"></a>00298 vec mea = <a class="code" href="classbdm_1_1epdf.html#0c8d13d788094605b3cdab20a6e22b75" title="return expected value">mean</a>();
220	<a name="l00299"></a>00299 vec std = sqrt(<a class="code" href="classbdm_1_1epdf.html#4543d40ccfc5d549b0c831a57a3ed84c" title="return expected variance (not covariance!)">variance</a>());
221	<a name="l00300"></a>00300 lb = mea - 2 * std;
222	<a name="l00301"></a>00301 ub = mea + 2 * std;
223	<a name="l00302"></a>00302 };
224	<a name="l00304"></a>00304
225	<a name="l00310"></a>00310
226	<a name="l00312"></a><a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff">00312</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff" title="Name its rv.">set_rv</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv0) {<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a> = rv0; } <span class="comment">//it_assert_debug(isnamed(),""); };</span>
227	<a name="l00314"></a><a class="code" href="classbdm_1_1epdf.html#c4b863ff84c7a4882fb3ad18556027f9">00314</a> <span class="comment"></span> <span class="keywordtype">bool</span> <a class="code" href="classbdm_1_1epdf.html#c4b863ff84c7a4882fb3ad18556027f9" title="True if rv is assigned.">isnamed</a>()<span class="keyword"> const </span>{<span class="keywordtype">bool</span> b = (<a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a> == <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>._dsize()); <span class="keywordflow">return</span> b;}
228	<a name="l00316"></a><a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21">00316</a> <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21" title="Return name (fails when isnamed is false).">_rv</a>()<span class="keyword"> const </span>{it_assert_debug(<a class="code" href="classbdm_1_1epdf.html#c4b863ff84c7a4882fb3ad18556027f9" title="True if rv is assigned.">isnamed</a>(), <span class="stringliteral">""</span>); <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Description of the random variable.">rv</a>;}
229	<a name="l00318"></a>00318
230	<a name="l00321"></a>00321
231	<a name="l00323"></a><a class="code" href="classbdm_1_1epdf.html#7083a65f7b7a0d0d13b2c516bd2ec29c">00323</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1epdf.html#7083a65f7b7a0d0d13b2c516bd2ec29c" title="Size of the random variable.">dimension</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>;}
232	<a name="l00331"></a><a class="code" href="classbdm_1_1epdf.html#fa9efc98a065b31c630fe18cc418f725">00331</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1epdf.html#fa9efc98a065b31c630fe18cc418f725">from_setting</a>(<span class="keyword">const</span> Setting &<span class="keyword">set</span>){
233	<a name="l00332"></a>00332 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>* r = UI::build<RV>(<span class="keyword">set</span>,<span class="stringliteral">"rv"</span>);
234	<a name="l00333"></a>00333 <span class="keywordflow">if</span> (r){
235	<a name="l00334"></a>00334 <a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff" title="Name its rv.">set_rv</a>(*r);
236	<a name="l00335"></a>00335 <span class="keyword">delete</span> r;
237	<a name="l00336"></a>00336 }
238	<a name="l00337"></a>00337 }
239	<a name="l00338"></a>00338
240	<a name="l00339"></a>00339 };
241	<a name="l00340"></a>00340
242	<a name="l00341"></a>00341
243	<a name="l00343"></a>00343 <span class="comment">//TODO Samplecond can be generalized</span>
244	<a name="l00344"></a>00344
245	<a name="l00345"></a><a class="code" href="classbdm_1_1mpdf.html">00345</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
246	<a name="l00346"></a>00346 {
247	<a name="l00347"></a>00347 <span class="keyword">protected</span>:
248	<a name="l00349"></a><a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6">00349</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>;
249	<a name="l00351"></a><a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288">00351</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>;
250	<a name="l00352"></a>00352
251	<a name="l00353"></a>00353 <span class="keyword">private</span>:
252	<a name="l00355"></a>00355 <a class="code" href="classbdm_1_1shared__ptr.html">shared_ptr<epdf></a> shep;
253	<a name="l00356"></a>00356
254	<a name="l00357"></a>00357 <span class="keyword">public</span>:
255	<a name="l00360"></a>00360
256	<a name="l00361"></a>00361 <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>():<a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>(0), <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>() { }
257	<a name="l00362"></a>00362
258	<a name="l00363"></a>00363 <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> &m):<a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>(m.<a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>), <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>(m.<a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>), shep(m.shep) { }
259	<a name="l00365"></a>00365
260	<a name="l00368"></a>00368
261	<a name="l00370"></a>00370 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#b024f2e636e5507c755de4a6df1f1e99" title="Returns a sample from the density conditioned on cond, .">samplecond</a>(<span class="keyword">const</span> vec &cond);
262	<a name="l00371"></a>00371
263	<a name="l00373"></a>00373 <span class="keyword">virtual</span> mat <a class="code" href="classbdm_1_1mpdf.html#6e497d501331057309d2fddf63259d1a" title="Returns.">samplecond_m</a>(<span class="keyword">const</span> vec &cond, <span class="keywordtype">int</span> N);
264	<a name="l00374"></a>00374
265	<a name="l00376"></a><a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530">00376</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a>(<span class="keyword">const</span> vec &cond) {it_error(<span class="stringliteral">"Not implemented"</span>);};
266	<a name="l00377"></a>00377
267	<a name="l00379"></a>00379 <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1mpdf.html#cfcbb897bf42015a7979b17d891fd9b5" title="Shortcut for conditioning and evaluation of the internal epdf. In some cases, this...">evallogcond</a>(<span class="keyword">const</span> vec &dt, <span class="keyword">const</span> vec &cond);
268	<a name="l00380"></a>00380
269	<a name="l00382"></a>00382 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#f97c61ecfb1f44ef6f08dee7bc15eab3" title="Matrix version of evallogcond.">evallogcond_m</a>(<span class="keyword">const</span> mat &Dt, <span class="keyword">const</span> vec &cond);
270	<a name="l00383"></a>00383
271	<a name="l00385"></a>00385 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#f97c61ecfb1f44ef6f08dee7bc15eab3" title="Matrix version of evallogcond.">evallogcond_m</a>(<span class="keyword">const</span> Array<vec> &Dt, <span class="keyword">const</span> vec &cond);
272	<a name="l00386"></a>00386
273	<a name="l00389"></a>00389
274	<a name="l00390"></a>00390 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> _rv() { <span class="keywordflow">return</span> shep->_rv(); }
275	<a name="l00391"></a>00391 RV _rvc() { it_assert_debug(isnamed(), <span class="stringliteral">""</span>); <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>; }
276	<a name="l00392"></a>00392 <span class="keywordtype">int</span> dimension() { <span class="keywordflow">return</span> shep->dimension(); }
277	<a name="l00393"></a>00393 <span class="keywordtype">int</span> dimensionc() { <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>; }
278	<a name="l00394"></a>00394
279	<a name="l00395"></a>00395 epdf *e() { <span class="keywordflow">return</span> shep.get(); }
280	<a name="l00396"></a>00396
281	<a name="l00397"></a>00397 <span class="keywordtype">void</span> set_ep(shared_ptr<epdf> ep) { shep = ep; }
282	<a name="l00398"></a>00398
283	<a name="l00407"></a>00407 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mpdf.html#30a0bceb948e287d497c578e371b1667">from_setting</a>(<span class="keyword">const</span> Setting &<span class="keyword">set</span>);
284	<a name="l00409"></a>00409
285	<a name="l00412"></a>00412 <span class="keywordtype">void</span> set_rvc(<span class="keyword">const</span> RV &rvc0) { <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a> = rvc0; }
286	<a name="l00413"></a>00413 <span class="keywordtype">void</span> set_rv(<span class="keyword">const</span> RV &rv0) { shep->set_rv(rv0); }
287	<a name="l00414"></a>00414 <span class="keywordtype">bool</span> isnamed() { <span class="keywordflow">return</span> (shep->isnamed()) && (<a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a> == <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>._dsize()); }
288	<a name="l00416"></a>00416 };
289	<a name="l00417"></a>00417
290	<a name="l00443"></a><a class="code" href="classbdm_1_1datalink.html">00443</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1datalink.html" title="DataLink is a connection between two data vectors Up and Down.">datalink</a>
291	<a name="l00444"></a>00444 {
292	<a name="l00445"></a>00445 <span class="keyword">protected</span>:
293	<a name="l00447"></a><a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c">00447</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a>;
294	<a name="l00448"></a>00448
295	<a name="l00450"></a><a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045">00450</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a>;
296	<a name="l00451"></a>00451
297	<a name="l00453"></a><a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb">00453</a> ivec <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>;
298	<a name="l00454"></a>00454
299	<a name="l00455"></a>00455 <span class="keyword">public</span>:
300	<a name="l00457"></a><a class="code" href="classbdm_1_1datalink.html#b8b5d5074fc1310a37c3713d58b900e2">00457</a> <a class="code" href="classbdm_1_1datalink.html#b8b5d5074fc1310a37c3713d58b900e2" title="Constructor.">datalink</a>():<a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a>(0), <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a>(0) { }
301	<a name="l00458"></a>00458 <a class="code" href="classbdm_1_1datalink.html#b8b5d5074fc1310a37c3713d58b900e2" title="Constructor.">datalink</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up) { <a class="code" href="classbdm_1_1datalink.html#fb8337aac8329eb0ca96bfb3ff73ec67" title="set connection, rv must be fully present in rv_up">set_connection</a>(rv, rv_up); }
302	<a name="l00459"></a>00459
303	<a name="l00461"></a><a class="code" href="classbdm_1_1datalink.html#fb8337aac8329eb0ca96bfb3ff73ec67">00461</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink.html#fb8337aac8329eb0ca96bfb3ff73ec67" title="set connection, rv must be fully present in rv_up">set_connection</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up) {
304	<a name="l00462"></a>00462 <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> = rv._dsize();
305	<a name="l00463"></a>00463 <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a> = rv_up._dsize();
306	<a name="l00464"></a>00464 <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a> = rv.<a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a>(rv_up);
307	<a name="l00465"></a>00465
308	<a name="l00466"></a>00466 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>.length() == <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a>, <span class="stringliteral">"rv is not fully in rv_up"</span>);
309	<a name="l00467"></a>00467 }
310	<a name="l00468"></a>00468
311	<a name="l00470"></a><a class="code" href="classbdm_1_1datalink.html#724486a37002e7a3a73366c3af9b82ad">00470</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink.html#724486a37002e7a3a73366c3af9b82ad" title="set connection using indices">set_connection</a>(<span class="keywordtype">int</span> ds, <span class="keywordtype">int</span> us, <span class="keyword">const</span> ivec &upind) {
312	<a name="l00471"></a>00471 <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> = ds;
313	<a name="l00472"></a>00472 <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a> = us;
314	<a name="l00473"></a>00473 <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a> = upind;
315	<a name="l00474"></a>00474
316	<a name="l00475"></a>00475 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>.length() == <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a>, <span class="stringliteral">"rv is not fully in rv_up"</span>);
317	<a name="l00476"></a>00476 }
318	<a name="l00477"></a>00477
319	<a name="l00479"></a><a class="code" href="classbdm_1_1datalink.html#5b5368c3899eed2f670d09142c2548e5">00479</a> vec <a class="code" href="classbdm_1_1datalink.html#5b5368c3899eed2f670d09142c2548e5" title="Get val for myself from val of &quot;Up&quot;.">pushdown</a>(<span class="keyword">const</span> vec &val_up) {
320	<a name="l00480"></a>00480 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a> == val_up.length(), <span class="stringliteral">"Wrong val_up"</span>);
321	<a name="l00481"></a>00481 <span class="keywordflow">return</span> get_vec(val_up, <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>);
322	<a name="l00482"></a>00482 }
323	<a name="l00483"></a>00483
324	<a name="l00485"></a><a class="code" href="classbdm_1_1datalink.html#610e966f2edb1fcd6a5b28267c571793">00485</a> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink.html#610e966f2edb1fcd6a5b28267c571793" title="Fill val of &quot;Up&quot; by my pieces.">pushup</a>(vec &val_up, <span class="keyword">const</span> vec &val) {
325	<a name="l00486"></a>00486 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> == val.length(), <span class="stringliteral">"Wrong val"</span>);
326	<a name="l00487"></a>00487 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a> == val_up.length(), <span class="stringliteral">"Wrong val_up"</span>);
327	<a name="l00488"></a>00488 set_subvector(val_up, <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>, val);
328	<a name="l00489"></a>00489 }
329	<a name="l00490"></a>00490 };
330	<a name="l00491"></a>00491
331	<a name="l00493"></a><a class="code" href="classbdm_1_1datalink__m2e.html">00493</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1datalink__m2e.html" title="Data link with a condition.">datalink_m2e</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1datalink.html" title="DataLink is a connection between two data vectors Up and Down.">datalink</a>
332	<a name="l00494"></a>00494 {
333	<a name="l00495"></a>00495 <span class="keyword">protected</span>:
334	<a name="l00497"></a><a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9">00497</a> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>;
335	<a name="l00498"></a>00498
336	<a name="l00500"></a><a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783">00500</a> ivec <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indices of the upper val">v2c_up</a>;
337	<a name="l00501"></a>00501
338	<a name="l00503"></a><a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc">00503</a> ivec <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, indices of the local cond">v2c_lo</a>;
339	<a name="l00504"></a>00504
340	<a name="l00505"></a>00505 <span class="keyword">public</span>:
341	<a name="l00507"></a><a class="code" href="classbdm_1_1datalink__m2e.html#5ad951e9cf850c9bfe4e492e91a50af6">00507</a> <a class="code" href="classbdm_1_1datalink__m2e.html#5ad951e9cf850c9bfe4e492e91a50af6" title="Constructor.">datalink_m2e</a>():<a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>(0) { }
342	<a name="l00508"></a>00508
343	<a name="l00509"></a>00509 <span class="keywordtype">void</span> set_connection(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rvc, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up) {
344	<a name="l00510"></a>00510 <a class="code" href="classbdm_1_1datalink.html#fb8337aac8329eb0ca96bfb3ff73ec67" title="set connection, rv must be fully present in rv_up">datalink::set_connection</a>(rv, rv_up);
345	<a name="l00511"></a>00511 <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> = rvc._dsize();
346	<a name="l00512"></a>00512 <span class="comment">//establish v2c connection</span>
347	<a name="l00513"></a>00513 rvc.<a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a>(rv_up, <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, indices of the local cond">v2c_lo</a>, <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indices of the upper val">v2c_up</a>);
348	<a name="l00514"></a>00514 }
349	<a name="l00515"></a>00515
350	<a name="l00517"></a><a class="code" href="classbdm_1_1datalink__m2e.html#7ba850822e3cf9a4b079afc7e7c01703">00517</a> vec <a class="code" href="classbdm_1_1datalink__m2e.html#7ba850822e3cf9a4b079afc7e7c01703" title="Construct condition.">get_cond</a>(<span class="keyword">const</span> vec &val_up) {
351	<a name="l00518"></a>00518 vec tmp(<a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>);
352	<a name="l00519"></a>00519 set_subvector(tmp, <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, indices of the local cond">v2c_lo</a>, val_up(<a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indices of the upper val">v2c_up</a>));
353	<a name="l00520"></a>00520 <span class="keywordflow">return</span> tmp;
354	<a name="l00521"></a>00521 }
355	<a name="l00522"></a>00522
356	<a name="l00523"></a>00523 <span class="keywordtype">void</span> pushup_cond(vec &val_up, <span class="keyword">const</span> vec &val, <span class="keyword">const</span> vec &cond) {
357	<a name="l00524"></a>00524 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> == val.length(), <span class="stringliteral">"Wrong val"</span>);
358	<a name="l00525"></a>00525 it_assert_debug(<a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &quot;Up&quot; should be.">upsize</a> == val_up.length(), <span class="stringliteral">"Wrong val_up"</span>);
359	<a name="l00526"></a>00526 set_subvector(val_up, <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indices of the upper val">v2v_up</a>, val);
360	<a name="l00527"></a>00527 set_subvector(val_up, <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indices of the upper val">v2c_up</a>, cond);
361	<a name="l00528"></a>00528 }
362	<a name="l00529"></a>00529 };
363	<a name="l00530"></a>00530
364	<a name="l00533"></a><a class="code" href="classbdm_1_1datalink__m2m.html">00533</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1datalink__m2m.html">datalink_m2m</a>: <span class="keyword">public</span> <a class="code" href="classbdm_1_1datalink__m2e.html" title="Data link with a condition.">datalink_m2e</a>
365	<a name="l00534"></a>00534 {
366	<a name="l00535"></a>00535 <span class="keyword">protected</span>:
367	<a name="l00537"></a><a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d">00537</a> ivec <a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indices of the upper cond">c2c_up</a>;
368	<a name="l00539"></a><a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5">00539</a> ivec <a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indices of the local cond">c2c_lo</a>;
369	<a name="l00540"></a>00540
370	<a name="l00541"></a>00541 <span class="keyword">public</span>:
371	<a name="l00543"></a><a class="code" href="classbdm_1_1datalink__m2m.html#3baed32cb66f0e2f731032cec0a4d287">00543</a> <a class="code" href="classbdm_1_1datalink__m2m.html#3baed32cb66f0e2f731032cec0a4d287" title="Constructor.">datalink_m2m</a>() {};
372	<a name="l00544"></a>00544 <span class="keywordtype">void</span> set_connection(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rvc, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv_up, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rvc_up) {
373	<a name="l00545"></a>00545 datalink_m2e::set_connection(rv, rvc, rv_up);
374	<a name="l00546"></a>00546 <span class="comment">//establish c2c connection</span>
375	<a name="l00547"></a>00547 rvc.<a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a>(rvc_up, <a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indices of the local cond">c2c_lo</a>, <a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indices of the upper cond">c2c_up</a>);
376	<a name="l00548"></a>00548 it_assert_debug(<a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indices of the local cond">c2c_lo</a>.length() + <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, indices of the local cond">v2c_lo</a>.length() == <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>, <span class="stringliteral">"cond is not fully given"</span>);
377	<a name="l00549"></a>00549 }
378	<a name="l00550"></a>00550
379	<a name="l00552"></a><a class="code" href="classbdm_1_1datalink__m2m.html#911fb7fcfa80d4e31db8e0cc03141a90">00552</a> vec <a class="code" href="classbdm_1_1datalink__m2m.html#911fb7fcfa80d4e31db8e0cc03141a90" title="Get cond for myself from val and cond of &quot;Up&quot;.">get_cond</a>(<span class="keyword">const</span> vec &val_up, <span class="keyword">const</span> vec &cond_up) {
380	<a name="l00553"></a>00553 vec tmp(<a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>);
381	<a name="l00554"></a>00554 set_subvector(tmp, <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, indices of the local cond">v2c_lo</a>, val_up(<a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indices of the upper val">v2c_up</a>));
382	<a name="l00555"></a>00555 set_subvector(tmp, <a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indices of the local cond">c2c_lo</a>, cond_up(<a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indices of the upper cond">c2c_up</a>));
383	<a name="l00556"></a>00556 <span class="keywordflow">return</span> tmp;
384	<a name="l00557"></a>00557 }
385	<a name="l00559"></a>00559
386	<a name="l00560"></a>00560 };
387	<a name="l00561"></a>00561
388	<a name="l00567"></a><a class="code" href="classbdm_1_1logger.html">00567</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
389	<a name="l00568"></a>00568 {
390	<a name="l00569"></a>00569 <span class="keyword">protected</span>:
391	<a name="l00571"></a><a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f">00571</a> Array<RV> <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>;
392	<a name="l00573"></a><a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02">00573</a> Array<string> <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>;
393	<a name="l00574"></a>00574 <span class="keyword">public</span>:
394	<a name="l00576"></a><a class="code" href="classbdm_1_1logger.html#f9c2fe7db5ddbda0098a4a5bd2c29fbf">00576</a> <a class="code" href="classbdm_1_1logger.html#f9c2fe7db5ddbda0098a4a5bd2c29fbf" title="Default constructor.">logger</a>() : <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>(0), <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>(0) {}
395	<a name="l00577"></a>00577
396	<a name="l00580"></a><a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">00580</a> <span class="keyword">virtual</span> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &rv, <span class="keywordtype">string</span> prefix = <span class="stringliteral">""</span>) {
397	<a name="l00581"></a>00581 <span class="keywordtype">int</span> id;
398	<a name="l00582"></a>00582 <span class="keywordflow">if</span> (rv._dsize() > 0) {
399	<a name="l00583"></a>00583 <span class="keywordtype">id</span> = <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>.length();
400	<a name="l00584"></a>00584 <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a> = concat(<a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>, prefix); <span class="comment">// diff</span>
401	<a name="l00585"></a>00585 <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>.set_length(<span class="keywordtype">id</span> + 1, <span class="keyword">true</span>);
402	<a name="l00586"></a>00586 <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>(<span class="keywordtype">id</span>) = rv;
403	<a name="l00587"></a>00587 }
404	<a name="l00588"></a>00588 <span class="keywordflow">else</span> { <span class="keywordtype">id</span> = -1;}
405	<a name="l00589"></a>00589 <span class="keywordflow">return</span> id; <span class="comment">// identifier of the last entry</span>
406	<a name="l00590"></a>00590 }
407	<a name="l00591"></a>00591
408	<a name="l00593"></a>00593 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<span class="keywordtype">int</span> <span class="keywordtype">id</span>, <span class="keyword">const</span> vec &v) = 0;
409	<a name="l00595"></a>00595 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<span class="keywordtype">int</span> <span class="keywordtype">id</span>, <span class="keyword">const</span> <span class="keywordtype">double</span> &d) = 0;
410	<a name="l00596"></a>00596
411	<a name="l00598"></a>00598 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1logger.html#4f31c03578a7b9b30d8b1f8c338a37e8" title="Shifts storage position for another time step.">step</a>() = 0;
412	<a name="l00599"></a>00599
413	<a name="l00601"></a><a class="code" href="classbdm_1_1logger.html#0e3fc7d13e23a8f71fc7c8a5ea68559d">00601</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1logger.html#0e3fc7d13e23a8f71fc7c8a5ea68559d" title="Finalize storing information.">finalize</a>() {};
414	<a name="l00602"></a>00602
415	<a name="l00604"></a><a class="code" href="classbdm_1_1logger.html#cd6291a73491f755dc7bc88d022eb2fc">00604</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1logger.html#cd6291a73491f755dc7bc88d022eb2fc" title="Initialize the storage.">init</a>() {};
416	<a name="l00605"></a>00605
417	<a name="l00606"></a>00606 };
418	<a name="l00607"></a>00607
419	<a name="l00611"></a><a class="code" href="classbdm_1_1mepdf.html">00611</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1mepdf.html" title="Unconditional mpdf, allows using epdf in the role of mpdf.">mepdf</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> {
420	<a name="l00612"></a>00612 <span class="keyword">public</span>:
421	<a name="l00614"></a><a class="code" href="classbdm_1_1mepdf.html#4a2b01e7bb510c6b8946757cae6c9155">00614</a> <a class="code" href="classbdm_1_1mepdf.html#4a2b01e7bb510c6b8946757cae6c9155" title="Default constructor.">mepdf</a>() { }
422	<a name="l00615"></a>00615
423	<a name="l00616"></a>00616 <a class="code" href="classbdm_1_1mepdf.html#4a2b01e7bb510c6b8946757cae6c9155" title="Default constructor.">mepdf</a>(<a class="code" href="classbdm_1_1shared__ptr.html">shared_ptr<epdf></a> em) {
424	<a name="l00617"></a>00617 set_ep(em);
425	<a name="l00618"></a>00618 <a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a> = 0;
426	<a name="l00619"></a>00619 }
427	<a name="l00620"></a>00620
428	<a name="l00622"></a>00622 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mepdf.html#b2b898058169e7ed9a6dc7cd4e1c147c" title="empty">condition</a>(<span class="keyword">const</span> vec &cond);
429	<a name="l00623"></a>00623
430	<a name="l00631"></a>00631 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mepdf.html#27ef18a3e3b1b5de02739d1324bf25fb">from_setting</a>(<span class="keyword">const</span> Setting &<span class="keyword">set</span>);
431	<a name="l00632"></a>00632 };
432	<a name="l00633"></a>00633 <a class="code" href="user__info_8h.html#4f9de2f17e844047726487b99def99c6" title="Macro for registration of class into map of user-infos, registered class is scriptable...">UIREGISTER</a>(mepdf);
433	<a name="l00634"></a>00634
434	<a name="l00638"></a><a class="code" href="classbdm_1_1compositepdf.html">00638</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1compositepdf.html" title="Chain rule of pdfs - abstract part common for mprod and merger. this abstract class...">compositepdf</a> {
435	<a name="l00639"></a>00639 <span class="keyword">protected</span>:
436	<a name="l00641"></a><a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a">00641</a> Array<mpdf*> <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>;
437	<a name="l00642"></a>00642 <span class="keywordtype">bool</span> owning_mpdfs;
438	<a name="l00643"></a>00643 <span class="keyword">public</span>:
439	<a name="l00644"></a>00644 <a class="code" href="classbdm_1_1compositepdf.html" title="Chain rule of pdfs - abstract part common for mprod and merger. this abstract class...">compositepdf</a>():<a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>(0){};
440	<a name="l00645"></a>00645 <a class="code" href="classbdm_1_1compositepdf.html" title="Chain rule of pdfs - abstract part common for mprod and merger. this abstract class...">compositepdf</a>(Array<mpdf*> A0, <span class="keywordtype">bool</span> own=<span class="keyword">false</span>){set_elements(A0,own);};
441	<a name="l00646"></a>00646 <span class="keywordtype">void</span> set_elements(Array<mpdf*> A0, <span class="keywordtype">bool</span> own=<span class="keyword">false</span>) {<a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>=A0;owning_mpdfs=own;};
442	<a name="l00648"></a>00648 RV <a class="code" href="classbdm_1_1compositepdf.html#c73d39acc4378eee6a63155c3517d3c9" title="find common rv, flag">getrv</a>(<span class="keywordtype">bool</span> checkoverlap = <span class="keyword">false</span>);
443	<a name="l00650"></a>00650 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1compositepdf.html#c71fc2c51f49d797e61f479f543c75ce" title="common rvc of all mpdfs is written to rvc">setrvc</a>(<span class="keyword">const</span> RV &rv, RV &rvc);
444	<a name="l00651"></a>00651 ~compositepdf(){<span class="keywordflow">if</span> (owning_mpdfs) <span class="keywordflow">for</span>(<span class="keywordtype">int</span> i=0;i<<a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>.length();i++){<span class="keyword">delete</span> <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>(i);}};
445	<a name="l00652"></a>00652 };
446	<a name="l00653"></a>00653
447	<a name="l00661"></a><a class="code" href="classbdm_1_1DS.html">00661</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1DS.html" title="Abstract class for discrete-time sources of data.">DS</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
448	<a name="l00662"></a>00662 {
449	<a name="l00663"></a>00663 <span class="keyword">protected</span>:
450	<a name="l00664"></a>00664 <span class="keywordtype">int</span> dtsize;
451	<a name="l00665"></a>00665 <span class="keywordtype">int</span> utsize;
452	<a name="l00667"></a><a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec">00667</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>;
453	<a name="l00669"></a><a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702">00669</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>; <span class="comment">//</span>
454	<a name="l00671"></a><a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867">00671</a> <span class="comment"></span> <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867" title="Remember its own index in Logger L.">L_dt</a>, L_ut;
455	<a name="l00672"></a>00672 <span class="keyword">public</span>:
456	<a name="l00674"></a><a class="code" href="classbdm_1_1DS.html#bedd82e08ae6686b4b1732a122c2201a">00674</a> <a class="code" href="classbdm_1_1DS.html#bedd82e08ae6686b4b1732a122c2201a" title="default constructors">DS</a>() : <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>(), <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>() {};
457	<a name="l00676"></a><a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba">00676</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba" title="Returns full vector of observed data=[output, input].">getdata</a>(vec &dt) {it_error(<span class="stringliteral">"abstract class"</span>);};
458	<a name="l00678"></a><a class="code" href="classbdm_1_1DS.html#a0bf451c7fea73ff4150ff9e83da4832">00678</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#a0bf451c7fea73ff4150ff9e83da4832" title="Returns data records at indeces.">getdata</a>(vec &dt, <span class="keyword">const</span> ivec &indeces) {it_error(<span class="stringliteral">"abstract class"</span>);};
459	<a name="l00680"></a><a class="code" href="classbdm_1_1DS.html#f5baeee8528e5e4c5c7a344b2e3f192a">00680</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#f5baeee8528e5e4c5c7a344b2e3f192a" title="Accepts action variable and schedule it for application.">write</a>(vec &ut) {it_error(<span class="stringliteral">"abstract class"</span>);};
460	<a name="l00682"></a><a class="code" href="classbdm_1_1DS.html#6bced13f84f6bf04449728817e049e1d">00682</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#6bced13f84f6bf04449728817e049e1d" title="Accepts action variables at specific indeces.">write</a>(vec &ut, <span class="keyword">const</span> ivec &indeces) {it_error(<span class="stringliteral">"abstract class"</span>);};
461	<a name="l00683"></a>00683
462	<a name="l00685"></a>00685 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#9a14c64d881e838ddc3a5617ee8a973b" title="Moves from to , i.e. perfroms the actions and reads response of the system.">step</a>() = 0;
463	<a name="l00686"></a>00686
464	<a name="l00688"></a><a class="code" href="classbdm_1_1DS.html#3c9e417d85052396fe74c164857ee67e">00688</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#3c9e417d85052396fe74c164857ee67e" title="Register DS for logging into logger L.">log_add</a>(<a class="code" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a> &L) {
465	<a name="l00689"></a>00689 it_assert_debug(dtsize == <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>._dsize(), <span class="stringliteral">""</span>);
466	<a name="l00690"></a>00690 it_assert_debug(utsize == <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>._dsize(), <span class="stringliteral">""</span>);
467	<a name="l00691"></a>00691
468	<a name="l00692"></a>00692 <a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867" title="Remember its own index in Logger L.">L_dt</a> = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>, <span class="stringliteral">""</span>);
469	<a name="l00693"></a>00693 L_ut = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>, <span class="stringliteral">""</span>);
470	<a name="l00694"></a>00694 }
471	<a name="l00696"></a><a class="code" href="classbdm_1_1DS.html#e03d2c93cddd5c3a03107ddd3ae999cf">00696</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#e03d2c93cddd5c3a03107ddd3ae999cf" title="Register DS for logging into logger L.">logit</a>(<a class="code" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a> &L) {
472	<a name="l00697"></a>00697 vec tmp(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>._dsize() + <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>._dsize());
473	<a name="l00698"></a>00698 <a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba" title="Returns full vector of observed data=[output, input].">getdata</a>(tmp);
474	<a name="l00699"></a>00699 <span class="comment">// d is first in getdata</span>
475	<a name="l00700"></a>00700 L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867" title="Remember its own index in Logger L.">L_dt</a>, tmp.left(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>._dsize()));
476	<a name="l00701"></a>00701 <span class="comment">// u follows after d in getdata</span>
477	<a name="l00702"></a>00702 L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(L_ut, tmp.mid(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>._dsize(), <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>._dsize()));
478	<a name="l00703"></a>00703 }
479	<a name="l00705"></a><a class="code" href="classbdm_1_1DS.html#ed2113851294ed96abe60fc3f05aa660">00705</a> <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1DS.html#ed2113851294ed96abe60fc3f05aa660" title="access function">_drv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> concat(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>, <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>);}
480	<a name="l00707"></a><a class="code" href="classbdm_1_1DS.html#a6a9a375da64c82fc227bda3418a1f33">00707</a> <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classbdm_1_1DS.html#a6a9a375da64c82fc227bda3418a1f33" title="access function">_urv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>;}
481	<a name="l00709"></a><a class="code" href="classbdm_1_1DS.html#821297dfd839d68691324f024e57252f">00709</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1DS.html#821297dfd839d68691324f024e57252f" title="set random rvariables">set_drv</a>(<span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &drv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &urv) { <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a> = drv; <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a> = urv;}
482	<a name="l00710"></a>00710 };
483	<a name="l00711"></a>00711
484	<a name="l00733"></a><a class="code" href="classbdm_1_1BM.html">00733</a> <span class="keyword">class </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> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1root.html" title="Root class of BDM objects.">root</a>
485	<a name="l00734"></a>00734 {
486	<a name="l00735"></a>00735 <span class="keyword">protected</span>:
487	<a name="l00737"></a><a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed">00737</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a>;
488	<a name="l00739"></a><a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a">00739</a> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>;
489	<a name="l00741"></a><a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee">00741</a> <span class="keywordtype">bool</span> <a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>;
490	<a name="l00742"></a>00742 <span class="keyword">public</span>:
491	<a name="l00745"></a>00745
492	<a name="l00746"></a>00746 <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of a system, i.e. all uncertainty is modeled by probabilities.">BM</a>() : <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>(0), <a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>(true), <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(4), <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(4) {
493	<a name="l00747"></a>00747 <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a> = -1;<span class="comment">/empty IDs/</span>
494	<a name="l00748"></a>00748 <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a> = 0;
495	<a name="l00749"></a>00749 <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(0) = 1;<span class="comment">/log only mean/</span>
496	<a name="l00750"></a>00750 };
497	<a name="l00751"></a>00751 <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of a system, i.e. all uncertainty is modeled by probabilities.">BM</a>(<span class="keyword">const</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> &B) : <a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a>(B.<a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a>), <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>(B.<a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>), <a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>(B.<a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a>) {}
498	<a name="l00754"></a><a class="code" href="classbdm_1_1BM.html#40f01563651c1e4fe926189e41a6b8ab">00754</a> <span class="keyword">virtual</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>* <a class="code" href="classbdm_1_1BM.html#40f01563651c1e4fe926189e41a6b8ab">_copy_</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> NULL;};
499	<a name="l00756"></a>00756
500	<a name="l00759"></a>00759
501	<a name="l00763"></a>00763 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#60b1779a577367c369a932cabd3a6188" title="Incremental Bayes rule.">bayes</a>(<span class="keyword">const</span> vec &dt) = 0;
502	<a name="l00765"></a>00765 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#1dee3fddaf021e62d925289660a707dc" title="Batch Bayes rule (columns of Dt are observations).">bayesB</a>(<span class="keyword">const</span> mat &Dt);
503	<a name="l00768"></a><a class="code" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">00768</a> <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">logpred</a>(<span class="keyword">const</span> vec &dt)<span class="keyword"> const </span>{it_error(<span class="stringliteral">"Not implemented"</span>); <span class="keywordflow">return</span> 0.0;}
504	<a name="l00770"></a><a class="code" href="classbdm_1_1BM.html#0e8ebe61fb14990abe1254bd3dda5fae">00770</a> vec <a class="code" href="classbdm_1_1BM.html#0e8ebe61fb14990abe1254bd3dda5fae" title="Matrix version of logpred.">logpred_m</a>(<span class="keyword">const</span> mat &dt)<span class="keyword"> const </span>{vec tmp(dt.cols()); <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i < dt.cols(); i++) {tmp(i) = <a class="code" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">logpred</a>(dt.get_col(i));} <span class="keywordflow">return</span> tmp;}
505	<a name="l00771"></a>00771
506	<a name="l00773"></a><a class="code" href="classbdm_1_1BM.html#688d7a2aced1e06aa1c468d73a9e5eba">00773</a> <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>* <a class="code" href="classbdm_1_1BM.html#688d7a2aced1e06aa1c468d73a9e5eba" title="Constructs a predictive density .">epredictor</a>()<span class="keyword"> const </span>{it_error(<span class="stringliteral">"Not implemented"</span>); <span class="keywordflow">return</span> NULL;};
507	<a name="l00775"></a><a class="code" href="classbdm_1_1BM.html#598b25e3f3d96a5bc00a5faeb5b3c912">00775</a> <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a>* <a class="code" href="classbdm_1_1BM.html#598b25e3f3d96a5bc00a5faeb5b3c912" title="Constructs a conditional density 1-step ahead predictor.">predictor</a>()<span class="keyword"> const </span>{it_error(<span class="stringliteral">"Not implemented"</span>); <span class="keywordflow">return</span> NULL;};
508	<a name="l00777"></a>00777
509	<a name="l00782"></a>00782
510	<a name="l00784"></a><a class="code" href="classbdm_1_1BM.html#de95c590d650b2db423d53992f0906cb">00784</a> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> rvc;
511	<a name="l00786"></a><a class="code" href="classbdm_1_1BM.html#b5e8e016be16e4029c789063cb3f3b33">00786</a> <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>& <a class="code" href="classbdm_1_1BM.html#b5e8e016be16e4029c789063cb3f3b33" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> rvc;}
512	<a name="l00787"></a>00787
513	<a name="l00789"></a><a class="code" href="classbdm_1_1BM.html#6799f4b16a6a59ed58b1d0d6e17116f4">00789</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#6799f4b16a6a59ed58b1d0d6e17116f4" title="Substitute val for rvc.">condition</a>(<span class="keyword">const</span> vec &val) {it_error(<span class="stringliteral">"Not implemented!"</span>);};
514	<a name="l00790"></a>00790
515	<a name="l00792"></a>00792
516	<a name="l00793"></a>00793
517	<a name="l00796"></a>00796
518	<a name="l00797"></a>00797 <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>& _drv()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a>;}
519	<a name="l00798"></a>00798 <span class="keywordtype">void</span> set_drv(<span class="keyword">const</span> RV &rv) {<a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a> = rv;}
520	<a name="l00799"></a>00799 <span class="keywordtype">void</span> set_rv(<span class="keyword">const</span> RV &rv) {<span class="keyword">const_cast<</span>epdf&<span class="keyword">></span>(posterior()).set_rv(rv);}
521	<a name="l00800"></a>00800 <span class="keywordtype">double</span> _ll()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>;}
522	<a name="l00801"></a>00801 <span class="keywordtype">void</span> set_evalll(<span class="keywordtype">bool</span> evl0) {<a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee" title="If true, the filter will compute likelihood of the data record and store it in ll...">evalll</a> = evl0;}
523	<a name="l00802"></a>00802 <span class="keyword">virtual</span> <span class="keyword">const</span> epdf& posterior() <span class="keyword">const</span> = 0;
524	<a name="l00803"></a>00803 <span class="keyword">virtual</span> <span class="keyword">const</span> epdf* _e() <span class="keyword">const</span> = 0;
525	<a name="l00805"></a>00805
526	<a name="l00808"></a>00808
527	<a name="l00810"></a><a class="code" href="classbdm_1_1BM.html#adb539fa2585c8a2a91c54e5aebf330e">00810</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#adb539fa2585c8a2a91c54e5aebf330e" title="Set boolean options from a string, recognized are: &quot;logbounds,logll&quot;.">set_options</a>(<span class="keyword">const</span> <span class="keywordtype">string</span> &opt) {
528	<a name="l00811"></a>00811 <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(0) = 1;
529	<a name="l00812"></a>00812 <span class="keywordflow">if</span> (opt.find(<span class="stringliteral">"logbounds"</span>) != string::npos) {<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(1) = 1; <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(2) = 1;}
530	<a name="l00813"></a>00813 <span class="keywordflow">if</span> (opt.find(<span class="stringliteral">"logll"</span>) != string::npos) {<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(3) = 1;}
531	<a name="l00814"></a>00814 }
532	<a name="l00816"></a><a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca">00816</a> ivec <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>;
533	<a name="l00817"></a>00817
534	<a name="l00819"></a><a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207">00819</a> ivec <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>;
535	<a name="l00821"></a><a class="code" href="classbdm_1_1BM.html#2298ca6af9b13a78a4c9e18ab43f1827">00821</a> <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#2298ca6af9b13a78a4c9e18ab43f1827" title="Add all logged variables to a logger.">log_add</a>(<a class="code" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a> &L, <span class="keyword">const</span> <span class="keywordtype">string</span> &name = <span class="stringliteral">""</span>) {
536	<a name="l00822"></a>00822 <span class="comment">// internal</span>
537	<a name="l00823"></a>00823 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> r;
538	<a name="l00824"></a>00824 <span class="keywordflow">if</span> (posterior().isnamed()) {r = posterior().<a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21" title="Return name (fails when isnamed is false).">_rv</a>();}
539	<a name="l00825"></a>00825 <span class="keywordflow">else</span> {r = <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>(<span class="stringliteral">"est"</span>, posterior().dimension());};
540	<a name="l00826"></a>00826
541	<a name="l00827"></a>00827 <span class="comment">// Add mean value</span>
542	<a name="l00828"></a>00828 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(0)) <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(0) = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(r, name + <span class="stringliteral">"mean_"</span>);
543	<a name="l00829"></a>00829 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(1)) <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(1) = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(r, name + <span class="stringliteral">"lb_"</span>);
544	<a name="l00830"></a>00830 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(2)) <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(2) = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(r, name + <span class="stringliteral">"ub_"</span>);
545	<a name="l00831"></a>00831 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(3)) <a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(3) = L.<a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">add</a>(<a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>(<span class="stringliteral">"ll"</span>, 1), name); <span class="comment">//TODO: "local" RV</span>
546	<a name="l00832"></a>00832 }
547	<a name="l00833"></a>00833 <span class="keyword">virtual</span> <span class="keywordtype">void</span> logit(<a class="code" href="classbdm_1_1logger.html" title="Class for storing results (and semi-results) of an experiment.">logger</a> &L) {
548	<a name="l00834"></a>00834 L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(0), posterior().mean());
549	<a name="l00835"></a>00835 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(1) \|\| <a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(2)) { <span class="comment">//if one of them is off, its LID==-1 and will not be stored</span>
550	<a name="l00836"></a>00836 vec ub, lb;
551	<a name="l00837"></a>00837 posterior().<a class="code" href="classbdm_1_1epdf.html#0a7705e80601957f7226d333a0a6d539" title="Lower and upper bounds of percentage % quantile, returns mean-2*sigma as default...">qbounds</a>(lb, ub);
552	<a name="l00838"></a>00838 L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(1), lb);
553	<a name="l00839"></a>00839 L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(2), ub);
554	<a name="l00840"></a>00840 }
555	<a name="l00841"></a>00841 <span class="keywordflow">if</span> (<a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207" title="Flags for logging - same size as LIDs, each entry correspond to the same in LIDs...">LFlags</a>(3)) L.<a class="code" href="classbdm_1_1logger.html#397c014b0aaf3581516e2a105d26a5c9" title="log this vector">logit</a>(<a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca" title="IDs of storages in loggers 4:[1=mean,2=lb,3=ub,4=ll].">LIDs</a>(3), <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>);
556	<a name="l00842"></a>00842 }
557	<a name="l00844"></a>00844 };
558	<a name="l00845"></a>00845
559	<a name="l00846"></a>00846
560	<a name="l00847"></a>00847 }; <span class="comment">//namespace</span>
561	<a name="l00848"></a>00848 <span class="preprocessor">#endif // BDMBASE_H</span>
562	</pre></div></div>
563	<hr size="1"><address style="text-align: right;"><small>Generated on Wed Aug 5 00:06:46 2009 for mixpp by
564	<a href="http://www.doxygen.org/index.html">
565	<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.9 </small></address>
566	</body>
567	</html>

Note: See TracBrowser for help on using the browser.

Context Navigation

root/library/doc/html/bdmbase_8h_source.html @ 472

Download in other formats: