<!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"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<title>CPPLapack: _dgbmatrix Class Reference</title>

<link href="tabs.css" rel="stylesheet" type="text/css"/>
<link href="doxygen.css" rel="stylesheet" type="text/css" />

<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
  $(document).ready(function() { searchBox.OnSelectItem(0); });
</script>

</head>
<body>
<div id="top"><!-- do not remove this div! -->


<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  
  
  <td style="padding-left: 0.5em;">
   <div id="projectname">CPPLapack
   
   </div>
   
  </td>
  
  
  
 </tr>
 </tbody>
</table>
</div>

<!-- Generated by Doxygen 1.7.6.1 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
  <div id="navrow1" class="tabs">
    <ul class="tablist">
      <li><a href="index.html"><span>Main&#160;Page</span></a></li>
      <li><a href="pages.html"><span>Related&#160;Pages</span></a></li>
      <li class="current"><a href="annotated.html"><span>Classes</span></a></li>
      <li><a href="files.html"><span>Files</span></a></li>
      <li><a href="dirs.html"><span>Directories</span></a></li>
      <li>
        <div id="MSearchBox" class="MSearchBoxInactive">
        <span class="left">
          <img id="MSearchSelect" src="search/mag_sel.png"
               onmouseover="return searchBox.OnSearchSelectShow()"
               onmouseout="return searchBox.OnSearchSelectHide()"
               alt=""/>
          <input type="text" id="MSearchField" value="Search" accesskey="S"
               onfocus="searchBox.OnSearchFieldFocus(true)" 
               onblur="searchBox.OnSearchFieldFocus(false)" 
               onkeyup="searchBox.OnSearchFieldChange(event)"/>
          </span><span class="right">
            <a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
          </span>
        </div>
      </li>
    </ul>
  </div>
  <div id="navrow2" class="tabs2">
    <ul class="tablist">
      <li><a href="annotated.html"><span>Class&#160;List</span></a></li>
      <li><a href="functions.html"><span>Class&#160;Members</span></a></li>
    </ul>
  </div>
</div>
<div class="header">
  <div class="summary">
<a href="#pub-methods">Public Member Functions</a> &#124;
<a href="#pub-attribs">Public Attributes</a> &#124;
<a href="#friends">Friends</a>  </div>
  <div class="headertitle">
<div class="title">_dgbmatrix Class Reference</div>  </div>
</div><!--header-->
<div class="contents">
<!-- doxytag: class="_dgbmatrix" -->
<p>(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class  
 <a href="class__dgbmatrix.html#details">More...</a></p>

<p><code>#include &lt;<a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>&gt;</code></p>

<p><a href="class__dgbmatrix-members.html">List of all members.</a></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">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a175eaab594a9f9f1d88330f3e511d8aa">_dgbmatrix</a> ()</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a54b165895369496cb5c263b55ace48d5">_dgbmatrix</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a3d238d72484bdc8ab3051d6e27710f8b">~_dgbmatrix</a> ()</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__zgbmatrix.html">_zgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a300545260e3715527e79173b1291e883">to_zgbmatrix</a> () const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a982423ea885f6583208d05625c31b03e">to_dgematrix</a> () const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top">double &amp;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a0bb6b797e67d9b5aef5c25b6f3de03e2">operator()</a> (const long &amp;, const long &amp;) const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#af59336a13c233cf8af0c54d170aa0f22">write</a> (const char *) const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a736c8efcb6a6f9fd458090aa468ee2ce">nullify</a> () const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a> () const </td></tr>
<tr><td colspan="2"><h2><a name="pub-attribs"></a>
Public Attributes</h2></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">long&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a">m</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">matrix row size  <a href="#a587a07ecbbfc6bbdb0e05574ea6b088a"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">long&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971">n</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">matrix column size  <a href="#a95aebc7e2d4ea66b1b7a08e981df6971"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">long&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c">kl</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">lower band width  <a href="#aeea4fdf76630ca1b51f93c5bed64662c"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">long&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea">ku</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">upper band width  <a href="#a116dadc8ed5246ee0b4cd0a9fffc53ea"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">double *&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb">array</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">1D array to store matrix data  <a href="#a3d7680e0014aa00249f76a86e13de9eb"></a><br/></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">double **&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d">darray</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">array of pointers of column head addresses  <a href="#ae9c142fd6903cf00df8e1d5363e9848d"></a><br/></td></tr>
<tr><td colspan="2"><h2><a name="friends"></a>
Friends</h2></td></tr>
<tr><td class="memItemLeft" align="right" valign="top">std::ostream &amp;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#abc9cd92dfb30bf77e6e98d7c910228d8">operator&lt;&lt;</a> (std::ostream &amp;, const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#abc0db5fc9798b89672e2e103a07b097b">t</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a45ae3ef397cc029be013792ce75471c6">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a8feabd36261c12eea77cadbc33916efe">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ae969cf19039e99763b4ef7e5bce4e7e6">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a46e30066362c28d9acea84681280212c">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ab15b39bacb422a2428295c7bc76994e1">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aef7ec311adfe90b7c26c0e92f8aebdd1">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a3e98bdfeb5d939881272054ad728a175">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aadaa4a751fb7153ca4a18679e1b235b7">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a4c110f12272b94d105a7b40b9c32626c">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a1899129fe441186353f81ad1dbcf6a35">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a0eaee3557a07d4103ed7654583b9845e">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aff572701e2dee5c2206c81798d72a6dd">operator+</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aeb9361c4eec76c4964258446cd6bc24c">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a6212a92d6cd140914d20b82fd8620856">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a01bc3a5f90fcbf23476c5427f3ae508c">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a4f919fea9e763b1afb1852ee711137c5">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aae7d46d5ff9708363386b1a041fce336">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a209b6000adc349cff070ee25ba1d8941">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aaf20d6c83d831430cc819a3b670393aa">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a6039af37d6b1e3ac7444f8b7475f5a02">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ad8a2c801718b3570fae41e3591886095">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a828552e4f2601fb88f6f62493a8b7e2d">operator-</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dcovector.html">_dcovector</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a0f435627957fb3957ce083a40e4a2df0">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdcovector.html">dcovector</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dcovector.html">_dcovector</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a929a817aaf554b47204b2fff2058af16">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dcovector.html">_dcovector</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a66c6987c5521e76467a4c149545aeb2f">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a7d45270990d9f629335bfe565ae81876">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ad08c9cddd6aa158109e227ba907ba9ce">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#aadf4f8e6254cd0f581d2fd1ad034191f">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a2ea7426af09460cc63d8127658c47909">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a6338b568730f34d9ebe9a69871b1b797">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ab255cc79e2cd0ee188e1598f6decf672">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#acee43528df478ebfaa53a2fc314e19a0">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#ae416997168ff0b923f4fd006fe034908">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#af6e2cae9fdd9b58983b2e7dead4b3a76">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a64ab19005827a80e63b387d5393ee9f9">operator*</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const double &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#af951deea3cb54784b21f1d26c876aa1e">operator/</a> (const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;, const double &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dgbmatrix.html#a7d1eaacff74c4ea410209cd9fb84d68c">operator*</a> (const double &amp;, const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;)</td></tr>
</table>
<hr/><a name="details" id="details"></a><h2>Detailed Description</h2>
<div class="textblock"><p>(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>
</div><hr/><h2>Constructor &amp; Destructor Documentation</h2>
<a class="anchor" id="a175eaab594a9f9f1d88330f3e511d8aa"></a><!-- doxytag: member="_dgbmatrix::_dgbmatrix" ref="a175eaab594a9f9f1d88330f3e511d8aa" args="()" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html#a175eaab594a9f9f1d88330f3e511d8aa">_dgbmatrix::_dgbmatrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td><code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p><a class="el" href="class__dgbmatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class">_dgbmatrix</a> constructor </p>

<p>Definition at line <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-constructor_8hpp_source.html">_dgbmatrix-constructor.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> =0;
  <a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> =0;
  <a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> =0;
  <a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a> =0;
  <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a> =NULL;
  <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a> =NULL;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a54b165895369496cb5c263b55ace48d5"></a><!-- doxytag: member="_dgbmatrix::_dgbmatrix" ref="a54b165895369496cb5c263b55ace48d5" args="(const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html#a175eaab594a9f9f1d88330f3e511d8aa">_dgbmatrix::_dgbmatrix</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em></td><td>)</td>
          <td><code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p><a class="el" href="class__dgbmatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class">_dgbmatrix</a> copy constructor </p>

<p>Definition at line <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00015">15</a> of file <a class="el" href="__dgbmatrix-constructor_8hpp_source.html">_dgbmatrix-constructor.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>, and <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> =mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>;
  <a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> =mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>;
  <a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> =mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>;
  <a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a> =mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>;
  <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a> =mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>;
  <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a> =mat.<a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>;
  
  mat.<a class="code" href="class__dgbmatrix.html#a736c8efcb6a6f9fd458090aa468ee2ce">nullify</a>();
}
</pre></div>
</div>
</div>
<a class="anchor" id="a3d238d72484bdc8ab3051d6e27710f8b"></a><!-- doxytag: member="_dgbmatrix::~_dgbmatrix" ref="a3d238d72484bdc8ab3051d6e27710f8b" args="()" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html#a3d238d72484bdc8ab3051d6e27710f8b">_dgbmatrix::~_dgbmatrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td><code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p><a class="el" href="class__dgbmatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision General Band Matrix Class">_dgbmatrix</a> destructor </p>

<p>Definition at line <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00033">33</a> of file <a class="el" href="__dgbmatrix-constructor_8hpp_source.html">_dgbmatrix-constructor.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keyword">delete</span>[] <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>;
  <span class="keyword">delete</span>[] <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>;
}
</pre></div>
</div>
</div>
<hr/><h2>Member Function Documentation</h2>
<a class="anchor" id="a300545260e3715527e79173b1291e883"></a><!-- doxytag: member="_dgbmatrix::to_zgbmatrix" ref="a300545260e3715527e79173b1291e883" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__zgbmatrix.html">_zgbmatrix</a> <a class="el" href="class__dgbmatrix.html#a300545260e3715527e79173b1291e883">_dgbmatrix::to_zgbmatrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>cast to <a class="el" href="class__zgbmatrix.html" title="(DO NOT USE) Smart-temporary Complex Double-precision General Band Matrix Class">_zgbmatrix</a> </p>

<p>Definition at line <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-cast_8hpp_source.html">_dgbmatrix-cast.hpp</a>.</p>

<p>References <a class="el" href="dcovector-misc_8hpp_source.html#l00134">_()</a>, <a class="el" href="zgbmatrix_8hpp_source.html#l00013">zgbmatrix::array</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">destroy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classzgbmatrix.html" title="Complex Double-precision General Band Matrix Class.">zgbmatrix</a> newmat(<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>);
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;(<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)*<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    newmat.array[<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>] =comple(<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>[<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>],0.0);
  }
  
  <a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a982423ea885f6583208d05625c31b03e"></a><!-- doxytag: member="_dgbmatrix::to_dgematrix" ref="a982423ea885f6583208d05625c31b03e" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> <a class="el" href="class__dgbmatrix.html#a982423ea885f6583208d05625c31b03e">_dgbmatrix::to_dgematrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>convert to <a class="el" href="class__dgematrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class">_dgematrix</a> </p>

<p>Definition at line <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">21</a> of file <a class="el" href="__dgbmatrix-cast_8hpp_source.html">_dgbmatrix-cast.hpp</a>.</p>

<p>References <a class="el" href="dcovector-misc_8hpp_source.html#l00134">_()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">destroy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a>(<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>).zero() );
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j) =(*this)(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  <a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a0bb6b797e67d9b5aef5c25b6f3de03e2"></a><!-- doxytag: member="_dgbmatrix::operator()" ref="a0bb6b797e67d9b5aef5c25b6f3de03e2" args="(const long &amp;, const long &amp;) const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double &amp; _dgbmatrix::operator() </td>
          <td>(</td>
          <td class="paramtype">const long &amp;&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const long &amp;&#160;</td>
          <td class="paramname"><em>j</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>operator() for const object </p>

<p>Definition at line <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-io_8hpp_source.html">_dgbmatrix-io.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>.</p>

<p>Referenced by <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>( <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;0 || j&lt;0 || <a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>&lt;=<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a> || n&lt;=j || i-j&gt;<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> || j-<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&gt;<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a> ){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;The required component is out of the matrix size.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a> &lt;&lt; <span class="stringliteral">&quot;,&quot;</span> &lt;&lt; j &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="comment">//return array[ku+i+(kl+ku)*j];</span>
  <span class="keywordflow">return</span> <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>-j+<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>];
}
</pre></div>
</div>
</div>
<a class="anchor" id="af59336a13c233cf8af0c54d170aa0f22"></a><!-- doxytag: member="_dgbmatrix::write" ref="af59336a13c233cf8af0c54d170aa0f22" args="(const char *) const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dgbmatrix.html#af59336a13c233cf8af0c54d170aa0f22">_dgbmatrix::write</a> </td>
          <td>(</td>
          <td class="paramtype">const char *&#160;</td>
          <td class="paramname"><em>filename</em></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>Definition at line <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">42</a> of file <a class="el" href="__dgbmatrix-io_8hpp_source.html">_dgbmatrix-io.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">destroy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  std::ofstream ofs(filename, std::ios::trunc);
  ofs.setf(std::cout.flags());
  ofs.precision(std::cout.precision());
  ofs.width(std::cout.width());
  ofs.fill(std::cout.fill());
  
  ofs &lt;&lt; <span class="stringliteral">&quot;#dgbmatrix&quot;</span> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; <a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; <a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; <a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; <a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a> &lt;&lt; std::endl;
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      ofs &lt;&lt; <a class="code" href="class__dgbmatrix.html#a0bb6b797e67d9b5aef5c25b6f3de03e2">operator()</a>(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j) &lt;&lt; <span class="stringliteral">&quot; &quot;</span>;
    }
    ofs &lt;&lt; std::endl;
  }
  
  ofs.close();
  <a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
}
</pre></div>
</div>
</div>
<a class="anchor" id="a736c8efcb6a6f9fd458090aa468ee2ce"></a><!-- doxytag: member="_dgbmatrix::nullify" ref="a736c8efcb6a6f9fd458090aa468ee2ce" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dgbmatrix.html#a736c8efcb6a6f9fd458090aa468ee2ce">_dgbmatrix::nullify</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>nullify all the matrix data </p>

<p>Definition at line <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-misc_8hpp_source.html">_dgbmatrix-misc.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">kl</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">ku</a>, <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">m</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">n</a>.</p>

<p>Referenced by <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00015">_dgbmatrix()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, and <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>=0;
  <a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>=0;
  <a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>=0;
  <a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>=0;
  <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>=NULL;
  <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>=NULL;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a170fdd995797f5d394863105e658df9e"></a><!-- doxytag: member="_dgbmatrix::destroy" ref="a170fdd995797f5d394863105e658df9e" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">_dgbmatrix::destroy</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>destroy all the matrix data </p>

<p>Definition at line <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">15</a> of file <a class="el" href="__dgbmatrix-misc_8hpp_source.html">_dgbmatrix-misc.hpp</a>.</p>

<p>References <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">array</a>, and <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">darray</a>.</p>

<p>Referenced by <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00099">dgbmatrix::operator*=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00048">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00016">dgbmatrix::operator+=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">operator-()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00058">dgbmatrix::operator-=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">operator&lt;&lt;()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">t()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keyword">delete</span> [] <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>;
  <span class="keyword">delete</span> [] <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>;
  <a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>=NULL;
  <a class="code" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d" title="array of pointers of column head addresses">darray</a>=NULL;
}
</pre></div>
</div>
</div>
<hr/><h2>Friends And Related Function Documentation</h2>
<a class="anchor" id="abc9cd92dfb30bf77e6e98d7c910228d8"></a><!-- doxytag: member="_dgbmatrix::operator&lt;&lt;" ref="abc9cd92dfb30bf77e6e98d7c910228d8" args="(std::ostream &amp;, const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">std::ostream&amp; operator&lt;&lt; </td>
          <td>(</td>
          <td class="paramtype">std::ostream &amp;&#160;</td>
          <td class="paramname"><em>s</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>Definition at line <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">23</a> of file <a class="el" href="__dgbmatrix-io_8hpp_source.html">_dgbmatrix-io.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>; j++){
      <span class="keywordflow">if</span>( <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-j&gt;mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> || j-<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&gt;mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a> ){ s &lt;&lt; <span class="stringliteral">&quot; x&quot;</span>; }
      <span class="keywordflow">else</span>{ s &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; mat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j); }
    }
    s &lt;&lt; std::endl;
  }
  
  mat.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> s;
}
</pre></div>
</div>
</div>
<a class="anchor" id="abc0db5fc9798b89672e2e103a07b097b"></a><!-- doxytag: member="_dgbmatrix::t" ref="abc0db5fc9798b89672e2e103a07b097b" args="(const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> <a class="el" href="class__dgbmatrix.html#abc0db5fc9798b89672e2e103a07b097b">t</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em></td><td>)</td>
          <td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>return its transposed dgbmatrix </p>

<p>Definition at line <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-calc_8hpp_source.html">_dgbmatrix-calc.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat(mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>, mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>);
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-newmat.kl); j&lt;std::min(newmat.n,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+newmat.ku+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j) =mat(j,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>);
    }
  }
  
  mat.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a024e4680944cd918565399ed987d4b20"></a><!-- doxytag: member="_dgbmatrix::i" ref="a024e4680944cd918565399ed987d4b20" args="(const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> <a class="el" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em></td><td>)</td>
          <td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>return its inverse matrix </p>

<p>Definition at line <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">19</a> of file <a class="el" href="__dgbmatrix-calc_8hpp_source.html">_dgbmatrix-calc.hpp</a>.</p>

<p>Referenced by <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;This matrix is not square and has no inverse matrix.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> mat_cp(mat);
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> mat_inv(mat_cp.m,mat_cp.n);
  mat_inv.identity();
  mat_cp.dgbsv(mat_inv);

  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(mat_inv);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a45ae3ef397cc029be013792ce75471c6"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a45ae3ef397cc029be013792ce75471c6" args="(const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a>&amp; operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em></td><td>)</td>
          <td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>+_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-unary_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-unary_8hpp_source.html">_dgbmatrix-unary.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a8feabd36261c12eea77cadbc33916efe"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a8feabd36261c12eea77cadbc33916efe" args="(const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em></td><td>)</td>
          <td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>-_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-unary_8hpp_source.html#l00010">10</a> of file <a class="el" href="__dgbmatrix-unary_8hpp_source.html">_dgbmatrix-unary.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;(mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)*mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>[<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>]=-mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>[<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>];
  }
  
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ae969cf19039e99763b4ef7e5bce4e7e6"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="ae969cf19039e99763b4ef7e5bce4e7e6" args="(const _dgbmatrix &amp;, const dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html">_dgbmatrix-dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matB);
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a46e30066362c28d9acea84681280212c"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a46e30066362c28d9acea84681280212c" args="(const _dgbmatrix &amp;, const _dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+_dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html">_dgbmatrix-_dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> matB;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ab15b39bacb422a2428295c7bc76994e1"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="ab15b39bacb422a2428295c7bc76994e1" args="(const _dgbmatrix &amp;, const dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html">_dgbmatrix-dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matB.<a class="code" href="classdsymatrix.html#a45ff312f6c694841edfd30809e09ba7b">to_dgematrix</a>());
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j) +=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="aef7ec311adfe90b7c26c0e92f8aebdd1"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="aef7ec311adfe90b7c26c0e92f8aebdd1" args="(const _dgbmatrix &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html">_dgbmatrix-_dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a3e98bdfeb5d939881272054ad728a175"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a3e98bdfeb5d939881272054ad728a175" args="(const _dgbmatrix &amp;, const dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html">_dgbmatrix-dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>&gt;=matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a> &amp;&amp; matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>&gt;=matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
        matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    <span class="keywordflow">return</span> matA;
  }
  
  <span class="keywordflow">else</span>{
    <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,std::max(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>),std::max(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>,matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>));
    newmat.zero();
    
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
  }
}
</pre></div>
</div>
</div>
<a class="anchor" id="aadaa4a751fb7153ca4a18679e1b235b7"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="aadaa4a751fb7153ca4a18679e1b235b7" args="(const _dgbmatrix &amp;, const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html">_dgbmatrix-_dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>&gt;matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> &amp;&amp; matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>&gt;matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matB.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> matA;
  }

  <span class="keywordflow">else</span> <span class="keywordflow">if</span>(matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>&gt;matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> &amp;&amp; matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>&gt;matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> matB;
  }
  
  <span class="keywordflow">else</span>{
    <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,std::max(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>),std::max(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>,matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>));
    newmat.zero();
    
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    matB.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
  }
}
</pre></div>
</div>
</div>
<a class="anchor" id="a4c110f12272b94d105a7b40b9c32626c"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a4c110f12272b94d105a7b40b9c32626c" args="(const _dgbmatrix &amp;, const dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html">_dgbmatrix-dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matB.<a class="code" href="classdgsmatrix.html#ad4251f53be60d1fea88f3cfaf6f3e409">to_dgematrix</a>() );
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a1899129fe441186353f81ad1dbcf6a35"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a1899129fe441186353f81ad1dbcf6a35" args="(const _dgbmatrix &amp;, const _dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix+_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html">_dgbmatrix-_dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matB.<a class="code" href="class__dgsmatrix.html#a5cc312d49cf097ebee56c1550988f463">to_dgematrix</a>() );
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dgsmatrix.html#aa2851712fa4f99d8d608815ef274a362">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a0eaee3557a07d4103ed7654583b9845e"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="a0eaee3557a07d4103ed7654583b9845e" args="(const _dgbmatrix &amp;, const dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="aff572701e2dee5c2206c81798d72a6dd"></a><!-- doxytag: member="_dgbmatrix::operator+" ref="aff572701e2dee5c2206c81798d72a6dd" args="(const _dgbmatrix &amp;, const _dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="aeb9361c4eec76c4964258446cd6bc24c"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="aeb9361c4eec76c4964258446cd6bc24c" args="(const _dgbmatrix &amp;, const dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html#l00027">27</a> of file <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html">_dgbmatrix-dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(-matB);
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a6212a92d6cd140914d20b82fd8620856"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a6212a92d6cd140914d20b82fd8620856" args="(const _dgbmatrix &amp;, const _dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-_dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html#l00026">26</a> of file <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html">_dgbmatrix-_dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j) =matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> matB;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a01bc3a5f90fcbf23476c5427f3ae508c"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a01bc3a5f90fcbf23476c5427f3ae508c" args="(const _dgbmatrix &amp;, const dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html#l00027">27</a> of file <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html">_dgbmatrix-dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( (-matB).<a class="code" href="class__dgbmatrix.html#a982423ea885f6583208d05625c31b03e">to_dgematrix</a>() );
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a4f919fea9e763b1afb1852ee711137c5"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a4f919fea9e763b1afb1852ee711137c5" args="(const _dgbmatrix &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00031">31</a> of file <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html">_dgbmatrix-_dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;newmat.n; j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)=-matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="aae7d46d5ff9708363386b1a041fce336"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="aae7d46d5ff9708363386b1a041fce336" args="(const _dgbmatrix &amp;, const dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html#l00044">44</a> of file <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html">_dgbmatrix-dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a subtraction.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) - (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>&gt;=matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a> &amp;&amp; matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>&gt;=matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
        matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    <span class="keywordflow">return</span> matA;
  }
  
  <span class="keywordflow">else</span>{
    <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,std::max(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>),std::max(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>,matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>));
    newmat.zero();
    
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
  }
}
</pre></div>
</div>
</div>
<a class="anchor" id="a209b6000adc349cff070ee25ba1d8941"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a209b6000adc349cff070ee25ba1d8941" args="(const _dgbmatrix &amp;, const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">57</a> of file <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html">_dgbmatrix-_dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a subtraction.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) - (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>&gt;matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a> &amp;&amp; matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>&gt;matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matB.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> matA;
  }
  
  <span class="keywordflow">else</span>{
    <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,std::max(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>),std::max(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>,matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>));
    newmat.zero();
    
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matB(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
      }
    }
    
    matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    matB.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
    <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
  }
}
</pre></div>
</div>
</div>
<a class="anchor" id="aaf20d6c83d831430cc819a3b670393aa"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="aaf20d6c83d831430cc819a3b670393aa" args="(const _dgbmatrix &amp;, const dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html#l00027">27</a> of file <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html">_dgbmatrix-dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( (-matB).<a class="code" href="class__dgbmatrix.html#a982423ea885f6583208d05625c31b03e">to_dgematrix</a>() );
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a6039af37d6b1e3ac7444f8b7475f5a02"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a6039af37d6b1e3ac7444f8b7475f5a02" args="(const _dgbmatrix &amp;, const _dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix-_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html#l00028">28</a> of file <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html">_dgbmatrix-_dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> || matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>!=matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a summation.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) + (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( (-matB).<a class="code" href="class__dgbmatrix.html#a982423ea885f6583208d05625c31b03e">to_dgematrix</a>() );
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); j&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)-=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dgsmatrix.html#aa2851712fa4f99d8d608815ef274a362">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="ad8a2c801718b3570fae41e3591886095"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="ad8a2c801718b3570fae41e3591886095" args="(const _dgbmatrix &amp;, const dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="a828552e4f2601fb88f6f62493a8b7e2d"></a><!-- doxytag: member="_dgbmatrix::operator&#45;" ref="a828552e4f2601fb88f6f62493a8b7e2d" args="(const _dgbmatrix &amp;, const _dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="a0f435627957fb3957ce083a40e4a2df0"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a0f435627957fb3957ce083a40e4a2df0" args="(const _dgbmatrix &amp;, const dcovector &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dcovector.html">_dcovector</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdcovector.html">dcovector</a> &amp;&#160;</td>
          <td class="paramname"><em>vec</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*dcovector operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dcovector_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-dcovector_8hpp_source.html">_dgbmatrix-dcovector.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=vec.<a class="code" href="classdcovector.html#a08e64e5ecb5121100787d0f5d48f0af4" title="vector size">l</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These matrix and vector can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; vec.<a class="code" href="classdcovector.html#a08e64e5ecb5121100787d0f5d48f0af4" title="vector size">l</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdcovector.html" title="Real Double-precision Column Vector Class.">dcovector</a> newvec(mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>);
  dgbmv_( <span class="charliteral">&#39;n&#39;</span>, mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>, mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>, 1.0, mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>,
          mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1, vec.<a class="code" href="classdcovector.html#a9b188bf160d8f7e24ea1d41deaceaf6d" title="1D array to store vector data">array</a>, 1, 0.0, newvec.array, 1 );
  
  mat.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newvec);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a929a817aaf554b47204b2fff2058af16"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a929a817aaf554b47204b2fff2058af16" args="(const _dgbmatrix &amp;, const _dcovector &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dcovector.html">_dcovector</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dcovector.html">_dcovector</a> &amp;&#160;</td>
          <td class="paramname"><em>vec</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*_dcovector operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html">_dgbmatrix-_dcovector.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=vec.<a class="code" href="class__dcovector.html#a2db3a577da8bc07533cf2da09fddaa58" title="vector size">l</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These matrix and vector can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; vec.<a class="code" href="class__dcovector.html#a2db3a577da8bc07533cf2da09fddaa58" title="vector size">l</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdcovector.html" title="Real Double-precision Column Vector Class.">dcovector</a> newvec(mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>);
  dgbmv_( <span class="charliteral">&#39;n&#39;</span>, mat.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>, mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>, 1.0, mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>,
          mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1, vec.<a class="code" href="class__dcovector.html#a27693f26cc6e9e25c349c7fbf7cc2b67" title="1D array to store vector data">array</a>, 1, 0.0, newvec.array, 1 );
  
  mat.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  vec.<a class="code" href="class__dcovector.html#a3e16f52298eac6c97958a9cbc6388bc4">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newvec);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a66c6987c5521e76467a4c149545aeb2f"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a66c6987c5521e76467a4c149545aeb2f" args="(const _dgbmatrix &amp;, const dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgematrix.html">dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html#l00051">51</a> of file <a class="el" href="__dgbmatrix-dgematrix_8hpp_source.html">_dgbmatrix-dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;newmat.n; j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); k&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a7d45270990d9f629335bfe565ae81876"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a7d45270990d9f629335bfe565ae81876" args="(const _dgbmatrix &amp;, const _dgematrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgematrix.html">_dgematrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*_dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html#l00049">49</a> of file <a class="el" href="__dgbmatrix-__dgematrix_8hpp_source.html">_dgbmatrix-_dgematrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;newmat.n; j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); k&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dgematrix.html#ac6d270d6c2d1d73afdfb2029eca920d9">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="ad08c9cddd6aa158109e227ba907ba9ce"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="ad08c9cddd6aa158109e227ba907ba9ce" args="(const _dgbmatrix &amp;, const dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdsymatrix.html">dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*dgematrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html#l00051">51</a> of file <a class="el" href="__dgbmatrix-dsymatrix_8hpp_source.html">_dgbmatrix-dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;newmat.n; j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); k&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="aadf4f8e6254cd0f581d2fd1ad034191f"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="aadf4f8e6254cd0f581d2fd1ad034191f" args="(const _dgbmatrix &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00059">59</a> of file <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html">_dgbmatrix-_dsymatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;newmat.n; j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); k&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+=matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a2ea7426af09460cc63d8127658c47909"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a2ea7426af09460cc63d8127658c47909" args="(const _dgbmatrix &amp;, const dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgbmatrix.html">dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html#l00084">84</a> of file <a class="el" href="__dgbmatrix-dgbmatrix_8hpp_source.html">_dgbmatrix-dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>, std::min(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>,matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>-1), std::min(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>,matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a>-1) );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-newmat.kl); j&lt;std::min(newmat.n,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+newmat.ku+1); j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max( std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>), std::max(<span class="keywordtype">long</span>(0),j-matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>) );
          k&lt; std::min( std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1), std::min(matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>,j+matB.<a class="code" href="classdgbmatrix.html#aeeb3e79431677905a4ece19c6616edd4" title="lower band width">kl</a>+1) ); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+= matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a6338b568730f34d9ebe9a69871b1b797"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a6338b568730f34d9ebe9a69871b1b797" args="(const _dgbmatrix &amp;, const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00100">100</a> of file <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html">_dgbmatrix-_dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgbmatrix.html" title="Real Double-precision General Band Matrix Class.">dgbmatrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, std::min(matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>,matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>-1), std::min(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>,matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>-1) );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=0; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;newmat.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>; <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-newmat.kl); j&lt;std::min(newmat.n,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+newmat.ku+1); j++){
      <span class="keywordflow">for</span>(<span class="keywordtype">long</span> k=std::max( std::max(<span class="keywordtype">long</span>(0),<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>-matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>), std::max(<span class="keywordtype">long</span>(0),j-matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>) );
          k&lt; std::min( std::min(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>,<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>+matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1), std::min(matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,j+matB.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+1) ); k++){
        newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,j)+= matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="ab255cc79e2cd0ee188e1598f6decf672"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="ab255cc79e2cd0ee188e1598f6decf672" args="(const _dgbmatrix &amp;, const dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdgsmatrix.html">dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html#l00051">51</a> of file <a class="el" href="__dgbmatrix-dgsmatrix_8hpp_source.html">_dgbmatrix-dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</a> );
  newmat.zero();
  <span class="keywordflow">for</span>(std::vector&lt;dcomponent&gt;::const_iterator it=matB.<a class="code" href="classdgsmatrix.html#a9c4f940eeedf7a7869cf03ae4a3558c6" title="matrix data">data</a>.begin(); it!=matB.<a class="code" href="classdgsmatrix.html#a9c4f940eeedf7a7869cf03ae4a3558c6" title="matrix data">data</a>.end(); it++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=std::max(<span class="keywordtype">long</span>(0),<span class="keywordtype">long</span>(it-&gt;i)-(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)); <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,<span class="keywordtype">long</span>(it-&gt;i)+matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,it-&gt;j) += matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,it-&gt;<a class="code" href="class__dgematrix.html#ab4e0cbe78518b4af9da89406bb2b9d19">i</a>)*it-&gt;v;
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="acee43528df478ebfaa53a2fc314e19a0"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="acee43528df478ebfaa53a2fc314e19a0" args="(const _dgbmatrix &amp;, const _dgsmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matA</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgsmatrix.html">_dgsmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>matB</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html#l00053">53</a> of file <a class="el" href="__dgbmatrix-__dgsmatrix_8hpp_source.html">_dgbmatrix-_dgsmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  <span class="keywordflow">if</span>(matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a>){
    ERROR_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;These two matrises can not make a product.&quot;</span> &lt;&lt; std::endl
              &lt;&lt; <span class="stringliteral">&quot;Your input was (&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;) * (&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;).&quot;</span> &lt;&lt; std::endl;
    exit(1);
  }
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>, matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</a> );
  newmat.zero();
  <span class="keywordflow">for</span>(std::vector&lt;dcomponent&gt;::const_iterator it=matB.<a class="code" href="class__dgsmatrix.html#a65ae5bc0001de27f9fe1fb13870b4a30" title="matrix data">data</a>.begin(); it!=matB.<a class="code" href="class__dgsmatrix.html#a65ae5bc0001de27f9fe1fb13870b4a30" title="matrix data">data</a>.end(); it++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>=std::max(<span class="keywordtype">long</span>(0),<span class="keywordtype">long</span>(it-&gt;i)-(matA.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)); <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>&lt;std::min(matA.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>,<span class="keywordtype">long</span>(it-&gt;i)+matA.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>); <a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>++){
      newmat(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,it-&gt;j) += matA(<a class="code" href="class__dgbmatrix.html#a024e4680944cd918565399ed987d4b20">i</a>,it-&gt;<a class="code" href="class__dgematrix.html#ab4e0cbe78518b4af9da89406bb2b9d19">i</a>)*it-&gt;v;
    }
  }
  
  matA.<a class="code" href="class__dgbmatrix.html#a170fdd995797f5d394863105e658df9e">destroy</a>();
  matB.<a class="code" href="class__dgsmatrix.html#aa2851712fa4f99d8d608815ef274a362">destroy</a>();
  <span class="keywordflow">return</span> <a class="code" href="dcovector-misc_8hpp.html#a297beda6502b0d92c99576d2fb9548cb">_</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="ae416997168ff0b923f4fd006fe034908"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="ae416997168ff0b923f4fd006fe034908" args="(const _dgbmatrix &amp;, const dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="classdssmatrix.html">dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="af6e2cae9fdd9b58983b2e7dead4b3a76"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="af6e2cae9fdd9b58983b2e7dead4b3a76" args="(const _dgbmatrix &amp;, const _dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgematrix.html">_dgematrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname">, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dssmatrix.html">_dssmatrix</a> &amp;&#160;</td>
          <td class="paramname">&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="a64ab19005827a80e63b387d5393ee9f9"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a64ab19005827a80e63b387d5393ee9f9" args="(const _dgbmatrix &amp;, const double &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const double &amp;&#160;</td>
          <td class="paramname"><em>d</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix*double operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dgbmatrix-double_8hpp_source.html">_dgbmatrix-double.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  dscal_((mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)*mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, d, mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>, 1);
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="af951deea3cb54784b21f1d26c876aa1e"></a><!-- doxytag: member="_dgbmatrix::operator/" ref="af951deea3cb54784b21f1d26c876aa1e" args="(const _dgbmatrix &amp;, const double &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator/ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const double &amp;&#160;</td>
          <td class="paramname"><em>d</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>_dgbmatrix/double operator </p>

<p>Definition at line <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00011">11</a> of file <a class="el" href="__dgbmatrix-double_8hpp_source.html">_dgbmatrix-double.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  dscal_((mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)*mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, 1./d, mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>, 1);
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a7d1eaacff74c4ea410209cd9fb84d68c"></a><!-- doxytag: member="_dgbmatrix::operator*" ref="a7d1eaacff74c4ea410209cd9fb84d68c" args="(const double &amp;, const _dgbmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const double &amp;&#160;</td>
          <td class="paramname"><em>d</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">const <a class="el" href="class__dgbmatrix.html">_dgbmatrix</a> &amp;&#160;</td>
          <td class="paramname"><em>mat</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>double*_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="double-__dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="double-__dgbmatrix_8hpp_source.html">double-_dgbmatrix.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  dscal_((mat.<a class="code" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c" title="lower band width">kl</a>+mat.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>+1)*mat.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a>, d, mat.<a class="code" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb" title="1D array to store matrix data">array</a>, 1);
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<hr/><h2>Member Data Documentation</h2>
<a class="anchor" id="a587a07ecbbfc6bbdb0e05574ea6b088a"></a><!-- doxytag: member="_dgbmatrix::m" ref="a587a07ecbbfc6bbdb0e05574ea6b088a" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long <a class="el" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a">_dgbmatrix::m</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>matrix row size </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00009">9</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00099">dgbmatrix::operator*=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00048">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00016">dgbmatrix::operator+=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">operator-()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00058">dgbmatrix::operator-=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">operator&lt;&lt;()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">t()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>

</div>
</div>
<a class="anchor" id="a95aebc7e2d4ea66b1b7a08e981df6971"></a><!-- doxytag: member="_dgbmatrix::n" ref="a95aebc7e2d4ea66b1b7a08e981df6971" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long <a class="el" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971">_dgbmatrix::n</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>matrix column size </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00010">10</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00019">i()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00099">dgbmatrix::operator*=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00048">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00016">dgbmatrix::operator+=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">operator-()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00058">dgbmatrix::operator-=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00011">operator/()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">operator&lt;&lt;()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">t()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>

</div>
</div>
<a class="anchor" id="aeea4fdf76630ca1b51f93c5bed64662c"></a><!-- doxytag: member="_dgbmatrix::kl" ref="aeea4fdf76630ca1b51f93c5bed64662c" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long <a class="el" href="class__dgbmatrix.html#aeea4fdf76630ca1b51f93c5bed64662c">_dgbmatrix::kl</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>lower band width </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00011">11</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00099">dgbmatrix::operator*=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00048">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00016">dgbmatrix::operator+=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">operator-()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00058">dgbmatrix::operator-=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00011">operator/()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">operator&lt;&lt;()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">t()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>

</div>
</div>
<a class="anchor" id="a116dadc8ed5246ee0b4cd0a9fffc53ea"></a><!-- doxytag: member="_dgbmatrix::ku" ref="a116dadc8ed5246ee0b4cd0a9fffc53ea" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long <a class="el" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea">_dgbmatrix::ku</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>upper band width </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00012">12</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00099">dgbmatrix::operator*=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00048">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00016">dgbmatrix::operator+=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dgbmatrix_8hpp_source.html#l00057">operator-()</a>, <a class="el" href="dgbmatrix-__dgbmatrix_8hpp_source.html#l00058">dgbmatrix::operator-=()</a>, <a class="el" href="dgematrix-__dgbmatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00011">operator/()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00023">operator&lt;&lt;()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, <a class="el" href="__dgbmatrix-calc_8hpp_source.html#l00003">t()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00021">to_dgematrix()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00042">write()</a>.</p>

</div>
</div>
<a class="anchor" id="a3d7680e0014aa00249f76a86e13de9eb"></a><!-- doxytag: member="_dgbmatrix::array" ref="a3d7680e0014aa00249f76a86e13de9eb" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double* <a class="el" href="class__dgbmatrix.html#a3d7680e0014aa00249f76a86e13de9eb">_dgbmatrix::array</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>1D array to store matrix data </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00013">13</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">destroy()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-__dcovector_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="__dgbmatrix-unary_8hpp_source.html#l00010">operator-()</a>, <a class="el" href="__dgbmatrix-double_8hpp_source.html#l00011">operator/()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, <a class="el" href="__dgbmatrix-cast_8hpp_source.html#l00003">to_zgbmatrix()</a>, and <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00033">~_dgbmatrix()</a>.</p>

</div>
</div>
<a class="anchor" id="ae9c142fd6903cf00df8e1d5363e9848d"></a><!-- doxytag: member="_dgbmatrix::darray" ref="ae9c142fd6903cf00df8e1d5363e9848d" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double** <a class="el" href="class__dgbmatrix.html#ae9c142fd6903cf00df8e1d5363e9848d">_dgbmatrix::darray</a><code> [mutable]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

<p>array of pointers of column head addresses </p>

<p>Definition at line <a class="el" href="__dgbmatrix_8hpp_source.html#l00014">14</a> of file <a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00166">_()</a>, <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00003">_dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00015">destroy()</a>, <a class="el" href="dgbmatrix-constructor_8hpp_source.html#l00033">dgbmatrix::dgbmatrix()</a>, <a class="el" href="__dgbmatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgbmatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="dgbmatrix-misc_8hpp_source.html#l00068">dgbmatrix::shallow_copy()</a>, and <a class="el" href="__dgbmatrix-constructor_8hpp_source.html#l00033">~_dgbmatrix()</a>.</p>

</div>
</div>
<hr/>The documentation for this class was generated from the following files:<ul>
<li><a class="el" href="__dgbmatrix_8hpp_source.html">_dgbmatrix.hpp</a></li>
<li><a class="el" href="__dgbmatrix-cast_8hpp_source.html">_dgbmatrix-cast.hpp</a></li>
<li><a class="el" href="__dgbmatrix-constructor_8hpp_source.html">_dgbmatrix-constructor.hpp</a></li>
<li><a class="el" href="__dgbmatrix-io_8hpp_source.html">_dgbmatrix-io.hpp</a></li>
<li><a class="el" href="__dgbmatrix-misc_8hpp_source.html">_dgbmatrix-misc.hpp</a></li>
</ul>
</div><!-- contents -->
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
     onmouseover="return searchBox.OnSearchSelectShow()"
     onmouseout="return searchBox.OnSearchSelectHide()"
     onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark">&#160;</span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark">&#160;</span>Classes</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark">&#160;</span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark">&#160;</span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark">&#160;</span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark">&#160;</span>Friends</a></div>

<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0" 
        name="MSearchResults" id="MSearchResults">
</iframe>
</div>



<hr class="footer"/><address class="footer"><small>
Generated by &#160;<a href="http://www.doxygen.org/index.html">
<img class="footer" src="doxygen.png" alt="doxygen"/>
</a> 1.7.6.1
</small></address>

</body>
</html>
