29 | | <a name="l00014"></a>00014 <span class="keyword">using</span> std::cout; |
30 | | <a name="l00015"></a>00015 <span class="keyword">using</span> std::endl; |
31 | | <a name="l00016"></a>00016 |
32 | | <a name="l00017"></a>00017 <span class="keyword">namespace </span>itpp { |
33 | | <a name="l00018"></a>00018 Array<int> to_Arr(<span class="keyword">const</span> ivec &indices); |
34 | | <a name="l00019"></a>00019 |
35 | | <a name="l00024"></a><a class="code" href="classitpp_1_1Gamma__RNG.html">00024</a> <span class="keyword">class </span><a class="code" href="classitpp_1_1Gamma__RNG.html" title="Gamma distribution.">Gamma_RNG</a> { |
36 | | <a name="l00025"></a>00025 <span class="keyword">public</span>: |
37 | | <a name="l00027"></a>00027 <a class="code" href="classitpp_1_1Gamma__RNG.html#e132860f825f427d63652edd02c25b54" title="constructor. Set lambda.">Gamma_RNG</a>(<span class="keywordtype">double</span> a=1.0, <span class="keywordtype">double</span> b=1.0); |
38 | | <a name="l00029"></a><a class="code" href="classitpp_1_1Gamma__RNG.html#dfaae19411e39aa87e1f72e409b6babe">00029</a> <span class="keywordtype">void</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#dfaae19411e39aa87e1f72e409b6babe" title="Set lambda.">setup</a>(<span class="keywordtype">double</span> a0, <span class="keywordtype">double</span> b0) { alpha=a0; beta=b0;} |
39 | | <a name="l00031"></a>00031 <span class="keywordtype">double</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#8e666953249038f262cece2d3ed952d6" title="get lambda">get_setup</a>() <span class="keyword">const</span>; |
40 | | <a name="l00033"></a><a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e">00033</a> <span class="keywordtype">double</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator()</a>() { <span class="keywordflow">return</span> sample(); } |
41 | | <a name="l00035"></a>00035 vec <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator()</a>(<span class="keywordtype">int</span> n); |
42 | | <a name="l00037"></a>00037 mat <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator()</a>(<span class="keywordtype">int</span> h, <span class="keywordtype">int</span> w); |
43 | | <a name="l00038"></a>00038 <span class="keyword">protected</span>: |
44 | | <a name="l00039"></a>00039 <span class="keyword">private</span>: |
45 | | <a name="l00041"></a>00041 <span class="keywordtype">double</span> sample(); |
46 | | <a name="l00043"></a>00043 <span class="keywordtype">double</span> alpha; |
47 | | <a name="l00045"></a>00045 <span class="keywordtype">double</span> beta; |
48 | | <a name="l00047"></a>00047 Random_Generator RNG; |
49 | | <a name="l00048"></a>00048 Normal_RNG NRNG; |
50 | | <a name="l00050"></a>00050 <span class="keyword">inline</span> <span class="keywordtype">double</span> exp_rand(){<span class="keywordflow">return</span> -std::log(RNG.random_01());} |
51 | | <a name="l00051"></a>00051 <span class="keyword">inline</span> <span class="keywordtype">double</span> unif_rand(){<span class="keywordflow">return</span> RNG.random_01();} |
52 | | <a name="l00052"></a>00052 <span class="keyword">inline</span> <span class="keywordtype">double</span> norm_rand(){<span class="keywordflow">return</span> NRNG.sample();} |
53 | | <a name="l00053"></a>00053 |
54 | | <a name="l00054"></a>00054 }; |
55 | | <a name="l00055"></a>00055 |
56 | | <a name="l00056"></a>00056 <span class="keywordtype">bool</span> qr(<span class="keyword">const</span> mat &A, mat &R); |
57 | | <a name="l00057"></a>00057 |
58 | | <a name="l00058"></a>00058 } |
59 | | </pre></div><hr size="1"><address style="text-align: right;"><small>Generated on Fri Apr 18 11:15:14 2008 for mixpp by |
| 30 | <a name="l00014"></a>00014 <span class="preprocessor">#ifndef ITEX_H</span> |
| 31 | <a name="l00015"></a>00015 <span class="preprocessor"></span><span class="preprocessor">#define ITEX_H </span> |
| 32 | <a name="l00016"></a>00016 <span class="preprocessor"></span> |
| 33 | <a name="l00017"></a>00017 <span class="keyword">using</span> std::cout; |
| 34 | <a name="l00018"></a>00018 <span class="keyword">using</span> std::endl; |
| 35 | <a name="l00019"></a>00019 |
| 36 | <a name="l00020"></a>00020 <span class="keyword">namespace </span>itpp { |
| 37 | <a name="l00021"></a>00021 Array<int> to_Arr ( <span class="keyword">const</span> ivec &indices ); |
| 38 | <a name="l00022"></a>00022 ivec linspace ( <span class="keywordtype">int</span> from, <span class="keywordtype">int</span> to ); |
| 39 | <a name="l00023"></a>00023 |
| 40 | <a name="l00028"></a><a class="code" href="classitpp_1_1Gamma__RNG.html">00028</a> <span class="keyword">class </span><a class="code" href="classitpp_1_1Gamma__RNG.html" title="Gamma distribution.">Gamma_RNG</a> { |
| 41 | <a name="l00029"></a>00029 <span class="keyword">public</span>: |
| 42 | <a name="l00031"></a>00031 <a class="code" href="classitpp_1_1Gamma__RNG.html#e132860f825f427d63652edd02c25b54" title="constructor. Set lambda.">Gamma_RNG</a> ( <span class="keywordtype">double</span> a=1.0, <span class="keywordtype">double</span> b=1.0 ); |
| 43 | <a name="l00033"></a><a class="code" href="classitpp_1_1Gamma__RNG.html#dfaae19411e39aa87e1f72e409b6babe">00033</a> <span class="keywordtype">void</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#dfaae19411e39aa87e1f72e409b6babe" title="Set lambda.">setup</a> ( <span class="keywordtype">double</span> a0, <span class="keywordtype">double</span> b0 ) { alpha=a0; beta=b0;} |
| 44 | <a name="l00035"></a>00035 <span class="keywordtype">double</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#8e666953249038f262cece2d3ed952d6" title="get lambda">get_setup</a>() <span class="keyword">const</span>; |
| 45 | <a name="l00037"></a><a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e">00037</a> <span class="keywordtype">double</span> <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator() </a>() { <span class="keywordflow">return</span> sample(); } |
| 46 | <a name="l00039"></a>00039 vec <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator() </a>( <span class="keywordtype">int</span> n ); |
| 47 | <a name="l00041"></a>00041 mat <a class="code" href="classitpp_1_1Gamma__RNG.html#98e2e65456bca33fca602b2bd607235e" title="Get one sample.">operator() </a>( <span class="keywordtype">int</span> h, <span class="keywordtype">int</span> w ); |
| 48 | <a name="l00042"></a>00042 <span class="keyword">protected</span>: |
| 49 | <a name="l00043"></a>00043 <span class="keyword">private</span>: |
| 50 | <a name="l00045"></a>00045 <span class="keywordtype">double</span> sample(); |
| 51 | <a name="l00047"></a>00047 <span class="keywordtype">double</span> alpha; |
| 52 | <a name="l00049"></a>00049 <span class="keywordtype">double</span> beta; |
| 53 | <a name="l00051"></a>00051 Random_Generator RNG; |
| 54 | <a name="l00052"></a>00052 Normal_RNG NRNG; |
| 55 | <a name="l00054"></a>00054 <span class="keyword">inline</span> <span class="keywordtype">double</span> exp_rand() {<span class="keywordflow">return</span> -std::log ( RNG.random_01() );} |
| 56 | <a name="l00055"></a>00055 <span class="keyword">inline</span> <span class="keywordtype">double</span> unif_rand() {<span class="keywordflow">return</span> RNG.random_01();} |
| 57 | <a name="l00056"></a>00056 <span class="keyword">inline</span> <span class="keywordtype">double</span> norm_rand() {<span class="keywordflow">return</span> NRNG.sample();} |
| 58 | <a name="l00057"></a>00057 |
| 59 | <a name="l00058"></a>00058 }; |
| 60 | <a name="l00059"></a>00059 |
| 61 | <a name="l00060"></a>00060 <span class="keywordtype">bool</span> qr ( <span class="keyword">const</span> mat &A, mat &R ); |
| 62 | <a name="l00061"></a>00061 |
| 63 | <a name="l00062"></a>00062 } |
| 64 | <a name="l00063"></a>00063 <span class="preprocessor">#endif //ITEX_H</span> |
| 65 | </pre></div></div> |
| 66 | <hr size="1"><address style="text-align: right;"><small>Generated on Tue Apr 29 20:46:24 2008 for mixpp by |