<!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>linbox: examples/smithvalence.C</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.7.4 -->
<script type="text/javascript">
function hasClass(ele,cls) {
return ele.className.match(new RegExp('(\\s|^)'+cls+'(\\s|$)'));
}
function addClass(ele,cls) {
if (!this.hasClass(ele,cls)) ele.className += " "+cls;
}
function removeClass(ele,cls) {
if (hasClass(ele,cls)) {
var reg = new RegExp('(\\s|^)'+cls+'(\\s|$)');
ele.className=ele.className.replace(reg,' ');
}
}
function toggleVisibility(linkObj) {
var base = linkObj.getAttribute('id');
var summary = document.getElementById(base + '-summary');
var content = document.getElementById(base + '-content');
var trigger = document.getElementById(base + '-trigger');
if ( hasClass(linkObj,'closed') ) {
summary.style.display = 'none';
content.style.display = 'block';
trigger.src = 'open.png';
removeClass(linkObj,'closed');
addClass(linkObj,'opened');
} else if ( hasClass(linkObj,'opened') ) {
summary.style.display = 'block';
content.style.display = 'none';
trigger.src = 'closed.png';
removeClass(linkObj,'opened');
addClass(linkObj,'closed');
}
return false;
}
</script>
<div id="top">
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td style="padding-left: 0.5em;">
<div id="projectname">linbox</div>
</td>
</tr>
</tbody>
</table>
</div>
<div id="navrow1" class="tabs">
<ul class="tablist">
<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="modules.html"><span>Modules</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li><a href="files.html"><span>Files</span></a></li>
<li><a href="dirs.html"><span>Directories</span></a></li>
<li><a href="examples.html"><span>Examples</span></a></li>
</ul>
</div>
</div>
<div class="header">
<div class="headertitle">
<div class="title">examples/smithvalence.C</div> </div>
</div>
<div class="contents">
<p>Valence of sparse matrix over Z or Zp.</p>
<div class="fragment"><pre class="fragment"><span class="comment">/* -*- mode: C++; tab-width: 8; indent-tabs-mode: t; c-basic-offset: 8 -*- */</span>
<span class="comment">// vim:sts=8:sw=8:ts=8:noet:sr:cino=>s,f0,{0,g0,(0,\:0,t0,+0,=s</span>
<span class="comment">/*</span>
<span class="comment"> * examples/smithvalence.C</span>
<span class="comment"> *</span>
<span class="comment"> * Copyright (C) 2010 J-G Dumas</span>
<span class="comment"> *</span>
<span class="comment"> * This file is part of LinBox.</span>
<span class="comment"> *</span>
<span class="comment"> * LinBox is free software: you can redistribute it and/or modify</span>
<span class="comment"> * it under the terms of the GNU Lesser General Public License as</span>
<span class="comment"> * published by the Free Software Foundation, either version 2 of</span>
<span class="comment"> * the License, or (at your option) any later version.</span>
<span class="comment"> *</span>
<span class="comment"> * LinBox is distributed in the hope that it will be useful,</span>
<span class="comment"> * but WITHOUT ANY WARRANTY; without even the implied warranty of</span>
<span class="comment"> * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span>
<span class="comment"> * GNU Lesser General Public License for more details.</span>
<span class="comment"> *</span>
<span class="comment"> * You should have received a copy of the GNU Lesser General Public</span>
<span class="comment"> * License along with LinBox. If not, see</span>
<span class="comment"> * <http://www.gnu.org/licenses/>.</span>
<span class="comment"> */</span>
<span class="comment">//#include "linbox-config.h"</span>
<span class="preprocessor">#include <iostream></span>
<span class="preprocessor">#include "linbox/field/gf2.h"</span>
<span class="preprocessor">#include "<a class="code" href="modular-double_8h.html" title="Standard representation of Z/mZ over double .">linbox/field/modular-double.h</a>"</span>
<span class="preprocessor">#include "linbox/field/givaro-zpz.h"</span>
<span class="preprocessor">#include "linbox/field/field-traits.h"</span>
<span class="preprocessor">#include "linbox/blackbox/transpose.h"</span>
<span class="preprocessor">#include "linbox/blackbox/compose.h"</span>
<span class="preprocessor">#include "linbox/blackbox/sparse.h"</span>
<span class="preprocessor">#include "linbox/solutions/rank.h"</span>
<span class="preprocessor">#include "linbox/solutions/valence.h"</span>
<span class="preprocessor">#include "linbox/algorithms/smith-form-sparseelim-local.h"</span>
<span class="preprocessor">#include "linbox/util/matrix-stream.h"</span>
<span class="preprocessor">#include "<a class="code" href="integer_8h.html" title="This is a representation of arbitrary integers.">linbox/integer.h</a>"</span>
<span class="preprocessor">#include <givaro/givintnumtheo.h></span>
<span class="keyword">using</span> <a name="a0"></a><a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">LinBox::integer</a>;
<span class="keyword">template</span><<span class="keyword">class</span> Field>
<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& TempLRank(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& r, <span class="keywordtype">char</span> * filename, <span class="keyword">const</span> <a name="_a1"></a><a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>& F)
{
std::ifstream input(filename);
<a name="_a2"></a><a class="code" href="class_lin_box_1_1_matrix_stream.html" title="MatrixStream.">LinBox::MatrixStream< Field ></a> msf( F, input );
<a name="_a3"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html" title="vector of sparse rows.">LinBox::SparseMatrix<Field,typename LinBox::Vector<Field>::SparseSeq</a>> FA(msf);
input.close();
::Givaro::Timer tim; tim.start();
<a name="a4"></a><a class="code" href="namespace_lin_box.html#ac9442bdcae8ef5fc90e5a255630bffa7" title="Rank of A.">LinBox::rankin</a>(r, FA);
tim.stop();
F.write(std::cout << <span class="stringliteral">"Rank over "</span>) << <span class="stringliteral">" is "</span> << r << <span class="charliteral">' '</span> << tim << std::endl;
<span class="keywordflow">return</span> r;
}
<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& TempLRank(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& r, <span class="keywordtype">char</span> * filename, <span class="keyword">const</span> LinBox::GF2& F2)
{
std::ifstream input(filename);
<a name="_a5"></a><a class="code" href="class_lin_box_1_1_zero_one.html" title="Time and space efficient representation of sparse {0,1}-matrices.">LinBox::ZeroOne<LinBox::GF2></a> A;
A.<a name="a6"></a><a class="code" href="class_lin_box_1_1_zero_one.html#a7f1e7ece6c12a9a3feb091f90c7c9ac1" title="Read the matrix from a stream in the JGD's SMS format.">read</a>(input);
input.close();
::Givaro::Timer tim; tim.start();
<a class="code" href="namespace_lin_box.html#ac9442bdcae8ef5fc90e5a255630bffa7" title="Rank of A.">LinBox::rankin</a>(r, A, <a class="code" href="struct_lin_box_1_1_method.html#abc071c5d851893868ab7a0c0d8ee34cb" title="Method::SparseElimination : no doc.">LinBox::Method::SparseElimination</a>() );
tim.stop();
F2.write(std::cout << <span class="stringliteral">"Rank over "</span>) << <span class="stringliteral">" is "</span> << r << <span class="charliteral">' '</span> << tim << std::endl;
<span class="keywordflow">return</span> r;
}
<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& LRank(<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span>& r, <span class="keywordtype">char</span> * filename, <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> p)
{
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> maxmod16; LinBox::FieldTraits<LinBox::GivaroZpz< ::Givaro::Std16> >::maxModulus(maxmod16);
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> maxmod32; LinBox::FieldTraits<LinBox::GivaroZpz< ::Givaro::Std32> >::maxModulus(maxmod32);
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> maxmod53; LinBox::FieldTraits<LinBox::Modular<double> >::maxModulus(maxmod53);
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> maxmod64; LinBox::FieldTraits<LinBox::GivaroZpz< ::Givaro::Std64> >::maxModulus(maxmod64);
<span class="keywordflow">if</span> (p == 2) {
LinBox::GF2 F2;
<span class="keywordflow">return</span> TempLRank(r, filename, F2);
}
<span class="keywordflow">else</span> <span class="keywordflow">if</span> (p <= maxmod16) {
<span class="keyword">typedef</span> <a name="_a7"></a><a class="code" href="class_lin_box_1_1_givaro_zpz.html">LinBox::GivaroZpz< ::Givaro::Std16></a> <a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>;
Field F(p);
<span class="keywordflow">return</span> TempLRank(r, filename, F);
}
<span class="keywordflow">else</span> <span class="keywordflow">if</span> (p <= maxmod32) {
<span class="keyword">typedef</span> <a name="_a8"></a><a class="code" href="class_lin_box_1_1_givaro_zpz.html">LinBox::GivaroZpz< ::Givaro::Std32></a> <a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>;
Field F(p);
<span class="keywordflow">return</span> TempLRank(r, filename, F);
}
<span class="keywordflow">else</span> <span class="keywordflow">if</span> (p <= maxmod53) {
<span class="keyword">typedef</span> <a name="_a9"></a><a class="code" href="class_lin_box_1_1_modular_3_01double_01_4.html" title="Standard representation of .">LinBox::Modular<double></a> <a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>;
Field F(p);
<span class="keywordflow">return</span> TempLRank(r, filename, F);
}
<span class="keywordflow">else</span> <span class="keywordflow">if</span> (p <= maxmod64) {
<span class="keyword">typedef</span> <a name="_a10"></a><a class="code" href="class_lin_box_1_1_givaro_zpz.html" title="wrapper of Givaro's ZpzDom.">LinBox::GivaroZpz< ::Givaro::Std64></a> <a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>;
Field F(p);
<span class="keywordflow">return</span> TempLRank(r, filename, F);
}
<span class="keywordflow">else</span> {
<span class="keyword">typedef</span> <a class="code" href="class_lin_box_1_1_givaro_zpz.html" title="wrapper of Givaro's ZpzDom.">LinBox::GivaroZpz<integer></a> <a class="code" href="class_lin_box_1_1_modular_3_01uint32__t_01_4.html" title="Specialization of class Modular for uint32_t element type.">Field</a>;
Field F(p);
<span class="keywordflow">return</span> TempLRank(r, filename, F);
}
<span class="keywordflow">return</span> r;
}
std::vector<size_t>& PRank(std::vector<size_t>& ranks, <span class="keywordtype">char</span> * filename, <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> p, <span class="keywordtype">size_t</span> e, <span class="keywordtype">size_t</span> intr)
{
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> maxmod;
LinBox::FieldTraits<LinBox::GivaroZpz< ::Givaro::Std64> >::maxModulus(maxmod);
<span class="keywordflow">if</span> (p <= maxmod) {
<span class="keyword">typedef</span> <a class="code" href="class_lin_box_1_1_givaro_zpz.html" title="wrapper of Givaro's ZpzDom.">LinBox::GivaroZpz< ::Givaro::Std64></a> Ring;
int64_t lp(p);
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> q = pow(p,e); int64_t lq(q);
<span class="keywordflow">if</span> (q > <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a>(lq)) {
std::cerr << <span class="stringliteral">"Sorry power rank mod large composite not yet implemented"</span> << std::endl;
q = p;
<span class="keywordflow">do</span> {
q *= p; lq = (int64_t)q;
} <span class="keywordflow">while</span> (q == <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a>(lq));
q/=p; lq = (int64_t)q;
std::cerr << <span class="stringliteral">"Trying: "</span> << lq << std::endl;
}
Ring F(lq);
std::ifstream input(filename);
<a class="code" href="class_lin_box_1_1_matrix_stream.html" title="MatrixStream.">LinBox::MatrixStream<Ring></a> ms( F, input );
<a class="code" href="class_lin_box_1_1_sparse_matrix.html" title="vector of sparse rows.">LinBox::SparseMatrix<Ring, LinBox::Vector<Ring>::SparseSeq</a> > A (ms);
input.close();
<a name="_a11"></a><a class="code" href="class_lin_box_1_1_power_gauss_domain.html" title="Repository of functions for rank modulo a prime power by elimination on sparse matrices.">LinBox::PowerGaussDomain< Ring ></a> PGD( F );
PGD.prime_power_rankin( lq, lp, ranks, A, A.<a name="a12"></a>rowdim(), A.<a name="a13"></a>coldim(), std::vector<size_t>());
F.write(std::cout << <span class="stringliteral">"Ranks over "</span>) << <span class="stringliteral">" are "</span> ;
<span class="keywordflow">for</span>(std::vector<size_t>::const_iterator rit=ranks.begin(); rit != ranks.end(); ++rit)
std::cout << *rit << <span class="charliteral">' '</span>;
std::cout << std::endl;
}
<span class="keywordflow">else</span> {
std::cerr << <span class="stringliteral">"Sorry power rank mod large composite not yet implemented"</span> << std::endl;
std::cerr << <span class="stringliteral">"Assuming integer rank"</span> << std::endl;
ranks.resize(0); ranks.push_back(intr);
}
<span class="keywordflow">return</span> ranks;
}
<span class="keyword">using namespace </span>LinBox;
<span class="keywordtype">int</span> main (<span class="keywordtype">int</span> argc, <span class="keywordtype">char</span> **argv)
{
<span class="comment">// commentator.setMaxDetailLevel (-1);</span>
<span class="comment">// commentator.setMaxDepth (-1);</span>
<span class="comment">// commentator.setReportStream (std::cerr);</span>
<span class="keywordflow">if</span> (argc < 2 || argc > 4) {
std::cerr << <span class="stringliteral">"Usage: valence <matrix-file-in-supported-format> [-ata|-aat|valence] [coprime]"</span> << std::endl;
<span class="keywordflow">return</span> -1;
}
std::ifstream input (argv[1]);
<span class="keywordflow">if</span> (!input) { std::cerr << <span class="stringliteral">"Error opening matrix file "</span> << argv[1] << std::endl; <span class="keywordflow">return</span> -1; }
<a name="_a14"></a><a class="code" href="class_lin_box_1_1_p_i_d__integer.html" title="Domain for integer operations.">PID_integer</a> ZZ;
<a class="code" href="class_lin_box_1_1_matrix_stream.html" title="MatrixStream.">MatrixStream< PID_integer ></a> ms( ZZ, input );
<span class="keyword">typedef</span> <a class="code" href="class_lin_box_1_1_sparse_matrix.html" title="vector of sparse rows.">SparseMatrix<PID_integer></a> Blackbox;
Blackbox A (ms);
input.close();
std::cout << <span class="stringliteral">"A is "</span> << A.rowdim() << <span class="stringliteral">" by "</span> << A.coldim() << std::endl;
PID_integer::Element val_A;
::Givaro::Timer chrono; chrono.start();
<span class="keywordflow">if</span> (argc >= 3) {
<a name="_a15"></a><a class="code" href="class_lin_box_1_1_transpose.html" title="transpose matrix without copying.">Transpose<Blackbox></a> T(&A);
<span class="keywordflow">if</span> (strcmp(argv[2],<span class="stringliteral">"-ata"</span>) == 0) {
<a name="_a16"></a><a class="code" href="class_lin_box_1_1_compose.html" title="Blackbox of a product: , i.e .">Compose< Transpose<Blackbox></a>, Blackbox > C (&T, &A);
std::cout << <span class="stringliteral">"A^T A is "</span> << C.rowdim() << <span class="stringliteral">" by "</span> << C.coldim() << std::endl;
valence(val_A, C);
}
<span class="keywordflow">else</span> <span class="keywordflow">if</span> (strcmp(argv[2],<span class="stringliteral">"-aat"</span>) == 0) {
<a class="code" href="class_lin_box_1_1_compose.html" title="Blackbox of a product: , i.e .">Compose< Blackbox, Transpose<Blackbox></a> > C (&A, &T);
std::cout << <span class="stringliteral">"A A^T is "</span> << C.rowdim() << <span class="stringliteral">" by "</span> << C.coldim() << std::endl;
valence(val_A, C);
}
<span class="keywordflow">else</span> {
std::cout << <span class="stringliteral">"Suppose primes are contained in "</span> << argv[2] << std::endl;
val_A = <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a>(argv[2]);
}
}
<span class="keywordflow">else</span> {
<span class="keywordflow">if</span> (A.rowdim() != A.coldim()) {
std::cerr << <span class="stringliteral">"Valence works only on square matrices, try either to change the dimension in the matrix file, or to compute the valence of A A^T or A^T A, via the -aat or -ata options."</span> << std::endl;
exit(0);
}
<span class="keywordflow">else</span>
valence (val_A, A);
}
std::cout << <span class="stringliteral">"Valence is "</span> << val_A << std::endl;
std::vector<integer> Moduli;
std::vector<size_t> exponents;
::Givaro::IntFactorDom<> FTD;
<span class="keyword">typedef</span> std::pair<integer,unsigned long> PairIntRk;
std::vector< PairIntRk > smith;
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> coprimeV=2;
<span class="keywordflow">if</span> (argc >= 4) {
coprimeV = <a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a>(argv[3]);
}
<span class="keywordflow">while</span> ( gcd(val_A,coprimeV) > 1 ) {
FTD.nextprimein(coprimeV);
}
<span class="keywordflow">if</span> (argc >= 4) {
std::cout << <span class="stringliteral">"Suppose "</span> << argv[3] << <span class="stringliteral">" is coprime with Smith form"</span> << std::endl;
}
std::cout << <span class="stringliteral">"integer rank: "</span> << std::endl;
<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> coprimeR; LRank(coprimeR, argv[1], coprimeV);
smith.push_back(PairIntRk(coprimeV, coprimeR));
<span class="comment">// std::cerr << "Rank mod " << coprimeV << " is " << coprimeR << std::endl;</span>
std::cout << <span class="stringliteral">"Some factors (5000 factoring loop bound): "</span>;
FTD.set(Moduli, exponents, val_A, 5000);
std::vector<size_t>::const_iterator eit=exponents.begin();
<span class="keywordflow">for</span>(std::vector<integer>::const_iterator mit=Moduli.begin();
mit != Moduli.end(); ++mit,++eit)
std::cout << *mit << <span class="charliteral">'^'</span> << *eit << <span class="charliteral">' '</span>;
std::cout << std::endl;
std::vector<integer> SmithDiagonal(coprimeR,<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a>(1));
<span class="keywordflow">for</span>(std::vector<integer>::const_iterator mit=Moduli.begin();
mit != Moduli.end(); ++mit) {
<span class="keywordtype">unsigned</span> <span class="keywordtype">long</span> r; LRank(r, argv[1], *mit);
<span class="comment">// std::cerr << "Rank mod " << *mit << " is " << r << std::endl;</span>
smith.push_back(PairIntRk(*mit, r));
<span class="keywordflow">for</span>(<span class="keywordtype">size_t</span> i=r; i < coprimeR; ++i)
SmithDiagonal[i] *= *mit;
}
eit=exponents.begin();
std::vector<PairIntRk>::const_iterator sit=smith.begin();
<span class="keywordflow">for</span>( ++sit; sit != smith.end(); ++sit, ++eit) {
<span class="keywordflow">if</span> (sit->second != coprimeR) {
std::vector<size_t> ranks;
ranks.push_back(sit->second);
<span class="keywordflow">if</span> (*eit > 1) {
PRank(ranks, argv[1], sit->first, *eit, coprimeR);
}
<span class="keywordflow">else</span> {
PRank(ranks, argv[1], sit->first, 2, coprimeR);
}
<span class="keywordflow">if</span> (ranks.size() == 1) ranks.push_back(coprimeR);
<span class="keywordflow">for</span>(<span class="keywordtype">size_t</span> expo = (*eit)<<1; ranks.back() < coprimeR; expo<<=1) {
PRank(ranks, argv[1], sit->first, expo, coprimeR);
<span class="keywordflow">if</span> (ranks.size() < expo) {
std::cerr << <span class="stringliteral">"Larger prime power not yet implemented"</span> << std::endl;
<span class="keywordflow">break</span>;
}
}
std::vector<size_t>::const_iterator rit=ranks.begin();
<span class="comment">// unsigned long modrank = *rit;</span>
<span class="keywordflow">for</span>(++rit; rit!= ranks.end(); ++rit) {
<span class="keywordflow">if</span> ((*rit)>= coprimeR) <span class="keywordflow">break</span>;
<span class="keywordflow">for</span>(<span class="keywordtype">size_t</span> i=(*rit); i < coprimeR; ++i)
SmithDiagonal[i] *= sit->first;
<span class="comment">// modrank = *rit;</span>
}
}
}
<a class="code" href="group__integers.html#gad62eceb96963b157a2357aba991f6d6e" title="Integers in LinBox.">integer</a> si=1;
<span class="keywordtype">size_t</span> num=0;
<span class="keywordflow">for</span>( std::vector<integer>::const_iterator dit=SmithDiagonal.begin();
dit != SmithDiagonal.end(); ++dit) {
<span class="keywordflow">if</span> (*dit == si) ++num;
<span class="keywordflow">else</span> {
std::cerr << <span class="charliteral">'['</span> << si << <span class="charliteral">','</span> << num << <span class="stringliteral">"] "</span>;
num=1;
si = *dit;
}
}
std::cerr << <span class="charliteral">'['</span> << si << <span class="charliteral">','</span> << num << <span class="stringliteral">"] "</span> << std::endl;
chrono.stop();
std::cerr << chrono << std::endl;
<span class="keywordflow">return</span> 0;
}
</pre></div> </div>
</div>
<hr class="footer"/><address class="footer"><small>Generated on Tue Aug 30 2011 for linbox by 
<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.4 </small></address>
</body>
</html>
