<!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: RationalSolver< Ring, Field, RandomPrime, DixonTraits > Class Template Reference</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 class="current"><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 id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="annotated.html"><span>Data Structures</span></a></li>
<li><a href="hierarchy.html"><span>Class Hierarchy</span></a></li>
<li><a href="functions.html"><span>Data Fields</span></a></li>
</ul>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="namespace_lin_box.html">LinBox</a> </li>
<li class="navelem"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html">RationalSolver< Ring, Field, RandomPrime, DixonTraits ></a> </li>
</ul>
</div>
</div>
<div class="header">
<div class="summary">
<a href="#pub-methods">Public Member Functions</a> </div>
<div class="headertitle">
<div class="title">RationalSolver< Ring, Field, RandomPrime, DixonTraits > Class Template Reference</div> </div>
</div>
<div class="contents">
<!-- doxytag: class="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >" -->
<p>partial specialization of p-adic based solver with Dixon algorithm.
<a href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#details">More...</a></p>
<p><code>#include <rational-solver.h></code></p>
<table class="memberdecls">
<tr><td colspan="2"><h2><a name="pub-methods"></a>
Public Member Functions</h2></td></tr>
<tr><td class="memItemLeft" align="right" valign="top"> </td><td class="memItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#af99cee064d4c6582f2eee75fb1512471">RationalSolver</a> (const Ring &r=Ring(), const RandomPrime &rp=RandomPrime(DEFAULT_PRIMESIZE))</td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Constructor. <a href="#af99cee064d4c6582f2eee75fb1512471"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top"> </td><td class="memItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#ae9718565677866b3fce935bd0f59821d">RationalSolver</a> (const Prime &p, const Ring &r=Ring(), const RandomPrime &rp=RandomPrime(DEFAULT_PRIMESIZE))</td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Constructor, trying the prime p first. <a href="#ae9718565677866b3fce935bd0f59821d"></a><br/></td></tr>
<tr><td class="memTemplParams" colspan="2">template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#a71a1d5168a3bbad0055ae8d786a43cc7">solve</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, const bool s=false, const int maxPrimes=DEFAULT_MAXPRIMES, const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> level=SL_DEFAULT) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Solve a linear system <code>Ax=b</code> over quotient field of a ring. <a href="#a71a1d5168a3bbad0055ae8d786a43cc7"></a><br/></td></tr>
<tr><td class="memTemplParams" colspan="2"><a class="anchor" id="a21d12adbd6c9fc8cb7435eadd1efa11c"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::solve" ref="a21d12adbd6c9fc8cb7435eadd1efa11c" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, const int maxPrimes, const SolverLevel level=SL_DEFAULT) const " -->
template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#a21d12adbd6c9fc8cb7435eadd1efa11c">solve</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, const int maxPrimes, const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> level=SL_DEFAULT) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">overload so that the bool 'oldMatrix' argument is not accidentally set to true <br/></td></tr>
<tr><td class="memTemplParams" colspan="2">template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#abeb85fb68910ed9443426ab5f431fd2a">solveNonsingular</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, bool s=false, int maxPrimes=DEFAULT_MAXPRIMES) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Solve a nonsingular, square linear system <code>Ax=b</code> over quotient field of a ring. <a href="#abeb85fb68910ed9443426ab5f431fd2a"></a><br/></td></tr>
<tr><td class="memTemplParams" colspan="2">template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#a588490938a24f50392eeeb0ac9680964">solveSingular</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, int maxPrimes=DEFAULT_MAXPRIMES, const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> level=SL_DEFAULT) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Solve a general rectangular linear system <code>Ax=b</code> over quotient field of a ring. <a href="#a588490938a24f50392eeeb0ac9680964"></a><br/></td></tr>
<tr><td class="memTemplParams" colspan="2">template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#a8c8ead8738c125c3fc186cb11ddb8f54">findRandomSolution</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, int maxPrimes=DEFAULT_MAXPRIMES, const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> level=SL_DEFAULT) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Find a random solution of the general linear system <code>Ax=b</code> over quotient field of a ring. <a href="#a8c8ead8738c125c3fc186cb11ddb8f54"></a><br/></td></tr>
<tr><td class="memTemplParams" colspan="2">template<class IMatrix , class Vector1 , class Vector2 > </td></tr>
<tr><td class="memTemplItemLeft" align="right" valign="top"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> </td><td class="memTemplItemRight" valign="bottom"><a class="el" href="class_lin_box_1_1_rational_solver_3_01_ring_00_01_field_00_01_random_prime_00_01_dixon_traits_01_4.html#a2ee8508ff30ab03a0b2d816e72a291be">monolithicSolve</a> (Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, bool makeMinDenomCert, bool randomSolution, int maxPrimes=DEFAULT_MAXPRIMES, const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> level=SL_DEFAULT) const </td></tr>
<tr><td class="mdescLeft"> </td><td class="mdescRight">Big solving routine to perform random solving and certificate generation. <a href="#a2ee8508ff30ab03a0b2d816e72a291be"></a><br/></td></tr>
</table>
<hr/><a name="details" id="details"></a><h2>Detailed Description</h2>
<div class="textblock"><h3>template<class Ring, class Field, class RandomPrime><br/>
class LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits ></h3>
<p>partial specialization of p-adic based solver with Dixon algorithm. </p>
<p>See the following reference for details on this algorithm: </p>
<dl class="bib"><dt><b><a class="el" href="bib.html#_bib000011">Bibliography:</a></b></dt><dd></dd></dl>
<p>- John D. Dixon <em>Exact Solution of linear equations using p-adic expansions</em>. Numerische Mathematik, volume 40, pages 137-141, 1982. </p>
</div><hr/><h2>Constructor & Destructor Documentation</h2>
<a class="anchor" id="af99cee064d4c6582f2eee75fb1512471"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::RationalSolver" ref="af99cee064d4c6582f2eee75fb1512471" args="(const Ring &r=Ring(), const RandomPrime &rp=RandomPrime(DEFAULT_PRIMESIZE))" -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="class_lin_box_1_1_rational_solver.html">RationalSolver</a> </td>
<td>(</td>
<td class="paramtype">const Ring & </td>
<td class="paramname"><em>r</em> = <code>Ring()</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const RandomPrime & </td>
<td class="paramname"><em>rp</em> = <code>RandomPrime(DEFAULT_PRIMESIZE)</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td><code> [inline]</code></td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Constructor. </p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">r</td><td>a Ring, set by default </td></tr>
<tr><td class="paramname">rp</td><td>a RandomPrime generator, set by default </td></tr>
</table>
</dd>
</dl>
</div>
</div>
<a class="anchor" id="ae9718565677866b3fce935bd0f59821d"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::RationalSolver" ref="ae9718565677866b3fce935bd0f59821d" args="(const Prime &p, const Ring &r=Ring(), const RandomPrime &rp=RandomPrime(DEFAULT_PRIMESIZE))" -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="class_lin_box_1_1_rational_solver.html">RationalSolver</a> </td>
<td>(</td>
<td class="paramtype">const Prime & </td>
<td class="paramname"><em>p</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Ring & </td>
<td class="paramname"><em>r</em> = <code>Ring()</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const RandomPrime & </td>
<td class="paramname"><em>rp</em> = <code>RandomPrime(DEFAULT_PRIMESIZE)</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td><code> [inline]</code></td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Constructor, trying the prime p first. </p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">p</td><td>a Prime </td></tr>
<tr><td class="paramname">r</td><td>a Ring, set by default </td></tr>
<tr><td class="paramname">rp</td><td>a RandomPrime generator, set by default </td></tr>
</table>
</dd>
</dl>
</div>
</div>
<hr/><h2>Member Function Documentation</h2>
<a class="anchor" id="a71a1d5168a3bbad0055ae8d786a43cc7"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::solve" ref="a71a1d5168a3bbad0055ae8d786a43cc7" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, const bool s=false, const int maxPrimes=DEFAULT_MAXPRIMES, const SolverLevel level=SL_DEFAULT) const " -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> solve </td>
<td>(</td>
<td class="paramtype">Vector1 & </td>
<td class="paramname"><em>num</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">Integer & </td>
<td class="paramname"><em>den</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const IMatrix & </td>
<td class="paramname"><em>A</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Vector2 & </td>
<td class="paramname"><em>b</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const bool </td>
<td class="paramname"><em>s</em> = <code>false</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const int </td>
<td class="paramname"><em>maxPrimes</em> = <code>DEFAULT_MAXPRIMES</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> </td>
<td class="paramname"><em>level</em> = <code>SL_DEFAULT</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td> const</td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Solve a linear system <code>Ax=b</code> over quotient field of a ring. </p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">num</td><td>Vector of numerators of the solution </td></tr>
<tr><td class="paramname">den</td><td>The common denominator. 1/den * num is the rational solution of <code>Ax=b</code>. </td></tr>
<tr><td class="paramname">A</td><td>Matrix of linear system </td></tr>
<tr><td class="paramname">b</td><td>Right-hand side of system </td></tr>
<tr><td class="paramname">s</td><td></td></tr>
<tr><td class="paramname">maxPrimes</td><td>maximum number of moduli to try </td></tr>
<tr><td class="paramname">level</td><td>level of certification to be used</td></tr>
</table>
</dd>
</dl>
<dl class="return"><dt><b>Returns:</b></dt><dd>status of solution. if <code></code>(return != SS_FAILED), and <code></code>(level >= SL_LASVEGAS), solution is guaranteed correct. <code>SS_FAILED</code> - all primes used were bad <code>SS_OK</code> - solution found. <code>SS_INCONSISTENT</code> - system appreared inconsistent. certificate is in <code>lastCertificate</code> if <code></code>(level >= SL_CERTIFIED) </dd></dl>
</div>
</div>
<a class="anchor" id="abeb85fb68910ed9443426ab5f431fd2a"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::solveNonsingular" ref="abeb85fb68910ed9443426ab5f431fd2a" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, bool s=false, int maxPrimes=DEFAULT_MAXPRIMES) const " -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> solveNonsingular </td>
<td>(</td>
<td class="paramtype">Vector1 & </td>
<td class="paramname"><em>num</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">Integer & </td>
<td class="paramname"><em>den</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const IMatrix & </td>
<td class="paramname"><em>A</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Vector2 & </td>
<td class="paramname"><em>b</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">bool </td>
<td class="paramname"><em>s</em> = <code>false</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">int </td>
<td class="paramname"><em>maxPrimes</em> = <code>DEFAULT_MAXPRIMES</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td> const</td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Solve a nonsingular, square linear system <code>Ax=b</code> over quotient field of a ring. </p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">num</td><td>Vector of numerators of the solution </td></tr>
<tr><td class="paramname">den</td><td>The common denominator. <code>1/den * num</code> is the rational solution of <code>Ax = b</code> </td></tr>
<tr><td class="paramname">A</td><td>Matrix of linear system (it must be square) </td></tr>
<tr><td class="paramname">b</td><td>Right-hand side of system </td></tr>
<tr><td class="paramname">s</td><td>unused </td></tr>
<tr><td class="paramname">maxPrimes</td><td>maximum number of moduli to try</td></tr>
</table>
</dd>
</dl>
<dl class="return"><dt><b>Returns:</b></dt><dd>status of solution :<ul>
<li><code>SS_FAILED</code> all primes used were bad;</li>
<li><code>SS_OK</code> solution found, guaranteed correct;</li>
<li><code>SS_SINGULAR</code> system appreared singular mod all primes.</li>
</ul>
</dd></dl>
</div>
</div>
<a class="anchor" id="a588490938a24f50392eeeb0ac9680964"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::solveSingular" ref="a588490938a24f50392eeeb0ac9680964" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, int maxPrimes=DEFAULT_MAXPRIMES, const SolverLevel level=SL_DEFAULT) const " -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> solveSingular </td>
<td>(</td>
<td class="paramtype">Vector1 & </td>
<td class="paramname"><em>num</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">Integer & </td>
<td class="paramname"><em>den</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const IMatrix & </td>
<td class="paramname"><em>A</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Vector2 & </td>
<td class="paramname"><em>b</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">int </td>
<td class="paramname"><em>maxPrimes</em> = <code>DEFAULT_MAXPRIMES</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> </td>
<td class="paramname"><em>level</em> = <code>SL_DEFAULT</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td> const</td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Solve a general rectangular linear system <code>Ax=b</code> over quotient field of a ring. </p>
<p>If A is known to be square and nonsingular, calling solveNonsingular is more efficient.</p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">num</td><td>Vector of numerators of the solution </td></tr>
<tr><td class="paramname">den</td><td>The common denominator. <code>1/den * num</code> is the rational solution of <code>Ax = b</code> </td></tr>
<tr><td class="paramname">A</td><td>Matrix of linear system </td></tr>
<tr><td class="paramname">b</td><td>Right-hand side of system </td></tr>
<tr><td class="paramname">maxPrimes</td><td>maximum number of moduli to try </td></tr>
<tr><td class="paramname">level</td><td>level of certification to be used</td></tr>
</table>
</dd>
</dl>
<dl class="return"><dt><b>Returns:</b></dt><dd>status of solution. if <code>(return != SS_FAILED)</code>, and <code>(level >= SL_LASVEGAS)</code>, solution is guaranteed correct.<ul>
<li><code>SS_FAILED</code> all primes used were bad</li>
<li><code>SS_OK</code> solution found.</li>
<li><code>SS_INCONSISTENT</code> system appreared inconsistent. certificate is in <code>lastCertificate</code> if <code>(level >= SL_CERTIFIED)</code> </li>
</ul>
</dd></dl>
</div>
</div>
<a class="anchor" id="a8c8ead8738c125c3fc186cb11ddb8f54"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::findRandomSolution" ref="a8c8ead8738c125c3fc186cb11ddb8f54" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, int maxPrimes=DEFAULT_MAXPRIMES, const SolverLevel level=SL_DEFAULT) const " -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> findRandomSolution </td>
<td>(</td>
<td class="paramtype">Vector1 & </td>
<td class="paramname"><em>num</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">Integer & </td>
<td class="paramname"><em>den</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const IMatrix & </td>
<td class="paramname"><em>A</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Vector2 & </td>
<td class="paramname"><em>b</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">int </td>
<td class="paramname"><em>maxPrimes</em> = <code>DEFAULT_MAXPRIMES</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> </td>
<td class="paramname"><em>level</em> = <code>SL_DEFAULT</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td> const</td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Find a random solution of the general linear system <code>Ax=b</code> over quotient field of a ring. </p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">num</td><td>Vector of numerators of the solution </td></tr>
<tr><td class="paramname">den</td><td>The common denominator. <code>1/den * num</code> is the rational solution of <code>Ax = b</code>. </td></tr>
<tr><td class="paramname">A</td><td>Matrix of linear system </td></tr>
<tr><td class="paramname">b</td><td>Right-hand side of system </td></tr>
<tr><td class="paramname">maxPrimes</td><td>maximum number of moduli to try </td></tr>
<tr><td class="paramname">level</td><td>level of certification to be used</td></tr>
</table>
</dd>
</dl>
<dl class="return"><dt><b>Returns:</b></dt><dd>status of solution. if <code>(return != SS_FAILED)</code>, and <code>(level >= SL_LASVEGAS)</code>, solution is guaranteed correct.<ul>
<li><code>SS_FAILED</code> all primes used were bad</li>
<li><code>SS_OK</code> solution found.</li>
<li><code>SS_INCONSISTENT</code> system appreared inconsistent. certificate is in lastCertificate if <code>(level >= SL_CERTIFIED)</code> </li>
</ul>
</dd></dl>
</div>
</div>
<a class="anchor" id="a2ee8508ff30ab03a0b2d816e72a291be"></a><!-- doxytag: member="LinBox::RationalSolver< Ring, Field, RandomPrime, DixonTraits >::monolithicSolve" ref="a2ee8508ff30ab03a0b2d816e72a291be" args="(Vector1 &num, Integer &den, const IMatrix &A, const Vector2 &b, bool makeMinDenomCert, bool randomSolution, int maxPrimes=DEFAULT_MAXPRIMES, const SolverLevel level=SL_DEFAULT) const " -->
<div class="memitem">
<div class="memproto">
<table class="memname">
<tr>
<td class="memname"><a class="el" href="group__padic.html#gabf95d1fb4fcd12ca89717f0be14212a1">SolverReturnStatus</a> monolithicSolve </td>
<td>(</td>
<td class="paramtype">Vector1 & </td>
<td class="paramname"><em>num</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">Integer & </td>
<td class="paramname"><em>den</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const IMatrix & </td>
<td class="paramname"><em>A</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const Vector2 & </td>
<td class="paramname"><em>b</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">bool </td>
<td class="paramname"><em>makeMinDenomCert</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">bool </td>
<td class="paramname"><em>randomSolution</em>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">int </td>
<td class="paramname"><em>maxPrimes</em> = <code>DEFAULT_MAXPRIMES</code>, </td>
</tr>
<tr>
<td class="paramkey"></td>
<td></td>
<td class="paramtype">const <a class="el" href="group__padic.html#gaee2ff986111a9d28c71343fbdc651a9f">SolverLevel</a> </td>
<td class="paramname"><em>level</em> = <code>SL_DEFAULT</code> </td>
</tr>
<tr>
<td></td>
<td>)</td>
<td></td><td> const</td>
</tr>
</table>
</div>
<div class="memdoc">
<p>Big solving routine to perform random solving and certificate generation. </p>
<p>Same arguments and return as findRandomSolution, except</p>
<dl><dt><b>Parameters:</b></dt><dd>
<table class="params">
<tr><td class="paramname">num</td><td>Vector of numerators of the solution </td></tr>
<tr><td class="paramname">den</td><td>The common denominator. <code>1/den * num</code> is the rational solution of <code>Ax = b</code> </td></tr>
<tr><td class="paramname">A</td><td></td></tr>
<tr><td class="paramname">b</td><td></td></tr>
<tr><td class="paramname">randomSolution</td><td>parameter to determine whether to randomize or not (since solveSingular calls this function as well) </td></tr>
<tr><td class="paramname">makeMinDenomCert</td><td>determines whether a partial certificate for the minimal denominator of a rational solution is made </td></tr>
<tr><td class="paramname">maxPrimes</td><td></td></tr>
<tr><td class="paramname">level</td><td>When <code>(randomSolution == true && makeMinDenomCert == true)</code>,<ul>
<li>If <code>(level == SL_MONTECARLO)</code> this function has the same effect as calling findRandomSolution.</li>
<li>If <code>(level >= SL_LASVEGAS && return == SS_OK)</code>, <code>lastCertifiedDenFactor</code> contains a certified factor of the min-solution's denominator.</li>
<li>If <code>(level >= SL_CERTIFIED && return == SS_OK)</code>, <code>lastZBNumer</code> and <code>lastCertificate</code> are updated as well. </li>
</ul>
</td></tr>
</table>
</dd>
</dl>
</div>
</div>
<hr/>The documentation for this class was generated from the following files:<ul>
<li><a class="el" href="rational-solver_8h.html">rational-solver.h</a></li>
<li>rational-solver.inl</li>
</ul>
</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>
