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

Revision 353, 111.0 kB (checked in by smidl, 16 years ago)

doc

RevLine 
[3]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">
[261]3<title>mixpp: libBM.h Source File</title>
[290]4<link href="tabs.css" rel="stylesheet" type="text/css">
[3]5<link href="doxygen.css" rel="stylesheet" type="text/css">
6</head><body>
[290]7<!-- Generated by Doxygen 1.5.8 -->
[271]8<script type="text/javascript">
9<!--
10function 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}
22function 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}
46window.onload = initDynSections;
47-->
48</script>
[91]49<div class="navigation" id="top">
50  <div class="tabs">
51    <ul>
[271]52      <li><a href="main.html"><span>Main&nbsp;Page</span></a></li>
[91]53      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
[219]54      <li><a href="modules.html"><span>Modules</span></a></li>
[290]55      <li><a href="annotated.html"><span>Classes</span></a></li>
[91]56      <li class="current"><a href="files.html"><span>Files</span></a></li>
57    </ul>
58  </div>
[290]59  <div class="tabs">
60    <ul>
61      <li><a href="files.html"><span>File&nbsp;List</span></a></li>
62      <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
63    </ul>
64  </div>
[261]65<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
[271]66<a name="l00013"></a>00013 <span class="preprocessor">#ifndef BM_H</span>
67<a name="l00014"></a>00014 <span class="preprocessor"></span><span class="preprocessor">#define BM_H</span>
68<a name="l00015"></a>00015 <span class="preprocessor"></span>
69<a name="l00016"></a>00016
70<a name="l00017"></a>00017 <span class="preprocessor">#include "../itpp_ext.h"</span>
[353]71<a name="l00018"></a>00018 <span class="preprocessor">#include "../libconfig/libconfig.h++"</span>
72<a name="l00019"></a>00019 <span class="preprocessor">#include &lt;map&gt;</span>
73<a name="l00020"></a>00020
74<a name="l00021"></a>00021 <span class="keyword">using namespace </span>libconfig;
75<a name="l00022"></a>00022
76<a name="l00023"></a>00023 <span class="keyword">using namespace </span>itpp;
77<a name="l00024"></a>00024 <span class="keyword">using namespace </span>std;
78<a name="l00025"></a>00025
79<a name="l00026"></a>00026 <span class="keyword">namespace </span>bdm {
80<a name="l00027"></a>00027
81<a name="l00029"></a>00029
82<a name="l00030"></a><a class="code" href="classbdm_1_1bdmroot.html">00030</a> <span class="keyword">class </span><a class="code" href="classbdm_1_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
83<a name="l00031"></a>00031 <span class="keyword">public</span>:
84<a name="l00033"></a><a class="code" href="classbdm_1_1bdmroot.html#9c191bccc52a133d32d84edbc1a19c49">00033</a>         <span class="keyword">virtual</span> <a class="code" href="classbdm_1_1bdmroot.html#9c191bccc52a133d32d84edbc1a19c49" title="make sure this is a virtual object">~bdmroot</a>()
85<a name="l00034"></a>00034         {
86<a name="l00035"></a>00035         }
87<a name="l00036"></a>00036
88<a name="l00038"></a><a class="code" href="classbdm_1_1bdmroot.html#3d95e9de79ca1d52b06f1b8f80660f09">00038</a>         <span class="keyword">virtual</span> <span class="keywordtype">string</span> <a class="code" href="classbdm_1_1bdmroot.html#3d95e9de79ca1d52b06f1b8f80660f09" title="This method returns a basic info about the current instance.">ToString</a>()
89<a name="l00039"></a>00039         {
90<a name="l00040"></a>00040                 <span class="keywordflow">return</span> <span class="stringliteral">""</span>;
91<a name="l00041"></a>00041         }
92<a name="l00042"></a>00042
93<a name="l00044"></a><a class="code" href="classbdm_1_1bdmroot.html#d4602efd3ad9fb9e4532426e4f488df8">00044</a>         <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1bdmroot.html#d4602efd3ad9fb9e4532426e4f488df8" title="This method arrange instance properties according the data stored in the Setting...">from_setting</a>( <span class="keyword">const</span> Setting &amp;root )
94<a name="l00045"></a>00045         {
95<a name="l00046"></a>00046         }
[271]96<a name="l00047"></a>00047
[353]97<a name="l00049"></a><a class="code" href="classbdm_1_1bdmroot.html#a80210bbe266bf26f100f342460a313f">00049</a>         <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1bdmroot.html#a80210bbe266bf26f100f342460a313f" title="This method save all the instance properties into the Setting structure.">to_setting</a>( Setting &amp;root )
98<a name="l00050"></a>00050         {       
99<a name="l00051"></a>00051         }
100<a name="l00052"></a>00052 };
101<a name="l00053"></a>00053
102<a name="l00054"></a>00054 <span class="keyword">typedef</span> std::map&lt;string, int&gt; RVmap;
103<a name="l00055"></a>00055 <span class="keyword">extern</span> ivec RV_SIZES;
104<a name="l00056"></a>00056 <span class="keyword">extern</span> <a class="code" href="classitpp_1_1Array.html">Array&lt;string&gt;</a> RV_NAMES;
105<a name="l00057"></a>00057
106<a name="l00059"></a><a class="code" href="classbdm_1_1str.html">00059</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> {
107<a name="l00060"></a>00060 <span class="keyword">public</span>:
108<a name="l00062"></a><a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0">00062</a>         ivec <a class="code" href="classbdm_1_1str.html#64675a73fc7324177175f19b3baeb0b0" title="vector id ids (non-unique!)">ids</a>;
109<a name="l00064"></a><a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a">00064</a>         ivec <a class="code" href="classbdm_1_1str.html#905dc0f240e40432c103708c4a34f80a" title="vector of times">times</a>;
110<a name="l00066"></a><a class="code" href="classbdm_1_1str.html#7c6b9dccc6814d7bd1e8e2208a0d3f3f">00066</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 ) {
111<a name="l00067"></a>00067                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( times0.length() ==ids0.length(),<span class="stringliteral">"Incompatible input"</span> );
112<a name="l00068"></a>00068         };
113<a name="l00069"></a>00069 };
114<a name="l00070"></a>00070
115<a name="l00109"></a><a class="code" href="classbdm_1_1RV.html">00109</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> {
116<a name="l00110"></a>00110 <span class="keyword">protected</span>:
117<a name="l00112"></a><a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb">00112</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#0cfe477965243a89eb87a06b024e2fcb" title="size of the data vector">dsize</a>;
118<a name="l00114"></a><a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385">00114</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1RV.html#b11bfdceb831ed30d33b852dd4a6b385" title="number of individual rvs">len</a>;
119<a name="l00116"></a><a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b">00116</a>         ivec <a class="code" href="classbdm_1_1RV.html#0a9a4769ab5a322495bea0cd8091500b" title="Vector of unique IDs.">ids</a>;
120<a name="l00118"></a><a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896">00118</a>         ivec <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a>;
[271]121<a name="l00119"></a>00119
[353]122<a name="l00120"></a>00120 <span class="keyword">private</span>:
123<a name="l00122"></a>00122         <span class="keywordtype">void</span> init ( <a class="code" href="classitpp_1_1Array.html">Array&lt;std::string&gt;</a> in_names, ivec in_sizes, ivec in_times );
124<a name="l00123"></a>00123         <span class="keywordtype">int</span> init ( <span class="keyword">const</span>  <span class="keywordtype">string</span> &amp;name, <span class="keywordtype">int</span> size );
125<a name="l00124"></a>00124 <span class="keyword">public</span>:
126<a name="l00127"></a>00127
127<a name="l00129"></a><a class="code" href="classbdm_1_1RV.html#3d93da01a64dab31074b25b7d93e5eae">00129</a>         <a class="code" href="classbdm_1_1RV.html#3d93da01a64dab31074b25b7d93e5eae" title="Full constructor.">RV</a> ( <a class="code" href="classitpp_1_1Array.html">Array&lt;std::string&gt;</a> in_names, ivec in_sizes, ivec in_times ) {init ( in_names,in_sizes,in_times );};
128<a name="l00131"></a><a class="code" href="classbdm_1_1RV.html#5bd1bd9ec51b91578ecc83a6b2398a20">00131</a>         <a class="code" href="classbdm_1_1RV.html#5bd1bd9ec51b91578ecc83a6b2398a20" title="Constructor with times=0.">RV</a> ( <a class="code" href="classitpp_1_1Array.html">Array&lt;std::string&gt;</a> in_names, ivec in_sizes ) {init ( in_names,in_sizes,<a class="code" href="group__specmat.html#g9df9aa3553b6a5c07924954b85466cec" title="A Int (rows,cols)-matrix of zeros.">zeros_i</a> ( in_names.<a class="code" href="classitpp_1_1Array.html#53629018f146084de506c33c69fa19a3" title="Returns the number of data elements in the array object.">length</a>() ) );};
129<a name="l00133"></a><a class="code" href="classbdm_1_1RV.html#a222dfef569e8f123180d40fde97c1cd">00133</a>         <a class="code" href="classbdm_1_1RV.html#a222dfef569e8f123180d40fde97c1cd" title="Constructor with sizes=1, times=0.">RV</a> ( <a class="code" href="classitpp_1_1Array.html">Array&lt;std::string&gt;</a> in_names ) {init ( in_names,<a class="code" href="group__specmat.html#g7c2c02190b4aa222cc692641cc3a73c7" title="A Int (rows,cols)-matrix of ones.">ones_i</a> ( in_names.<a class="code" href="classitpp_1_1Array.html#53629018f146084de506c33c69fa19a3" title="Returns the number of data elements in the array object.">length</a>() ),<a class="code" href="group__specmat.html#g9df9aa3553b6a5c07924954b85466cec" title="A Int (rows,cols)-matrix of zeros.">zeros_i</a> ( in_names.<a class="code" href="classitpp_1_1Array.html#53629018f146084de506c33c69fa19a3" title="Returns the number of data elements in the array object.">length</a>() ) );}
130<a name="l00135"></a><a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412">00135</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 ) {};
131<a name="l00137"></a>00137         <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 );
132<a name="l00139"></a>00139
133<a name="l00142"></a>00142
134<a name="l00144"></a>00144         <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 );
135<a name="l00145"></a>00145         <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>;} ;
136<a name="l00147"></a>00147         <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>;
137<a name="l00148"></a>00148         ivec cumsizes() <span class="keyword">const</span>;
138<a name="l00149"></a>00149         <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>;} ;
139<a name="l00150"></a>00150         <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 );};
140<a name="l00151"></a>00151         <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 ) );};
141<a name="l00152"></a>00152         <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 );};
142<a name="l00153"></a>00153         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 ) );};
143<a name="l00154"></a>00154         <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;};
144<a name="l00156"></a>00156
145<a name="l00157"></a>00157         <span class="comment">//TODO why not inline and later??</span>
146<a name="l00158"></a>00158
147<a name="l00161"></a>00161
148<a name="l00163"></a>00163         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> &amp;rv2 ) <span class="keyword">const</span>;
149<a name="l00165"></a>00165         <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> &amp;rv2 ) <span class="keyword">const</span>;
150<a name="l00167"></a>00167         <span class="keywordtype">bool</span> <span class="keyword">add</span> ( <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html#e319ca900261a5005a9b4b353ec16412" title="Constructor of empty RV.">RV</a> &amp;rv2 );
151<a name="l00169"></a>00169         <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> &amp;rv2 ) <span class="keyword">const</span>;
152<a name="l00171"></a>00171         <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 indeces ind.">subselect</a> ( <span class="keyword">const</span> ivec &amp;ind ) <span class="keyword">const</span>;
153<a name="l00173"></a><a class="code" href="classbdm_1_1RV.html#6f3daf0f857fff53a1c0a445fe160966">00173</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 indeces ind.">operator() </a>( <span class="keyword">const</span> ivec &amp;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 indeces ind.">subselect</a> ( ind );};
154<a name="l00175"></a><a class="code" href="classbdm_1_1RV.html#452fb5f7ab1eafe17932bf9f458d0415">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_1RV.html#6f3daf0f857fff53a1c0a445fe160966" title="Select only variables at indeces ind.">operator() </a>( <span class="keywordtype">int</span> di1, <span class="keywordtype">int</span> di2 )<span class="keyword"> const </span>{
155<a name="l00176"></a>00176                 ivec sz=cumsizes();
156<a name="l00177"></a>00177                 <span class="keywordtype">int</span> i1=0;
157<a name="l00178"></a>00178                 <span class="keywordflow">while</span> ( sz ( i1 ) &lt;di1 ) i1++;
158<a name="l00179"></a>00179                 <span class="keywordtype">int</span> i2=i1;
159<a name="l00180"></a>00180                 <span class="keywordflow">while</span> ( sz ( i2 ) &lt;di2 ) i2++;
160<a name="l00181"></a>00181                 <span class="keywordflow">return</span> <a class="code" href="classbdm_1_1RV.html#c0c04ab99ded62e1e05bcefb151ebeb9" title="Select only variables at indeces ind.">subselect</a> ( <a class="code" href="group__specmat.html#ge5632baa525e231b015f3d45f3503cac" title="Linspace (works in the same way as the matlab version).">linspace</a> ( i1,i2 ) );
161<a name="l00182"></a>00182         };
162<a name="l00184"></a>00184         <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 );
163<a name="l00186"></a>00186
[271]164<a name="l00189"></a>00189
[353]165<a name="l00191"></a>00191         <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>;
166<a name="l00194"></a>00194         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>;
167<a name="l00197"></a>00197         <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>;
168<a name="l00199"></a><a class="code" href="classbdm_1_1RV.html#9a5c5aeeaf39299d07763377b05a8241">00199</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> <a class="code" href="group__protocol.html#g27dde5be5aac902c0de838b4b872cd2f" title="ADD DOCUMENTATION HERE.">min</a> ( <a class="code" href="classbdm_1_1RV.html#ef6703854820bbf61736738f1bac6896" title="Vector of shifts from current time.">times</a> );};
[271]169<a name="l00201"></a>00201
[353]170<a name="l00202"></a>00202 };
171<a name="l00203"></a>00203
[270]172<a name="l00204"></a>00204
[353]173<a name="l00206"></a>00206 RV concat ( <span class="keyword">const</span> RV &amp;rv1, <span class="keyword">const</span> RV &amp;rv2 );
174<a name="l00207"></a>00207
175<a name="l00209"></a>00209 <span class="keyword">extern</span> RV RV0;
[271]176<a name="l00210"></a>00210
177<a name="l00212"></a>00212
[353]178<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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
[271]179<a name="l00214"></a>00214 <span class="keyword">protected</span>:
[353]180<a name="l00216"></a><a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62">00216</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1fnc.html#52156cb4a52a62d51fc7455985797a62" title="Length of the output vector.">dimy</a>;
181<a name="l00217"></a>00217 <span class="keyword">public</span>:
182<a name="l00219"></a><a class="code" href="classbdm_1_1fnc.html#65e9fb11f6cc7884018d8f229bcabbb7">00219</a>         <a class="code" href="classbdm_1_1fnc.html#65e9fb11f6cc7884018d8f229bcabbb7" title="default constructor">fnc</a> ( ) {};
183<a name="l00221"></a><a class="code" href="classbdm_1_1fnc.html#6277b11d7fffc7ef8a2fa3e84ae5bad4">00221</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 ) {
184<a name="l00222"></a>00222                 <span class="keywordflow">return</span> vec ( 0 );
185<a name="l00223"></a>00223         };
186<a name="l00224"></a>00224
187<a name="l00226"></a><a class="code" href="classbdm_1_1fnc.html#0786e40fade2663a70d654c1dda5d73e">00226</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 ) {};
188<a name="l00227"></a>00227
189<a name="l00229"></a><a class="code" href="classbdm_1_1fnc.html#083832294da9d1e40804158b979c4341">00229</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>;}
190<a name="l00230"></a>00230 };
191<a name="l00231"></a>00231
192<a name="l00232"></a>00232 <span class="keyword">class </span>mpdf;
193<a name="l00233"></a>00233
194<a name="l00235"></a>00235
195<a name="l00236"></a><a class="code" href="classbdm_1_1epdf.html">00236</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> {
196<a name="l00237"></a>00237 <span class="keyword">protected</span>:
197<a name="l00239"></a><a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce">00239</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1epdf.html#16adac20ec7fe07e1ea0b27d917788ce" title="dimension of the random variable">dim</a>;
198<a name="l00241"></a><a class="code" href="classbdm_1_1epdf.html#62c5b8ff71d9ebe6cd58d3c342eb1dc8">00241</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>;
199<a name="l00242"></a>00242
200<a name="l00243"></a>00243 <span class="keyword">public</span>:
201<a name="l00255"></a>00255         <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> ( ) {};
202<a name="l00256"></a>00256         <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> &amp;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> ) {};
203<a name="l00257"></a>00257         epdf ( <span class="keyword">const</span> RV &amp;rv0 ) {<a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff" title="Name its rv.">set_rv</a> ( rv0 );};
204<a name="l00258"></a>00258         <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;}
205<a name="l00260"></a>00260
206<a name="l00263"></a>00263
207<a name="l00265"></a><a class="code" href="classbdm_1_1epdf.html#75ba57a89681a8efec488f5796c7cc1b">00265</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"not implemneted"</span> );<span class="keywordflow">return</span> vec ( 0 );};
208<a name="l00267"></a>00267         <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>;
209<a name="l00269"></a><a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249">00269</a>         <span class="keyword">virtual</span> <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1epdf.html#deab266d63c236c277538867d5c3f249" title="Compute log-probability of argument val.">evallog</a> ( <span class="keyword">const</span> vec &amp;val )<span class="keyword"> const </span>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"not implemneted"</span> );<span class="keywordflow">return</span> 0.0;};
210<a name="l00271"></a><a class="code" href="classbdm_1_1epdf.html#34956d4dd3176eeb5937cf48a1546b62">00271</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>{
211<a name="l00272"></a>00272                 vec x ( Val.cols() );
212<a name="l00273"></a>00273                 <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#deab266d63c236c277538867d5c3f249" title="Compute log-probability of argument val.">evallog</a> ( Val.get_col ( i ) ) ;}
213<a name="l00274"></a>00274                 <span class="keywordflow">return</span> x;
214<a name="l00275"></a>00275         }
215<a name="l00277"></a><a class="code" href="classbdm_1_1epdf.html#e584eac5579c1b6384947ecf66166c77">00277</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="Description of the random variable.">rv</a> )<span class="keyword"> const  </span>{<a class="code" href="group__errorhandlingfunc.html#g365842195965b61b7ab3ac7938c58b77" title="Display a warning message.">it_warning</a> ( <span class="stringliteral">"Not implemented"</span> ); <span class="keywordflow">return</span> NULL;}
216<a name="l00279"></a><a class="code" href="classbdm_1_1epdf.html#3fb2ece54f720b62ad325e61214fa0a1">00279</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="Description of the random variable.">rv</a> )<span class="keyword"> const </span>{<a class="code" href="group__errorhandlingfunc.html#g365842195965b61b7ab3ac7938c58b77" title="Display a warning message.">it_warning</a> ( <span class="stringliteral">"Not implemented"</span> ); <span class="keywordflow">return</span> NULL;}
217<a name="l00281"></a><a class="code" href="classbdm_1_1epdf.html#0c8d13d788094605b3cdab20a6e22b75">00281</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"not implemneted"</span> );<span class="keywordflow">return</span> vec ( 0 );};
218<a name="l00283"></a><a class="code" href="classbdm_1_1epdf.html#4543d40ccfc5d549b0c831a57a3ed84c">00283</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"not implemneted"</span> );<span class="keywordflow">return</span> vec ( 0 );};
219<a name="l00285"></a><a class="code" href="classbdm_1_1epdf.html#0a7705e80601957f7226d333a0a6d539">00285</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 &amp;lb, vec &amp;ub, <span class="keywordtype">double</span> percentage=0.95 )<span class="keyword"> const </span>{
220<a name="l00286"></a>00286                 vec mea=<a class="code" href="classbdm_1_1epdf.html#0c8d13d788094605b3cdab20a6e22b75" title="return expected value">mean</a>(); vec std=<a class="code" href="group__miscfunc.html#g20af7c97287e8275db8c2b0f650310ac" title="Square root of the elements.">sqrt</a> ( <a class="code" href="classbdm_1_1epdf.html#4543d40ccfc5d549b0c831a57a3ed84c" title="return expected variance (not covariance!)">variance</a>() );
221<a name="l00287"></a>00287                 lb = mea-2*std; ub=mea+2*std;
222<a name="l00288"></a>00288         };
[287]223<a name="l00290"></a>00290
[353]224<a name="l00296"></a>00296
225<a name="l00298"></a><a class="code" href="classbdm_1_1epdf.html#f423e28448dbb69ef4905295ec8de8ff">00298</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> &amp;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>
226<a name="l00300"></a><a class="code" href="classbdm_1_1epdf.html#c4b863ff84c7a4882fb3ad18556027f9">00300</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>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() );<span class="keywordflow">return</span> b;}
227<a name="l00302"></a><a class="code" href="classbdm_1_1epdf.html#a4ab378d5e004c3ff3e2d4e64f7bba21">00302</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="Return name (fails when isnamed is false).">_rv</a>()<span class="keyword"> const </span>{<a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <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>;}
228<a name="l00304"></a>00304
229<a name="l00307"></a>00307
230<a name="l00309"></a><a class="code" href="classbdm_1_1epdf.html#7083a65f7b7a0d0d13b2c516bd2ec29c">00309</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>;}
[287]231<a name="l00311"></a>00311
[353]232<a name="l00312"></a>00312 };
233<a name="l00313"></a>00313
[287]234<a name="l00314"></a>00314
[353]235<a name="l00316"></a>00316 <span class="comment">//TODO Samplecond can be generalized</span>
236<a name="l00317"></a>00317
237<a name="l00318"></a><a class="code" href="classbdm_1_1mpdf.html">00318</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> {
238<a name="l00319"></a>00319 <span class="keyword">protected</span>:
239<a name="l00321"></a><a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6">00321</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1mpdf.html#7c1900976ff13dbc09c9729b3bbff9e6" title="dimension of the condition">dimc</a>;
240<a name="l00323"></a><a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288">00323</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>;
241<a name="l00325"></a><a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0">00325</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>;
242<a name="l00326"></a>00326 <span class="keyword">public</span>:
243<a name="l00329"></a>00329
244<a name="l00330"></a>00330         <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> ( ) {};
245<a name="l00332"></a><a class="code" href="classbdm_1_1mpdf.html#7089fcbd918f201fdf4256119e31f173">00332</a>         <a class="code" href="classbdm_1_1mpdf.html#7089fcbd918f201fdf4256119e31f173" title="copy constructor does not set pointer ep - has to be done in offsprings!">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> &amp;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> ) {};
246<a name="l00334"></a>00334
247<a name="l00337"></a>00337
248<a name="l00339"></a><a class="code" href="classbdm_1_1mpdf.html#f0c1db6fcbb3aae2dd6123884457a367">00339</a>         <span class="keyword">virtual</span> vec <a class="code" href="classbdm_1_1mpdf.html#f0c1db6fcbb3aae2dd6123884457a367" title="Returns a sample from the density conditioned on cond, .">samplecond</a> ( <span class="keyword">const</span> vec &amp;cond ) {
249<a name="l00340"></a>00340                 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 );
250<a name="l00341"></a>00341                 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#75ba57a89681a8efec488f5796c7cc1b" title="Returns a sample,  from density .">sample</a>();
251<a name="l00342"></a>00342                 <span class="keywordflow">return</span> temp;
252<a name="l00343"></a>00343         };
253<a name="l00345"></a><a class="code" href="classbdm_1_1mpdf.html#afe4185b26baeb03688202e254d3b005">00345</a>         <span class="keyword">virtual</span> mat <a class="code" href="classbdm_1_1mpdf.html#afe4185b26baeb03688202e254d3b005" title="Returns.">samplecond_m</a> ( <span class="keyword">const</span> vec &amp;cond, <span class="keywordtype">int</span> N ) {
254<a name="l00346"></a>00346                 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 );
255<a name="l00347"></a>00347                 mat 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#7083a65f7b7a0d0d13b2c516bd2ec29c" title="Size of the random variable.">dimension</a>(),N ); vec 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#7083a65f7b7a0d0d13b2c516bd2ec29c" title="Size of the random variable.">dimension</a>() );
256<a name="l00348"></a>00348                 <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#75ba57a89681a8efec488f5796c7cc1b" title="Returns a sample,  from density .">sample</a>() ;temp.set_col ( i, smp );}
257<a name="l00349"></a>00349                 <span class="keywordflow">return</span> temp;
258<a name="l00350"></a>00350         };
259<a name="l00352"></a><a class="code" href="classbdm_1_1mpdf.html#db94784a9aacf74c65c5a014f6743530">00352</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 ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"Not implemented"</span> );};
260<a name="l00353"></a>00353
261<a name="l00355"></a><a class="code" href="classbdm_1_1mpdf.html#6336a8a72462e2a56a3989a220f18b1b">00355</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 ) {
262<a name="l00356"></a>00356                 <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#deab266d63c236c277538867d5c3f249" title="Compute log-probability of argument val.">evallog</a> ( dt );          <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( std::isfinite ( tmp ),<span class="stringliteral">"Infinite value"</span> ); <span class="keywordflow">return</span> tmp;
263<a name="l00357"></a>00357         };
264<a name="l00358"></a>00358
265<a name="l00360"></a><a class="code" href="classbdm_1_1mpdf.html#0b0ed1ed663071bb7cf4a1349eb94fcb">00360</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 );};
266<a name="l00361"></a>00361
267<a name="l00364"></a>00364
268<a name="l00365"></a>00365         <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> _rv() {<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#a4ab378d5e004c3ff3e2d4e64f7bba21" title="Return name (fails when isnamed is false).">_rv</a>();}
269<a name="l00366"></a>00366         RV _rvc() {<a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( 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>;}
270<a name="l00367"></a>00367         <span class="keywordtype">int</span> dimension() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;dimension();}
271<a name="l00368"></a>00368         <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>;}
272<a name="l00369"></a>00369         epdf&amp; _epdf() {<span class="keywordflow">return</span> *<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>;}
273<a name="l00370"></a>00370         epdf* _e() {<span class="keywordflow">return</span> <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>;}
274<a name="l00372"></a>00372
275<a name="l00375"></a>00375         <span class="keywordtype">void</span> set_rvc ( <span class="keyword">const</span> RV &amp;rvc0 ) {<a class="code" href="classbdm_1_1mpdf.html#5a5f08950daa08b85b01ddf4e1c36288" title="random variable in condition">rvc</a>=rvc0;}
276<a name="l00376"></a>00376         <span class="keywordtype">void</span> set_rv ( <span class="keyword">const</span> RV &amp;rv0 ) {<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;set_rv ( rv0 );}
277<a name="l00377"></a>00377         <span class="keywordtype">bool</span> isnamed() {<span class="keywordflow">return</span> ( <a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>-&gt;isnamed() ) &amp;&amp; ( <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() );}
278<a name="l00379"></a>00379 };
279<a name="l00380"></a>00380
280<a name="l00406"></a><a class="code" href="classbdm_1_1datalink.html">00406</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> {
281<a name="l00407"></a>00407 <span class="keyword">protected</span>:
282<a name="l00409"></a><a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c">00409</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a>;
283<a name="l00411"></a><a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045">00411</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a>;
284<a name="l00413"></a><a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb">00413</a>         ivec <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a>;
285<a name="l00414"></a>00414 <span class="keyword">public</span>:
286<a name="l00416"></a><a class="code" href="classbdm_1_1datalink.html#b8b5d5074fc1310a37c3713d58b900e2">00416</a>         <a class="code" href="classbdm_1_1datalink.html#b8b5d5074fc1310a37c3713d58b900e2" title="Constructor.">datalink</a> () {};
287<a name="l00417"></a>00417         <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> &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;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 );};
288<a name="l00419"></a><a class="code" href="classbdm_1_1datalink.html#fb8337aac8329eb0ca96bfb3ff73ec67">00419</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> &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;rv_up ) {
289<a name="l00420"></a>00420                 <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> = rv.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>();
290<a name="l00421"></a>00421                 <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a> = rv_up.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>();
291<a name="l00422"></a>00422                 <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a>= ( rv.<a class="code" href="classbdm_1_1RV.html#cbebdb5e0d30101a6eb63550ef701c55">dataind</a> ( rv_up ) );
[290]292<a name="l00423"></a>00423
[353]293<a name="l00424"></a>00424                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces 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> );
294<a name="l00425"></a>00425         }
295<a name="l00427"></a><a class="code" href="classbdm_1_1datalink.html#724486a37002e7a3a73366c3af9b82ad">00427</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink.html#724486a37002e7a3a73366c3af9b82ad" title="set connection using indeces">set_connection</a> ( <span class="keywordtype">int</span> ds, <span class="keywordtype">int</span> us, <span class="keyword">const</span> ivec &amp;upind ) {
296<a name="l00428"></a>00428                 <a class="code" href="classbdm_1_1datalink.html#bc8ff50395fc1066a1e486dea18d987c" title="Remember how long val should be.">downsize</a> = ds;
297<a name="l00429"></a>00429                 <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a> = us;
298<a name="l00430"></a>00430                 <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a>= upind;
299<a name="l00431"></a>00431
300<a name="l00432"></a>00432                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces 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> );
301<a name="l00433"></a>00433         }
302<a name="l00435"></a><a class="code" href="classbdm_1_1datalink.html#5b5368c3899eed2f670d09142c2548e5">00435</a>         vec <a class="code" href="classbdm_1_1datalink.html#5b5368c3899eed2f670d09142c2548e5" title="Get val for myself from val of &amp;quot;Up&amp;quot;.">pushdown</a> ( <span class="keyword">const</span> vec &amp;val_up ) {
303<a name="l00436"></a>00436                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> );
304<a name="l00437"></a>00437                 <span class="keywordflow">return</span> get_vec ( val_up,<a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a> );
305<a name="l00438"></a>00438         }
306<a name="l00440"></a><a class="code" href="classbdm_1_1datalink.html#610e966f2edb1fcd6a5b28267c571793">00440</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink.html#610e966f2edb1fcd6a5b28267c571793" title="Fill val of &amp;quot;Up&amp;quot; by my pieces.">pushup</a> ( vec &amp;val_up, <span class="keyword">const</span> vec &amp;val ) {
307<a name="l00441"></a>00441                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <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> );
308<a name="l00442"></a>00442                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> );
309<a name="l00443"></a>00443                 set_subvector ( val_up, <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a>, val );
310<a name="l00444"></a>00444         }
311<a name="l00445"></a>00445 };
312<a name="l00446"></a>00446
313<a name="l00448"></a><a class="code" href="classbdm_1_1datalink__m2e.html">00448</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> <a class="code" href="classbdm_1_1datalink.html" title="DataLink is a connection between two data vectors Up and Down.">datalink</a> {
314<a name="l00449"></a>00449 <span class="keyword">protected</span>:
315<a name="l00451"></a><a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9">00451</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>;
316<a name="l00453"></a><a class="code" href="classbdm_1_1datalink__m2e.html#43fe3e2d7aa04134053ec15793722783">00453</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>;
317<a name="l00455"></a><a class="code" href="classbdm_1_1datalink__m2e.html#9e0048d87bda78e9b8d38761440c7abc">00455</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>;
318<a name="l00456"></a>00456
319<a name="l00457"></a>00457 <span class="keyword">public</span>:
320<a name="l00458"></a>00458         <a class="code" href="classbdm_1_1datalink__m2e.html" title="data link between">datalink_m2e</a>() {};
321<a name="l00460"></a><a class="code" href="classbdm_1_1datalink__m2e.html#add918bf877c48714cff51ee4ad989b0">00460</a>         <span class="keywordtype">void</span> <a class="code" href="classbdm_1_1datalink__m2e.html#add918bf877c48714cff51ee4ad989b0" title="Constructor.">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> &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 ) {
322<a name="l00461"></a>00461                 <a class="code" href="classbdm_1_1datalink__m2e.html#add918bf877c48714cff51ee4ad989b0" title="Constructor.">datalink::set_connection</a> ( rv,rv_up );
323<a name="l00462"></a>00462                 <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a>=  rvc.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>();
324<a name="l00463"></a>00463                 <span class="comment">//establish v2c connection</span>
325<a name="l00464"></a>00464                 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> );
326<a name="l00465"></a>00465         }
327<a name="l00467"></a><a class="code" href="classbdm_1_1datalink__m2e.html#7ba850822e3cf9a4b079afc7e7c01703">00467</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 ) {
328<a name="l00468"></a>00468                 vec tmp ( <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> );
329<a name="l00469"></a>00469                 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> ) );
330<a name="l00470"></a>00470                 <span class="keywordflow">return</span> tmp;
331<a name="l00471"></a>00471         }
332<a name="l00472"></a>00472         <span class="keywordtype">void</span> pushup_cond ( vec &amp;val_up, <span class="keyword">const</span> vec &amp;val, <span class="keyword">const</span> vec &amp;cond ) {
333<a name="l00473"></a>00473                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <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> );
334<a name="l00474"></a>00474                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <a class="code" href="classbdm_1_1datalink.html#0f0d0af6848114d96108d3a6af25a045" title="Remember how long val of &amp;quot;Up&amp;quot; should be.">upsize</a>==val_up.length(),<span class="stringliteral">"Wrong val_up"</span> );
335<a name="l00475"></a>00475                 set_subvector ( val_up, <a class="code" href="classbdm_1_1datalink.html#240b37d793c6bdc32824ed0a0541b1eb" title="val-to-val link, indeces of the upper val">v2v_up</a>, val );
336<a name="l00476"></a>00476                 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 );
337<a name="l00477"></a>00477         }
338<a name="l00478"></a>00478 };
339<a name="l00481"></a><a class="code" href="classbdm_1_1datalink__m2m.html">00481</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> {
340<a name="l00482"></a>00482 <span class="keyword">protected</span>:
341<a name="l00484"></a><a class="code" href="classbdm_1_1datalink__m2m.html#d64dda6f92654628276bbcebb7600d0d">00484</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>;
342<a name="l00486"></a><a class="code" href="classbdm_1_1datalink__m2m.html#0ea05c47de4538b7042ff878e6e309a5">00486</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>;
343<a name="l00487"></a>00487 <span class="keyword">public</span>:
344<a name="l00489"></a><a class="code" href="classbdm_1_1datalink__m2m.html#3baed32cb66f0e2f731032cec0a4d287">00489</a>         <a class="code" href="classbdm_1_1datalink__m2m.html#3baed32cb66f0e2f731032cec0a4d287" title="Constructor.">datalink_m2m</a>() {};
345<a name="l00490"></a>00490         <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> &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 ) {
346<a name="l00491"></a>00491                 <a class="code" href="classbdm_1_1datalink__m2e.html#add918bf877c48714cff51ee4ad989b0" title="Constructor.">datalink_m2e::set_connection</a> ( rv, rvc, rv_up );
347<a name="l00492"></a>00492                 <span class="comment">//establish c2c connection</span>
348<a name="l00493"></a>00493                 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> );
349<a name="l00494"></a>00494                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( <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> );
350<a name="l00495"></a>00495         }
351<a name="l00497"></a><a class="code" href="classbdm_1_1datalink__m2m.html#911fb7fcfa80d4e31db8e0cc03141a90">00497</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 ) {
352<a name="l00498"></a>00498                 vec tmp ( <a class="code" href="classbdm_1_1datalink__m2e.html#3fb26525e0bfa7932a37c5cd404801f9" title="Remember how long cond should be.">condsize</a> );
353<a name="l00499"></a>00499                 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> ) );
354<a name="l00500"></a>00500                 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> ) );
355<a name="l00501"></a>00501                 <span class="keywordflow">return</span> tmp;
356<a name="l00502"></a>00502         }
357<a name="l00504"></a>00504
358<a name="l00505"></a>00505 };
359<a name="l00506"></a>00506
360<a name="l00512"></a><a class="code" href="classbdm_1_1logger.html">00512</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
361<a name="l00513"></a>00513 <span class="keyword">protected</span>:
362<a name="l00515"></a><a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f">00515</a>         <a class="code" href="classitpp_1_1Array.html">Array&lt;RV&gt;</a> <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>;
363<a name="l00517"></a><a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02">00517</a>         <a class="code" href="classitpp_1_1Array.html">Array&lt;string&gt;</a> <a class="code" href="classbdm_1_1logger.html#68ca1f05cf3d414ba0f39211977b8d02" title="Names of logged quantities, e.g. names of algorithm variants.">names</a>;
364<a name="l00518"></a>00518 <span class="keyword">public</span>:
365<a name="l00520"></a><a class="code" href="classbdm_1_1logger.html#f9c2fe7db5ddbda0098a4a5bd2c29fbf">00520</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 ) {}
[290]366<a name="l00521"></a>00521
[353]367<a name="l00524"></a><a class="code" href="classbdm_1_1logger.html#47aa3b098a81d02a8ffaea2c9fa09ac7">00524</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> prefix=<span class="stringliteral">""</span> ) {
368<a name="l00525"></a>00525                 <span class="keywordtype">int</span> id;
369<a name="l00526"></a>00526                 <span class="keywordflow">if</span> ( rv.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() &gt;0 ) {
370<a name="l00527"></a>00527                         <span class="keywordtype">id</span>=<a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>.<a class="code" href="classitpp_1_1Array.html#53629018f146084de506c33c69fa19a3" title="Returns the number of data elements in the array object.">length</a>();
371<a name="l00528"></a>00528                         <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>
372<a name="l00529"></a>00529                         <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a>.<a class="code" href="classitpp_1_1Array.html#3709407656d26aae2c5d3e844316d4e0" title="Resizing an Array&amp;lt;T&amp;gt;.">set_length</a> ( <span class="keywordtype">id</span>+1,<span class="keyword">true</span> );
373<a name="l00530"></a>00530                         <a class="code" href="classbdm_1_1logger.html#c54566a4554089443a20b2a49ed1a61f" title="RVs of all logged variables.">entries</a> ( <span class="keywordtype">id</span> ) = rv;
374<a name="l00531"></a>00531                 }
375<a name="l00532"></a>00532                 <span class="keywordflow">else</span> { <span class="keywordtype">id</span> =-1;}
376<a name="l00533"></a>00533                 <span class="keywordflow">return</span> id; <span class="comment">// identifier of the last entry</span>
377<a name="l00534"></a>00534         }
378<a name="l00535"></a>00535
379<a name="l00537"></a>00537         <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;
380<a name="l00539"></a>00539         <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> &amp;d ) =0;
381<a name="l00540"></a>00540
382<a name="l00542"></a>00542         <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;
383<a name="l00543"></a>00543
384<a name="l00545"></a><a class="code" href="classbdm_1_1logger.html#0e3fc7d13e23a8f71fc7c8a5ea68559d">00545</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>() {};
385<a name="l00546"></a>00546
386<a name="l00548"></a><a class="code" href="classbdm_1_1logger.html#cd6291a73491f755dc7bc88d022eb2fc">00548</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>() {};
387<a name="l00549"></a>00549
388<a name="l00550"></a>00550 };
389<a name="l00551"></a>00551
390<a name="l00555"></a><a class="code" href="classbdm_1_1mepdf.html">00555</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> {
391<a name="l00556"></a>00556 <span class="keyword">public</span>:
392<a name="l00558"></a><a class="code" href="classbdm_1_1mepdf.html#8c320c7f878fb665e93aec1499ef53ff">00558</a>         <a class="code" href="classbdm_1_1mepdf.html#8c320c7f878fb665e93aec1499ef53ff" title="Default constructor.">mepdf</a> ( <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> ( ) {<a class="code" href="classbdm_1_1mpdf.html#5eea43c56d38e4441bfb30270db949c0" title="pointer to internal epdf">ep</a>= em ;};
393<a name="l00559"></a>00559         <a class="code" href="classbdm_1_1mepdf.html#8c320c7f878fb665e93aec1499ef53ff" 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> ( ) {<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 );};
394<a name="l00560"></a><a class="code" href="classbdm_1_1mepdf.html#b2b898058169e7ed9a6dc7cd4e1c147c">00560</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 ) {}
395<a name="l00561"></a>00561 };
396<a name="l00562"></a>00562
397<a name="l00565"></a><a class="code" href="classbdm_1_1compositepdf.html">00565</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> {
398<a name="l00566"></a>00566 <span class="keyword">protected</span>:
399<a name="l00568"></a><a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a">00568</a>         <span class="keywordtype">int</span> <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a>;
400<a name="l00570"></a><a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a">00570</a>         <a class="code" href="classitpp_1_1Array.html">Array&lt;mpdf*&gt;</a> <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a>;
[302]401<a name="l00571"></a>00571 <span class="keyword">public</span>:
[353]402<a name="l00572"></a>00572         <a class="code" href="classbdm_1_1compositepdf.html" title="Abstract composition of pdfs, will be used for specific classes this abstract class...">compositepdf</a> ( <a class="code" href="classitpp_1_1Array.html">Array&lt;mpdf*&gt;</a> A0 ) : <a class="code" href="classbdm_1_1compositepdf.html#7ca3cd8ed92a7154f5b9ff13b1e9d52a" title="Number of mpdfs in the composite.">n</a> ( A0.<a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>() ), <a class="code" href="classbdm_1_1compositepdf.html#23faf2debc4dde10836393b8c665914a" title="Elements of composition.">mpdfs</a> ( A0 ) {};
403<a name="l00574"></a>00574         <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> );
404<a name="l00576"></a>00576         <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 );
405<a name="l00577"></a>00577 };
406<a name="l00578"></a>00578
407<a name="l00586"></a><a class="code" href="classbdm_1_1DS.html">00586</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> {
408<a name="l00587"></a>00587 <span class="keyword">protected</span>:
409<a name="l00588"></a>00588         <span class="keywordtype">int</span> dtsize;
410<a name="l00589"></a>00589         <span class="keywordtype">int</span> utsize;
411<a name="l00591"></a><a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec">00591</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>;
412<a name="l00593"></a><a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702">00593</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>
413<a name="l00595"></a><a class="code" href="classbdm_1_1DS.html#b4b7bdb9d98a11ab587aa0ad803bb867">00595</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;
414<a name="l00596"></a>00596 <span class="keyword">public</span>:
415<a name="l00598"></a><a class="code" href="classbdm_1_1DS.html#bedd82e08ae6686b4b1732a122c2201a">00598</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> ( ) {};
416<a name="l00600"></a><a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba">00600</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 ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"abstract class"</span> );};
417<a name="l00602"></a><a class="code" href="classbdm_1_1DS.html#a0bf451c7fea73ff4150ff9e83da4832">00602</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 &amp;dt, <span class="keyword">const</span> ivec &amp;indeces ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"abstract class"</span> );};
418<a name="l00604"></a><a class="code" href="classbdm_1_1DS.html#f5baeee8528e5e4c5c7a344b2e3f192a">00604</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 ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"abstract class"</span> );};
419<a name="l00606"></a><a class="code" href="classbdm_1_1DS.html#6bced13f84f6bf04449728817e049e1d">00606</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 &amp;ut, <span class="keyword">const</span> ivec &amp;indeces ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"abstract class"</span> );};
420<a name="l00607"></a>00607
421<a name="l00609"></a>00609         <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;
422<a name="l00610"></a>00610
423<a name="l00612"></a><a class="code" href="classbdm_1_1DS.html#3c9e417d85052396fe74c164857ee67e">00612</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> &amp;L ) {
424<a name="l00613"></a>00613                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( dtsize==<a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>(),<span class="stringliteral">""</span> );
425<a name="l00614"></a>00614                 <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a> ( utsize==<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>(),<span class="stringliteral">""</span> );
426<a name="l00615"></a>00615
427<a name="l00616"></a>00616                 <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> );
428<a name="l00617"></a>00617                 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> );
429<a name="l00618"></a>00618         }
430<a name="l00620"></a><a class="code" href="classbdm_1_1DS.html#e03d2c93cddd5c3a03107ddd3ae999cf">00620</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> &amp;L ) {
431<a name="l00621"></a>00621                 vec tmp ( <a class="code" href="classbdm_1_1DS.html#d8db125ac35f3b8eeb773fc9a4cb45ec" title="Description of data returned by getdata().">Drv</a>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() +<a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() );
432<a name="l00622"></a>00622                 <a class="code" href="classbdm_1_1DS.html#c1ec862dec7164fee3284c6f9a51b0ba" title="Returns full vector of observed data=[output, input].">getdata</a> ( tmp );
433<a name="l00623"></a>00623                 <span class="comment">// d is first in getdata</span>
434<a name="l00624"></a>00624                 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>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() ) );
435<a name="l00625"></a>00625                 <span class="comment">// u follows after d in getdata</span>
436<a name="l00626"></a>00626                 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>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>(), <a class="code" href="classbdm_1_1DS.html#9ac5970128fd39f4f6226723cd1de702" title="Description of data witten by by write().">Urv</a>.<a class="code" href="classbdm_1_1RV.html#de30156104f61d86c94f758861418089">_dsize</a>() ) );
437<a name="l00627"></a>00627         }
438<a name="l00629"></a><a class="code" href="classbdm_1_1DS.html#ed2113851294ed96abe60fc3f05aa660">00629</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> );}
439<a name="l00631"></a><a class="code" href="classbdm_1_1DS.html#a6a9a375da64c82fc227bda3418a1f33">00631</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_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>;}
440<a name="l00633"></a><a class="code" href="classbdm_1_1DS.html#821297dfd839d68691324f024e57252f">00633</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> &amp;drv, <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> &amp;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;}
441<a name="l00634"></a>00634 };
442<a name="l00635"></a>00635
443<a name="l00657"></a><a class="code" href="classbdm_1_1BM.html">00657</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_1bdmroot.html" title="Root class of BDM objects.">bdmroot</a> {
444<a name="l00658"></a>00658 <span class="keyword">protected</span>:
445<a name="l00660"></a><a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed">00660</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>;
446<a name="l00662"></a><a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a">00662</a>         <span class="keywordtype">double</span> <a class="code" href="classbdm_1_1BM.html#4064b6559d962633e4372b12f4cd204a" title="Logarithm of marginalized data likelihood.">ll</a>;
447<a name="l00664"></a><a class="code" href="classbdm_1_1BM.html#faff0ad12556fe7dc0e2807d4fd938ee">00664</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>;
448<a name="l00665"></a>00665 <span class="keyword">public</span>:
449<a name="l00668"></a>00668
450<a name="l00669"></a>00669         <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) {
451<a name="l00670"></a>00670                 <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> <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#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>};
452<a name="l00671"></a>00671         <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> &amp;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> ) {}
453<a name="l00674"></a><a class="code" href="classbdm_1_1BM.html#40f01563651c1e4fe926189e41a6b8ab">00674</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;};
454<a name="l00676"></a>00676
[302]455<a name="l00679"></a>00679
[353]456<a name="l00683"></a>00683         <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;
457<a name="l00685"></a>00685         <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 );
458<a name="l00688"></a><a class="code" href="classbdm_1_1BM.html#50257e0c1e5b5c73153ea6e716ad8ae0">00688</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> 0.0;}
459<a name="l00690"></a><a class="code" href="classbdm_1_1BM.html#0e8ebe61fb14990abe1254bd3dda5fae">00690</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;}
460<a name="l00691"></a>00691
461<a name="l00693"></a><a class="code" href="classbdm_1_1BM.html#688d7a2aced1e06aa1c468d73a9e5eba">00693</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> NULL;};
462<a name="l00695"></a><a class="code" href="classbdm_1_1BM.html#598b25e3f3d96a5bc00a5faeb5b3c912">00695</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>{<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"Not implemented"</span> );<span class="keywordflow">return</span> NULL;};
[302]463<a name="l00697"></a>00697
[353]464<a name="l00702"></a>00702
465<a name="l00704"></a><a class="code" href="classbdm_1_1BM.html#de95c590d650b2db423d53992f0906cb">00704</a>         <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> rvc;
466<a name="l00706"></a><a class="code" href="classbdm_1_1BM.html#b5e8e016be16e4029c789063cb3f3b33">00706</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#b5e8e016be16e4029c789063cb3f3b33" title="access function">_rvc</a>()<span class="keyword"> const </span>{<span class="keywordflow">return</span> rvc;}
[302]467<a name="l00707"></a>00707
[353]468<a name="l00709"></a><a class="code" href="classbdm_1_1BM.html#6799f4b16a6a59ed58b1d0d6e17116f4">00709</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 &amp;val ) {<a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a> ( <span class="stringliteral">"Not implemented!"</span> );};
469<a name="l00710"></a>00710
470<a name="l00712"></a>00712
471<a name="l00713"></a>00713
[302]472<a name="l00716"></a>00716
[353]473<a name="l00717"></a>00717         <span class="keyword">const</span> <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a>&amp; _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>;}
474<a name="l00718"></a>00718         <span class="keywordtype">void</span> set_drv ( <span class="keyword">const</span> RV &amp;rv ) {<a class="code" href="classbdm_1_1BM.html#c400357e37d27a4834b2b1d9211009ed" title="Random variable of the data (optional).">drv</a>=rv;}
475<a name="l00719"></a>00719         <span class="keywordtype">void</span> set_rv ( <span class="keyword">const</span> RV &amp;rv ) {<span class="keyword">const_cast&lt;</span>epdf&amp;<span class="keyword">&gt;</span> ( posterior() ).set_rv ( rv );}
476<a name="l00720"></a>00720         <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>;}
477<a name="l00721"></a>00721         <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;}
478<a name="l00722"></a>00722         <span class="keyword">virtual</span> <span class="keyword">const</span> epdf&amp; posterior() <span class="keyword">const</span> =0;
479<a name="l00723"></a>00723         <span class="keyword">virtual</span> <span class="keyword">const</span> epdf* _e() <span class="keyword">const</span> =0;
480<a name="l00725"></a>00725
481<a name="l00728"></a>00728
482<a name="l00730"></a><a class="code" href="classbdm_1_1BM.html#adb539fa2585c8a2a91c54e5aebf330e">00730</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: &amp;quot;logbounds,logll&amp;quot;.">set_options</a> ( <span class="keyword">const</span> <span class="keywordtype">string</span> &amp;opt ) {
483<a name="l00731"></a>00731                 <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;
484<a name="l00732"></a>00732                 <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;}
485<a name="l00733"></a>00733                 <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;}
486<a name="l00734"></a>00734         }
487<a name="l00736"></a><a class="code" href="classbdm_1_1BM.html#109c1a626a69031658e3a44e9e500cca">00736</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>;
488<a name="l00737"></a>00737
489<a name="l00739"></a><a class="code" href="classbdm_1_1BM.html#398b7e7a162cf37d0022d6126495f207">00739</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>;
490<a name="l00741"></a><a class="code" href="classbdm_1_1BM.html#2298ca6af9b13a78a4c9e18ab43f1827">00741</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> &amp;L, <span class="keyword">const</span> <span class="keywordtype">string</span> &amp;name=<span class="stringliteral">""</span> ) {
491<a name="l00742"></a>00742                 <span class="comment">// internal</span>
492<a name="l00743"></a>00743                 <a class="code" href="classbdm_1_1RV.html" title="Class representing variables, most often random variables.">RV</a> r;
493<a name="l00744"></a>00744                 <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>();}
494<a name="l00745"></a>00745                 <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() );};
495<a name="l00746"></a>00746
496<a name="l00747"></a>00747                 <span class="comment">// Add mean value</span>
497<a name="l00748"></a>00748                 <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> );
498<a name="l00749"></a>00749                 <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> );
499<a name="l00750"></a>00750                 <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> );
500<a name="l00751"></a>00751                 <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>
501<a name="l00752"></a>00752         }
502<a name="l00753"></a>00753         <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> &amp;L ) {
503<a name="l00754"></a>00754                 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().<a class="code" href="group__statistics.html#gc370b4cfdc6baa037b7442ae6e2c4b5c" title="The mean value.">mean</a>() );
504<a name="l00755"></a>00755                 <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>
505<a name="l00756"></a>00756                         vec ub,lb;
506<a name="l00757"></a>00757                         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 );
507<a name="l00758"></a>00758                         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 );
508<a name="l00759"></a>00759                         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 );
509<a name="l00760"></a>00760                 }
510<a name="l00761"></a>00761                 <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> );
511<a name="l00762"></a>00762         }
512<a name="l00764"></a>00764 };
513<a name="l00765"></a>00765
514<a name="l00766"></a>00766
515<a name="l00767"></a>00767 }; <span class="comment">//namespace</span>
516<a name="l00768"></a>00768 <span class="preprocessor">#endif // BM_H</span>
[91]517</pre></div></div>
[353]518<hr size="1"><address style="text-align: right;"><small>Generated on Tue Jun 2 10:02:11 2009 for mixpp by&nbsp;
[3]519<a href="http://www.doxygen.org/index.html">
[290]520<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.8 </small></address>
[3]521</body>
522</html>
Note: See TracBrowser for help on using the browser.