<?xml version="1.0" encoding="iso-8859-1" ?> 
<?xml-stylesheet type="text/css" href="hsl_ma42.css"?> 
<?xml-stylesheet type="text/xsl" href="pmathml.xsl"?> 
<!--http://www.w3.org/Math/XSL/pmathml.xsl--> 
<html  
xmlns="http://www.w3.org/1999/xhtml"  
><head>
   <title>HSL_MA42 &#x2014; Sparse unsymmetric system: out-of-core frontal method</title> 
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1" /> 
<meta name="generator" content="TeX4ht (http://www.cse.ohio-state.edu/~gurari/TeX4ht/)" /> 
<meta name="originator" content="TeX4ht (http://www.cse.ohio-state.edu/~gurari/TeX4ht/)" /> 
<!-- html,pmathml,xhtml,mozilla --> 
<meta name="src" content="hsl_ma42.tex" /> 
<meta name="date" content="2011-02-08 15:26:00" /> 
<link rel="stylesheet" type="text/css" href="hsl_ma42.css" /> 
<link rel="stylesheet" href="/css/hsl.css" type="text/css" /><link rel="stylesheet" href="/css/catalog.css" type="text/css" /></head><body 
>
   <div id="body"> 
<div id="header">  
<div class="w100"> 
<div class="stfc-logo"> 
<a accesskey="1" tabindex="2" href="http://www.stfc.ac.uk"> 
<img alt="STFC Home Page" src="/images/stfc-logo.jpg"/> </a> 
<a class="hide_tab" accesskey="S" tabindex="1" href="#pgContent">Skip navigation</a> 
</div> 
<div class="clear"></div> 
</div> 
<div class="w100">  
<div class="RCUK"> 
<a href="http://framework.rcuk.ac.uk/" target="_blank" title="RCUK website - opens in a new window">Excellence with Impact - part of Research Councils UK</a> 
</div> 
<div class="header-links"> 
<a href="index.html">HSL Home</a> : <a href="/contact.html">Contact HSL</a> 
</div> 
<div class="clear"></div> 
</div> 
</div> 
  <div id="menu"> 
<ul> 
<li> 
<a href="http://www.stfc.ac.uk">STFC home</a> 
<span class="mTags"><br /> 
The Science and Technology Facilities Council 
</span> 
</li> 
<li> 
<a href="/index.html">HSL home</a> 
<span class="mTags"><br /> 
About the HSL Mathematical Software Library 
</span> 
</li> 
<li> 
<a href="/catalogue/">HSL 2011</a> 
<span class="mTags"><br /> 
Our current release 
</span> 
</li> 
<li> 
<a href="/archive/index.html"> 
HSL archive</a> 
<span class="mTags"><br /> 
Older packages available at no cost 
</span> 
</li> 
<li> 
<a href="/catalogue/"> 
Academic access</a> 
<span class="mTags"><br /> 
Free access for teaching and research 
</span> 
</li> 
<li> 
<a href="http://www.cse.scitech.ac.uk/nag/index.shtml"> 
The Numerical Analysis Group</a> 
<span class="mTags"><br /> 
About the research group behind HSL 
</span> 
</li> 
<!--<li> 
<a href="http://galahad.rl.ac.uk/">GALAHAD</a> 
<span class="mTags"><br /> 
Our optimization software library 
</span> 
</li>--> 
<li> 
<a href="/contact.html">Contact us</a> 
<span class="mTags"><br /> 
Our contact details 
</span> 
</li> 
</ul> 
</div> 
 <div id="bodytext"> 
<a name="pgContent"></a> 
 <div> 
<form action="http://www.hsl.rl.ac.uk/search.html" id="cse-search-box"> 
<div> 
<input type="hidden" name="cx" value="011781803477211733479:s9roibysody" /> 
<input type="hidden" name="cof" value="FORID:11" /> 
<input type="hidden" name="ie" value="UTF-8" /> 
<input type="text" name="q" size="31" /> 
<input type="submit" name="sa" value="Search" /> 
</div> 
 
</form> 
<script type="text/javascript" src="http://www.google.com/cse/brand?form=cse-search-box&amp;lang=en"></script> 
</div> 
 <div class="download"> 
 <h2>Version 1.2.0</h2> 
 <p class="vdate"></p> <p class="doclink"><a href="/specs/hsl_ma42.pdf">User documentation</a></p> <p>Code Download</p> 
 <ul> 
  
  
 <li><a href="/download/HSL_MA42/1.2.0/s">Single</a></li> 

<li><a href="/download/HSL_MA42/1.2.0/d">Double</a></li> 
  
  
  </ul> 
</div> 

   <h3 class="likesectionHead"><a 
 id="x1-1000"></a>HSL_MA42 Sparse unsymmetric system: out-of-core frontal method</h3>
<!--l. 27--><p class="noindent" >The module <span 
class="cmtt-10">HSL</span><span 
class="cmtt-10">_MA42 </span>solves one or more sets of sparse linear equations,
<!--l. 29--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><mstyle mathvariant="bold"><mi 
>A</mi><mi 
>x</mi></mstyle> <mo 
class="MathClass-rel">=</mo> <mstyle mathvariant="bold"><mi 
>b</mi></mstyle></math> or
<!--l. 29--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><msup><mrow 
><mstyle mathvariant="bold"><mi 
>A</mi></mstyle></mrow><mrow 
><mi 
>T</mi></mrow></msup 
><mstyle mathvariant="bold"><mi 
>x</mi></mstyle> <mo 
class="MathClass-rel">=</mo> <mstyle mathvariant="bold"><mi 
>b</mi></mstyle></math>, by the
frontal method, optionally using direct-access files for the matrix factors. Use is made of
high level BLAS kernels. The code has low in-core memory requirements. The matrix
<!--l. 32--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><mstyle mathvariant="bold"><mi 
>A</mi></mstyle></math> may
be input by the user in either of the following ways:
</p><!--l. 36--><p class="indent" >
     </p><dl class="description"><dt class="description">
 </dt><dd 
class="description">(i) by elements in a finite-element calculation,
     </dd><dt class="description">
 </dt><dd 
class="description">(ii) by equations (matrix rows).</dd></dl>
<!--l. 41--><p class="noindent" >In both cases, the coefficient matrix and right-hand side(s) are of the form
</p><!--l. 44--><p class="indent" >
<!--tex4ht:inline--></p><!--l. 44--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="block" >
                 <mstyle mathvariant="bold"><mi 
>A</mi></mstyle> <mo 
class="MathClass-rel">=</mo><munderover accentunder="false" accent="false"><mrow  
><mo mathsize="big" 
> &#x2211;</mo>
  </mrow><mrow 
><mi 
>k</mi><mo 
class="MathClass-rel">=</mo><mn>1</mn></mrow><mrow 
><mi 
>m</mi></mrow></munderover 
><msup><mrow 
><mstyle mathvariant="bold"><mi 
>A</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
><mo 
class="MathClass-punc">,</mo><mspace width="2em" class="qquad"/><mstyle mathvariant="bold"><mi 
>b</mi></mstyle> <mo 
class="MathClass-rel">=</mo><munderover accentunder="false" accent="false"><mrow  
><mo mathsize="big" 
> &#x2211;</mo>
</mrow><mrow 
><mi 
>k</mi><mo 
class="MathClass-rel">=</mo><mn>1</mn></mrow><mrow 
><mi 
>m</mi></mrow></munderover 
><msup><mrow 
><mstyle mathvariant="bold"><mi 
>b</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
><mo 
class="MathClass-punc">.</mo>
</math>
<!--l. 45--><p class="nopar" >
</p><!--l. 47--><p class="indent" >   In case (i), the summation is over finite elements. The matrix
<!--l. 47--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><msup><mrow 
><mstyle mathvariant="bold"><mi 
>A</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
></math> is
nonzero only in those rows and columns which correspond to variables in the
<!--l. 49--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><mi 
>k</mi></math>-th element.
The vector <!--l. 49--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><msup><mrow 
><mstyle mathvariant="bold"><mi 
>b</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
></math>
is nonzero only in those rows which correspond to variables in element
<!--l. 51--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><mi 
>k</mi></math>.

</p><!--l. 53--><p class="indent" >   In case (ii), the summation is over equations and
<!--l. 53--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><msup><mrow 
><mstyle mathvariant="bold"><mi 
>A</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
></math> and
<!--l. 54--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><msup><mrow 
><mstyle mathvariant="bold"><mi 
>b</mi></mstyle></mrow><mrow 
><mrow ><mo 
class="MathClass-open">(</mo><mrow><mi 
>k</mi></mrow><mo 
class="MathClass-close">)</mo></mrow></mrow></msup 
></math> are nonzero
only in row <!--l. 54--><math 
 xmlns="http://www.w3.org/1998/Math/MathML" display="inline" ><mi 
>k</mi></math>.
</p>
    </div> <!-- bodytext -->  <div class="clear"></div> 
 
<div class="footer"> 
<div class="rcuk-footer"> 
<a href="http://www.rcuk.ac.uk" > 
<img src="/images/rcuk-footer.gif" alt="Click to visit RCUK web site"/> 
 
</a> 
</div> 
 
<div class="footer-links"> 
 
<a href="/site_map.html" accesskey="3">SITE-MAP</a> :  
<a href="http://www.stfc.ac.uk/About/Accessibility.aspx" accesskey="0">ACCESSIBILITY</a> :  
<a href="http://www.stfc.ac.uk/privacy.aspx">PRIVACY POLICY</a> :  
<a href="http://www.stfc.ac.uk/access.aspx">ACCESS TO INFORMATION</a> :  
<a href="http://www.stfc.ac.uk/disclaimer.aspx" accesskey="5">TERMS OF USE</a> : 
<a accesskey="9" href="mailto:hsl@stfc.ac.uk?subject=Feedback From Website">WEBMASTER</a> 
 
<span class="hide">|<a accesskey="7" tabindex="7" href="http://www.stfc.ac.uk/About/pol/sfirst/contents.aspx">Complaints Procedure</a></span> 
<span class="hide">|<a accesskey="2" tabindex="8" href="http://www.stfc.ac.uk/PMC/News.aspx">News</a></span> 
</div> 
</div> 
 <div class="clear"></div> 
</div> <!-- body --> 
  
</body> 
</html> 




