root/doc/html/libBM_8h-source.html @ 264

Revision 264, 82.7 kB (checked in by smidl, 16 years ago)

doc

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: libBM.h Source File</title>
4<link href="doxygen.css" rel="stylesheet" type="text/css">
5<link href="tabs.css" rel="stylesheet" type="text/css">
6</head><body>
7<!-- Generated by Doxygen 1.5.6 -->
8<div class="navigation" id="top">
9  <div class="tabs">
10    <ul>
11      <li><a href="index.html"><span>Main&nbsp;Page</span></a></li>
12      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
13      <li><a href="modules.html"><span>Modules</span></a></li>
14      <li><a href="namespaces.html"><span>Namespaces</span></a></li>
15      <li><a href="classes.html"><span>Classes</span></a></li>
16      <li class="current"><a href="files.html"><span>Files</span></a></li>
17    </ul>
18  </div>
19<h1>libBM.h</h1><a href="libBM_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
20<a name="l00017"></a>00017 <span class="preprocessor">#ifndef BM_H</span>
21<a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define BM_H</span>
22<a name="l00019"></a>00019 <span class="preprocessor"></span>
23<a name="l00020"></a>00020
24<a name="l00021"></a>00021 <span class="preprocessor">#include "../itpp_ext.h"</span>
25<a name="l00022"></a>00022 <span class="comment">//#include &lt;std&gt;</span>
26<a name="l00023"></a>00023
27<a name="l00024"></a>00024 <span class="keyword">namespace </span>bdm {
28<a name="l00025"></a>00025         <span class="keyword">using namespace </span>itpp;
29<a name="l00026"></a>00026         <span class="keyword">using</span> std::string;
30<a name="l00027"></a>00027
31<a name="l00029"></a><a class="code" href="classbdm_1_1bdmroot.html">00029</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
32<a name="l00030"></a>00030                 <span class="keyword">virtual</span> <span class="keywordtype">void</span> print() {}
33<a name="l00031"></a>00031         };
34<a name="l00032"></a>00032
35<a name="l00034"></a><a class="code" href="classbdm_1_1str.html">00034</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1str.html" title="Structure of RV (used internally), i.e. expanded RVs.">str</a> {
36<a name="l00035"></a>00035         <span class="keyword">public</span>:
37<a name="l00037"></a><a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0">00037</a>                 ivec <a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0" title="vector id ids (non-unique!)">ids</a>;
38<a name="l00039"></a><a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a">00039</a>                 ivec <a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a" title="vector of times">times</a>;
39<a name="l00041"></a><a class="code" href="classbdm_1_1str.html#7c6b9dccc6814d7bd1e8e2208a0d3f3f">00041</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 ) {
40<a name="l00042"></a>00042                         it_assert_debug ( times0.length() ==ids0.length(),<span class="stringliteral">"Incompatible input"</span> );
41<a name="l00043"></a>00043                 };
42<a name="l00044"></a>00044         };
43<a name="l00045"></a>00045
44<a name="l00052"></a><a class="code" href="classbdm_1_1RV.html">00052</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
45<a name="l00053"></a>00053         <span class="keyword">protected</span>:
46<a name="l00055"></a><a class="code" href="classbdm_1_1RV.html#38203096d950671def44de8cf1c6bb2e">00055</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#38203096d950671def44de8cf1c6bb2e" title="size = sum of sizes">tsize</a>;
47<a name="l00057"></a><a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385">00057</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="len = number of individual rvs">len</a>;
48<a name="l00059"></a><a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b">00059</a>                 ivec <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>;
49<a name="l00061"></a><a class="code" href="classbdm_1_1RV.html#0112d25045f6f8fabdc9e821029928ee">00061</a>                 ivec <a class="code" href="classbdm_1_1RV.html#0112d25045f6f8fabdc9e821029928ee" title="Vector of sizes.">sizes</a>;
50<a name="l00063"></a><a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896">00063</a>                 ivec <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>;
51<a name="l00065"></a><a class="code" href="classbdm_1_1RV.html#9a87c46918c3c0db0bdca854a48a44aa">00065</a>                 Array&lt;std::string&gt; <a class="code" href="classbdm_1_1RV.html#9a87c46918c3c0db0bdca854a48a44aa" title="Array of names.">names</a>;
52<a name="l00066"></a>00066
53<a name="l00067"></a>00067         <span class="keyword">private</span>:
54<a name="l00069"></a>00069                 <span class="keywordtype">void</span> init ( ivec in_ids, Array&lt;std::string&gt; in_names, ivec in_sizes, ivec in_times );
55<a name="l00070"></a>00070         <span class="keyword">public</span>:
56<a name="l00072"></a>00072                 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> ( Array&lt;std::string&gt; in_names, ivec in_sizes, ivec in_times );
57<a name="l00074"></a>00074                 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> ( Array&lt;std::string&gt; in_names, ivec in_sizes );
58<a name="l00076"></a>00076                 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> ( Array&lt;std::string&gt; in_names );
59<a name="l00078"></a>00078                 <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> ();
60<a name="l00079"></a>00079
61<a name="l00081"></a>00081                 <span class="keyword">friend</span> std::ostream &amp;<a class="code" href="classbdm_1_1RV.html#cad6f15535d10437f6690cb4492176a8" title="Printing output e.g. for debugging.">operator&lt;&lt; </a>( std::ostream &amp;os, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv );
62<a name="l00082"></a>00082
63<a name="l00084"></a><a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67">00084</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#38203096d950671def44de8cf1c6bb2e" title="size = sum of sizes">tsize</a>;} ;
64<a name="l00086"></a><a class="code" href="classbdm_1_1RV.html#e9ec8c3e756651ff352ab5e3d3acda4b">00086</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#e9ec8c3e756651ff352ab5e3d3acda4b" title="Return length (number of entries) of the RV.">length</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="len = number of individual rvs">len</a>;} ;
65<a name="l00087"></a>00087
66<a name="l00088"></a>00088                 <span class="comment">//TODO why not inline and later??</span>
67<a name="l00089"></a>00089
68<a name="l00091"></a>00091                 ivec <a class="code" href="classbdm_1_1RV.html#1367055b2608ad9d50d01a2da6cec869" title="Find indexes of self in another rv,.">findself</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv2 ) <span class="keyword">const</span>;
69<a name="l00093"></a>00093                 <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" title="Class representing variables, most often random variables.">RV</a> &amp;rv2 ) <span class="keyword">const</span>;
70<a name="l00095"></a>00095                 <span class="keywordtype">bool</span> <span class="keyword">add</span> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv2 );
71<a name="l00097"></a>00097                 <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#aec44dabdf0a6d90fbae95e1356eda39" title="Subtract another variable from the current one.">subt</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv2 ) <span class="keyword">const</span>;
72<a name="l00099"></a>00099                 <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#c0c04ab99ded62e1e05bcefb151ebeb9" title="Select only variables at indeces ind.">subselect</a> ( <span class="keyword">const</span> ivec &amp;ind ) <span class="keyword">const</span>;
73<a name="l00101"></a>00101                 <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 indeces ind.">operator() </a>( <span class="keyword">const</span> ivec &amp;ind ) <span class="keyword">const</span>;
74<a name="l00103"></a>00103                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#bd77cb0899f7fba1d0e34396137b7f1f" title="Shift time shifted by delta.">t</a> ( <span class="keywordtype">int</span> delta );
75<a name="l00105"></a>00105                 <a class="code" href="classbdm_1_1str.html" title="Structure of RV (used internally), i.e. expanded RVs.">str</a> <a class="code" href="classbdm_1_1RV.html#8bdd7f4072749c7cc864b78a418ca388" title="generate str from rv, by expanding sizes">tostr</a>() <span class="keyword">const</span>;
76<a name="l00108"></a>00108                 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> &amp;crv ) <span class="keyword">const</span>;
77<a name="l00111"></a>00111                 <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> &amp;rv2, ivec &amp;selfi, ivec &amp;rv2i ) <span class="keyword">const</span>;
78<a name="l00113"></a><a class="code" href="classbdm_1_1RV.html#8bfda300b769a05598c55cb94fe61f52">00113</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#8bfda300b769a05598c55cb94fe61f52" title="Minimum time-offset.">mint</a> () {<span class="keywordflow">return</span> min ( <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a> );};
79<a name="l00114"></a>00114
80<a name="l00116"></a><a class="code" href="classbdm_1_1RV.html#ae0c6e729851a43c82a5b1c199474ef2">00116</a>                 Array&lt;std::string&gt;&amp; <a class="code" href="classbdm_1_1RV.html#ae0c6e729851a43c82a5b1c199474ef2" title="access function">_names</a>() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#9a87c46918c3c0db0bdca854a48a44aa" title="Array of names.">names</a>;};
81<a name="l00117"></a>00117
82<a name="l00119"></a><a class="code" href="classbdm_1_1RV.html#172390b8525ef5aac684228241fff7b5">00119</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#172390b8525ef5aac684228241fff7b5" title="access function">id</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a> ( at );};
83<a name="l00121"></a><a class="code" href="classbdm_1_1RV.html#8e64d02005a02aef72f214286977d8de">00121</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#8e64d02005a02aef72f214286977d8de" title="access function">size</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#0112d25045f6f8fabdc9e821029928ee" title="Vector of sizes.">sizes</a> ( at );};
84<a name="l00123"></a><a class="code" href="classbdm_1_1RV.html#aff4de8185750b506492623847520518">00123</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#aff4de8185750b506492623847520518" title="access function">time</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a> ( at );};
85<a name="l00125"></a><a class="code" href="classbdm_1_1RV.html#d05af87595cfb7fa2e0fadeae5597fa8">00125</a>                 std::string <a class="code" href="classbdm_1_1RV.html#d05af87595cfb7fa2e0fadeae5597fa8" title="access function">name</a> ( <span class="keywordtype">int</span> at ) {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#9a87c46918c3c0db0bdca854a48a44aa" title="Array of names.">names</a> ( at );};
86<a name="l00126"></a>00126
87<a name="l00128"></a><a class="code" href="classbdm_1_1RV.html#602cea0aecebc5146d81657df5127543">00128</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#602cea0aecebc5146d81657df5127543" title="access function">set_id</a> ( <span class="keywordtype">int</span> at, <span class="keywordtype">int</span> id0 ) {<a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a> ( at ) =id0;};
88<a name="l00130"></a><a class="code" href="classbdm_1_1RV.html#0b2c9e73ff66847c3644ebc3eb559a03">00130</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#0b2c9e73ff66847c3644ebc3eb559a03" title="access function">set_size</a> ( <span class="keywordtype">int</span> at, <span class="keywordtype">int</span> size0 ) {<a class="code" href="classbdm_1_1RV.html#0112d25045f6f8fabdc9e821029928ee" title="Vector of sizes.">sizes</a> ( at ) =size0; <a class="code" href="classbdm_1_1RV.html#38203096d950671def44de8cf1c6bb2e" title="size = sum of sizes">tsize</a>=sum ( <a class="code" href="classbdm_1_1RV.html#0112d25045f6f8fabdc9e821029928ee" title="Vector of sizes.">sizes</a> );};
89<a name="l00132"></a><a class="code" href="classbdm_1_1RV.html#fb2892d12dc978549e9491aed327dfd4">00132</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#fb2892d12dc978549e9491aed327dfd4" title="access function">set_time</a> ( <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;};
90<a name="l00133"></a>00133
91<a name="l00135"></a>00135                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1RV.html#176e9c0c35b51fee0ad8813a39883a63" title="Assign unused ids to this rv.">newids</a>();
92<a name="l00136"></a>00136         };
93<a name="l00137"></a>00137
94<a name="l00139"></a>00139         RV <a class="code" href="namespacebdm.html#b9016687c0e874ca5cdcf75ae28811aa" title="Concat two random variables.">concat</a> ( <span class="keyword">const</span> RV &amp;rv1, <span class="keyword">const</span> RV &amp;rv2 );
95<a name="l00140"></a>00140
96<a name="l00142"></a>00142         <span class="keyword">extern</span> RV <a class="code" href="namespacebdm.html#bdd4760f14ec7d6fea083aa2170f6236" title="Default empty RV that can be used as default argument.">RV0</a>;
97<a name="l00143"></a>00143
98<a name="l00145"></a>00145
99<a name="l00146"></a><a class="code" href="classbdm_1_1fnc.html">00146</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
100<a name="l00147"></a>00147         <span class="keyword">protected</span>:
101<a name="l00149"></a><a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62">00149</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62" title="Length of the output vector.">dimy</a>;
102<a name="l00150"></a>00150         <span class="keyword">public</span>:
103<a name="l00152"></a><a class="code" href="classbdm_1_1fnc.html#d3667c337c88a1c51f672688d8070409">00152</a>                 <a class="code" href="classbdm_1_1fnc.html#d3667c337c88a1c51f672688d8070409" title="default constructor">fnc</a> ( <span class="keywordtype">int</span> dy ) :<a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62" title="Length of the output vector.">dimy</a> ( dy ) {};
104<a name="l00154"></a><a class="code" href="classbdm_1_1fnc.html#6277b11d7fffc7ef8a2fa3e84ae5bad4">00154</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 &amp;cond ) {
105<a name="l00155"></a>00155                         <span class="keywordflow">return</span> vec ( 0 );
106<a name="l00156"></a>00156                 };
107<a name="l00157"></a>00157
108<a name="l00159"></a><a class="code" href="classbdm_1_1fnc.html#0786e40fade2663a70d654c1dda5d73e">00159</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 &amp;val ) {};
109<a name="l00160"></a>00160
110<a name="l00162"></a><a class="code" href="classbdm_1_1fnc.html#a2277a400fc9f4d6c0bf24dc7156183f">00162</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1fnc.html#a2277a400fc9f4d6c0bf24dc7156183f" title="access function">_dimy</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>;}
111<a name="l00163"></a>00163
112<a name="l00165"></a><a class="code" href="classbdm_1_1fnc.html#66455bb4154dcf871964e30a637386d1">00165</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1fnc.html#66455bb4154dcf871964e30a637386d1" title="Destructor for future use;.">~fnc</a>() {};
113<a name="l00166"></a>00166         };
114<a name="l00167"></a>00167
115<a name="l00168"></a>00168         <span class="keyword">class </span>mpdf;
116<a name="l00169"></a>00169
117<a name="l00171"></a>00171
118<a name="l00172"></a><a class="code" href="classbdm_1_1epdf.html">00172</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
119<a name="l00173"></a>00173         <span class="keyword">protected</span>:
120<a name="l00175"></a><a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8">00175</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="Identified of the random variable.">rv</a>;
121<a name="l00176"></a>00176         <span class="keyword">public</span>:
122<a name="l00178"></a><a class="code" href="classbdm_1_1epdf.html#992f3de956fb3439c62225156aafb6f6">00178</a>                 <a class="code" href="classbdm_1_1epdf.html#992f3de956fb3439c62225156aafb6f6" title="default constructor">epdf</a>() :<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a> ( ) {};
123<a name="l00179"></a>00179
124<a name="l00181"></a><a class="code" href="classbdm_1_1epdf.html#fdfed7aa18f2870951888b1dde977415">00181</a>                 <a class="code" href="classbdm_1_1epdf.html#992f3de956fb3439c62225156aafb6f6" title="default constructor">epdf</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv0 ) :<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a> ( rv0 ) {};
125<a name="l00182"></a>00182
126<a name="l00183"></a>00183 <span class="comment">//      //! Returns the required moment of the epdf</span>
127<a name="l00184"></a>00184 <span class="comment">//      virtual vec moment ( const int order = 1 );</span>
128<a name="l00185"></a>00185
129<a name="l00187"></a>00187                 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10" title="Returns a sample,  from density .">sample</a> () <span class="keyword">const</span> =0;
130<a name="l00189"></a>00189                 <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>;
131<a name="l00190"></a>00190
132<a name="l00192"></a>00192                 <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1epdf.html#a8e39e2ff5e5cf5471bc159530d3b2d3" title="Compute log-probability of argument val.">evallog</a> ( <span class="keyword">const</span> vec &amp;val ) <span class="keyword">const</span> =0;
133<a name="l00193"></a>00193
134<a name="l00195"></a><a class="code" href="classbdm_1_1epdf.html#34956d4dd3176eeb5937cf48a1546b62">00195</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 &amp;Val )<span class="keyword"> const </span>{
135<a name="l00196"></a>00196                         vec x ( Val.cols() );
136<a name="l00197"></a>00197                         <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;Val.cols();i++ ) {x ( i ) =<a class="code" href="classbdm_1_1epdf.html#a8e39e2ff5e5cf5471bc159530d3b2d3" title="Compute log-probability of argument val.">evallog</a> ( Val.get_col ( i ) ) ;}
137<a name="l00198"></a>00198                         <span class="keywordflow">return</span> x;
138<a name="l00199"></a>00199                 }
139<a name="l00201"></a><a class="code" href="classbdm_1_1epdf.html#e584eac5579c1b6384947ecf66166c77">00201</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> &amp;<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified 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;}
140<a name="l00203"></a><a class="code" href="classbdm_1_1epdf.html#3fb2ece54f720b62ad325e61214fa0a1">00203</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> &amp;<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified 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;}
141<a name="l00204"></a>00204
142<a name="l00206"></a>00206                 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#565300f603ca00c989e16389d9eb42f5" title="return expected value">mean</a>() <span class="keyword">const</span> =0;
143<a name="l00207"></a>00207
144<a name="l00209"></a>00209                 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1epdf.html#45c2b64faa4bf526799e4154a6ec6097" title="return expected variance (not covariance!)">variance</a>() <span class="keyword">const</span> = 0;
145<a name="l00210"></a>00210
146<a name="l00212"></a><a class="code" href="classbdm_1_1epdf.html#a4fa8637e7307cb051c400c2ce0ca0c8">00212</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1epdf.html#a4fa8637e7307cb051c400c2ce0ca0c8" title="Destructor for future use;.">~epdf</a>() {};
147<a name="l00214"></a><a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21">00214</a>                 <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>&amp; <a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21" title="access function, possibly dangerous!">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>;}
148<a name="l00216"></a><a class="code" href="classbdm_1_1epdf.html#62e88cbce0ce77a8692f5e15d76e805f">00216</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1epdf.html#62e88cbce0ce77a8692f5e15d76e805f" title="modifier function - useful when copying epdfs">_renewrv</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;in_rv ) {<a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8" title="Identified of the random variable.">rv</a>=in_rv;}
149<a name="l00218"></a>00218         };
150<a name="l00219"></a>00219
151<a name="l00220"></a>00220
152<a name="l00222"></a>00222 <span class="comment">//TODO Samplecond can be generalized</span>
153<a name="l00223"></a>00223
154<a name="l00224"></a><a class="code" href="classbdm_1_1mpdf.html">00224</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
155<a name="l00225"></a>00225         <span class="keyword">protected</span>:
156<a name="l00227"></a><a class="code" href="classbdm_1_1mpdf.html#9bcfb45435d30983f436d41c298cbb51">00227</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#9bcfb45435d30983f436d41c298cbb51" title="modeled random variable">rv</a>;
157<a name="l00229"></a><a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288">00229</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>;
158<a name="l00231"></a><a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0">00231</a>                 <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_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>;
159<a name="l00232"></a>00232         <span class="keyword">public</span>:
160<a name="l00233"></a>00233
161<a name="l00235"></a><a class="code" href="classbdm_1_1mpdf.html#e4848a428d8ef0549c6e4a9ed386d9f2">00235</a>                 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#e4848a428d8ef0549c6e4a9ed386d9f2" title="Returns a sample from the density conditioned on cond, .">samplecond</a> ( <span class="keyword">const</span> vec &amp;cond, <span class="keywordtype">double</span> &amp;ll ) {
162<a name="l00236"></a>00236                         this-&gt;<a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );
163<a name="l00237"></a>00237                         vec temp= <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10" title="Returns a sample,  from density .">sample</a>();
164<a name="l00238"></a>00238                         ll=<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#a8e39e2ff5e5cf5471bc159530d3b2d3" title="Compute log-probability of argument val.">evallog</a> ( temp );<span class="keywordflow">return</span> temp;
165<a name="l00239"></a>00239                 };
166<a name="l00241"></a><a class="code" href="classbdm_1_1mpdf.html#ee26963a637b2ea1fb1933652981e652">00241</a>                 <span class="keyword">virtual</span> mat <a class="code" href="classbdm_1_1mpdf.html#ee26963a637b2ea1fb1933652981e652" title="Returns.">samplecond_m</a> ( <span class="keyword">const</span> vec &amp;cond, vec &amp;ll, <span class="keywordtype">int</span> N ) {
167<a name="l00242"></a>00242                         this-&gt;<a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );
168<a name="l00243"></a>00243                         mat temp ( <a class="code" href="classbdm_1_1mpdf.html#9bcfb45435d30983f436d41c298cbb51" title="modeled random variable">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>(),N ); vec smp ( <a class="code" href="classbdm_1_1mpdf.html#9bcfb45435d30983f436d41c298cbb51" title="modeled random variable">rv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>() );
169<a name="l00244"></a>00244                         <span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;N;i++ ) {smp=<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#09721eb148bc49d6ec2f1956039f3a10" title="Returns a sample,  from density .">sample</a>() ;temp.set_col ( i, smp );ll ( i ) =<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#a8e39e2ff5e5cf5471bc159530d3b2d3" title="Compute log-probability of argument val.">evallog</a> ( smp );}
170<a name="l00245"></a>00245                         <span class="keywordflow">return</span> temp;
171<a name="l00246"></a>00246                 };
172<a name="l00248"></a><a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530">00248</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 &amp;cond ) {it_error ( <span class="stringliteral">"Not implemented"</span> );};
173<a name="l00249"></a>00249
174<a name="l00251"></a><a class="code" href="classbdm_1_1mpdf.html#6336a8a72462e2a56a3989a220f18b1b">00251</a>                 <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1mpdf.html#6336a8a72462e2a56a3989a220f18b1b" title="Shortcut for conditioning and evaluation of the internal epdf. In some cases, this...">evallogcond</a> ( <span class="keyword">const</span> vec &amp;dt, <span class="keyword">const</span> vec &amp;cond ) {
175<a name="l00252"></a>00252                         <span class="keywordtype">double</span> tmp; this-&gt;<a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );tmp = <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#a8e39e2ff5e5cf5471bc159530d3b2d3" title="Compute log-probability of argument val.">evallog</a> ( dt );          it_assert_debug ( std::isfinite ( tmp ),<span class="stringliteral">"Infinite value"</span> ); <span class="keywordflow">return</span> tmp;
176<a name="l00253"></a>00253                 };
177<a name="l00254"></a>00254
178<a name="l00256"></a><a class="code" href="classbdm_1_1mpdf.html#0b0ed1ed663071bb7cf4a1349eb94fcb">00256</a>                 <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#0b0ed1ed663071bb7cf4a1349eb94fcb" title="Matrix version of evallogcond.">evallogcond_m</a> ( <span class="keyword">const</span> mat &amp;Dt, <span class="keyword">const</span> vec &amp;cond ) {this-&gt;<a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( cond );<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;<a class="code" href="classbdm_1_1epdf.html#34956d4dd3176eeb5937cf48a1546b62" title="Compute log-probability of multiple values argument val.">evallog_m</a> ( Dt );};
179<a name="l00257"></a>00257
180<a name="l00259"></a><a class="code" href="classbdm_1_1mpdf.html#2d83cf2a5ba40cab11dfe158ced9ee9a">00259</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1mpdf.html#2d83cf2a5ba40cab11dfe158ced9ee9a" title="Destructor for future use;.">~mpdf</a>() {};
181<a name="l00260"></a>00260
182<a name="l00262"></a><a class="code" href="classbdm_1_1mpdf.html#5fd998aa9e223231ab20f29762e594f7">00262</a>                 <a class="code" href="classbdm_1_1mpdf.html#5fd998aa9e223231ab20f29762e594f7" title="Default constructor.">mpdf</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv0, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rvc0 ) :<a class="code" href="classbdm_1_1mpdf.html#9bcfb45435d30983f436d41c298cbb51" title="modeled random variable">rv</a> ( rv0 ),<a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a> ( rvc0 ) {};
183<a name="l00264"></a><a class="code" href="classbdm_1_1mpdf.html#b3aba7311038bf990d706a64cab60cf8">00264</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#b3aba7311038bf990d706a64cab60cf8" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>;}
184<a name="l00266"></a><a class="code" href="classbdm_1_1mpdf.html#222d5280e309c5a053ba73841e98c151">00266</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#222d5280e309c5a053ba73841e98c151" title="access function">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#9bcfb45435d30983f436d41c298cbb51" title="modeled random variable">rv</a>;}
185<a name="l00268"></a><a class="code" href="classbdm_1_1mpdf.html#1892fe3933488942253679f068e9e7f6">00268</a>                 <a class="code" href="classbdm_1_1epdf.html" title="Probability density function with numerical statistics, e.g. posterior density.">epdf</a>&amp; <a class="code" href="classbdm_1_1mpdf.html#1892fe3933488942253679f068e9e7f6" title="access function">_epdf</a>() {<span class="keywordflow">return</span> *<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>;}
186<a name="l00270"></a><a class="code" href="classbdm_1_1mpdf.html#05e843fd11c410a99dad2b88c55aca80">00270</a>                 <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_1mpdf.html#05e843fd11c410a99dad2b88c55aca80" title="access function">_e</a>() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>;}
187<a name="l00271"></a>00271         };
188<a name="l00272"></a>00272
189<a name="l00304"></a><a class="code" href="classbdm_1_1logger.html">00304</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1logger.html" title="DataLink is a connection between two data vectors Up and Down.">logger</a> : <span class="keyword">public</span> <a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
190<a name="l00305"></a>00305         <span class="keyword">protected</span>:
191<a name="l00307"></a><a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f">00307</a>                 Array&lt;RV&gt; <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>;
192<a name="l00309"></a><a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02">00309</a>                 Array&lt;string&gt; <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>;
193<a name="l00310"></a>00310         <span class="keyword">public</span>:
194<a name="l00312"></a><a class="code" href="classbdm_1_1logger.html#f9c2fe7db5ddbda0098a4a5bd2c29fbf">00312</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 ) {}
195<a name="l00313"></a>00313
196<a name="l00315"></a><a class="code" href="classbdm_1_1logger.html#cfb4807e20b3d00f3d33e55ada528dc1">00315</a>                 <span class="keyword">virtual</span> <span class="keywordtype">int</span> <span class="keyword">add</span> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv, <span class="keywordtype">string</span> name=<span class="stringliteral">""</span> ) {
197<a name="l00316"></a>00316                         <span class="keywordtype">int</span> <span class="keywordtype">id</span>=<a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>.length();
198<a name="l00317"></a>00317                         <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>=<a class="code" href="namespacebdm.html#b9016687c0e874ca5cdcf75ae28811aa" title="Concat two random variables.">concat</a> ( <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>, name ); <span class="comment">// diff</span>
199<a name="l00318"></a>00318                         <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> );
200<a name="l00319"></a>00319                         <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a> ( <span class="keywordtype">id</span> ) = rv;
201<a name="l00320"></a>00320                         <span class="keywordflow">return</span> id; <span class="comment">// identifier of the last entry</span>
202<a name="l00321"></a>00321                 }
203<a name="l00322"></a>00322
204<a name="l00324"></a>00324                 <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 &amp;v ) =0;
205<a name="l00325"></a>00325
206<a name="l00327"></a>00327                 <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;
207<a name="l00328"></a>00328
208<a name="l00330"></a><a class="code" href="classbdm_1_1logger.html#0e3fc7d13e23a8f71fc7c8a5ea68559d">00330</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>() {};
209<a name="l00331"></a>00331
210<a name="l00333"></a><a class="code" href="classbdm_1_1logger.html#cd6291a73491f755dc7bc88d022eb2fc">00333</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>() {};
211<a name="l00334"></a>00334
212<a name="l00336"></a><a class="code" href="classbdm_1_1logger.html#f2ade56238a5d408e65a77c53f516448">00336</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1logger.html#f2ade56238a5d408e65a77c53f516448" title="for future use">~logger</a>() {};
213<a name="l00337"></a>00337         };
214<a name="l00338"></a>00338
215<a name="l00339"></a>00339         <span class="keyword">class </span>datalink_e2e {
216<a name="l00340"></a>00340         <span class="keyword">protected</span>:
217<a name="l00342"></a>00342                 <span class="keywordtype">int</span> valsize;
218<a name="l00344"></a>00344                 <span class="keywordtype">int</span> valupsize;
219<a name="l00346"></a>00346                 ivec v2v_up;
220<a name="l00347"></a>00347         <span class="keyword">public</span>:
221<a name="l00349"></a>00349                 datalink_e2e ( <span class="keyword">const</span> RV &amp;rv, <span class="keyword">const</span> RV &amp;rv_up ) :
222<a name="l00350"></a>00350                                 valsize ( rv.count() ), valupsize ( rv_up.count() ), v2v_up ( rv.dataind ( rv_up ) )  {
223<a name="l00351"></a>00351                         it_assert_debug ( v2v_up.length() ==valsize,<span class="stringliteral">"rv is not fully in rv_up"</span> );
224<a name="l00352"></a>00352                 }
225<a name="l00354"></a>00354                 vec get_val ( <span class="keyword">const</span> vec &amp;val_up ) {it_assert_debug ( valupsize==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> ); <span class="keywordflow">return</span> get_vec ( val_up,v2v_up );}
226<a name="l00356"></a>00356                 <span class="keywordtype">void</span> fill_val ( vec &amp;val_up, <span class="keyword">const</span> vec &amp;val ) {
227<a name="l00357"></a>00357                         it_assert_debug ( valsize==val.length(),<span class="stringliteral">"Wrong val"</span> );
228<a name="l00358"></a>00358                         it_assert_debug ( valupsize==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> );
229<a name="l00359"></a>00359                         set_subvector ( val_up, v2v_up, val );
230<a name="l00360"></a>00360                 }
231<a name="l00361"></a>00361         };
232<a name="l00362"></a>00362
233<a name="l00364"></a><a class="code" href="classbdm_1_1datalink__m2e.html">00364</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1datalink__m2e.html" title="data link between">datalink_m2e</a>: <span class="keyword">public</span> datalink_e2e {
234<a name="l00365"></a>00365         <span class="keyword">protected</span>:
235<a name="l00367"></a><a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9">00367</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>;
236<a name="l00369"></a><a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783">00369</a>                 ivec <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a>;
237<a name="l00371"></a><a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc">00371</a>                 ivec <a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, ideces of the local cond">v2c_lo</a>;
238<a name="l00372"></a>00372
239<a name="l00373"></a>00373         <span class="keyword">public</span>:
240<a name="l00375"></a><a class="code" href="classbdm_1_1datalink__m2e.html#31bb5b6cd04b4d99a02fc8d1ad5ce202">00375</a>                 <a class="code" href="classbdm_1_1datalink__m2e.html#31bb5b6cd04b4d99a02fc8d1ad5ce202" title="Constructor.">datalink_m2e</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv,  <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rvc, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv_up ) :
241<a name="l00376"></a>00376                                 datalink_e2e ( rv,rv_up ), <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> ( rvc.count() ) {
242<a name="l00377"></a>00377                         <span class="comment">//establish v2c connection</span>
243<a name="l00378"></a>00378                         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, ideces of the local cond">v2c_lo</a>, <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a> );
244<a name="l00379"></a>00379                 }
245<a name="l00381"></a><a class="code" href="classbdm_1_1datalink__m2e.html#7ba850822e3cf9a4b079afc7e7c01703">00381</a>                 vec <a class="code" href="classbdm_1_1datalink__m2e.html#7ba850822e3cf9a4b079afc7e7c01703" title="Construct condition.">get_cond</a> ( <span class="keyword">const</span> vec &amp;val_up ) {
246<a name="l00382"></a>00382                         vec tmp ( <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> );
247<a name="l00383"></a>00383                         set_subvector ( tmp,<a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, ideces 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, indeces of the upper val">v2c_up</a> ) );
248<a name="l00384"></a>00384                         <span class="keywordflow">return</span> tmp;
249<a name="l00385"></a>00385                 }
250<a name="l00386"></a>00386                 <span class="keywordtype">void</span> fill_val_cond ( vec &amp;val_up, <span class="keyword">const</span> vec &amp;val, <span class="keyword">const</span> vec &amp;cond ) {
251<a name="l00387"></a>00387                         it_assert_debug ( valsize==val.length(),<span class="stringliteral">"Wrong val"</span> );
252<a name="l00388"></a>00388                         it_assert_debug ( valupsize==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> );
253<a name="l00389"></a>00389                         set_subvector ( val_up, v2v_up, val );
254<a name="l00390"></a>00390                         set_subvector ( val_up, <a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783" title="upper_val-to-local_cond link, indeces of the upper val">v2c_up</a>, cond );
255<a name="l00391"></a>00391                 }
256<a name="l00392"></a>00392         };
257<a name="l00395"></a><a class="code" href="classbdm_1_1datalink__m2m.html">00395</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 between">datalink_m2e</a> {
258<a name="l00396"></a>00396         <span class="keyword">protected</span>:
259<a name="l00398"></a><a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d">00398</a>                 ivec <a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indeces of the upper cond">c2c_up</a>;
260<a name="l00400"></a><a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5">00400</a>                 ivec <a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>;
261<a name="l00401"></a>00401         <span class="keyword">public</span>:
262<a name="l00403"></a><a class="code" href="classbdm_1_1datalink__m2m.html#540c4b60686fdf1a7915300ae11bc41b">00403</a>                 <a class="code" href="classbdm_1_1datalink__m2m.html#540c4b60686fdf1a7915300ae11bc41b" title="Constructor.">datalink_m2m</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rvc, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;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> &amp;rvc_up ) :
263<a name="l00404"></a>00404                                 <a class="code" href="classbdm_1_1datalink__m2e.html" title="data link between">datalink_m2e</a> ( rv, rvc, rv_up ) {
264<a name="l00405"></a>00405                         <span class="comment">//establish c2c connection</span>
265<a name="l00406"></a>00406                         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, indeces of the local cond">c2c_lo</a>, <a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indeces of the upper cond">c2c_up</a> );
266<a name="l00407"></a>00407                         it_assert_debug ( <a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indeces 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, ideces 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> );
267<a name="l00408"></a>00408                 }
268<a name="l00410"></a><a class="code" href="classbdm_1_1datalink__m2m.html#911fb7fcfa80d4e31db8e0cc03141a90">00410</a>                 vec <a class="code" href="classbdm_1_1datalink__m2m.html#911fb7fcfa80d4e31db8e0cc03141a90" title="Get cond for myself from val and cond of &amp;quot;Up&amp;quot;.">get_cond</a> ( <span class="keyword">const</span> vec &amp;val_up, <span class="keyword">const</span> vec &amp;cond_up ) {
269<a name="l00411"></a>00411                         vec tmp ( <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> );
270<a name="l00412"></a>00412                         set_subvector ( tmp,<a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc" title="upper_val-to-local_cond link, ideces 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, indeces of the upper val">v2c_up</a> ) );
271<a name="l00413"></a>00413                         set_subvector ( tmp,<a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5" title="cond-to-cond link, indeces of the local cond">c2c_lo</a>,cond_up ( <a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d" title="cond-to-cond link, indeces of the upper cond">c2c_up</a> ) );
272<a name="l00414"></a>00414                         <span class="keywordflow">return</span> tmp;
273<a name="l00415"></a>00415                 }
274<a name="l00417"></a>00417
275<a name="l00418"></a>00418         };
276<a name="l00419"></a>00419
277<a name="l00423"></a><a class="code" href="classbdm_1_1mepdf.html">00423</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> {
278<a name="l00424"></a>00424         <span class="keyword">public</span>:
279<a name="l00426"></a><a class="code" href="classbdm_1_1mepdf.html#27cc32d1a810516c25bda51ddf886376">00426</a>                 <a class="code" href="classbdm_1_1mepdf.html#27cc32d1a810516c25bda51ddf886376" title="Default constructor.">mepdf</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>* em ) :<a class="code" href="classbdm_1_1mpdf.html" title="Conditional probability density, e.g. modeling some dependencies.">mpdf</a> ( em-&gt;<a class="code" href="classbdm_1_1mpdf.html#222d5280e309c5a053ba73841e98c151" title="access function">_rv</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#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>=<span class="keyword">const_cast&lt;</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">&gt;</span> ( em );};
280<a name="l00427"></a><a class="code" href="classbdm_1_1mepdf.html#b2b898058169e7ed9a6dc7cd4e1c147c">00427</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1mepdf.html#b2b898058169e7ed9a6dc7cd4e1c147c" title="Update ep so that it represents this mpdf conditioned on rvc = cond.">condition</a> ( <span class="keyword">const</span> vec &amp;cond ) {}
281<a name="l00428"></a>00428         };
282<a name="l00429"></a>00429
283<a name="l00432"></a><a class="code" href="classbdm_1_1compositepdf.html">00432</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1compositepdf.html" title="Abstract composition of pdfs, will be used for specific classes this abstract class...">compositepdf</a> {
284<a name="l00433"></a>00433         <span class="keyword">protected</span>:
285<a name="l00435"></a><a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a">00435</a>                 <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a>;
286<a name="l00437"></a><a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a">00437</a>                 Array&lt;mpdf*&gt; <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>;
287<a name="l00438"></a>00438         <span class="keyword">public</span>:
288<a name="l00439"></a>00439                 <a class="code" href="classbdm_1_1compositepdf.html" title="Abstract composition of pdfs, will be used for specific classes this abstract class...">compositepdf</a> ( Array&lt;mpdf*&gt; A0 ) : <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ( A0.length() ), <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( A0 ) {};
289<a name="l00441"></a>00441                 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> <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> );
290<a name="l00443"></a>00443                 <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> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv, <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rvc );
291<a name="l00444"></a>00444         };
292<a name="l00445"></a>00445
293<a name="l00453"></a><a class="code" href="classbdm_1_1DS.html">00453</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
294<a name="l00454"></a>00454         <span class="keyword">protected</span>:
295<a name="l00456"></a><a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec">00456</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="Observed variables, returned by getdata().">Drv</a>;
296<a name="l00458"></a><a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702">00458</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="Action variables, accepted by write().">Urv</a>; <span class="comment">//</span>
297<a name="l00460"></a><a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867">00460</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;
298<a name="l00461"></a>00461         <span class="keyword">public</span>:
299<a name="l00462"></a>00462                 <a class="code" href="classbdm_1_1DS.html" title="Abstract class for discrete-time sources of data.">DS</a>() :<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Observed variables, returned by getdata().">Drv</a> ( <a class="code" href="namespacebdm.html#bdd4760f14ec7d6fea083aa2170f6236" title="Default empty RV that can be used as default argument.">RV0</a> ),<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Action variables, accepted by write().">Urv</a> ( <a class="code" href="namespacebdm.html#bdd4760f14ec7d6fea083aa2170f6236" title="Default empty RV that can be used as default argument.">RV0</a> ) {};
300<a name="l00463"></a>00463                 <a class="code" href="classbdm_1_1DS.html" title="Abstract class for discrete-time sources of data.">DS</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;Drv0, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;Urv0 ) :<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Observed variables, returned by getdata().">Drv</a> ( Drv0 ),<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Action variables, accepted by write().">Urv</a> ( Urv0 ) {};
301<a name="l00465"></a><a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba">00465</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 &amp;dt ) {it_error ( <span class="stringliteral">"abstract class"</span> );};
302<a name="l00467"></a><a class="code" href="classbdm_1_1DS.html#a0bf451c7fea73ff4150ff9e83da4832">00467</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 &amp;dt, <span class="keyword">const</span> ivec &amp;indeces ) {it_error ( <span class="stringliteral">"abstract class"</span> );};
303<a name="l00469"></a><a class="code" href="classbdm_1_1DS.html#f5baeee8528e5e4c5c7a344b2e3f192a">00469</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 &amp;ut ) {it_error ( <span class="stringliteral">"abstract class"</span> );};
304<a name="l00471"></a><a class="code" href="classbdm_1_1DS.html#6bced13f84f6bf04449728817e049e1d">00471</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 &amp;ut, <span class="keyword">const</span> ivec &amp;indeces ) {it_error ( <span class="stringliteral">"abstract class"</span> );};
305<a name="l00472"></a>00472
306<a name="l00474"></a>00474                 <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;
307<a name="l00475"></a>00475
308<a name="l00477"></a><a class="code" href="classbdm_1_1DS.html#3c9e417d85052396fe74c164857ee67e">00477</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="DataLink is a connection between two data vectors Up and Down.">logger</a> &amp;L ) {
309<a name="l00478"></a>00478                         <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#cfb4807e20b3d00f3d33e55ada528dc1" title="returns an identifier which will be later needed for calling the log() function">add</a> ( <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Observed variables, returned by getdata().">Drv</a>,<span class="stringliteral">""</span> );
310<a name="l00479"></a>00479                         L_ut=L.<a class="code" href="classbdm_1_1logger.html#cfb4807e20b3d00f3d33e55ada528dc1" title="returns an identifier which will be later needed for calling the log() function">add</a> ( <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Action variables, accepted by write().">Urv</a>,<span class="stringliteral">""</span> );
311<a name="l00480"></a>00480                 }
312<a name="l00482"></a><a class="code" href="classbdm_1_1DS.html#e03d2c93cddd5c3a03107ddd3ae999cf">00482</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="DataLink is a connection between two data vectors Up and Down.">logger</a> &amp;L ) {
313<a name="l00483"></a>00483                         vec tmp(<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Observed variables, returned by getdata().">Drv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>()+<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Action variables, accepted by write().">Urv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>());
314<a name="l00484"></a>00484                         <a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba" title="Returns full vector of observed data=[output, input].">getdata</a>(tmp);
315<a name="l00485"></a>00485                         <span class="comment">// d is first in getdata</span>
316<a name="l00486"></a>00486                         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="Observed variables, returned by getdata().">Drv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>() ) );
317<a name="l00487"></a>00487                         <span class="comment">// u follows after d in getdata</span>
318<a name="l00488"></a>00488                         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="Observed variables, returned by getdata().">Drv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>(), <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Action variables, accepted by write().">Urv</a>.<a class="code" href="classbdm_1_1RV.html#2174751a00ce19f941edd2c1a861be67" title="Return number of scalars in the RV.">count</a>() ) );
319<a name="l00489"></a>00489                 }
320<a name="l00490"></a>00490         };
321<a name="l00491"></a>00491
322<a name="l00496"></a><a class="code" href="classbdm_1_1BM.html">00496</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a> :<span class="keyword">public</span> <a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
323<a name="l00497"></a>00497         <span class="keyword">protected</span>:
324<a name="l00499"></a><a class="code" href="classbdm_1_1BM.html#18d6db4af8ee42077741d9e3618153ca">00499</a>                 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> rv;
325<a name="l00501"></a><a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a">00501</a>                 <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>;
326<a name="l00503"></a><a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee">00503</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>;
327<a name="l00504"></a>00504         <span class="keyword">public</span>:
328<a name="l00505"></a>00505
329<a name="l00507"></a><a class="code" href="classbdm_1_1BM.html#a33000bb124c32e58742f056ecab4293">00507</a>                 <a class="code" href="classbdm_1_1BM.html#a33000bb124c32e58742f056ecab4293" title="Default constructor.">BM</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv0, <span class="keywordtype">double</span> ll0=0,<span class="keywordtype">bool</span> evalll0=<span class="keyword">true</span> ) :rv ( rv0 ), <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a> ( ll0 ),<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> ( evalll0 ) {<span class="comment">//Fixme: test rv</span>
330<a name="l00508"></a>00508                 };
331<a name="l00510"></a><a class="code" href="classbdm_1_1BM.html#241b3701190ff1f729fe873a2eef0055">00510</a>                 <a class="code" href="classbdm_1_1BM.html#a33000bb124c32e58742f056ecab4293" title="Default constructor.">BM</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a> &amp;B ) : rv ( B.rv ), <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> ) {}
332<a name="l00511"></a>00511
333<a name="l00515"></a>00515                 <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 &amp;dt ) = 0;
334<a name="l00517"></a>00517                 <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 &amp;Dt );
335<a name="l00519"></a>00519                 <span class="keyword">virtual</span> <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>&amp; <a class="code" href="classbdm_1_1BM.html#963258c4c2dd05be001003b19aceefef" title="Returns a reference to the epdf representing posterior density on parameters.">_epdf</a>() <span class="keyword">const</span> =0;
336<a name="l00520"></a>00520
337<a name="l00522"></a>00522                 <span class="keyword">virtual</span> <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>* <a class="code" href="classbdm_1_1BM.html#4ed0f8b880e606316ae800f3a011c3a6" title="Returns a pointer to the epdf representing posterior density on parameters. Use with...">_e</a>() <span class="keyword">const</span> =0;
338<a name="l00523"></a>00523
339<a name="l00526"></a><a class="code" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">00526</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 &amp;dt )<span class="keyword"> const</span>{it_error ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> 0.0;}
340<a name="l00528"></a><a class="code" href="classbdm_1_1BM.html#0e8ebe61fb14990abe1254bd3dda5fae">00528</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 &amp;dt )<span class="keyword"> const</span>{vec tmp ( dt.cols() );<span class="keywordflow">for</span> ( <span class="keywordtype">int</span> i=0;i&lt;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;}
341<a name="l00529"></a>00529
342<a name="l00531"></a><a class="code" href="classbdm_1_1BM.html#710e7d69c0d8791fb41a7cd4683cca2c">00531</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#710e7d69c0d8791fb41a7cd4683cca2c" title="Constructs a predictive density (marginal density on data).">predictor</a> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv )<span class="keyword"> const </span>{it_error ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> NULL;};
343<a name="l00532"></a>00532
344<a name="l00534"></a><a class="code" href="classbdm_1_1BM.html#374281cad3caa16fa30609bb4ac16505">00534</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1BM.html#374281cad3caa16fa30609bb4ac16505" title="Destructor for future use;.">~BM</a>() {};
345<a name="l00536"></a><a class="code" href="classbdm_1_1BM.html#40a3c891996391e3135518053a917793">00536</a>                 <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>&amp; <a class="code" href="classbdm_1_1BM.html#40a3c891996391e3135518053a917793" title="access function">_rv</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> rv;}
346<a name="l00538"></a><a class="code" href="classbdm_1_1BM.html#5be65d37dedfe33a3671e7065f523a70">00538</a>                 <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1BM.html#5be65d37dedfe33a3671e7065f523a70" title="access function">_ll</a>()<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>;}
347<a name="l00540"></a><a class="code" href="classbdm_1_1BM.html#236b3abbcc93594fc97cd86d82c1a83f">00540</a>                 <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BM.html#236b3abbcc93594fc97cd86d82c1a83f" title="access function">set_evalll</a> ( <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;}
348<a name="l00541"></a>00541
349<a name="l00544"></a><a class="code" href="classbdm_1_1BM.html#3efb3098172f1f67564a312fe732473e">00544</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1BM.html" title="Bayesian Model of the world, i.e. all uncertainty is modeled by probabilities.">BM</a>* <a class="code" href="classbdm_1_1BM.html#3efb3098172f1f67564a312fe732473e">_copy_</a> ( <span class="keywordtype">bool</span> changerv=<span class="keyword">false</span> ) {it_error ( <span class="stringliteral">"function _copy_ not implemented for this BM"</span> ); <span class="keywordflow">return</span> NULL;};
350<a name="l00545"></a>00545         };
351<a name="l00546"></a>00546
352<a name="l00556"></a><a class="code" href="classbdm_1_1BMcond.html">00556</a>         <span class="keyword">class </span><a class="code" href="classbdm_1_1BMcond.html" title="Conditional Bayesian Filter.">BMcond</a> :<span class="keyword">public</span> <a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
353<a name="l00557"></a>00557         <span class="keyword">protected</span>:
354<a name="l00559"></a><a class="code" href="classbdm_1_1BMcond.html#9a12750776d977408aada06a70093297">00559</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_1BMcond.html#9a12750776d977408aada06a70093297" title="Identificator of the conditioning variable.">rvc</a>;
355<a name="l00560"></a>00560         <span class="keyword">public</span>:
356<a name="l00562"></a>00562                 <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1BMcond.html#a67bac171f89ef18d0191f8f74189f25" title="Substitute val for rvc.">condition</a> ( <span class="keyword">const</span> vec &amp;val ) =0;
357<a name="l00564"></a><a class="code" href="classbdm_1_1BMcond.html#dde99f401cbf2e8e6a926eaa764b457d">00564</a>                 <a class="code" href="classbdm_1_1BMcond.html#dde99f401cbf2e8e6a926eaa764b457d" title="Default constructor.">BMcond</a> ( <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;rv0 ) :<a class="code" href="classbdm_1_1BMcond.html#9a12750776d977408aada06a70093297" title="Identificator of the conditioning variable.">rvc</a> ( rv0 ) {};
358<a name="l00566"></a><a class="code" href="classbdm_1_1BMcond.html#7d68d46cfcc36f9cc07a15a5403b5a1a">00566</a>                 <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1BMcond.html#7d68d46cfcc36f9cc07a15a5403b5a1a" title="Destructor for future use.">~BMcond</a>() {};
359<a name="l00568"></a><a class="code" href="classbdm_1_1BMcond.html#7506910f93250b44fea505ec4ffb19dc">00568</a>                 <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>&amp; <a class="code" href="classbdm_1_1BMcond.html#7506910f93250b44fea505ec4ffb19dc" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1BMcond.html#9a12750776d977408aada06a70093297" title="Identificator of the conditioning variable.">rvc</a>;}
360<a name="l00569"></a>00569         };
361<a name="l00570"></a>00570
362<a name="l00571"></a>00571 }; <span class="comment">//namespace</span>
363<a name="l00573"></a>00573 <span class="comment"></span><span class="preprocessor">#endif // BM_H</span>
364</pre></div></div>
365<hr size="1"><address style="text-align: right;"><small>Generated on Fri Feb 6 19:49:45 2009 for mixpp by&nbsp;
366<a href="http://www.doxygen.org/index.html">
367<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.6 </small></address>
368</body>
369</html>
Note: See TracBrowser for help on using the browser.