<!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/solve.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/solve.C</div> </div>
</div>
<div class="contents">
<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/solve.C</span>
<span class="comment"> *</span>
<span class="comment"> * Copyright (C) 2005, 2010 J-G Dumas, D. Saunders, P. Giorgi</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="preprocessor">#include <iostream></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/blackbox/sparse.h"</span>
<span class="preprocessor">#include "linbox/solutions/solve.h"</span>
<span class="preprocessor">#include "linbox/util/matrix-stream.h"</span>
<span class="preprocessor">#include "<a class="code" href="methods_8h.html" title="NO DOC.">linbox/solutions/methods.h</a>"</span>
<span class="keyword">using namespace </span>LinBox;
<span class="keyword">using namespace </span>std;
<span class="keywordtype">int</span> main (<span class="keywordtype">int</span> argc, <span class="keywordtype">char</span> **argv)
{
commentator.<a name="a0"></a>setMaxDetailLevel (-1);
commentator.<a name="a1"></a>setMaxDepth (-1);
commentator.<a name="a2"></a>setReportStream (std::cerr);
<span class="keywordflow">if</span> (argc < 2 || argc > 4) {
cerr << <span class="stringliteral">"Usage: solve <matrix-file-in-supported-format> [<dense-vector-file>] [<p>]"</span> << endl;
<span class="keywordflow">return</span> 0;
}
srand48( BaseTimer::seed() );
std::ifstream input (argv[1]);
<span class="keywordflow">if</span> (!input) { cerr << <span class="stringliteral">"Error opening matrix file "</span> << argv[1] << endl; <span class="keywordflow">return</span> -1; }
std::ifstream invect;
<span class="keywordtype">bool</span> createB = <span class="keyword">false</span>;
<span class="keywordtype">int</span> ModComp = 0;
<span class="keywordflow">if</span> (argc == 2) {
createB = <span class="keyword">true</span>;
ModComp = 0;
}
<span class="keywordflow">if</span> (argc == 3) {
invect.open (argv[2], std::ifstream::in);
<span class="keywordflow">if</span> (!invect) {
createB = <span class="keyword">true</span>;
ModComp = 2;
}
<span class="keywordflow">else</span> {
createB = <span class="keyword">false</span>;
ModComp = 0;
}
}
<span class="keywordflow">if</span> (argc == 4) {
ModComp = 3;
invect.open (argv[2], std::ifstream::in);
<span class="keywordflow">if</span> (!invect) {
createB = <span class="keyword">true</span>;
}
<span class="keywordflow">else</span>
createB = <span class="keyword">false</span>;
}
<span class="keywordflow">if</span> (ModComp) {
<span class="keyword">typedef</span> <a name="_a3"></a><a class="code" href="class_lin_box_1_1_modular_3_01double_01_4.html" title="Standard representation of .">Modular<double></a> <a name="_a4"></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>;
<span class="keywordtype">double</span> q = atof(argv[ModComp]);
Field F(q);
<a name="_a5"></a><a class="code" href="class_lin_box_1_1_matrix_stream.html" title="MatrixStream.">MatrixStream< Field ></a> ms ( F, input );
<a name="_a6"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html" title="vector of sparse rows.">SparseMatrix<Field></a> A (ms); <span class="comment">// A.write(std::cout);</span>
cout << <span class="stringliteral">"A is "</span> << A.<a name="a7"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html#a0fee90e6b8ef7fcd8d76a7f3e1f268e6" title="Retreive row dimensions of Sparsemat matrix.">rowdim</a>() << <span class="stringliteral">" by "</span> << A.<a name="a8"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html#a32edb490d3597f5553152d14b102e227" title="Retreive column dimensions of Sparsemat matrix.">coldim</a>() << endl;
std::vector<Field::Element> X( A.<a class="code" href="class_lin_box_1_1_sparse_matrix.html#a32edb490d3597f5553152d14b102e227" title="Retreive column dimensions of Sparsemat matrix.">coldim</a>()),B(A.<a class="code" href="class_lin_box_1_1_sparse_matrix.html#a0fee90e6b8ef7fcd8d76a7f3e1f268e6" title="Retreive row dimensions of Sparsemat matrix.">rowdim</a>());
<span class="keywordflow">if</span> (createB) {
cerr << <span class="stringliteral">"Creating a random {-1,1} vector "</span> << endl;
std::vector<Field::Element> U( A.<a class="code" href="class_lin_box_1_1_sparse_matrix.html#a32edb490d3597f5553152d14b102e227" title="Retreive column dimensions of Sparsemat matrix.">coldim</a>() );
<span class="keywordflow">for</span>(std::vector<Field::Element>::iterator it=U.begin();
it != U.end(); ++it)
<span class="keywordflow">if</span> (drand48() <0.5)
F.<a name="a9"></a>init(*it,-1);
<span class="keywordflow">else</span>
F.init(*it,1);
A.<a name="a10"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html#aa37d57c030b1db0c198478268aa3cc15" title="Matrix-vector product .">apply</a>(B,U);
}
<span class="keywordflow">else</span> {
<span class="keywordflow">for</span>(std::vector<Field::Element>::iterator it=B.begin();
it != B.end(); ++it)
invect >> *it;
}
<span class="comment">// A.write(std::cout << "A: ") << std::endl;</span>
std::cout << <span class="stringliteral">"B is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=B.begin();it != B.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
Timer chrono;
<span class="comment">// Sparse Elimination</span>
chrono.clear();
chrono.start();
<a name="a11"></a><a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, A, B, <a class="code" href="struct_lin_box_1_1_method.html#abc071c5d851893868ab7a0c0d8ee34cb" title="Method::SparseElimination : no doc.">Method::SparseElimination</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Sparse Gauss) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=X.begin();it != X.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl<<std::endl;;
<span class="comment">// BlasElimination</span>
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, A, B, <a class="code" href="struct_lin_box_1_1_method.html#ae72ae4bbae70adf0514c807be2037706" title="Method::BlasElimination : no doc.">Method::BlasElimination</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(BlasElimination) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=X.begin();it != X.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl<< std::endl;
<span class="comment">// Wiedemann</span>
chrono.clear();
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, A, B, <a name="_a12"></a><a class="code" href="struct_lin_box_1_1_blackbox_specifier.html" title="BlackboxSpecifier.">Method::Blackbox</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Wiedemann) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=X.begin();it != X.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl<<std::endl;;
<span class="preprocessor">#if 0</span>
<span class="preprocessor"></span> <span class="comment">// Lanczos</span>
chrono.clear();
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, A, B, <a class="code" href="struct_lin_box_1_1_method.html#afd5f21592e0f3fae8a6197fb5d041e27" title="Method::Lanczos : no doc.">Method::Lanczos</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Lanczos) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=X.begin();it != X.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl<< std::endl;
<span class="comment">// Block Lanczos</span>
<a class="code" href="struct_lin_box_1_1_method.html#ac24179fa0ac59b75ab87d31275164bf0" title="Method::BlockLanczos : no doc.">Method::BlockLanczos</a> MBL;
MBL.preconditioner(Specifier::FULL_DIAGONAL);
chrono.clear();
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, A, B, MBL);
chrono.stop();
std::cout << <span class="stringliteral">"(Block Lanczos) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<Field::Element>::const_iterator it=X.begin();it != X.end(); ++it)
F.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl<< std::endl;
<span class="preprocessor">#endif</span>
<span class="preprocessor"></span>
}
<span class="keywordflow">else</span> {
<a name="_a13"></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 );
<a class="code" href="class_lin_box_1_1_sparse_matrix.html" title="vector of sparse rows.">SparseMatrix<PID_integer></a> A (ms);
PID_integer::Element d;
std::cout << <span class="stringliteral">"A is "</span> << A.<a name="a14"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html#a0fee90e6b8ef7fcd8d76a7f3e1f268e6" title="Retreive row dimensions of Sparsemat matrix.">rowdim</a>() << <span class="stringliteral">" by "</span> << A.<a name="a15"></a><a class="code" href="class_lin_box_1_1_sparse_matrix.html#a32edb490d3597f5553152d14b102e227" title="Retreive column dimensions of Sparsemat matrix.">coldim</a>() << std::endl;
std::vector<PID_integer::Element> X( A.<a class="code" href="class_lin_box_1_1_sparse_matrix.html#a32edb490d3597f5553152d14b102e227" title="Retreive column dimensions of Sparsemat matrix.">coldim</a>()),B(A.<a class="code" href="class_lin_box_1_1_sparse_matrix.html#a0fee90e6b8ef7fcd8d76a7f3e1f268e6" title="Retreive row dimensions of Sparsemat matrix.">rowdim</a>());
<span class="keywordflow">if</span> (createB) {
cerr << <span class="stringliteral">"Creating a random {-1,1} vector "</span> << endl;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::iterator it=B.begin();
it != B.end(); ++it)
<span class="keywordflow">if</span> (drand48() <0.5)
*it = -1;
<span class="keywordflow">else</span>
*it = 1;
}
<span class="keywordflow">else</span> {
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::iterator it=B.begin();
it != B.end(); ++it)
invect >> *it;
}
std::cout << <span class="stringliteral">"B is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=B.begin();
it != B.end(); ++it)
ZZ.<a name="a16"></a>write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"]"</span> << std::endl;
Timer chrono;
<span class="comment">// Wiedemann</span>
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, d, A, B, <a class="code" href="struct_lin_box_1_1_method.html#a95658bcef98e7d0868b2e35d7ce03392" title="Method::Wiedemann : no doc.">Method::Wiedemann</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Wiedemann) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=X.begin();it != X.end(); ++it)
ZZ.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"] / "</span>;
ZZ.write(std::cout, d) << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl;
<span class="comment">// BlasElimination</span>
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, d, A, B, <a class="code" href="struct_lin_box_1_1_method.html#ae72ae4bbae70adf0514c807be2037706" title="Method::BlasElimination : no doc.">Method::BlasElimination</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(BlasElimination) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=X.begin();it != X.end(); ++it)
ZZ.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"] / "</span>;
ZZ.write(std::cout, d)<< std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl;
<span class="comment">// Sparse Elimination</span>
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, d, A, B, <a class="code" href="struct_lin_box_1_1_method.html#abc071c5d851893868ab7a0c0d8ee34cb" title="Method::SparseElimination : no doc.">Method::SparseElimination</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(SparseElimination) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=X.begin();it != X.end(); ++it)
ZZ.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"] / "</span>;
ZZ.write(std::cout, d)<< std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl;
<span class="preprocessor">#if 0</span>
<span class="preprocessor"></span> <span class="comment">// Lanczos</span>
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, d, A, B, <a class="code" href="struct_lin_box_1_1_method.html#afd5f21592e0f3fae8a6197fb5d041e27" title="Method::Lanczos : no doc.">Method::Lanczos</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Lanczos) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=X.begin();it != X.end(); ++it)
ZZ.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"] / "</span>;
ZZ.write(std::cout, d) << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl;
<span class="comment">// Block Lanczos</span>
chrono.clear();
chrono.start();
<a class="code" href="namespace_lin_box.html#a353c2c55c4d1fd5b35e9b4e04c274982" title="Solve Ax=b over field F using Wiedemann's method, with inconsistency certificate.">solve</a> (X, d, A, B, <a class="code" href="struct_lin_box_1_1_method.html#ac24179fa0ac59b75ab87d31275164bf0" title="Method::BlockLanczos : no doc.">Method::BlockLanczos</a>());
chrono.stop();
std::cout << <span class="stringliteral">"(Block Lanczos) Solution is ["</span>;
<span class="keywordflow">for</span>(std::vector<PID_integer::Element>::const_iterator it=X.begin();it != X.end(); ++it)
ZZ.write(cout, *it) << <span class="stringliteral">" "</span>;
std::cout << <span class="stringliteral">"] / "</span>;
ZZ.write(std::cout, d) << std::endl;
std::cout << <span class="stringliteral">"CPU time (seconds): "</span> << chrono.usertime() << std::endl;
<span class="preprocessor">#endif</span>
<span class="preprocessor"></span> }
<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>
