root/library/doc/html/itpp__ext_8h_source.html @ 661

Revision 661, 10.2 kB (checked in by smidl, 15 years ago)

doc

Line 
1<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
2<html xmlns="http://www.w3.org/1999/xhtml">
3<head>
4<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
5<title>mixpp: itpp_ext.h Source File</title>
6<link href="tabs.css" rel="stylesheet" type="text/css"/>
7<link href="doxygen.css" rel="stylesheet" type="text/css"/>
8</head>
9<body>
10<!-- Generated by Doxygen 1.6.1 -->
11<script type="text/javascript">
12<!--
13function changeDisplayState (e){
14  var num=this.id.replace(/[^[0-9]/g,'');
15  var button=this.firstChild;
16  var sectionDiv=document.getElementById('dynsection'+num);
17  if (sectionDiv.style.display=='none'||sectionDiv.style.display==''){
18    sectionDiv.style.display='block';
19    button.src='open.gif';
20  }else{
21    sectionDiv.style.display='none';
22    button.src='closed.gif';
23  }
24}
25function initDynSections(){
26  var divs=document.getElementsByTagName('div');
27  var sectionCounter=1;
28  for(var i=0;i<divs.length-1;i++){
29    if(divs[i].className=='dynheader'&&divs[i+1].className=='dynsection'){
30      var header=divs[i];
31      var section=divs[i+1];
32      var button=header.firstChild;
33      if (button!='IMG'){
34        divs[i].insertBefore(document.createTextNode(' '),divs[i].firstChild);
35        button=document.createElement('img');
36        divs[i].insertBefore(button,divs[i].firstChild);
37      }
38      header.style.cursor='pointer';
39      header.onclick=changeDisplayState;
40      header.id='dynheader'+sectionCounter;
41      button.src='closed.gif';
42      section.id='dynsection'+sectionCounter;
43      section.style.display='none';
44      section.style.marginLeft='14px';
45      sectionCounter++;
46    }
47  }
48}
49window.onload = initDynSections;
50-->
51</script>
52<div class="navigation" id="top">
53  <div class="tabs">
54    <ul>
55      <li><a href="main.html"><span>Main&nbsp;Page</span></a></li>
56      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
57      <li><a href="annotated.html"><span>Classes</span></a></li>
58      <li class="current"><a href="files.html"><span>Files</span></a></li>
59    </ul>
60  </div>
61  <div class="tabs">
62    <ul>
63      <li><a href="files.html"><span>File&nbsp;List</span></a></li>
64      <li><a href="globals.html"><span>File&nbsp;Members</span></a></li>
65    </ul>
66  </div>
67<h1>itpp_ext.h</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">//</span>
68<a name="l00002"></a>00002 <span class="comment">// C++ Interface: itpp_ext</span>
69<a name="l00003"></a>00003 <span class="comment">//</span>
70<a name="l00004"></a>00004 <span class="comment">// Description:</span>
71<a name="l00005"></a>00005 <span class="comment">//</span>
72<a name="l00006"></a>00006 <span class="comment">//</span>
73<a name="l00007"></a>00007 <span class="comment">// Author: smidl &lt;smidl@utia.cas.cz&gt;, (C) 2008</span>
74<a name="l00008"></a>00008 <span class="comment">//</span>
75<a name="l00009"></a>00009 <span class="comment">// Copyright: See COPYING file that comes with this distribution</span>
76<a name="l00010"></a>00010 <span class="comment">//</span>
77<a name="l00011"></a>00011 <span class="comment">//</span>
78<a name="l00012"></a>00012 <span class="preprocessor">#include &lt;itpp/itbase.h&gt;</span>
79<a name="l00013"></a>00013
80<a name="l00014"></a>00014 <span class="preprocessor">#ifndef ITEX_H</span>
81<a name="l00015"></a>00015 <span class="preprocessor"></span><span class="preprocessor">#define ITEX_H</span>
82<a name="l00016"></a>00016 <span class="preprocessor"></span>
83<a name="l00017"></a>00017 <span class="keyword">using</span> std::cout;
84<a name="l00018"></a>00018 <span class="keyword">using</span> std::endl;
85<a name="l00019"></a>00019
86<a name="l00020"></a>00020 <span class="keyword">namespace </span>itpp {
87<a name="l00021"></a>00021 Array&lt;int&gt; to_Arr ( <span class="keyword">const</span> ivec &amp;indices );
88<a name="l00022"></a>00022 ivec linspace ( <span class="keywordtype">int</span> from, <span class="keywordtype">int</span> to );
89<a name="l00023"></a>00023
90<a name="l00024"></a>00024 vec get_vec ( <span class="keyword">const</span> vec &amp;v, <span class="keyword">const</span> ivec &amp;indexlist );
91<a name="l00025"></a>00025
92<a name="l00026"></a>00026 bvec operator&amp; ( <span class="keyword">const</span> bvec &amp;a, <span class="keyword">const</span> bvec &amp;b );
93<a name="l00027"></a>00027 bvec operator| ( <span class="keyword">const</span> bvec &amp;a, <span class="keyword">const</span> bvec &amp;b );
94<a name="l00028"></a>00028
95<a name="l00029"></a>00029 bvec operator&gt; ( <span class="keyword">const</span> vec &amp;t1, <span class="keyword">const</span> vec &amp;t2 );
96<a name="l00030"></a>00030
97<a name="l00031"></a>00031 bvec operator&lt; ( <span class="keyword">const</span> vec &amp;t1, <span class="keyword">const</span> vec &amp;t2 );
98<a name="l00032"></a>00032
99<a name="l00033"></a>00033 <span class="comment">// template&lt;class Num_T&gt;</span>
100<a name="l00034"></a>00034 <span class="comment">// void set_subvector(vec &amp;ov, ivec &amp;iv, const Vec&lt;Num_T&gt; &amp;v);</span>
101<a name="l00035"></a>00035
102<a name="l00036"></a>00036 <span class="keywordtype">void</span> set_subvector ( vec &amp;ov, <span class="keyword">const</span> ivec &amp;iv, <span class="keyword">const</span> vec &amp;v );
103<a name="l00037"></a>00037
104<a name="l00038"></a>00038 <span class="keyword">template</span>&lt;<span class="keyword">class</span> Num_T&gt; <span class="keyword">inline</span>
105<a name="l00039"></a>00039 <span class="keywordtype">void</span> set_col_part ( mat &amp;M, <span class="keywordtype">int</span> c, <span class="keyword">const</span> Vec&lt;Num_T&gt; &amp;v ) {
106<a name="l00040"></a>00040         copy_vector ( v.size(), v._data(), M._data() + c*M.rows() );
107<a name="l00041"></a>00041 }
108<a name="l00042"></a>00042
109<a name="l00043"></a>00043 <span class="keyword">const</span> <span class="keywordtype">double</span> inf = std::numeric_limits&lt;double&gt;::infinity();
110<a name="l00044"></a>00044
111<a name="l00045"></a>00045 <span class="preprocessor">#if 0</span>
112<a name="l00046"></a>00046 <span class="preprocessor"></span>
113<a name="l00050"></a>00050 <span class="keyword">class </span>Gamma_RNG {
114<a name="l00051"></a>00051 <span class="keyword">public</span>:
115<a name="l00053"></a>00053         Gamma_RNG ( <span class="keywordtype">double</span> a = 1.0, <span class="keywordtype">double</span> b = 1.0 );
116<a name="l00055"></a>00055         <span class="keywordtype">void</span> setup ( <span class="keywordtype">double</span> a0, <span class="keywordtype">double</span> b0 ) {
117<a name="l00056"></a>00056                 alpha = a0;
118<a name="l00057"></a>00057                 beta = b0;
119<a name="l00058"></a>00058         }
120<a name="l00060"></a>00060         <span class="keywordtype">double</span> get_setup() <span class="keyword">const</span>;
121<a name="l00062"></a>00062         <span class="keywordtype">double</span> operator() () {
122<a name="l00063"></a>00063                 <span class="keywordflow">return</span> sample();
123<a name="l00064"></a>00064         }
124<a name="l00066"></a>00066         vec operator() ( <span class="keywordtype">int</span> n );
125<a name="l00068"></a>00068         mat operator() ( <span class="keywordtype">int</span> h, <span class="keywordtype">int</span> w );
126<a name="l00069"></a>00069 <span class="keyword">protected</span>:
127<a name="l00070"></a>00070 <span class="keyword">private</span>:
128<a name="l00072"></a>00072         <span class="keywordtype">double</span> sample();
129<a name="l00074"></a>00074         <span class="keywordtype">double</span> alpha;
130<a name="l00076"></a>00076         <span class="keywordtype">double</span> beta;
131<a name="l00078"></a>00078         Random_Generator RNG;
132<a name="l00079"></a>00079         Normal_RNG NRNG;
133<a name="l00081"></a>00081         <span class="keyword">inline</span> <span class="keywordtype">double</span> exp_rand() {
134<a name="l00082"></a>00082                 <span class="keywordflow">return</span> -std::log ( RNG.random_01() );
135<a name="l00083"></a>00083         }
136<a name="l00084"></a>00084         <span class="keyword">inline</span> <span class="keywordtype">double</span> unif_rand() {
137<a name="l00085"></a>00085                 <span class="keywordflow">return</span> RNG.random_01();
138<a name="l00086"></a>00086         }
139<a name="l00087"></a>00087         <span class="keyword">inline</span> <span class="keywordtype">double</span> norm_rand() {
140<a name="l00088"></a>00088                 <span class="keywordflow">return</span> NRNG.sample();
141<a name="l00089"></a>00089         }
142<a name="l00090"></a>00090
143<a name="l00091"></a>00091 };
144<a name="l00092"></a>00092 <span class="keywordtype">bool</span> qr ( <span class="keyword">const</span> mat &amp;A, mat &amp;R );
145<a name="l00093"></a>00093 <span class="preprocessor">#endif</span>
146<a name="l00095"></a>00095 <span class="preprocessor">std::string num2str ( double d );</span>
147<a name="l00096"></a>00096 <span class="preprocessor"></span>
148<a name="l00098"></a>00098 std::string num2str ( <span class="keywordtype">int</span> i );
149<a name="l00099"></a>00099
150<a name="l00101"></a>00101 <span class="keywordtype">double</span> psi ( <span class="keywordtype">double</span> );
151<a name="l00102"></a>00102
152<a name="l00104"></a>00104 <span class="keywordtype">void</span> triu(mat &amp;A);
153<a name="l00105"></a>00105
154<a name="l00109"></a>00109 <span class="keywordtype">double</span> randun();
155<a name="l00111"></a>00111 vec randun(<span class="keywordtype">int</span> n);
156<a name="l00113"></a>00113 mat randun(<span class="keywordtype">int</span> n, <span class="keywordtype">int</span> m);
157<a name="l00114"></a>00114
158<a name="l00116"></a>00116 ivec unique(<span class="keyword">const</span> ivec &amp;in);
159<a name="l00117"></a>00117
160<a name="l00119"></a>00119 ivec unique_complement(<span class="keyword">const</span> ivec &amp;in, <span class="keyword">const</span> ivec &amp;base);
161<a name="l00120"></a>00120 }
162<a name="l00121"></a>00121
163<a name="l00122"></a>00122
164<a name="l00123"></a>00123 <span class="preprocessor">#endif //ITEX_H</span>
165</pre></div></div>
166<hr size="1"/><address style="text-align: right;"><small>Generated on Thu Oct 15 00:07:48 2009 for mixpp by&nbsp;
167<a href="http://www.doxygen.org/index.html">
168<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.6.1 </small></address>
169</body>
170</html>
Note: See TracBrowser for help on using the browser.