<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<title>ublas: io.hpp Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="doxygen.css" rel="stylesheet" type="text/css"/>
</head>
<body>
<!-- Generated by Doxygen 1.6.1 -->
<div class="navigation" id="top">
<div class="tabs">
<ul>
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="pages.html"><span>Related Pages</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
</ul>
</div>
<div class="tabs">
<ul>
<li><a href="files.html"><span>File List</span></a></li>
</ul>
</div>
<h1>io.hpp</h1><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">//</span>
<a name="l00002"></a>00002 <span class="comment">// Copyright (c) 2000-2010</span>
<a name="l00003"></a>00003 <span class="comment">// Joerg Walter, Mathias Koch, David Bellot</span>
<a name="l00004"></a>00004 <span class="comment">//</span>
<a name="l00005"></a>00005 <span class="comment">// Distributed under the Boost Software License, Version 1.0. (See</span>
<a name="l00006"></a>00006 <span class="comment">// accompanying file LICENSE_1_0.txt or copy at</span>
<a name="l00007"></a>00007 <span class="comment">// http://www.boost.org/LICENSE_1_0.txt)</span>
<a name="l00008"></a>00008 <span class="comment">//</span>
<a name="l00009"></a>00009 <span class="comment">// The authors gratefully acknowledge the support of</span>
<a name="l00010"></a>00010 <span class="comment">// GeNeSys mbH & Co. KG in producing this work.</span>
<a name="l00011"></a>00011 <span class="comment">//</span>
<a name="l00012"></a>00012
<a name="l00013"></a>00013 <span class="preprocessor">#ifndef _BOOST_UBLAS_IO_</span>
<a name="l00014"></a>00014 <span class="preprocessor"></span><span class="preprocessor">#define _BOOST_UBLAS_IO_</span>
<a name="l00015"></a>00015 <span class="preprocessor"></span>
<a name="l00016"></a>00016 <span class="comment">// Only forward definition required to define stream operations</span>
<a name="l00017"></a>00017 <span class="preprocessor">#include <iosfwd></span>
<a name="l00018"></a>00018 <span class="preprocessor">#include <sstream></span>
<a name="l00019"></a>00019 <span class="preprocessor">#include <boost/numeric/ublas/matrix_expression.hpp></span>
<a name="l00020"></a>00020
<a name="l00021"></a>00021
<a name="l00022"></a>00022 <span class="keyword">namespace </span>boost { <span class="keyword">namespace </span>numeric { <span class="keyword">namespace </span>ublas {
<a name="l00023"></a>00023
<a name="l00046"></a>00046 <span class="keyword">template</span><<span class="keyword">class</span> E, <span class="keyword">class</span> T, <span class="keyword">class</span> VE>
<a name="l00047"></a>00047 <span class="comment">// BOOST_UBLAS_INLINE This function seems to be big. So we do not let the compiler inline it.</span>
<a name="l00048"></a>00048 std::basic_ostream<E, T> &operator << (std::basic_ostream<E, T> &os,
<a name="l00049"></a>00049 <span class="keyword">const</span> vector_expression<VE> &v) {
<a name="l00050"></a>00050 <span class="keyword">typedef</span> <span class="keyword">typename</span> VE::size_type size_type;
<a name="l00051"></a>00051 size_type size = v ().size ();
<a name="l00052"></a>00052 std::basic_ostringstream<E, T, std::allocator<E> > s;
<a name="l00053"></a>00053 s.flags (os.flags ());
<a name="l00054"></a>00054 s.imbue (os.getloc ());
<a name="l00055"></a>00055 s.precision (os.precision ());
<a name="l00056"></a>00056 s << <span class="charliteral">'['</span> << size << <span class="stringliteral">"]("</span>;
<a name="l00057"></a>00057 <span class="keywordflow">if</span> (size > 0)
<a name="l00058"></a>00058 s << v () (0);
<a name="l00059"></a>00059 <span class="keywordflow">for</span> (size_type i = 1; i < size; ++ i)
<a name="l00060"></a>00060 s << <span class="charliteral">','</span> << v () (i);
<a name="l00061"></a>00061 s << <span class="charliteral">')'</span>;
<a name="l00062"></a>00062 <span class="keywordflow">return</span> os << s.str ().c_str ();
<a name="l00063"></a>00063 }
<a name="l00064"></a>00064
<a name="l00087"></a>00087 <span class="keyword">template</span><<span class="keyword">class</span> E, <span class="keyword">class</span> T, <span class="keyword">class</span> VT, <span class="keyword">class</span> VA>
<a name="l00088"></a>00088 <span class="comment">// BOOST_UBLAS_INLINE This function seems to be big. So we do not let the compiler inline it.</span>
<a name="l00089"></a>00089 std::basic_istream<E, T> &operator >> (std::basic_istream<E, T> &is,
<a name="l00090"></a>00090 vector<VT, VA> &v) {
<a name="l00091"></a>00091 <span class="keyword">typedef</span> <span class="keyword">typename</span> vector<VT, VA>::size_type size_type;
<a name="l00092"></a>00092 E ch;
<a name="l00093"></a>00093 size_type size;
<a name="l00094"></a>00094 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'['</span>) {
<a name="l00095"></a>00095 is.putback (ch);
<a name="l00096"></a>00096 is.setstate (std::ios_base::failbit);
<a name="l00097"></a>00097 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (is >> size >> ch && ch != <span class="charliteral">']'</span>) {
<a name="l00098"></a>00098 is.putback (ch);
<a name="l00099"></a>00099 is.setstate (std::ios_base::failbit);
<a name="l00100"></a>00100 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00101"></a>00101 vector<VT, VA> s (size);
<a name="l00102"></a>00102 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'('</span>) {
<a name="l00103"></a>00103 is.putback (ch);
<a name="l00104"></a>00104 is.setstate (std::ios_base::failbit);
<a name="l00105"></a>00105 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00106"></a>00106 <span class="keywordflow">for</span> (size_type i = 0; i < size; i ++) {
<a name="l00107"></a>00107 <span class="keywordflow">if</span> (is >> s (i) >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00108"></a>00108 is.putback (ch);
<a name="l00109"></a>00109 <span class="keywordflow">if</span> (i < size - 1)
<a name="l00110"></a>00110 is.setstate (std::ios_base::failbit);
<a name="l00111"></a>00111 <span class="keywordflow">break</span>;
<a name="l00112"></a>00112 }
<a name="l00113"></a>00113 }
<a name="l00114"></a>00114 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">')'</span>) {
<a name="l00115"></a>00115 is.putback (ch);
<a name="l00116"></a>00116 is.setstate (std::ios_base::failbit);
<a name="l00117"></a>00117 }
<a name="l00118"></a>00118 }
<a name="l00119"></a>00119 <span class="keywordflow">if</span> (! is.fail ())
<a name="l00120"></a>00120 v.swap (s);
<a name="l00121"></a>00121 }
<a name="l00122"></a>00122 <span class="keywordflow">return</span> is;
<a name="l00123"></a>00123 }
<a name="l00124"></a>00124
<a name="l00148"></a>00148 <span class="keyword">template</span><<span class="keyword">class</span> E, <span class="keyword">class</span> T, <span class="keyword">class</span> ME>
<a name="l00149"></a>00149 <span class="comment">// BOOST_UBLAS_INLINE This function seems to be big. So we do not let the compiler inline it.</span>
<a name="l00150"></a>00150 std::basic_ostream<E, T> &operator << (std::basic_ostream<E, T> &os,
<a name="l00151"></a>00151 <span class="keyword">const</span> matrix_expression<ME> &m) {
<a name="l00152"></a>00152 <span class="keyword">typedef</span> <span class="keyword">typename</span> ME::size_type size_type;
<a name="l00153"></a>00153 size_type size1 = m ().size1 ();
<a name="l00154"></a>00154 size_type size2 = m ().size2 ();
<a name="l00155"></a>00155 std::basic_ostringstream<E, T, std::allocator<E> > s;
<a name="l00156"></a>00156 s.flags (os.flags ());
<a name="l00157"></a>00157 s.imbue (os.getloc ());
<a name="l00158"></a>00158 s.precision (os.precision ());
<a name="l00159"></a>00159 s << <span class="charliteral">'['</span> << size1 << <span class="charliteral">','</span> << size2 << <span class="stringliteral">"]("</span>;
<a name="l00160"></a>00160 <span class="keywordflow">if</span> (size1 > 0) {
<a name="l00161"></a>00161 s << <span class="charliteral">'('</span> ;
<a name="l00162"></a>00162 <span class="keywordflow">if</span> (size2 > 0)
<a name="l00163"></a>00163 s << m () (0, 0);
<a name="l00164"></a>00164 <span class="keywordflow">for</span> (size_type j = 1; j < size2; ++ j)
<a name="l00165"></a>00165 s << <span class="charliteral">','</span> << m () (0, j);
<a name="l00166"></a>00166 s << <span class="charliteral">')'</span>;
<a name="l00167"></a>00167 }
<a name="l00168"></a>00168 <span class="keywordflow">for</span> (size_type i = 1; i < size1; ++ i) {
<a name="l00169"></a>00169 s << <span class="stringliteral">",("</span> ;
<a name="l00170"></a>00170 <span class="keywordflow">if</span> (size2 > 0)
<a name="l00171"></a>00171 s << m () (i, 0);
<a name="l00172"></a>00172 <span class="keywordflow">for</span> (size_type j = 1; j < size2; ++ j)
<a name="l00173"></a>00173 s << <span class="charliteral">','</span> << m () (i, j);
<a name="l00174"></a>00174 s << <span class="charliteral">')'</span>;
<a name="l00175"></a>00175 }
<a name="l00176"></a>00176 s << <span class="charliteral">')'</span>;
<a name="l00177"></a>00177 <span class="keywordflow">return</span> os << s.str ().c_str ();
<a name="l00178"></a>00178 }
<a name="l00179"></a>00179
<a name="l00199"></a>00199 <span class="keyword">template</span><<span class="keyword">class</span> E, <span class="keyword">class</span> T, <span class="keyword">class</span> MT, <span class="keyword">class</span> MF, <span class="keyword">class</span> MA>
<a name="l00200"></a>00200 <span class="comment">// BOOST_UBLAS_INLINE This function seems to be big. So we do not let the compiler inline it.</span>
<a name="l00201"></a>00201 std::basic_istream<E, T> &operator >> (std::basic_istream<E, T> &is,
<a name="l00202"></a>00202 matrix<MT, MF, MA> &m) {
<a name="l00203"></a>00203 <span class="keyword">typedef</span> <span class="keyword">typename</span> matrix<MT, MF, MA>::size_type size_type;
<a name="l00204"></a>00204 E ch;
<a name="l00205"></a>00205 size_type size1, size2;
<a name="l00206"></a>00206 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'['</span>) {
<a name="l00207"></a>00207 is.putback (ch);
<a name="l00208"></a>00208 is.setstate (std::ios_base::failbit);
<a name="l00209"></a>00209 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (is >> size1 >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00210"></a>00210 is.putback (ch);
<a name="l00211"></a>00211 is.setstate (std::ios_base::failbit);
<a name="l00212"></a>00212 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (is >> size2 >> ch && ch != <span class="charliteral">']'</span>) {
<a name="l00213"></a>00213 is.putback (ch);
<a name="l00214"></a>00214 is.setstate (std::ios_base::failbit);
<a name="l00215"></a>00215 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00216"></a>00216 matrix<MT, MF, MA> s (size1, size2);
<a name="l00217"></a>00217 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'('</span>) {
<a name="l00218"></a>00218 is.putback (ch);
<a name="l00219"></a>00219 is.setstate (std::ios_base::failbit);
<a name="l00220"></a>00220 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00221"></a>00221 <span class="keywordflow">for</span> (size_type i = 0; i < size1; i ++) {
<a name="l00222"></a>00222 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'('</span>) {
<a name="l00223"></a>00223 is.putback (ch);
<a name="l00224"></a>00224 is.setstate (std::ios_base::failbit);
<a name="l00225"></a>00225 <span class="keywordflow">break</span>;
<a name="l00226"></a>00226 }
<a name="l00227"></a>00227 <span class="keywordflow">for</span> (size_type j = 0; j < size2; j ++) {
<a name="l00228"></a>00228 <span class="keywordflow">if</span> (is >> s (i, j) >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00229"></a>00229 is.putback (ch);
<a name="l00230"></a>00230 <span class="keywordflow">if</span> (j < size2 - 1) {
<a name="l00231"></a>00231 is.setstate (std::ios_base::failbit);
<a name="l00232"></a>00232 <span class="keywordflow">break</span>;
<a name="l00233"></a>00233 }
<a name="l00234"></a>00234 }
<a name="l00235"></a>00235 }
<a name="l00236"></a>00236 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">')'</span>) {
<a name="l00237"></a>00237 is.putback (ch);
<a name="l00238"></a>00238 is.setstate (std::ios_base::failbit);
<a name="l00239"></a>00239 <span class="keywordflow">break</span>;
<a name="l00240"></a>00240 }
<a name="l00241"></a>00241 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00242"></a>00242 is.putback (ch);
<a name="l00243"></a>00243 <span class="keywordflow">if</span> (i < size1 - 1) {
<a name="l00244"></a>00244 is.setstate (std::ios_base::failbit);
<a name="l00245"></a>00245 <span class="keywordflow">break</span>;
<a name="l00246"></a>00246 }
<a name="l00247"></a>00247 }
<a name="l00248"></a>00248 }
<a name="l00249"></a>00249 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">')'</span>) {
<a name="l00250"></a>00250 is.putback (ch);
<a name="l00251"></a>00251 is.setstate (std::ios_base::failbit);
<a name="l00252"></a>00252 }
<a name="l00253"></a>00253 }
<a name="l00254"></a>00254 <span class="keywordflow">if</span> (! is.fail ())
<a name="l00255"></a>00255 m.swap (s);
<a name="l00256"></a>00256 }
<a name="l00257"></a>00257 <span class="keywordflow">return</span> is;
<a name="l00258"></a>00258 }
<a name="l00259"></a>00259
<a name="l00281"></a>00281 <span class="keyword">template</span><<span class="keyword">class</span> E, <span class="keyword">class</span> T, <span class="keyword">class</span> MT, <span class="keyword">class</span> MF1, <span class="keyword">class</span> MF2, <span class="keyword">class</span> MA>
<a name="l00282"></a>00282 <span class="comment">// BOOST_UBLAS_INLINE This function seems to be big. So we do not let the compiler inline it.</span>
<a name="l00283"></a>00283 std::basic_istream<E, T> &operator >> (std::basic_istream<E, T> &is,
<a name="l00284"></a>00284 symmetric_matrix<MT, MF1, MF2, MA> &m) {
<a name="l00285"></a>00285 <span class="keyword">typedef</span> <span class="keyword">typename</span> symmetric_matrix<MT, MF1, MF2, MA>::size_type size_type;
<a name="l00286"></a>00286 E ch;
<a name="l00287"></a>00287 size_type size1, size2;
<a name="l00288"></a>00288 MT value;
<a name="l00289"></a>00289 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'['</span>) {
<a name="l00290"></a>00290 is.putback (ch);
<a name="l00291"></a>00291 is.setstate (std::ios_base::failbit);
<a name="l00292"></a>00292 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (is >> size1 >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00293"></a>00293 is.putback (ch);
<a name="l00294"></a>00294 is.setstate (std::ios_base::failbit);
<a name="l00295"></a>00295 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (is >> size2 >> ch && (size2 != size1 || ch != <span class="charliteral">']'</span>)) { <span class="comment">// symmetric matrix must be square</span>
<a name="l00296"></a>00296 is.putback (ch);
<a name="l00297"></a>00297 is.setstate (std::ios_base::failbit);
<a name="l00298"></a>00298 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00299"></a>00299 symmetric_matrix<MT, MF1, MF2, MA> s (size1, size2);
<a name="l00300"></a>00300 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'('</span>) {
<a name="l00301"></a>00301 is.putback (ch);
<a name="l00302"></a>00302 is.setstate (std::ios_base::failbit);
<a name="l00303"></a>00303 } <span class="keywordflow">else</span> <span class="keywordflow">if</span> (! is.fail ()) {
<a name="l00304"></a>00304 <span class="keywordflow">for</span> (size_type i = 0; i < size1; i ++) {
<a name="l00305"></a>00305 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">'('</span>) {
<a name="l00306"></a>00306 is.putback (ch);
<a name="l00307"></a>00307 is.setstate (std::ios_base::failbit);
<a name="l00308"></a>00308 <span class="keywordflow">break</span>;
<a name="l00309"></a>00309 }
<a name="l00310"></a>00310 <span class="keywordflow">for</span> (size_type j = 0; j < size2; j ++) {
<a name="l00311"></a>00311 <span class="keywordflow">if</span> (is >> value >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00312"></a>00312 is.putback (ch);
<a name="l00313"></a>00313 <span class="keywordflow">if</span> (j < size2 - 1) {
<a name="l00314"></a>00314 is.setstate (std::ios_base::failbit);
<a name="l00315"></a>00315 <span class="keywordflow">break</span>;
<a name="l00316"></a>00316 }
<a name="l00317"></a>00317 }
<a name="l00318"></a>00318 <span class="keywordflow">if</span> (i <= j) {
<a name="l00319"></a>00319 <span class="comment">// this is the first time we read this element - set the value</span>
<a name="l00320"></a>00320 s(i,j) = value;
<a name="l00321"></a>00321 }
<a name="l00322"></a>00322 <span class="keywordflow">else</span> <span class="keywordflow">if</span> ( s(i,j) != value ) {
<a name="l00323"></a>00323 <span class="comment">// matrix is not symmetric</span>
<a name="l00324"></a>00324 is.setstate (std::ios_base::failbit);
<a name="l00325"></a>00325 <span class="keywordflow">break</span>;
<a name="l00326"></a>00326 }
<a name="l00327"></a>00327 }
<a name="l00328"></a>00328 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">')'</span>) {
<a name="l00329"></a>00329 is.putback (ch);
<a name="l00330"></a>00330 is.setstate (std::ios_base::failbit);
<a name="l00331"></a>00331 <span class="keywordflow">break</span>;
<a name="l00332"></a>00332 }
<a name="l00333"></a>00333 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">','</span>) {
<a name="l00334"></a>00334 is.putback (ch);
<a name="l00335"></a>00335 <span class="keywordflow">if</span> (i < size1 - 1) {
<a name="l00336"></a>00336 is.setstate (std::ios_base::failbit);
<a name="l00337"></a>00337 <span class="keywordflow">break</span>;
<a name="l00338"></a>00338 }
<a name="l00339"></a>00339 }
<a name="l00340"></a>00340 }
<a name="l00341"></a>00341 <span class="keywordflow">if</span> (is >> ch && ch != <span class="charliteral">')'</span>) {
<a name="l00342"></a>00342 is.putback (ch);
<a name="l00343"></a>00343 is.setstate (std::ios_base::failbit);
<a name="l00344"></a>00344 }
<a name="l00345"></a>00345 }
<a name="l00346"></a>00346 <span class="keywordflow">if</span> (! is.fail ())
<a name="l00347"></a>00347 m.swap (s);
<a name="l00348"></a>00348 }
<a name="l00349"></a>00349 <span class="keywordflow">return</span> is;
<a name="l00350"></a>00350 }
<a name="l00351"></a>00351
<a name="l00352"></a>00352
<a name="l00353"></a>00353 }}}
<a name="l00354"></a>00354
<a name="l00355"></a>00355 <span class="preprocessor">#endif</span>
</pre></div></div>
<hr size="1"/><address style="text-align: right;"><small>Generated on Sun Jul 4 20:30:51 2010 for ublas by
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.6.1 </small></address>
</body>
</html>
