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

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

doc

Line 
1<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
2<html><head><meta http-equiv="Content-Type" content="text/html;charset=UTF-8">
3<title>mixpp: matfunc.h Source File</title>
4<link href="tabs.css" rel="stylesheet" type="text/css">
5<link href="doxygen.css" rel="stylesheet" type="text/css">
6</head><body>
7<!-- Generated by Doxygen 1.5.8 -->
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>
49<div class="navigation" id="top">
50  <div class="tabs">
51    <ul>
52      <li><a href="main.html"><span>Main&nbsp;Page</span></a></li>
53      <li><a href="pages.html"><span>Related&nbsp;Pages</span></a></li>
54      <li><a href="modules.html"><span>Modules</span></a></li>
55      <li><a href="annotated.html"><span>Classes</span></a></li>
56      <li class="current"><a href="files.html"><span>Files</span></a></li>
57    </ul>
58  </div>
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>
65<h1>matfunc.h</h1><a href="matfunc_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
66<a name="l00030"></a>00030 <span class="preprocessor">#ifndef MATFUNC_H</span>
67<a name="l00031"></a>00031 <span class="preprocessor"></span><span class="preprocessor">#define MATFUNC_H</span>
68<a name="l00032"></a>00032 <span class="preprocessor"></span>
69<a name="l00033"></a>00033 <span class="preprocessor">#include &lt;<a class="code" href="mat_8h.html" title="Matrix Class Definitions.">itpp/base/mat.h</a>&gt;</span>
70<a name="l00034"></a>00034 <span class="preprocessor">#include &lt;<a class="code" href="log__exp_8h.html" title="Logarithmic and exponenential functions - header file.">itpp/base/math/log_exp.h</a>&gt;</span>
71<a name="l00035"></a>00035 <span class="preprocessor">#include &lt;<a class="code" href="elem__math_8h.html" title="Elementary mathematical functions - header file.">itpp/base/math/elem_math.h</a>&gt;</span>
72<a name="l00036"></a>00036 <span class="preprocessor">#include &lt;<a class="code" href="inv_8h.html" title="Definitions of matrix inversion routines.">itpp/base/algebra/inv.h</a>&gt;</span>
73<a name="l00037"></a>00037 <span class="preprocessor">#include &lt;<a class="code" href="svd_8h.html" title="Definitions of Singular Value Decompositions.">itpp/base/algebra/svd.h</a>&gt;</span>
74<a name="l00038"></a>00038
75<a name="l00039"></a>00039
76<a name="l00040"></a>00040 <span class="keyword">namespace </span>itpp
77<a name="l00041"></a>00041 {
78<a name="l00042"></a>00042
79<a name="l00047"></a>00047
80<a name="l00048"></a>00048
81<a name="l00050"></a>00050 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
82<a name="l00051"></a><a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d">00051</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v) { <span class="keywordflow">return</span> v.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>(); }
83<a name="l00052"></a>00052
84<a name="l00054"></a>00054 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
85<a name="l00055"></a><a class="code" href="group__matrix__functions.html#g3c1a2b0972c6a8e1215eb3f76d7c7512">00055</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g3c1a2b0972c6a8e1215eb3f76d7c7512" title="Length of vector.">size</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v) { <span class="keywordflow">return</span> v.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>(); }
86<a name="l00056"></a>00056
87<a name="l00058"></a>00058 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
88<a name="l00059"></a><a class="code" href="group__matrix__functions.html#g8b3dcc61484fce58ee4244d5731e18fa">00059</a> T <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v)
89<a name="l00060"></a>00060 {
90<a name="l00061"></a>00061   T M = 0;
91<a name="l00062"></a>00062
92<a name="l00063"></a>00063   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0;i &lt; v.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>();i++)
93<a name="l00064"></a>00064     M += v[i];
94<a name="l00065"></a>00065
95<a name="l00066"></a>00066   <span class="keywordflow">return</span> M;
96<a name="l00067"></a>00067 }
97<a name="l00068"></a>00068
98<a name="l00076"></a>00076 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
99<a name="l00077"></a><a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355">00077</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> dim = 1)
100<a name="l00078"></a>00078 {
101<a name="l00079"></a>00079   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((dim == 1) || (dim == 2), <span class="stringliteral">"sum: dimension need to be 1 or 2"</span>);
102<a name="l00080"></a>00080   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> out;
103<a name="l00081"></a>00081
104<a name="l00082"></a>00082   <span class="keywordflow">if</span> (dim == 1) {
105<a name="l00083"></a>00083     out.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), <span class="keyword">false</span>);
106<a name="l00084"></a>00084
107<a name="l00085"></a>00085     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); i++)
108<a name="l00086"></a>00086       out(i) = <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(m.<a class="code" href="classitpp_1_1Mat.html#095a4e54f39932b62633f7f4f252ab46" title="Get column c.">get_col</a>(i));
109<a name="l00087"></a>00087   }
110<a name="l00088"></a>00088   <span class="keywordflow">else</span> {
111<a name="l00089"></a>00089     out.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), <span class="keyword">false</span>);
112<a name="l00090"></a>00090
113<a name="l00091"></a>00091     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++)
114<a name="l00092"></a>00092       out(i) = <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(m.<a class="code" href="classitpp_1_1Mat.html#9e840908b256dd1353845356b9ffacb2" title="Get row r.">get_row</a>(i));
115<a name="l00093"></a>00093   }
116<a name="l00094"></a>00094
117<a name="l00095"></a>00095   <span class="keywordflow">return</span> out;
118<a name="l00096"></a>00096 }
119<a name="l00097"></a>00097
120<a name="l00098"></a>00098
121<a name="l00100"></a>00100 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
122<a name="l00101"></a><a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7">00101</a> T <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;X)
123<a name="l00102"></a>00102 {
124<a name="l00103"></a>00103   <span class="keyword">const</span> T * X_data = X.<a class="code" href="classitpp_1_1Mat.html#e88c1da88027f1e597cb5eb34823b56c" title="Access of the internal data structure (not recommended to use).">_data</a>();
125<a name="l00104"></a>00104   <span class="keyword">const</span> <span class="keywordtype">int</span> X_datasize = X.<a class="code" href="classitpp_1_1Mat.html#20119f290acfc60a8025a853ea00d366" title="Access to the internal data structure (not recommended to use).">_datasize</a>();
126<a name="l00105"></a>00105   T acc = 0;
127<a name="l00106"></a>00106
128<a name="l00107"></a>00107   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0;i &lt; X_datasize;i++)
129<a name="l00108"></a>00108     acc += X_data[i];
130<a name="l00109"></a>00109
131<a name="l00110"></a>00110   <span class="keywordflow">return</span> acc;
132<a name="l00111"></a>00111 }
133<a name="l00112"></a>00112
134<a name="l00113"></a>00113
135<a name="l00115"></a>00115 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
136<a name="l00116"></a><a class="code" href="group__matrix__functions.html#gc6a1b0e0af40f95bc636ce039c8f1a0f">00116</a> T <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v)
137<a name="l00117"></a>00117 {
138<a name="l00118"></a>00118   T M = 0;
139<a name="l00119"></a>00119
140<a name="l00120"></a>00120   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; v.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>(); i++)
141<a name="l00121"></a>00121     M += v[i] * v[i];
142<a name="l00122"></a>00122
143<a name="l00123"></a>00123   <span class="keywordflow">return</span> M;
144<a name="l00124"></a>00124 }
145<a name="l00125"></a>00125
146<a name="l00133"></a>00133 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
147<a name="l00134"></a><a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6">00134</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> dim = 1)
148<a name="l00135"></a>00135 {
149<a name="l00136"></a>00136   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((dim == 1) || (dim == 2), <span class="stringliteral">"sum_sqr: dimension need to be 1 or 2"</span>);
150<a name="l00137"></a>00137   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> out;
151<a name="l00138"></a>00138
152<a name="l00139"></a>00139   <span class="keywordflow">if</span> (dim == 1) {
153<a name="l00140"></a>00140     out.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), <span class="keyword">false</span>);
154<a name="l00141"></a>00141
155<a name="l00142"></a>00142     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); i++)
156<a name="l00143"></a>00143       out(i) = <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(m.<a class="code" href="classitpp_1_1Mat.html#095a4e54f39932b62633f7f4f252ab46" title="Get column c.">get_col</a>(i));
157<a name="l00144"></a>00144   }
158<a name="l00145"></a>00145   <span class="keywordflow">else</span> {
159<a name="l00146"></a>00146     out.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), <span class="keyword">false</span>);
160<a name="l00147"></a>00147
161<a name="l00148"></a>00148     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++)
162<a name="l00149"></a>00149       out(i) = <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(m.<a class="code" href="classitpp_1_1Mat.html#9e840908b256dd1353845356b9ffacb2" title="Get row r.">get_row</a>(i));
163<a name="l00150"></a>00150   }
164<a name="l00151"></a>00151
165<a name="l00152"></a>00152   <span class="keywordflow">return</span> out;
166<a name="l00153"></a>00153 }
167<a name="l00154"></a>00154
168<a name="l00156"></a>00156 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
169<a name="l00157"></a><a class="code" href="group__matrix__functions.html#g1609c616b581628be87a4faca7778f54">00157</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v)
170<a name="l00158"></a>00158 {
171<a name="l00159"></a>00159   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> out(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>());
172<a name="l00160"></a>00160
173<a name="l00161"></a>00161   out(0) = v(0);
174<a name="l00162"></a>00162   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 1; i &lt; v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(); i++)
175<a name="l00163"></a>00163     out(i) = out(i - 1) + v(i);
176<a name="l00164"></a>00164
177<a name="l00165"></a>00165   <span class="keywordflow">return</span> out;
178<a name="l00166"></a>00166 }
179<a name="l00167"></a>00167
180<a name="l00175"></a>00175 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
181<a name="l00176"></a><a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655">00176</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> dim = 1)
182<a name="l00177"></a>00177 {
183<a name="l00178"></a>00178   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((dim == 1) || (dim == 2), <span class="stringliteral">"cumsum: dimension need to be 1 or 2"</span>);
184<a name="l00179"></a>00179   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> out(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>());
185<a name="l00180"></a>00180
186<a name="l00181"></a>00181   <span class="keywordflow">if</span> (dim == 1) {
187<a name="l00182"></a>00182     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); i++)
188<a name="l00183"></a>00183       out.set_col(i, <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(m.<a class="code" href="classitpp_1_1Mat.html#095a4e54f39932b62633f7f4f252ab46" title="Get column c.">get_col</a>(i)));
189<a name="l00184"></a>00184   }
190<a name="l00185"></a>00185   <span class="keywordflow">else</span> {
191<a name="l00186"></a>00186     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++)
192<a name="l00187"></a>00187       out.set_row(i, <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(m.<a class="code" href="classitpp_1_1Mat.html#9e840908b256dd1353845356b9ffacb2" title="Get row r.">get_row</a>(i)));
193<a name="l00188"></a>00188   }
194<a name="l00189"></a>00189
195<a name="l00190"></a>00190   <span class="keywordflow">return</span> out;
196<a name="l00191"></a>00191 }
197<a name="l00192"></a>00192
198<a name="l00194"></a>00194 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
199<a name="l00195"></a><a class="code" href="group__matrix__functions.html#gff8a3812cb84a18f3b52c4bb116feb08">00195</a> T <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v)
200<a name="l00196"></a>00196 {
201<a name="l00197"></a>00197   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() &gt;= 1, <span class="stringliteral">"prod: size of vector should be at least 1"</span>);
202<a name="l00198"></a>00198   T out = v(0);
203<a name="l00199"></a>00199
204<a name="l00200"></a>00200   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 1; i &lt; v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(); i++)
205<a name="l00201"></a>00201     out *= v(i);
206<a name="l00202"></a>00202
207<a name="l00203"></a>00203   <span class="keywordflow">return</span> out;
208<a name="l00204"></a>00204 }
209<a name="l00205"></a>00205
210<a name="l00213"></a>00213 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
211<a name="l00214"></a><a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2">00214</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> dim = 1)
212<a name="l00215"></a>00215 {
213<a name="l00216"></a>00216   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((dim == 1) || (dim == 2), <span class="stringliteral">"prod: dimension need to be 1 or 2"</span>);
214<a name="l00217"></a>00217   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> out(m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>());
215<a name="l00218"></a>00218
216<a name="l00219"></a>00219   <span class="keywordflow">if</span> (dim == 1) {
217<a name="l00220"></a>00220     <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>() &gt;= 1) &amp;&amp; (m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() &gt;= 1),
218<a name="l00221"></a>00221               <span class="stringliteral">"prod: number of columns should be at least 1"</span>);
219<a name="l00222"></a>00222     out.set_size(m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), <span class="keyword">false</span>);
220<a name="l00223"></a>00223
221<a name="l00224"></a>00224     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); i++)
222<a name="l00225"></a>00225       out(i) = <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(m.<a class="code" href="classitpp_1_1Mat.html#095a4e54f39932b62633f7f4f252ab46" title="Get column c.">get_col</a>(i));
223<a name="l00226"></a>00226   }
224<a name="l00227"></a>00227   <span class="keywordflow">else</span> {
225<a name="l00228"></a>00228     <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>() &gt;= 1) &amp;&amp; (m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() &gt;= 1),
226<a name="l00229"></a>00229               <span class="stringliteral">"prod: number of rows should be at least 1"</span>);
227<a name="l00230"></a>00230     out.set_size(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), <span class="keyword">false</span>);
228<a name="l00231"></a>00231
229<a name="l00232"></a>00232     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++)
230<a name="l00233"></a>00233       out(i) = <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(m.<a class="code" href="classitpp_1_1Mat.html#9e840908b256dd1353845356b9ffacb2" title="Get row r.">get_row</a>(i));
231<a name="l00234"></a>00234   }
232<a name="l00235"></a>00235   <span class="keywordflow">return</span> out;
233<a name="l00236"></a>00236 }
234<a name="l00237"></a>00237
235<a name="l00239"></a>00239 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
236<a name="l00240"></a><a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e">00240</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v1, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v2)
237<a name="l00241"></a>00241 {
238<a name="l00242"></a>00242   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((v1.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == 3) &amp;&amp; (v2.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == 3),
239<a name="l00243"></a>00243             <span class="stringliteral">"cross: vectors should be of size 3"</span>);
240<a name="l00244"></a>00244
241<a name="l00245"></a>00245   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> r(3);
242<a name="l00246"></a>00246
243<a name="l00247"></a>00247   r(0) = v1(1) * v2(2) - v1(2) * v2(1);
244<a name="l00248"></a>00248   r(1) = v1(2) * v2(0) - v1(0) * v2(2);
245<a name="l00249"></a>00249   r(2) = v1(0) * v2(1) - v1(1) * v2(0);
246<a name="l00250"></a>00250
247<a name="l00251"></a>00251   <span class="keywordflow">return</span> r;
248<a name="l00252"></a>00252 }
249<a name="l00253"></a>00253
250<a name="l00254"></a>00254
251<a name="l00256"></a>00256 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
252<a name="l00257"></a><a class="code" href="group__matrix__functions.html#gb42e2269c4fb80730ac82ef7c85f9cd0">00257</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keywordtype">int</span> n)
253<a name="l00258"></a>00258 {
254<a name="l00259"></a>00259   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(n &gt;= v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(), <span class="stringliteral">"zero_pad() cannot shrink the vector!"</span>);
255<a name="l00260"></a>00260   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> v2(n);
256<a name="l00261"></a>00261   v2.set_subvector(0, v);
257<a name="l00262"></a>00262   <span class="keywordflow">if</span> (n &gt; v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>())
258<a name="l00263"></a>00263     v2.set_subvector(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(), n - 1, T(0));
259<a name="l00264"></a>00264
260<a name="l00265"></a>00265   <span class="keywordflow">return</span> v2;
261<a name="l00266"></a>00266 }
262<a name="l00267"></a>00267
263<a name="l00269"></a>00269 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
264<a name="l00270"></a><a class="code" href="group__matrix__functions.html#ga92be4bb938ac8c78f6313432210599c">00270</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v)
265<a name="l00271"></a>00271 {
266<a name="l00272"></a>00272   <span class="keywordtype">int</span> n = <a class="code" href="group__logexpfunc.html#g44da7c49dc0f990f661404a47e2bf9be" title="Calculate two to the power of x (2^x); x is integer.">pow2i</a>(<a class="code" href="group__logexpfunc.html#g4a774bd5418520f574eddf730333dada" title="Calculate the number of bits needed to represent a numer of levels saved in a vector...">levels2bits</a>(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>()));
267<a name="l00273"></a>00273
268<a name="l00274"></a>00274   <span class="keywordflow">return</span> (n == v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>()) ? v : <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(v, n);
269<a name="l00275"></a>00275 }
270<a name="l00276"></a>00276
271<a name="l00278"></a>00278 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
272<a name="l00279"></a><a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea">00279</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols)
273<a name="l00280"></a>00280 {
274<a name="l00281"></a>00281   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((rows &gt;= m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>()) &amp;&amp; (cols &gt;= m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>()),
275<a name="l00282"></a>00282             <span class="stringliteral">"zero_pad() cannot shrink the matrix!"</span>);
276<a name="l00283"></a>00283   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> m2(rows, cols);
277<a name="l00284"></a>00284   m2.<a class="code" href="classitpp_1_1Mat.html#73379f7980314d0dd64874fc4a92e8b1" title="This function is deprecated. Please use set_submatrix(int r, int c, const Mat&amp;lt;&amp;gt;...">set_submatrix</a>(0, 0, m);
278<a name="l00285"></a>00285   <span class="keywordflow">if</span> (cols &gt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>()) <span class="comment">// Zero</span>
279<a name="l00286"></a>00286     m2.<a class="code" href="classitpp_1_1Mat.html#73379f7980314d0dd64874fc4a92e8b1" title="This function is deprecated. Please use set_submatrix(int r, int c, const Mat&amp;lt;&amp;gt;...">set_submatrix</a>(0, m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() - 1, m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), cols - 1, T(0));
280<a name="l00287"></a>00287   <span class="keywordflow">if</span> (rows &gt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>()) <span class="comment">// Zero</span>
281<a name="l00288"></a>00288     m2.<a class="code" href="classitpp_1_1Mat.html#73379f7980314d0dd64874fc4a92e8b1" title="This function is deprecated. Please use set_submatrix(int r, int c, const Mat&amp;lt;&amp;gt;...">set_submatrix</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), rows - 1, 0, cols - 1, T(0));
282<a name="l00289"></a>00289
283<a name="l00290"></a>00290   <span class="keywordflow">return</span> m2;
284<a name="l00291"></a>00291 }
285<a name="l00292"></a>00292
286<a name="l00293"></a>00293
287<a name="l00296"></a>00296 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
288<a name="l00297"></a><a class="code" href="group__matrix__functions.html#gab6a4b173411de9290f0dc401c9dcec7">00297</a> T <a class="code" href="group__matrix__functions.html#gab6a4b173411de9290f0dc401c9dcec7">index_zero_pad</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keyword">const</span> <span class="keywordtype">int</span> index)
289<a name="l00298"></a>00298 {
290<a name="l00299"></a>00299   <span class="keywordflow">if</span> (index &gt;= 0 &amp;&amp; index &lt; v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>())
291<a name="l00300"></a>00300     <span class="keywordflow">return</span> v(index);
292<a name="l00301"></a>00301   <span class="keywordflow">else</span>
293<a name="l00302"></a>00302     <span class="keywordflow">return</span> T(0);
294<a name="l00303"></a>00303 }
295<a name="l00304"></a>00304
296<a name="l00305"></a>00305
297<a name="l00307"></a>00307 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
298<a name="l00308"></a><a class="code" href="group__matrix__functions.html#g22c587fbc615e295f2b881c2e5025006">00308</a> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g22c587fbc615e295f2b881c2e5025006" title="Transposition of the matrix m returning the transposed matrix in out.">transpose</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;out) { out = m.<a class="code" href="classitpp_1_1Mat.html#3be8603ee4668eb467860875a1c5462c" title="Matrix transpose.">T</a>(); }
299<a name="l00309"></a>00309
300<a name="l00311"></a>00311 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
301<a name="l00312"></a><a class="code" href="group__matrix__functions.html#g8c724070b69eb97eae2985e02088feca">00312</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g22c587fbc615e295f2b881c2e5025006" title="Transposition of the matrix m returning the transposed matrix in out.">transpose</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m) { <span class="keywordflow">return</span> m.<a class="code" href="classitpp_1_1Mat.html#3be8603ee4668eb467860875a1c5462c" title="Matrix transpose.">T</a>(); }
302<a name="l00313"></a>00313
303<a name="l00314"></a>00314
304<a name="l00317"></a>00317 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
305<a name="l00318"></a><a class="code" href="group__matrix__functions.html#ge71aa93cc176699065fc0e740a50387e">00318</a> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;out) { out = m.<a class="code" href="classitpp_1_1Mat.html#38c68bb449851504a4d5fcaa96fce1ad" title="Hermitian matrix transpose (conjugate transpose).">H</a>(); }
306<a name="l00319"></a>00319
307<a name="l00321"></a>00321 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
308<a name="l00322"></a><a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579">00322</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m) { <span class="keywordflow">return</span> m.<a class="code" href="classitpp_1_1Mat.html#38c68bb449851504a4d5fcaa96fce1ad" title="Hermitian matrix transpose (conjugate transpose).">H</a>(); }
309<a name="l00323"></a>00323
310<a name="l00324"></a>00324
311<a name="l00325"></a>00325
312<a name="l00335"></a>00335 <span class="keyword">template</span>&lt;<span class="keyword">class</span> Num_T&gt;
313<a name="l00336"></a><a class="code" href="group__matrix__functions.html#g49691bf233320a2dd8ae9929988f8f91">00336</a> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g49691bf233320a2dd8ae9929988f8f91" title="Returns true if matrix X is hermitian, false otherwise.">is_hermitian</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a>&amp; X)
314<a name="l00337"></a>00337 {
315<a name="l00338"></a>00338
316<a name="l00339"></a>00339   <span class="keywordflow">if</span> (X == X.<a class="code" href="classitpp_1_1Mat.html#38c68bb449851504a4d5fcaa96fce1ad" title="Hermitian matrix transpose (conjugate transpose).">H</a>())
317<a name="l00340"></a>00340     <span class="keywordflow">return</span> <span class="keyword">true</span>;
318<a name="l00341"></a>00341   <span class="keywordflow">else</span>
319<a name="l00342"></a>00342     <span class="keywordflow">return</span> <span class="keyword">false</span>;
320<a name="l00343"></a>00343 }
321<a name="l00344"></a>00344
322<a name="l00354"></a>00354 <span class="keyword">template</span>&lt;<span class="keyword">class</span> Num_T&gt;
323<a name="l00355"></a><a class="code" href="group__matrix__functions.html#g0e711a4a8db7d9a2ed5515a33dd65d4c">00355</a> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g0e711a4a8db7d9a2ed5515a33dd65d4c" title="Returns true if matrix X is unitary, false otherwise.">is_unitary</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a>&amp; X)
324<a name="l00356"></a>00356 {
325<a name="l00357"></a>00357
326<a name="l00358"></a>00358   <span class="keywordflow">if</span> (<a class="code" href="group__inverse.html#g125b083397cc9450b8a67a485b545bc8" title="Inverse of real square matrix.Calculate the inverse of the real matrix .">inv</a>(X) == X.<a class="code" href="classitpp_1_1Mat.html#38c68bb449851504a4d5fcaa96fce1ad" title="Hermitian matrix transpose (conjugate transpose).">H</a>())
327<a name="l00359"></a>00359     <span class="keywordflow">return</span> <span class="keyword">true</span>;
328<a name="l00360"></a>00360   <span class="keywordflow">else</span>
329<a name="l00361"></a>00361     <span class="keywordflow">return</span> <span class="keyword">false</span>;
330<a name="l00362"></a>00362 }
331<a name="l00363"></a>00363
332<a name="l00364"></a>00364
333<a name="l00373"></a>00373 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
334<a name="l00374"></a><a class="code" href="group__matrix__functions.html#g53cd7ee1e173c8a2e56d0b4c77e49c20">00374</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> repmat(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keywordtype">int</span> n)
335<a name="l00375"></a>00375 {
336<a name="l00376"></a>00376   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(n &gt; 0, <span class="stringliteral">"repmat(): Wrong repetition parameter"</span>);
337<a name="l00377"></a>00377   <span class="keywordtype">int</span> data_length = v.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>();
338<a name="l00378"></a>00378   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(data_length &gt; 0, <span class="stringliteral">"repmat(): Input vector can not be empty"</span>);
339<a name="l00379"></a>00379   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> assembly(data_length * n);
340<a name="l00380"></a>00380   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j = 0; j &lt; n; ++j) {
341<a name="l00381"></a>00381     assembly.<a class="code" href="classitpp_1_1Vec.html#fc812b3402a3d8d76ac670818cda5dc9" title="This function is deprecated. Please use set_subvector(i, v) instead.">set_subvector</a>(j * data_length, v);
342<a name="l00382"></a>00382   }
343<a name="l00383"></a>00383   <span class="keywordflow">return</span> assembly;
344<a name="l00384"></a>00384 }
345<a name="l00385"></a>00385
346<a name="l00386"></a>00386
347<a name="l00396"></a>00396 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
348<a name="l00397"></a><a class="code" href="group__matrix__functions.html#g608ecbac6ff6a1ce6d5e69c0dd789021">00397</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> repmat(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n)
349<a name="l00398"></a>00398 {
350<a name="l00399"></a>00399   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((m &gt; 0) &amp;&amp; (n &gt; 0), <span class="stringliteral">"repmat(): Wrong repetition parameters"</span>);
351<a name="l00400"></a>00400   <span class="keywordtype">int</span> data_rows = data.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>();
352<a name="l00401"></a>00401   <span class="keywordtype">int</span> data_cols = data.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
353<a name="l00402"></a>00402   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>((data_rows &gt; 0) &amp;&amp; (data_cols &gt; 0), <span class="stringliteral">"repmat(): Input matrix can "</span>
354<a name="l00403"></a>00403             <span class="stringliteral">"not be empty"</span>);
355<a name="l00404"></a>00404   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> assembly(data_rows*m, data_cols*n);
356<a name="l00405"></a>00405   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; m; ++i) {
357<a name="l00406"></a>00406     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j = 0; j &lt; n; ++j) {
358<a name="l00407"></a>00407       assembly.<a class="code" href="classitpp_1_1Mat.html#73379f7980314d0dd64874fc4a92e8b1" title="This function is deprecated. Please use set_submatrix(int r, int c, const Mat&amp;lt;&amp;gt;...">set_submatrix</a>(i*data_rows, j*data_cols, data);
359<a name="l00408"></a>00408     }
360<a name="l00409"></a>00409   }
361<a name="l00410"></a>00410   <span class="keywordflow">return</span> assembly;
362<a name="l00411"></a>00411 }
363<a name="l00412"></a>00412
364<a name="l00424"></a>00424 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt; <span class="keyword">inline</span>
365<a name="l00425"></a><a class="code" href="group__matrix__functions.html#g8d9541af4bd11193689d3f30b1cd5a74">00425</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> repmat(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose = <span class="keyword">false</span>)
366<a name="l00426"></a>00426 {
367<a name="l00427"></a>00427   <span class="keywordflow">return</span> repmat((<a class="code" href="group__matrix__functions.html#g22c587fbc615e295f2b881c2e5025006" title="Transposition of the matrix m returning the transposed matrix in out.">transpose</a> ? v.<a class="code" href="classitpp_1_1Vec.html#9931ef6889bfd119946a3af4740f0cc9" title="Matrix transpose. Converts to a matrix with the vector in the first row.">T</a>() : <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a>(v)), m, n);
368<a name="l00428"></a>00428 }
369<a name="l00429"></a>00429
370<a name="l00430"></a>00430
371<a name="l00442"></a>00442 <span class="keyword">template</span>&lt;<span class="keyword">class</span> Num_T&gt;
372<a name="l00443"></a><a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86">00443</a> <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a> <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a>&amp; X, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a>&amp; Y)
373<a name="l00444"></a>00444 {
374<a name="l00445"></a>00445   <a class="code" href="classitpp_1_1Mat.html" title="Matrix Class (Templated).">Mat&lt;Num_T&gt;</a> result(X.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() * Y.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), X.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>() * Y.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>());
375<a name="l00446"></a>00446
376<a name="l00447"></a>00447   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; X.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++)
377<a name="l00448"></a>00448     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j = 0; j &lt; X.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); j++)
378<a name="l00449"></a>00449       result.set_submatrix(i * Y.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), j * Y.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), X(i, j) * Y);
379<a name="l00450"></a>00450
380<a name="l00451"></a>00451   <span class="keywordflow">return</span> result;
381<a name="l00452"></a>00452 }
382<a name="l00453"></a>00453
383<a name="l00454"></a>00454
384<a name="l00467"></a>00467 cmat <a class="code" href="group__matrix__functions.html#gbb4749a685e730a2dc0ff55313d13a22" title="Square root of the real square matrix A.">sqrtm</a>(<span class="keyword">const</span> cmat&amp; A);
385<a name="l00468"></a>00468
386<a name="l00481"></a>00481 cmat <a class="code" href="group__matrix__functions.html#gbb4749a685e730a2dc0ff55313d13a22" title="Square root of the real square matrix A.">sqrtm</a>(<span class="keyword">const</span> mat&amp; A);
387<a name="l00482"></a>00482
388<a name="l00483"></a>00483
389<a name="l00492"></a>00492 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
390<a name="l00493"></a><a class="code" href="group__matrix__functions.html#g472a6d548c6d195bb349c812c8b2203c">00493</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">double</span> tol = -1.0)
391<a name="l00494"></a>00494 {
392<a name="l00495"></a>00495   <span class="keywordtype">int</span> rows = m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>();
393<a name="l00496"></a>00496   <span class="keywordtype">int</span> cols = m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
394<a name="l00497"></a>00497   <span class="keywordflow">if</span> ((rows == 0) || (cols == 0))
395<a name="l00498"></a>00498     <span class="keywordflow">return</span> 0;
396<a name="l00499"></a>00499
397<a name="l00500"></a>00500   vec sing_val = <a class="code" href="group__matrixdecomp.html#ga53763a60632139d71957ea7cf0080c9" title="Return singular values of a complex matrix A using SVD.">svd</a>(m);
398<a name="l00501"></a>00501
399<a name="l00502"></a>00502   <span class="keywordflow">if</span> (tol &lt; 0.0) { <span class="comment">// Calculate default tolerance</span>
400<a name="l00503"></a>00503     tol = eps * sing_val(0) * (rows &gt; cols ? rows : cols);
401<a name="l00504"></a>00504   }
402<a name="l00505"></a>00505
403<a name="l00506"></a>00506   <span class="comment">// Count number of nonzero singular values</span>
404<a name="l00507"></a>00507   <span class="keywordtype">int</span> r = 0;
405<a name="l00508"></a>00508   <span class="keywordflow">while</span> ((r &lt; sing_val.length()) &amp;&amp; (sing_val(r) &gt; tol)) {
406<a name="l00509"></a>00509     r++;
407<a name="l00510"></a>00510   }
408<a name="l00511"></a>00511
409<a name="l00512"></a>00512   <span class="keywordflow">return</span> r;
410<a name="l00513"></a>00513 }
411<a name="l00514"></a>00514
412<a name="l00516"></a>00516 <span class="keyword">template</span>&lt;&gt; <span class="keyword">inline</span>
413<a name="l00517"></a><a class="code" href="group__matrix__functions.html#g36e27d8fdd14e8d5b4a0653c101f1570">00517</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">double</span> tol)
414<a name="l00518"></a>00518 {
415<a name="l00519"></a>00519   <span class="keywordflow">return</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<a class="code" href="group__fixed.html#g36d7bf8746a81d073e4db8b0ad3bec86" title="Converts a fixmat to mat.">to_mat</a>(m), tol);
416<a name="l00520"></a>00520 }
417<a name="l00521"></a>00521
418<a name="l00523"></a>00523 <span class="keyword">template</span>&lt;&gt; <span class="keyword">inline</span>
419<a name="l00524"></a><a class="code" href="group__matrix__functions.html#g82519e6233d0082856d4acbf84353160">00524</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">double</span> tol)
420<a name="l00525"></a>00525 {
421<a name="l00526"></a>00526   <span class="keywordflow">return</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<a class="code" href="group__fixed.html#g36d7bf8746a81d073e4db8b0ad3bec86" title="Converts a fixmat to mat.">to_mat</a>(m), tol);
422<a name="l00527"></a>00527 }
423<a name="l00528"></a>00528
424<a name="l00530"></a>00530 <span class="keyword">template</span>&lt;&gt; <span class="keyword">inline</span>
425<a name="l00531"></a><a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d">00531</a> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g6adab7c28deb3fec407684458f0b719d" title="Specialisation of rank() function.">rank</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;, <span class="keywordtype">double</span>)
426<a name="l00532"></a>00532 {
427<a name="l00533"></a>00533   <a class="code" href="group__errorhandlingfunc.html#g22d38e98332f9edff88cc501463eedce" title="Abort unconditionally.">it_error</a>(<span class="stringliteral">"rank(bmat): Function not implemented for GF(2) algebra"</span>);
428<a name="l00534"></a>00534   <span class="keywordflow">return</span> 0;
429<a name="l00535"></a>00535 }
430<a name="l00536"></a>00536
431<a name="l00538"></a>00538
432<a name="l00539"></a>00539
433<a name="l00540"></a>00540
434<a name="l00541"></a>00541 <span class="comment">// -------------------- Diagonal matrix functions -------------------------</span>
435<a name="l00542"></a>00542
436<a name="l00545"></a>00545
437<a name="l00556"></a>00556 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
438<a name="l00557"></a><a class="code" href="group__diag.html#g1a5019e9703cd9f87c3914ab298714ff">00557</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keyword">const</span> <span class="keywordtype">int</span> K = 0)
439<a name="l00558"></a>00558 {
440<a name="l00559"></a>00559   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> m(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + <a class="code" href="group__fixed.html#ge0a38d51b2bf99a0ba07f2f8885e6adc" title="Absolute value.">std::abs</a>(K), v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + <a class="code" href="group__fixed.html#ge0a38d51b2bf99a0ba07f2f8885e6adc" title="Absolute value.">std::abs</a>(K));
441<a name="l00560"></a>00560   m = T(0);
442<a name="l00561"></a>00561   <span class="keywordflow">if</span> (K &gt; 0)
443<a name="l00562"></a>00562     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() - 1; i &gt;= 0; i--)
444<a name="l00563"></a>00563       m(i, i + K) = v(i);
445<a name="l00564"></a>00564   <span class="keywordflow">else</span>
446<a name="l00565"></a>00565     <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() - 1; i &gt;= 0; i--)
447<a name="l00566"></a>00566       m(i - K, i) = v(i);
448<a name="l00567"></a>00567
449<a name="l00568"></a>00568   <span class="keywordflow">return</span> m;
450<a name="l00569"></a>00569 }
451<a name="l00570"></a>00570
452<a name="l00580"></a>00580 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
453<a name="l00581"></a><a class="code" href="group__diag.html#g853c555c28b2c812f39962d272ae33ac">00581</a> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
454<a name="l00582"></a>00582 {
455<a name="l00583"></a>00583   m.<a class="code" href="classitpp_1_1Mat.html#1136d376109f2af8843593bccf9de88f" title="Set size of matrix. If copy = true then keep the data before resizing.">set_size</a>(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(), v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(), <span class="keyword">false</span>);
456<a name="l00584"></a>00584   m = T(0);
457<a name="l00585"></a>00585   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() - 1; i &gt;= 0; i--)
458<a name="l00586"></a>00586     m(i, i) = v(i);
459<a name="l00587"></a>00587 }
460<a name="l00588"></a>00588
461<a name="l00596"></a>00596 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
462<a name="l00597"></a><a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553">00597</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
463<a name="l00598"></a>00598 {
464<a name="l00599"></a>00599   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> t(<a class="code" href="group__protocol.html#g27dde5be5aac902c0de838b4b872cd2f" title="ADD DOCUMENTATION HERE.">std::min</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>()));
465<a name="l00600"></a>00600
466<a name="l00601"></a>00601   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; t.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>(); i++)
467<a name="l00602"></a>00602     t(i) = m(i, i);
468<a name="l00603"></a>00603
469<a name="l00604"></a>00604   <span class="keywordflow">return</span> t;
470<a name="l00605"></a>00605 }
471<a name="l00606"></a>00606
472<a name="l00616"></a>00616 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
473<a name="l00617"></a><a class="code" href="group__diag.html#g9f2b6719c7de2f116bc794e406c952cc">00617</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup)
474<a name="l00618"></a>00618 {
475<a name="l00619"></a>00619   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sup.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1, <span class="stringliteral">"bidiag()"</span>);
476<a name="l00620"></a>00620
477<a name="l00621"></a>00621   <span class="keywordtype">int</span> n = main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>();
478<a name="l00622"></a>00622   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> m(n, n);
479<a name="l00623"></a>00623   m = T(0);
480<a name="l00624"></a>00624   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
481<a name="l00625"></a>00625     m(i, i) = main(i);
482<a name="l00626"></a>00626     m(i, i + 1) = sup(i);
483<a name="l00627"></a>00627   }
484<a name="l00628"></a>00628   m(n - 1, n - 1) = main(n - 1);
485<a name="l00629"></a>00629
486<a name="l00630"></a>00630   <span class="keywordflow">return</span> m;
487<a name="l00631"></a>00631 }
488<a name="l00632"></a>00632
489<a name="l00642"></a>00642 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
490<a name="l00643"></a><a class="code" href="group__diag.html#g2d7f36dc45fb3098131111d172f5b08d">00643</a> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup, <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
491<a name="l00644"></a>00644 {
492<a name="l00645"></a>00645   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sup.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1, <span class="stringliteral">"bidiag()"</span>);
493<a name="l00646"></a>00646
494<a name="l00647"></a>00647   <span class="keywordtype">int</span> n = main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>();
495<a name="l00648"></a>00648   m.<a class="code" href="classitpp_1_1Mat.html#1136d376109f2af8843593bccf9de88f" title="Set size of matrix. If copy = true then keep the data before resizing.">set_size</a>(n, n);
496<a name="l00649"></a>00649   m = T(0);
497<a name="l00650"></a>00650   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
498<a name="l00651"></a>00651     m(i, i) = main(i);
499<a name="l00652"></a>00652     m(i, i + 1) = sup(i);
500<a name="l00653"></a>00653   }
501<a name="l00654"></a>00654   m(n - 1, n - 1) = main(n - 1);
502<a name="l00655"></a>00655 }
503<a name="l00656"></a>00656
504<a name="l00665"></a>00665 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
505<a name="l00666"></a><a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748">00666</a> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup)
506<a name="l00667"></a>00667 {
507<a name="l00668"></a>00668   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() == m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), <span class="stringliteral">"bidiag(): Matrix must be square!"</span>);
508<a name="l00669"></a>00669
509<a name="l00670"></a>00670   <span class="keywordtype">int</span> n = m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
510<a name="l00671"></a>00671   main.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(n);
511<a name="l00672"></a>00672   sup.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(n - 1);
512<a name="l00673"></a>00673   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
513<a name="l00674"></a>00674     main(i) = m(i, i);
514<a name="l00675"></a>00675     sup(i) = m(i, i + 1);
515<a name="l00676"></a>00676   }
516<a name="l00677"></a>00677   main(n - 1) = m(n - 1, n - 1);
517<a name="l00678"></a>00678 }
518<a name="l00679"></a>00679
519<a name="l00689"></a>00689 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
520<a name="l00690"></a><a class="code" href="group__diag.html#g0ff63d5f0db74b830fa918a1866c440b">00690</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sub)
521<a name="l00691"></a>00691 {
522<a name="l00692"></a>00692   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sup.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1 &amp;&amp; main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sub.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1, <span class="stringliteral">"bidiag()"</span>);
523<a name="l00693"></a>00693
524<a name="l00694"></a>00694   <span class="keywordtype">int</span> n = main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>();
525<a name="l00695"></a>00695   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> m(n, n);
526<a name="l00696"></a>00696   m = T(0);
527<a name="l00697"></a>00697   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
528<a name="l00698"></a>00698     m(i, i) = main(i);
529<a name="l00699"></a>00699     m(i, i + 1) = sup(i);
530<a name="l00700"></a>00700     m(i + 1, i) = sub(i);
531<a name="l00701"></a>00701   }
532<a name="l00702"></a>00702   m(n - 1, n - 1) = main(n - 1);
533<a name="l00703"></a>00703
534<a name="l00704"></a>00704   <span class="keywordflow">return</span> m;
535<a name="l00705"></a>00705 }
536<a name="l00706"></a>00706
537<a name="l00716"></a>00716 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
538<a name="l00717"></a><a class="code" href="group__diag.html#geb2b140f8d64b19d7a4462464c59cdca">00717</a> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup, <span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sub, <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
539<a name="l00718"></a>00718 {
540<a name="l00719"></a>00719   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sup.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1 &amp;&amp; main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == sub.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() + 1, <span class="stringliteral">"bidiag()"</span>);
541<a name="l00720"></a>00720
542<a name="l00721"></a>00721   <span class="keywordtype">int</span> n = main.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>();
543<a name="l00722"></a>00722   m.<a class="code" href="classitpp_1_1Mat.html#1136d376109f2af8843593bccf9de88f" title="Set size of matrix. If copy = true then keep the data before resizing.">set_size</a>(n, n);
544<a name="l00723"></a>00723   m = T(0);
545<a name="l00724"></a>00724   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
546<a name="l00725"></a>00725     m(i, i) = main(i);
547<a name="l00726"></a>00726     m(i, i + 1) = sup(i);
548<a name="l00727"></a>00727     m(i + 1, i) = sub(i);
549<a name="l00728"></a>00728   }
550<a name="l00729"></a>00729   m(n - 1, n - 1) = main(n - 1);
551<a name="l00730"></a>00730 }
552<a name="l00731"></a>00731
553<a name="l00740"></a>00740 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
554<a name="l00741"></a><a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e">00741</a> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;main, <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sup, <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;sub)
555<a name="l00742"></a>00742 {
556<a name="l00743"></a>00743   <a class="code" href="group__errorhandlingfunc.html#gd5c34b291e5018534fd2344486e2b5a1" title="Abort if t is not true.">it_assert</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>() == m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(), <span class="stringliteral">"tridiag(): Matrix must be square!"</span>);
557<a name="l00744"></a>00744
558<a name="l00745"></a>00745   <span class="keywordtype">int</span> n = m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
559<a name="l00746"></a>00746   main.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(n);
560<a name="l00747"></a>00747   sup.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(n - 1);
561<a name="l00748"></a>00748   sub.<a class="code" href="classitpp_1_1Vec.html#a4b02b7819e540bb3853c15e6480ada6" title="Set length of vector. if copy = true then keeping the old values.">set_size</a>(n - 1);
562<a name="l00749"></a>00749   <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i = 0; i &lt; n - 1; i++) {
563<a name="l00750"></a>00750     main(i) = m(i, i);
564<a name="l00751"></a>00751     sup(i) = m(i, i + 1);
565<a name="l00752"></a>00752     sub(i) = m(i + 1, i);
566<a name="l00753"></a>00753   }
567<a name="l00754"></a>00754   main(n - 1) = m(n - 1, n - 1);
568<a name="l00755"></a>00755 }
569<a name="l00756"></a>00756
570<a name="l00757"></a>00757
571<a name="l00761"></a>00761 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
572<a name="l00762"></a><a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb">00762</a> T <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
573<a name="l00763"></a>00763 {
574<a name="l00764"></a>00764   <span class="keywordflow">return</span> <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(m));
575<a name="l00765"></a>00765 }
576<a name="l00766"></a>00766
577<a name="l00768"></a>00768
578<a name="l00769"></a>00769
579<a name="l00770"></a>00770 <span class="comment">// ----------------- reshaping vectors and matrices ------------------------</span>
580<a name="l00771"></a>00771
581<a name="l00774"></a>00774
582<a name="l00776"></a>00776 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
583<a name="l00777"></a><a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8">00777</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;in)
584<a name="l00778"></a>00778 {
585<a name="l00779"></a>00779   <span class="keywordtype">int</span> i, s = in.<a class="code" href="classitpp_1_1Vec.html#44478ea086eda22a180102e77e352325" title="The size of the vector.">length</a>();
586<a name="l00780"></a>00780
587<a name="l00781"></a>00781   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> out(s);
588<a name="l00782"></a>00782   <span class="keywordflow">for</span> (i = 0;i &lt; s;i++)
589<a name="l00783"></a>00783     out[i] = in[s-1-i];
590<a name="l00784"></a>00784   <span class="keywordflow">return</span> out;
591<a name="l00785"></a>00785 }
592<a name="l00786"></a>00786
593<a name="l00788"></a>00788 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
594<a name="l00789"></a><a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793">00789</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
595<a name="l00790"></a>00790 {
596<a name="l00791"></a>00791   <span class="keywordtype">int</span> i, j, n = 0, r = m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), c = m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
597<a name="l00792"></a>00792   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> v(r * c);
598<a name="l00793"></a>00793
599<a name="l00794"></a>00794   <span class="keywordflow">for</span> (i = 0; i &lt; r; i++)
600<a name="l00795"></a>00795     <span class="keywordflow">for</span> (j = 0; j &lt; c; j++)
601<a name="l00796"></a>00796       v(n++) = m(i, j);
602<a name="l00797"></a>00797
603<a name="l00798"></a>00798   <span class="keywordflow">return</span> v;
604<a name="l00799"></a>00799 }
605<a name="l00800"></a>00800
606<a name="l00802"></a>00802 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
607<a name="l00803"></a><a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0">00803</a> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m)
608<a name="l00804"></a>00804 {
609<a name="l00805"></a>00805   <span class="keywordtype">int</span> i, j, n = 0, r = m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(), c = m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>();
610<a name="l00806"></a>00806   <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> v(r * c);
611<a name="l00807"></a>00807
612<a name="l00808"></a>00808   <span class="keywordflow">for</span> (j = 0; j &lt; c; j++)
613<a name="l00809"></a>00809     <span class="keywordflow">for</span> (i = 0; i &lt; r; i++)
614<a name="l00810"></a>00810       v(n++) = m(i, j);
615<a name="l00811"></a>00811
616<a name="l00812"></a>00812   <span class="keywordflow">return</span> v;
617<a name="l00813"></a>00813 }
618<a name="l00814"></a>00814
619<a name="l00821"></a>00821 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
620<a name="l00822"></a><a class="code" href="group__reshaping.html#gcf8bc686cea4d9b8d2ddd853ee0e3bc2">00822</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols)
621<a name="l00823"></a>00823 {
622<a name="l00824"></a>00824   <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a>(m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>()*m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>() == rows*cols, <span class="stringliteral">"Mat&lt;T&gt;::reshape: Sizes must match"</span>);
623<a name="l00825"></a>00825   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> temp(rows, cols);
624<a name="l00826"></a>00826   <span class="keywordtype">int</span> i, j, ii = 0, jj = 0;
625<a name="l00827"></a>00827   <span class="keywordflow">for</span> (j = 0; j &lt; m.<a class="code" href="classitpp_1_1Mat.html#49680ef454d1d75b737e3ce157f1ac7b" title="The number of columns.">cols</a>(); j++) {
626<a name="l00828"></a>00828     <span class="keywordflow">for</span> (i = 0; i &lt; m.<a class="code" href="classitpp_1_1Mat.html#f3028609c01e08fb07c5bf0bcd4c4e03" title="The number of rows.">rows</a>(); i++) {
627<a name="l00829"></a>00829       temp(ii++, jj) = m(i, j);
628<a name="l00830"></a>00830       <span class="keywordflow">if</span> (ii == rows) {
629<a name="l00831"></a>00831         jj++;
630<a name="l00832"></a>00832         ii = 0;
631<a name="l00833"></a>00833       }
632<a name="l00834"></a>00834     }
633<a name="l00835"></a>00835   }
634<a name="l00836"></a>00836   <span class="keywordflow">return</span> temp;
635<a name="l00837"></a>00837 }
636<a name="l00838"></a>00838
637<a name="l00845"></a>00845 <span class="keyword">template</span>&lt;<span class="keyword">class</span> T&gt;
638<a name="l00846"></a><a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb">00846</a> <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> <a class="code" href="classitpp_1_1Vec.html">Vec&lt;T&gt;</a> &amp;v, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols)
639<a name="l00847"></a>00847 {
640<a name="l00848"></a>00848   <a class="code" href="group__errorhandlingfunc.html#gb319550e696ee9d824d23c2a176bc3a6" title="Abort if t is not true and NDEBUG is not defined.">it_assert_debug</a>(v.<a class="code" href="classitpp_1_1Vec.html#a906c893cd6184a774e4da8a47217d6a" title="The size of the vector.">size</a>() == rows*cols, <span class="stringliteral">"Mat&lt;T&gt;::reshape: Sizes must match"</span>);
641<a name="l00849"></a>00849   <a class="code" href="classitpp_1_1Mat.html">Mat&lt;T&gt;</a> temp(rows, cols);
642<a name="l00850"></a>00850   <span class="keywordtype">int</span> i, j, ii = 0;
643<a name="l00851"></a>00851   <span class="keywordflow">for</span> (j = 0; j &lt; cols; j++) {
644<a name="l00852"></a>00852     <span class="keywordflow">for</span> (i = 0; i &lt; rows; i++) {
645<a name="l00853"></a>00853       temp(i, j) = v(ii++);
646<a name="l00854"></a>00854     }
647<a name="l00855"></a>00855   }
648<a name="l00856"></a>00856   <span class="keywordflow">return</span> temp;
649<a name="l00857"></a>00857 }
650<a name="l00858"></a>00858
651<a name="l00860"></a>00860
652<a name="l00861"></a>00861
653<a name="l00863"></a>00863 <span class="keywordtype">bool</span> all(<span class="keyword">const</span> bvec &amp;testvec);
654<a name="l00865"></a>00865 <span class="keywordtype">bool</span> any(<span class="keyword">const</span> bvec &amp;testvec);
655<a name="l00866"></a>00866
656<a name="l00868"></a>00868
657<a name="l00869"></a>00869 <span class="comment">// ----------------------------------------------------------------------</span>
658<a name="l00870"></a>00870 <span class="comment">// Instantiations</span>
659<a name="l00871"></a>00871 <span class="comment">// ----------------------------------------------------------------------</span>
660<a name="l00872"></a>00872
661<a name="l00873"></a>00873 <span class="preprocessor">#ifdef HAVE_EXTERN_TEMPLATE</span>
662<a name="l00874"></a>00874 <span class="preprocessor"></span>
663<a name="l00875"></a>00875 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> vec &amp;v);
664<a name="l00876"></a>00876 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> cvec &amp;v);
665<a name="l00877"></a>00877 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> svec &amp;v);
666<a name="l00878"></a>00878 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> ivec &amp;v);
667<a name="l00879"></a>00879 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g97bb6b35d6d4c683c7558842e4500c2d" title="Length of vector.">length</a>(<span class="keyword">const</span> bvec &amp;v);
668<a name="l00880"></a>00880
669<a name="l00881"></a>00881 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">double</span> <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> vec &amp;v);
670<a name="l00882"></a>00882 <span class="keyword">extern</span> <span class="keyword">template</span> std::complex&lt;double&gt; <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> cvec &amp;v);
671<a name="l00883"></a>00883 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">short</span> <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> svec &amp;v);
672<a name="l00884"></a>00884 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> ivec &amp;v);
673<a name="l00885"></a>00885 <span class="keyword">extern</span> <span class="keyword">template</span> bin <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> bvec &amp;v);
674<a name="l00886"></a>00886
675<a name="l00887"></a>00887 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">double</span> <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> vec &amp;v);
676<a name="l00888"></a>00888 <span class="keyword">extern</span> <span class="keyword">template</span> std::complex&lt;double&gt; <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> cvec &amp;v);
677<a name="l00889"></a>00889 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">short</span> <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> svec &amp;v);
678<a name="l00890"></a>00890 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> ivec &amp;v);
679<a name="l00891"></a>00891 <span class="keyword">extern</span> <span class="keyword">template</span> bin <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> bvec &amp;v);
680<a name="l00892"></a>00892
681<a name="l00893"></a>00893 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> vec &amp;v);
682<a name="l00894"></a>00894 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> cvec &amp;v);
683<a name="l00895"></a>00895 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> svec &amp;v);
684<a name="l00896"></a>00896 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> ivec &amp;v);
685<a name="l00897"></a>00897 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> bvec &amp;v);
686<a name="l00898"></a>00898
687<a name="l00899"></a>00899 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">double</span> <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> vec &amp;v);
688<a name="l00900"></a>00900 <span class="keyword">extern</span> <span class="keyword">template</span> std::complex&lt;double&gt; <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> cvec &amp;v);
689<a name="l00901"></a>00901 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">short</span> <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> svec &amp;v);
690<a name="l00902"></a>00902 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> ivec &amp;v);
691<a name="l00903"></a>00903 <span class="keyword">extern</span> <span class="keyword">template</span> bin <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> bvec &amp;v);
692<a name="l00904"></a>00904
693<a name="l00905"></a>00905 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> vec &amp;v1, <span class="keyword">const</span> vec &amp;v2);
694<a name="l00906"></a>00906 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> cvec &amp;v1, <span class="keyword">const</span> cvec &amp;v2);
695<a name="l00907"></a>00907 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> ivec &amp;v1, <span class="keyword">const</span> ivec &amp;v2);
696<a name="l00908"></a>00908 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> svec &amp;v1, <span class="keyword">const</span> svec &amp;v2);
697<a name="l00909"></a>00909 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#gc9576b3da83bec579bacbfd2cfd71b3e" title="Vector cross product. Vectors need to be of size 3.">cross</a>(<span class="keyword">const</span> bvec &amp;v1, <span class="keyword">const</span> bvec &amp;v2);
698<a name="l00910"></a>00910
699<a name="l00911"></a>00911 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> vec &amp;in);
700<a name="l00912"></a>00912 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> cvec &amp;in);
701<a name="l00913"></a>00913 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> svec &amp;in);
702<a name="l00914"></a>00914 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> ivec &amp;in);
703<a name="l00915"></a>00915 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__reshaping.html#ge53eb5ba319eb6a25649576deeb9bca8" title="Reverse the input vector.">reverse</a>(<span class="keyword">const</span> bvec &amp;in);
704<a name="l00916"></a>00916
705<a name="l00917"></a>00917 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> vec &amp;v, <span class="keywordtype">int</span> n);
706<a name="l00918"></a>00918 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> cvec &amp;v, <span class="keywordtype">int</span> n);
707<a name="l00919"></a>00919 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> ivec &amp;v, <span class="keywordtype">int</span> n);
708<a name="l00920"></a>00920 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> svec &amp;v, <span class="keywordtype">int</span> n);
709<a name="l00921"></a>00921 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> bvec &amp;v, <span class="keywordtype">int</span> n);
710<a name="l00922"></a>00922
711<a name="l00923"></a>00923 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> vec &amp;v);
712<a name="l00924"></a>00924 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> cvec &amp;v);
713<a name="l00925"></a>00925 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> ivec &amp;v);
714<a name="l00926"></a>00926 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> svec &amp;v);
715<a name="l00927"></a>00927 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> bvec &amp;v);
716<a name="l00928"></a>00928
717<a name="l00929"></a>00929 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> mat &amp;, <span class="keywordtype">int</span>, <span class="keywordtype">int</span>);
718<a name="l00930"></a>00930 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> cmat &amp;, <span class="keywordtype">int</span>, <span class="keywordtype">int</span>);
719<a name="l00931"></a>00931 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> imat &amp;, <span class="keywordtype">int</span>, <span class="keywordtype">int</span>);
720<a name="l00932"></a>00932 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> smat &amp;, <span class="keywordtype">int</span>, <span class="keywordtype">int</span>);
721<a name="l00933"></a>00933 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__matrix__functions.html#g70bc01fbb0018476a9ba6a74ed75d4ea" title="Zero-pad a matrix to size rows x cols.">zero_pad</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;, <span class="keywordtype">int</span>, <span class="keywordtype">int</span>);
722<a name="l00934"></a>00934
723<a name="l00935"></a>00935 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> mat &amp;m, <span class="keywordtype">int</span> dim);
724<a name="l00936"></a>00936 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> cmat &amp;m, <span class="keywordtype">int</span> dim);
725<a name="l00937"></a>00937 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">int</span> dim);
726<a name="l00938"></a>00938 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">int</span> dim);
727<a name="l00939"></a>00939 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#gfcf8f54c2f4fc257bbc3111700b08355" title="Sum of elements in the matrix m, either along columns or rows.">sum</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <span class="keywordtype">int</span> dim);
728<a name="l00940"></a>00940
729<a name="l00941"></a>00941 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">double</span> <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> mat &amp;X);
730<a name="l00942"></a>00942 <span class="keyword">extern</span> <span class="keyword">template</span> std::complex&lt;double&gt; <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> cmat &amp;X);
731<a name="l00943"></a>00943 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">short</span> <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> smat &amp;X);
732<a name="l00944"></a>00944 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> imat &amp;X);
733<a name="l00945"></a>00945 <span class="keyword">extern</span> <span class="keyword">template</span> bin <a class="code" href="group__matrix__functions.html#g81a7516f7fb453dcb623ad6d4c89a5b7" title="Sum of all elements in the given matrix. Fast version of sum(sum(X)).">sumsum</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;X);
734<a name="l00946"></a>00946
735<a name="l00947"></a>00947 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> mat &amp; m, <span class="keywordtype">int</span> dim);
736<a name="l00948"></a>00948 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> cmat &amp;m, <span class="keywordtype">int</span> dim);
737<a name="l00949"></a>00949 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">int</span> dim);
738<a name="l00950"></a>00950 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">int</span> dim);
739<a name="l00951"></a>00951 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#g4c81137a6b9b989a0c5a36972885c1e6" title="Sum of the square of elements in the matrix m.">sum_sqr</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <span class="keywordtype">int</span> dim);
740<a name="l00952"></a>00952
741<a name="l00953"></a>00953 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> mat &amp;m, <span class="keywordtype">int</span> dim);
742<a name="l00954"></a>00954 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> cmat &amp;m, <span class="keywordtype">int</span> dim);
743<a name="l00955"></a>00955 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">int</span> dim);
744<a name="l00956"></a>00956 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">int</span> dim);
745<a name="l00957"></a>00957 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__matrix__functions.html#g6b5efcaa7c61e429326f4e6e56f11655" title="Cumulative sum of elements in the matrix m.">cumsum</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <span class="keywordtype">int</span> dim);
746<a name="l00958"></a>00958
747<a name="l00959"></a>00959 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> mat &amp;m, <span class="keywordtype">int</span> dim);
748<a name="l00960"></a>00960 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> cmat &amp;v, <span class="keywordtype">int</span> dim);
749<a name="l00961"></a>00961 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">int</span> dim);
750<a name="l00962"></a>00962 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">int</span> dim);
751<a name="l00963"></a>00963 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__matrix__functions.html#g22fdeace9f20483008cf4dc6d471cbc2" title="Product of elements in the matrix m.">prod</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <span class="keywordtype">int</span> dim);
752<a name="l00964"></a>00964
753<a name="l00965"></a>00965 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> mat &amp;in);
754<a name="l00966"></a>00966 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> cmat &amp;in);
755<a name="l00967"></a>00967 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> vec &amp;in, mat &amp;m);
756<a name="l00968"></a>00968 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> cvec &amp;in, cmat &amp;m);
757<a name="l00969"></a>00969 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> vec &amp;v, <span class="keyword">const</span> <span class="keywordtype">int</span> K);
758<a name="l00970"></a>00970 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__diag.html#gb0938c932c1cbc19b2ce6c5ac1007553" title="Get the diagonal elements of the input matrix m.">diag</a>(<span class="keyword">const</span> cvec &amp;v, <span class="keyword">const</span> <span class="keywordtype">int</span> K);
759<a name="l00971"></a>00971
760<a name="l00972"></a>00972 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> vec &amp;, <span class="keyword">const</span> vec &amp;);
761<a name="l00973"></a>00973 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> cvec &amp;, <span class="keyword">const</span> cvec &amp;);
762<a name="l00974"></a>00974 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> vec &amp;, <span class="keyword">const</span> vec &amp;, mat &amp;);
763<a name="l00975"></a>00975 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> cvec &amp;, <span class="keyword">const</span> cvec &amp;, cmat &amp;);
764<a name="l00976"></a>00976 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> mat &amp;, vec &amp;, vec &amp;);
765<a name="l00977"></a>00977 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g6b51bb44b5a610731b34f7f0ed2f4748" title="Returns the main diagonal and the diagonal row above in the two output vectors main...">bidiag</a>(<span class="keyword">const</span> cmat &amp;, cvec &amp;, cvec &amp;);
766<a name="l00978"></a>00978
767<a name="l00979"></a>00979 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> vec &amp;main, <span class="keyword">const</span> vec &amp;, <span class="keyword">const</span> vec &amp;);
768<a name="l00980"></a>00980 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> cvec &amp;main, <span class="keyword">const</span> cvec &amp;, <span class="keyword">const</span> cvec &amp;);
769<a name="l00981"></a>00981 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> vec &amp;main, <span class="keyword">const</span> vec &amp;, <span class="keyword">const</span> vec &amp;, mat &amp;);
770<a name="l00982"></a>00982 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> cvec &amp;main, <span class="keyword">const</span> cvec &amp;, <span class="keyword">const</span> cvec &amp;, cmat &amp;);
771<a name="l00983"></a>00983 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> mat &amp;m, vec &amp;, vec &amp;, vec &amp;);
772<a name="l00984"></a>00984 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__diag.html#g16dfe03af37710450c28943a3fc4433e" title="Returns the main diagonal, the diagonal row above, and the diagonal row below int...">tridiag</a>(<span class="keyword">const</span> cmat &amp;m, cvec &amp;, cvec &amp;, cvec &amp;);
773<a name="l00985"></a>00985
774<a name="l00986"></a>00986 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">double</span> <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> mat &amp;in);
775<a name="l00987"></a>00987 <span class="keyword">extern</span> <span class="keyword">template</span> std::complex&lt;double&gt; <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> cmat &amp;in);
776<a name="l00988"></a>00988 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">short</span> <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> smat &amp;in);
777<a name="l00989"></a>00989 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">int</span> <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> imat &amp;in);
778<a name="l00990"></a>00990 <span class="keyword">extern</span> <span class="keyword">template</span> bin <a class="code" href="group__diag.html#gd10b393e16f1a426b6daed5b9b78c3fb" title="The trace of the matrix m, i.e. the sum of the diagonal elements.">trace</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;in);
779<a name="l00991"></a>00991
780<a name="l00992"></a>00992 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> transpose(<span class="keyword">const</span> mat &amp;m, mat &amp;out);
781<a name="l00993"></a>00993 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> transpose(<span class="keyword">const</span> cmat &amp;m, cmat &amp;out);
782<a name="l00994"></a>00994 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> transpose(<span class="keyword">const</span> smat &amp;m, smat &amp;out);
783<a name="l00995"></a>00995 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> transpose(<span class="keyword">const</span> imat &amp;m, imat &amp;out);
784<a name="l00996"></a>00996 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> transpose(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;out);
785<a name="l00997"></a>00997
786<a name="l00998"></a>00998 <span class="keyword">extern</span> <span class="keyword">template</span> mat transpose(<span class="keyword">const</span> mat &amp;m);
787<a name="l00999"></a>00999 <span class="keyword">extern</span> <span class="keyword">template</span> cmat transpose(<span class="keyword">const</span> cmat &amp;m);
788<a name="l01000"></a>01000 <span class="keyword">extern</span> <span class="keyword">template</span> smat transpose(<span class="keyword">const</span> smat &amp;m);
789<a name="l01001"></a>01001 <span class="keyword">extern</span> <span class="keyword">template</span> imat transpose(<span class="keyword">const</span> imat &amp;m);
790<a name="l01002"></a>01002 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> transpose(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m);
791<a name="l01003"></a>01003
792<a name="l01004"></a>01004 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> mat &amp;m, mat &amp;out);
793<a name="l01005"></a>01005 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> cmat &amp;m, cmat &amp;out);
794<a name="l01006"></a>01006 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> smat &amp;m, smat &amp;out);
795<a name="l01007"></a>01007 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> imat &amp;m, imat &amp;out);
796<a name="l01008"></a>01008 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">void</span> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;out);
797<a name="l01009"></a>01009
798<a name="l01010"></a>01010 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> mat &amp;m);
799<a name="l01011"></a>01011 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> cmat &amp;m);
800<a name="l01012"></a>01012 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> smat &amp;m);
801<a name="l01013"></a>01013 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> imat &amp;m);
802<a name="l01014"></a>01014 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__matrix__functions.html#g123e8bb976c29bd43d993fcd6e3f3579" title="Hermitian transpose (complex conjugate transpose) of the matrix m.">hermitian_transpose</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m);
803<a name="l01015"></a>01015
804<a name="l01016"></a>01016 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g49691bf233320a2dd8ae9929988f8f91" title="Returns true if matrix X is hermitian, false otherwise.">is_hermitian</a>(<span class="keyword">const</span> mat &amp;X);
805<a name="l01017"></a>01017 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g49691bf233320a2dd8ae9929988f8f91" title="Returns true if matrix X is hermitian, false otherwise.">is_hermitian</a>(<span class="keyword">const</span> cmat &amp;X);
806<a name="l01018"></a>01018
807<a name="l01019"></a>01019 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g0e711a4a8db7d9a2ed5515a33dd65d4c" title="Returns true if matrix X is unitary, false otherwise.">is_unitary</a>(<span class="keyword">const</span> mat &amp;X);
808<a name="l01020"></a>01020 <span class="keyword">extern</span> <span class="keyword">template</span> <span class="keywordtype">bool</span> <a class="code" href="group__matrix__functions.html#g0e711a4a8db7d9a2ed5515a33dd65d4c" title="Returns true if matrix X is unitary, false otherwise.">is_unitary</a>(<span class="keyword">const</span> cmat &amp;X);
809<a name="l01021"></a>01021
810<a name="l01022"></a>01022 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> mat &amp;m);
811<a name="l01023"></a>01023 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> cmat &amp;m);
812<a name="l01024"></a>01024 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> imat &amp;m);
813<a name="l01025"></a>01025 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> smat &amp;m);
814<a name="l01026"></a>01026 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__reshaping.html#g67069e107b618a919a243f0f21fa8793" title="Row vectorize the matrix [(0,0) (0,1) ... (N-1,N-2) (N-1,N-1)].">rvectorize</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m);
815<a name="l01027"></a>01027
816<a name="l01028"></a>01028 <span class="keyword">extern</span> <span class="keyword">template</span> vec <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> mat &amp;m);
817<a name="l01029"></a>01029 <span class="keyword">extern</span> <span class="keyword">template</span> cvec <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> cmat &amp;m);
818<a name="l01030"></a>01030 <span class="keyword">extern</span> <span class="keyword">template</span> ivec <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> imat &amp;m);
819<a name="l01031"></a>01031 <span class="keyword">extern</span> <span class="keyword">template</span> svec <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> smat &amp;m);
820<a name="l01032"></a>01032 <span class="keyword">extern</span> <span class="keyword">template</span> bvec <a class="code" href="group__reshaping.html#g7796ab7ca5b96ce924c4e3f798f3fbe0" title="Column vectorize the matrix [(0,0) (1,0) ... (N-2,N-1) (N-1,N-1)].">cvectorize</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m);
821<a name="l01033"></a>01033
822<a name="l01034"></a>01034 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> mat &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
823<a name="l01035"></a>01035 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> cmat &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
824<a name="l01036"></a>01036 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> imat &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
825<a name="l01037"></a>01037 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> smat &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
826<a name="l01038"></a>01038 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
827<a name="l01039"></a>01039
828<a name="l01040"></a>01040 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> vec &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
829<a name="l01041"></a>01041 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> cvec &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
830<a name="l01042"></a>01042 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> ivec &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
831<a name="l01043"></a>01043 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> svec &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
832<a name="l01044"></a>01044 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__reshaping.html#gcaa0a45900248d17792ec91b2ab94dfb" title="Reshape the vector into an rows*cols matrix.">reshape</a>(<span class="keyword">const</span> bvec &amp;m, <span class="keywordtype">int</span> rows, <span class="keywordtype">int</span> cols);
833<a name="l01045"></a>01045
834<a name="l01046"></a>01046 <span class="keyword">extern</span> <span class="keyword">template</span> mat <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> mat &amp;X, <span class="keyword">const</span> mat &amp;Y);
835<a name="l01047"></a>01047 <span class="keyword">extern</span> <span class="keyword">template</span> cmat <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> cmat &amp;X, <span class="keyword">const</span> cmat &amp;Y);
836<a name="l01048"></a>01048 <span class="keyword">extern</span> <span class="keyword">template</span> imat <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> imat &amp;X, <span class="keyword">const</span> imat &amp;Y);
837<a name="l01049"></a>01049 <span class="keyword">extern</span> <span class="keyword">template</span> smat <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> smat &amp;X, <span class="keyword">const</span> smat &amp;Y);
838<a name="l01050"></a>01050 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> <a class="code" href="group__matrix__functions.html#g3b60d55399d969b0194f80997d41ea86" title="Computes the Kronecker product of two matrices.">kron</a>(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;X, <span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;Y);
839<a name="l01051"></a>01051
840<a name="l01052"></a>01052 <span class="keyword">extern</span> <span class="keyword">template</span> vec repmat(<span class="keyword">const</span> vec &amp;v, <span class="keywordtype">int</span> n);
841<a name="l01053"></a>01053 <span class="keyword">extern</span> <span class="keyword">template</span> cvec repmat(<span class="keyword">const</span> cvec &amp;v, <span class="keywordtype">int</span> n);
842<a name="l01054"></a>01054 <span class="keyword">extern</span> <span class="keyword">template</span> ivec repmat(<span class="keyword">const</span> ivec &amp;v, <span class="keywordtype">int</span> n);
843<a name="l01055"></a>01055 <span class="keyword">extern</span> <span class="keyword">template</span> svec repmat(<span class="keyword">const</span> svec &amp;v, <span class="keywordtype">int</span> n);
844<a name="l01056"></a>01056 <span class="keyword">extern</span> <span class="keyword">template</span> bvec repmat(<span class="keyword">const</span> bvec &amp;v, <span class="keywordtype">int</span> n);
845<a name="l01057"></a>01057
846<a name="l01058"></a>01058 <span class="keyword">extern</span> <span class="keyword">template</span> mat repmat(<span class="keyword">const</span> vec &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose);
847<a name="l01059"></a>01059 <span class="keyword">extern</span> <span class="keyword">template</span> cmat repmat(<span class="keyword">const</span> cvec &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose);
848<a name="l01060"></a>01060 <span class="keyword">extern</span> <span class="keyword">template</span> imat repmat(<span class="keyword">const</span> ivec &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose);
849<a name="l01061"></a>01061 <span class="keyword">extern</span> <span class="keyword">template</span> smat repmat(<span class="keyword">const</span> svec &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose);
850<a name="l01062"></a>01062 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> repmat(<span class="keyword">const</span> bvec &amp;v, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n, <span class="keywordtype">bool</span> transpose);
851<a name="l01063"></a>01063
852<a name="l01064"></a>01064 <span class="keyword">extern</span> <span class="keyword">template</span> mat repmat(<span class="keyword">const</span> mat &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n);
853<a name="l01065"></a>01065 <span class="keyword">extern</span> <span class="keyword">template</span> cmat repmat(<span class="keyword">const</span> cmat &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n);
854<a name="l01066"></a>01066 <span class="keyword">extern</span> <span class="keyword">template</span> imat repmat(<span class="keyword">const</span> imat &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n);
855<a name="l01067"></a>01067 <span class="keyword">extern</span> <span class="keyword">template</span> smat repmat(<span class="keyword">const</span> smat &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n);
856<a name="l01068"></a>01068 <span class="keyword">extern</span> <span class="keyword">template</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> repmat(<span class="keyword">const</span> <a class="code" href="mat_8h.html#f90acd1af41bf2d1d8a4bb23662fff69" title="bin matrix">bmat</a> &amp;data, <span class="keywordtype">int</span> m, <span class="keywordtype">int</span> n);
857<a name="l01069"></a>01069
858<a name="l01070"></a>01070 <span class="preprocessor">#endif // HAVE_EXTERN_TEMPLATE</span>
859<a name="l01071"></a>01071 <span class="preprocessor"></span>
860<a name="l01073"></a>01073
861<a name="l01074"></a>01074 } <span class="comment">// namespace itpp</span>
862<a name="l01075"></a>01075
863<a name="l01076"></a>01076 <span class="preprocessor">#endif // #ifndef MATFUNC_H</span>
864</pre></div></div>
865<hr size="1"><address style="text-align: right;"><small>Generated on Tue Jun 2 10:02:12 2009 for mixpp by&nbsp;
866<a href="http://www.doxygen.org/index.html">
867<img src="doxygen.png" alt="doxygen" align="middle" border="0"></a> 1.5.8 </small></address>
868</body>
869</html>
Note: See TracBrowser for help on using the browser.