<!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: _dsymatrix 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">_dsymatrix Class Reference</div>  </div>
</div><!--header-->
<div class="contents">
<!-- doxytag: class="_dsymatrix" -->
<p>(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class  
 <a href="class__dsymatrix.html#details">More...</a></p>

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

<p><a href="class__dsymatrix-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__dsymatrix.html#a0923826d0991b6f0b239d918db29a6fc">_dsymatrix</a> ()</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#ad66cf7c617f4da42d8b1ee19ad103e17">_dsymatrix</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a9043f20cff6245ceeb267825c74e7f88">~_dsymatrix</a> ()</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__zhematrix.html">_zhematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a1bdcb3af835709870c90a234d617edc9">to_zhematrix</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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">to_dgematrix</a> () const </td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dssmatrix.html">_dssmatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#aef33c5a0a42724acd76c81dfd2318cc8">to_dssmatrix</a> (const double=DBL_MIN) 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__dsymatrix.html#a9ea016b607f516fa57b52e946ebdc2de">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__dsymatrix.html#a0c61eed301e45a7f3d3a678f3e9e800f">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__dsymatrix.html#a15a1bf572055917c5963ba65d83fa282">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__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</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__dsymatrix.html#abfd55ec47814e32e991ddd3c8de6b1f8">complete</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 const &amp;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f">m</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">matrix row size  <a href="#a0ae0bb2812e15cb1a0667fd6f6f6b44f"></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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b">n</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">matrix column size  <a href="#a8b56cdf252a15da3e211b558abd46c1b"></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__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a">array</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">1D array to store matrix data  <a href="#a656ab414ee51638aedc9bf4960fe428a"></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__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc">darray</a></td></tr>
<tr><td class="mdescLeft">&#160;</td><td class="mdescRight">array of pointers of column head addresses  <a href="#a9db2d341c99952e75adb7e183450c3dc"></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__dsymatrix.html#a077ea0a074855b18f11abcd19a77ed96">operator&lt;&lt;</a> (std::ostream &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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#ad6d0f9c8ec17d8237a26e98d0f10306c">t</a> (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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">void&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a0cff2ad3ba1ef5e9d56b0b42437225de">idamax</a> (long &amp;, long &amp;, const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">double&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#abec838aac77ffaf2de6539bc21425e28">damax</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top">const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a1791c142cdb35a6981bdc4a12f6f3455">operator+</a> (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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a916a82da86a46ee0d61fc281ece79714">operator-</a> (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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#ae2dfab27868426c693ebecd09b8c7489">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a4ece796a5cfa4b73fdbd0ca8c277f04f">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#acb40c53507cbcb1e7cd541fc0e55e442">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#ab49f1ed04f239df3b9b1d30f9e9ccc21">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#adae8f965901ba1791a5fc1c629668868">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a3d0715719a516fb07f42490224b49fb9">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a9aed430aafed12429123b8ffaa3547c7">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#ae66d01b7b77643ba8c1ef8c180c6eb1a">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a95c2875a85d6b74ea1fd25574bd94ad8">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a94eaaffd21f4cce63df3feb343947cc4">operator+</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#ae073dbc5e88525decfe52762e31e8b68">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a873eaf16f7c95d5097e38030e1a0170d">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#acb8b2c3a3aa1fa010ded175f5d7ef0bb">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a146298da44c82fd8f2753f3a965bcb3e">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a92948829a98575861e6556f62b7e9c62">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a126c457d3d7e16dd0c6af841df491a33">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a5e102121931a8b18b1e3b082b40ad646">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#af7068f5e29904707b5260316afa6152b">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#ad81c3bfb85f1ebd6f38bed5e746a3c19">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a4bb60887d2c29b980bf361dc35940800">operator-</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#af51cdd88dbb5453f102e096d05cc73a5">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a3f3f8a2d27f873d0842cc3599c0cd860">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#afda93bc54596d73f3921ca6572790814">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a977ef4982892149eb6bb4794e7f72a0f">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a8d2bf7e9fdf57f365e74ed38313ba8d0">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#aa28788252525cda33141d129fd6e6a09">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#afd3c21b75b16dfaed32f26f7332bd591">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dgematrix.html">_dgematrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a1e50366c55b277b6677bdcb87397ee5b">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a577fce177108b8fd791ac8a73813f132">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#ad9c9cafc9a48b0f44b9b39ccf5d149d8">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a3ea6f20cbc65b317fb261187bf3d8732">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html#a4b58e624b1d0fd3ceb8976cbbad920a1">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#aa1c618a4d03afb961c861c64aef29aaf">operator*</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;, const double &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a4acbdbaa8ba95f7ef2b461ac116ad4ae">operator/</a> (const <a class="el" href="class__dsymatrix.html">_dsymatrix</a> &amp;, const double &amp;)</td></tr>
<tr><td class="memItemLeft" align="right" valign="top"><a class="el" href="class__dsymatrix.html">_dsymatrix</a>&#160;</td><td class="memItemRight" valign="bottom"><a class="el" href="class__dsymatrix.html#a53c70446bb321088fa2e2a79eba8b58c">operator*</a> (const double &amp;, const <a class="el" href="class__dsymatrix.html">_dsymatrix</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 Symmetric Matrix Class </p>

<p>Definition at line <a class="el" href="__dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a>.</p>
</div><hr/><h2>Constructor &amp; Destructor Documentation</h2>
<a class="anchor" id="a0923826d0991b6f0b239d918db29a6fc"></a><!-- doxytag: member="_dsymatrix::_dsymatrix" ref="a0923826d0991b6f0b239d918db29a6fc" args="()" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html#a0923826d0991b6f0b239d918db29a6fc">_dsymatrix::_dsymatrix</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__dsymatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class">_dsymatrix</a> constructor without arguments </p>

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

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00011">array</a>, <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>
<div class="fragment"><pre class="fragment">  :<a class="code" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f" title="matrix row size">m</a>(<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>)
{VERBOSE_REPORT;
  <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> =0;
  <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a> =NULL;
  <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a> =NULL;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ad66cf7c617f4da42d8b1ee19ad103e17"></a><!-- doxytag: member="_dsymatrix::_dsymatrix" ref="ad66cf7c617f4da42d8b1ee19ad103e17" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html#a0923826d0991b6f0b239d918db29a6fc">_dsymatrix::_dsymatrix</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class">_dsymatrix</a> copy constructor </p>

<p>Definition at line <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00013">13</a> of file <a class="el" href="__dsymatrix-constructor_8hpp_source.html">_dsymatrix-constructor.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00011">array</a>, <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>, <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>, and <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00003">nullify()</a>.</p>
<div class="fragment"><pre class="fragment">  :<a class="code" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f" title="matrix row size">m</a>(<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>)
{VERBOSE_REPORT;
  <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> =mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>;
  <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a> =mat.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>;
  <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a> =mat.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>;
  
  mat.<a class="code" href="class__dsymatrix.html#a15a1bf572055917c5963ba65d83fa282">nullify</a>();
}
</pre></div>
</div>
</div>
<a class="anchor" id="a9043f20cff6245ceeb267825c74e7f88"></a><!-- doxytag: member="_dsymatrix::~_dsymatrix" ref="a9043f20cff6245ceeb267825c74e7f88" args="()" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html#a9043f20cff6245ceeb267825c74e7f88">_dsymatrix::~_dsymatrix</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td><code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>dsymatrix destructor </p>

<p>Definition at line <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00029">29</a> of file <a class="el" href="__dsymatrix-constructor_8hpp_source.html">_dsymatrix-constructor.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00011">array</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keyword">delete</span>[] <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>;
  <span class="keyword">delete</span>[] <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>;
}
</pre></div>
</div>
</div>
<hr/><h2>Member Function Documentation</h2>
<a class="anchor" id="a1bdcb3af835709870c90a234d617edc9"></a><!-- doxytag: member="_dsymatrix::to_zhematrix" ref="a1bdcb3af835709870c90a234d617edc9" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__zhematrix.html">_zhematrix</a> <a class="el" href="class__dsymatrix.html#a1bdcb3af835709870c90a234d617edc9">_dsymatrix::to_zhematrix</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__zhematrix.html" title="(DO NOT USE) Smart-temporary Complex Double-precision Hermitian Matrix Class">_zhematrix</a> </p>

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

<p>References <a class="el" href="dcovector-misc_8hpp_source.html#l00134">_()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00012">destroy()</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classzhematrix.html" title="Complex Double-precision Hermitian Matrix Class [l-type (UPLO=l) Strage].">zhematrix</a> newmat(<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) =comple((*<span class="keyword">this</span>)(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j),0.0);
    }
  }
  
  <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="acdf4ff1ebc021a28f63351bd6f6a368c"></a><!-- doxytag: member="_dsymatrix::to_dgematrix" ref="acdf4ff1ebc021a28f63351bd6f6a368c" 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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">_dsymatrix::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="__dsymatrix-cast_8hpp_source.html#l00022">22</a> of file <a class="el" href="__dsymatrix-cast_8hpp_source.html">_dsymatrix-cast.hpp</a>.</p>

<p>References <a class="el" href="dcovector-misc_8hpp_source.html#l00134">_()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00012">destroy()</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>

<p>Referenced by <a class="el" href="__dgsmatrix-__dsymatrix_8hpp_source.html#l00003">operator+()</a>, and <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html#l00026">operator-()</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="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>,<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) =(*this)(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  <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="aef33c5a0a42724acd76c81dfd2318cc8"></a><!-- doxytag: member="_dsymatrix::to_dssmatrix" ref="aef33c5a0a42724acd76c81dfd2318cc8" args="(const double=DBL_MIN) const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dssmatrix.html">_dssmatrix</a> <a class="el" href="class__dsymatrix.html#aef33c5a0a42724acd76c81dfd2318cc8">_dsymatrix::to_dssmatrix</a> </td>
          <td>(</td>
          <td class="paramtype">const double&#160;</td>
          <td class="paramname"> = <code>DBL_MIN</code></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">

</div>
</div>
<a class="anchor" id="a9ea016b607f516fa57b52e946ebdc2de"></a><!-- doxytag: member="_dsymatrix::operator()" ref="a9ea016b607f516fa57b52e946ebdc2de" args="(const long &amp;, const long &amp;) const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double &amp; _dsymatrix::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 object </p>

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

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>

<p>Referenced by <a class="el" href="__dsymatrix-io_8hpp_source.html#l00050">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;0 || j&lt;0 || <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a> || <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>&lt;=j ){
    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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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="keywordflow">if</span>( <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a> &gt;= j ){
    <span class="comment">//return array[i+n*j];</span>
    <span class="keywordflow">return</span> <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
  } <span class="keywordflow">else</span> {
    <span class="comment">//return array[j+n*i];</span>
    <span class="keywordflow">return</span> <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>][j];
  }
}
</pre></div>
</div>
</div>
<a class="anchor" id="a0c61eed301e45a7f3d3a678f3e9e800f"></a><!-- doxytag: member="_dsymatrix::write" ref="a0c61eed301e45a7f3d3a678f3e9e800f" args="(const char *) const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dsymatrix.html#a0c61eed301e45a7f3d3a678f3e9e800f">_dsymatrix::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="__dsymatrix-io_8hpp_source.html#l00050">50</a> of file <a class="el" href="__dsymatrix-io_8hpp_source.html">_dsymatrix-io.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00012">destroy()</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i</a>, <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>, and <a class="el" href="__dsymatrix-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;#dsymatrix&quot;</span> &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; std::endl;
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++ ){
      ofs &lt;&lt; <a class="code" href="class__dsymatrix.html#a9ea016b607f516fa57b52e946ebdc2de">operator()</a>(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) &lt;&lt; <span class="stringliteral">&quot; &quot;</span>;
    }
    ofs &lt;&lt; std::endl;
  }
  
  ofs.close();
  <a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
}
</pre></div>
</div>
</div>
<a class="anchor" id="a15a1bf572055917c5963ba65d83fa282"></a><!-- doxytag: member="_dsymatrix::nullify" ref="a15a1bf572055917c5963ba65d83fa282" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dsymatrix.html#a15a1bf572055917c5963ba65d83fa282">_dsymatrix::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="__dsymatrix-misc_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-misc_8hpp_source.html">_dsymatrix-misc.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00011">array</a>, <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>

<p>Referenced by <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00013">_dsymatrix()</a>, <a class="el" href="dsymatrix-constructor_8hpp_source.html#l00029">dsymatrix::dsymatrix()</a>, and <a class="el" href="dsymatrix-misc_8hpp_source.html#l00066">dsymatrix::shallow_copy()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>=0;
  <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>=NULL;
  <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>=NULL;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a09e5b4ca8579bef58e01e80c22e799b1"></a><!-- doxytag: member="_dsymatrix::destroy" ref="a09e5b4ca8579bef58e01e80c22e799b1" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">_dsymatrix::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="__dsymatrix-misc_8hpp_source.html#l00012">12</a> of file <a class="el" href="__dsymatrix-misc_8hpp_source.html">_dsymatrix-misc.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00011">array</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>.</p>

<p>Referenced by <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00059">operator*()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00049">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00015">dsymatrix::operator+=()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00031">operator-()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00038">dsymatrix::operator-=()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dsymatrix-io_8hpp_source.html#l00028">operator&lt;&lt;()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00022">to_dgematrix()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00003">to_zhematrix()</a>, and <a class="el" href="__dsymatrix-io_8hpp_source.html#l00050">write()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <span class="keyword">delete</span> [] <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>;
  <span class="keyword">delete</span> [] <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>;
  <a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>=NULL;
  <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>=NULL;
}
</pre></div>
</div>
</div>
<a class="anchor" id="abfd55ec47814e32e991ddd3c8de6b1f8"></a><!-- doxytag: member="_dsymatrix::complete" ref="abfd55ec47814e32e991ddd3c8de6b1f8" args="() const " -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dsymatrix.html#abfd55ec47814e32e991ddd3c8de6b1f8">_dsymatrix::complete</a> </td>
          <td>(</td>
          <td class="paramname"></td><td>)</td>
          <td> const<code> [inline]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>complete the upper-right components </p>

<p>Definition at line <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00022">22</a> of file <a class="el" href="__dsymatrix-misc_8hpp_source.html">_dsymatrix-misc.hpp</a>.</p>

<p>References <a class="el" href="__dsymatrix_8hpp_source.html#l00012">darray</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i</a>, and <a class="el" href="__dsymatrix_8hpp_source.html#l00010">n</a>.</p>

<p>Referenced by <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00049">operator*()</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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      <a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>][j] =<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
    }
  }
}
</pre></div>
</div>
</div>
<hr/><h2>Friends And Related Function Documentation</h2>
<a class="anchor" id="a077ea0a074855b18f11abcd19a77ed96"></a><!-- doxytag: member="_dsymatrix::operator&lt;&lt;" ref="a077ea0a074855b18f11abcd19a77ed96" args="(std::ostream &amp;, const dsymatrix &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="classdsymatrix.html">dsymatrix</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="dsymatrix-io_8hpp_source.html#l00077">77</a> of file <a class="el" href="dsymatrix-io_8hpp_source.html">dsymatrix-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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;mat.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      s &lt;&lt; <span class="stringliteral">&quot; &quot;</span> &lt;&lt; mat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) &lt;&lt; <span class="stringliteral">&quot; &quot;</span>;
    }
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>+1; j&lt;mat.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>; j++){
      s &lt;&lt; <span class="stringliteral">&quot;{&quot;</span> &lt;&lt; mat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) &lt;&lt; <span class="stringliteral">&quot;}&quot;</span>;
    }
    s &lt;&lt; std::endl;
  }
  <span class="keywordflow">return</span> s;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ad6d0f9c8ec17d8237a26e98d0f10306c"></a><!-- doxytag: member="_dsymatrix::t" ref="ad6d0f9c8ec17d8237a26e98d0f10306c" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> <a class="el" href="class__dsymatrix.html#ad6d0f9c8ec17d8237a26e98d0f10306c">t</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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 transposed <a class="el" href="class__dsymatrix.html" title="(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class">_dsymatrix</a> </p>

<p>Definition at line <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-calc_8hpp_source.html">_dsymatrix-calc.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
<span class="preprocessor">#ifdef  CPPL_DEBUG</span>
<span class="preprocessor"></span>  WARNING_REPORT;
  std::cerr &lt;&lt; <span class="stringliteral">&quot;This function call has no effect since the matrix is symmetric.&quot;</span> &lt;&lt; std::endl;
<span class="preprocessor">#endif//CPPL_DEBUG</span>
<span class="preprocessor"></span>  
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a0697d2173dfc389bf1312d7bf6a5bd3b"></a><!-- doxytag: member="_dsymatrix::i" ref="a0697d2173dfc389bf1312d7bf6a5bd3b" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> <a class="el" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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="__dsymatrix-calc_8hpp_source.html#l00015">15</a> of file <a class="el" href="__dsymatrix-calc_8hpp_source.html">_dsymatrix-calc.hpp</a>.</p>

<p>Referenced by <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00022">complete()</a>, <a class="el" href="__dsymatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00022">to_dgematrix()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00003">to_zhematrix()</a>, and <a class="el" href="__dsymatrix-io_8hpp_source.html#l00050">write()</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classdsymatrix.html" title="Real Double-precision Symmetric Matrix Class [l-type (UPLO=l) Strage].">dsymatrix</a> mat_cp(mat);
  <a class="code" href="classdsymatrix.html" title="Real Double-precision Symmetric Matrix Class [l-type (UPLO=l) Strage].">dsymatrix</a> mat_inv(mat_cp.n);
  mat_inv.identity();
  
  <span class="keywordtype">char</span> UPLO(<span class="charliteral">&#39;l&#39;</span>);
  <span class="keywordtype">long</span> NRHS(mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>), LDA(mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>), *IPIV(<span class="keyword">new</span> <span class="keywordtype">long</span>[mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>]), LDB(mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>), LWORK(-1), INFO(1);
  <span class="keywordtype">double</span> *WORK( <span class="keyword">new</span> <span class="keywordtype">double</span>[1] );
  dsysv_(UPLO, mat_cp.n, NRHS, mat_cp.array, LDA, IPIV, mat_inv.array, LDB, WORK, LWORK, INFO);
  LWORK = long(WORK[0]);
  <span class="keyword">delete</span> [] WORK;  WORK = <span class="keyword">new</span> <span class="keywordtype">double</span>[LWORK];
  dsysv_(UPLO, mat_cp.n, NRHS, mat_cp.array, LDA, IPIV, mat_inv.array, LDB, WORK, LWORK, INFO);
  <span class="keyword">delete</span> [] WORK; <span class="keyword">delete</span> [] IPIV;
  <span class="keywordflow">if</span>(INFO!=0){
    WARNING_REPORT;
    std::cerr &lt;&lt; <span class="stringliteral">&quot;Serious trouble happend. INFO = &quot;</span> &lt;&lt; INFO &lt;&lt; <span class="stringliteral">&quot;.&quot;</span> &lt;&lt; std::endl;
  }
  
  <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="a0cff2ad3ba1ef5e9d56b0b42437225de"></a><!-- doxytag: member="_dsymatrix::idamax" ref="a0cff2ad3ba1ef5e9d56b0b42437225de" args="(long &amp;, long &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">void <a class="el" href="class__dsymatrix.html#a0cff2ad3ba1ef5e9d56b0b42437225de">idamax</a> </td>
          <td>(</td>
          <td class="paramtype">long &amp;&#160;</td>
          <td class="paramname"><em>i</em>, </td>
        </tr>
        <tr>
          <td class="paramkey"></td>
          <td></td>
          <td class="paramtype">long &amp;&#160;</td>
          <td class="paramname"><em>j</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>mat</em>&#160;</td>
        </tr>
        <tr>
          <td></td>
          <td>)</td>
          <td></td><td><code> [friend]</code></td>
        </tr>
      </table>
</div>
<div class="memdoc">
<p>search the index of element having the largest absolute value in 0-based numbering system </p>

<p>Definition at line <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00044">44</a> of file <a class="el" href="__dsymatrix-calc_8hpp_source.html">_dsymatrix-calc.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classdsymatrix.html" title="Real Double-precision Symmetric Matrix Class [l-type (UPLO=l) Strage].">dsymatrix</a> newmat(mat);
  <a class="code" href="class__dsymatrix.html#a0cff2ad3ba1ef5e9d56b0b42437225de">idamax</a>(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>, j, newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="abec838aac77ffaf2de6539bc21425e28"></a><!-- doxytag: member="_dsymatrix::damax" ref="abec838aac77ffaf2de6539bc21425e28" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double <a class="el" href="class__dsymatrix.html#abec838aac77ffaf2de6539bc21425e28">damax</a> </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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 largest absolute value </p>

<p>Definition at line <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00052">52</a> of file <a class="el" href="__dsymatrix-calc_8hpp_source.html">_dsymatrix-calc.hpp</a>.</p>
<div class="fragment"><pre class="fragment">{VERBOSE_REPORT;
  <a class="code" href="classdsymatrix.html" title="Real Double-precision Symmetric Matrix Class [l-type (UPLO=l) Strage].">dsymatrix</a> newmat(mat);
  <span class="keywordflow">return</span> <a class="code" href="class__dsymatrix.html#abec838aac77ffaf2de6539bc21425e28">damax</a>(newmat);
}
</pre></div>
</div>
</div>
<a class="anchor" id="a1791c142cdb35a6981bdc4a12f6f3455"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a1791c142cdb35a6981bdc4a12f6f3455" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">const <a class="el" href="class__dsymatrix.html">_dsymatrix</a>&amp; operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>+_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-unary_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-unary_8hpp_source.html">_dsymatrix-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="a916a82da86a46ee0d61fc281ece79714"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a916a82da86a46ee0d61fc281ece79714" args="(const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>-_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-unary_8hpp_source.html#l00010">10</a> of file <a class="el" href="__dsymatrix-unary_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>*mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){ mat.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>[<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>]=-mat.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>[<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>]; }
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ae2dfab27868426c693ebecd09b8c7489"></a><!-- doxytag: member="_dsymatrix::operator+" ref="ae2dfab27868426c693ebecd09b8c7489" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix+dgematrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++ ){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)+=matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<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="a4ece796a5cfa4b73fdbd0ca8c277f04f"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a4ece796a5cfa4b73fdbd0ca8c277f04f" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dgematrix+dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++) {
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++) {
      matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)+=matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> matB;
}
</pre></div>
</div>
</div>
<a class="anchor" id="acb40c53507cbcb1e7cd541fc0e55e442"></a><!-- doxytag: member="_dsymatrix::operator+" ref="acb40c53507cbcb1e7cd541fc0e55e442" args="(const _dsymatrix &amp;, const dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix+dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      matA.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] +=matB.<a class="code" href="classdsymatrix.html#a786e420fb2551cffb147af789197a681" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
    }
  }
  
  <span class="keywordflow">return</span> matA;
}
</pre></div>
</div>
</div>
<a class="anchor" id="ab49f1ed04f239df3b9b1d30f9e9ccc21"></a><!-- doxytag: member="_dsymatrix::operator+" ref="ab49f1ed04f239df3b9b1d30f9e9ccc21" args="(const _dsymatrix &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix+_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      matA.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] +=matB.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
    }
  }

  matB.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> matA;
}
</pre></div>
</div>
</div>
<a class="anchor" id="adae8f965901ba1791a5fc1c629668868"></a><!-- doxytag: member="_dsymatrix::operator+" ref="adae8f965901ba1791a5fc1c629668868" args="(const _dsymatrix &amp;, 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> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix+dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) = newmat(j,<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>) = matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)+=matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<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="a3d0715719a516fb07f42490224b49fb9"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a3d0715719a516fb07f42490224b49fb9" args="(const _dsymatrix &amp;, 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> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix+_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) = newmat(j,<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>) = matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)+=matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a9aed430aafed12429123b8ffaa3547c7"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a9aed430aafed12429123b8ffaa3547c7" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix+dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a> || matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">to_dgematrix</a>() );
  <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++){
    newmat(it-&gt;i,it-&gt;j) += it-&gt;v;
  }
  
  matA.<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="ae66d01b7b77643ba8c1ef8c180c6eb1a"></a><!-- doxytag: member="_dsymatrix::operator+" ref="ae66d01b7b77643ba8c1ef8c180c6eb1a" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix+_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a> || matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">to_dgematrix</a>() );
  <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++){
    newmat(it-&gt;i,it-&gt;j) += it-&gt;v;
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a95c2875a85d6b74ea1fd25574bd94ad8"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a95c2875a85d6b74ea1fd25574bd94ad8" args="(const _dsymatrix &amp;, const dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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="a94eaaffd21f4cce63df3feb343947cc4"></a><!-- doxytag: member="_dsymatrix::operator+" ref="a94eaaffd21f4cce63df3feb343947cc4" args="(const _dsymatrix &amp;, const _dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator+ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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="ae073dbc5e88525decfe52762e31e8b68"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="ae073dbc5e88525decfe52762e31e8b68" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix-dgematrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html#l00027">27</a> of file <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++ ){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)+=matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<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="a873eaf16f7c95d5097e38030e1a0170d"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a873eaf16f7c95d5097e38030e1a0170d" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix-dgematrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html#l00026">26</a> of file <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++) {
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++) {
      matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) =matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)-matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> matB;
}
</pre></div>
</div>
</div>
<a class="anchor" id="acb8b2c3a3aa1fa010ded175f5d7ef0bb"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="acb8b2c3a3aa1fa010ded175f5d7ef0bb" args="(const _dsymatrix &amp;, const dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix-dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html#l00025">25</a> of file <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      matA.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] -=matB.<a class="code" href="classdsymatrix.html#a786e420fb2551cffb147af789197a681" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
    }
  }

  <span class="keywordflow">return</span> matA;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a146298da44c82fd8f2753f3a965bcb3e"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a146298da44c82fd8f2753f3a965bcb3e" args="(const _dsymatrix &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix-_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00026">26</a> of file <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      matA.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] -=matB.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>];
    }
  }
  
  matB.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">destroy</a>();
  <span class="keywordflow">return</span> matA;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a92948829a98575861e6556f62b7e9c62"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a92948829a98575861e6556f62b7e9c62" args="(const _dsymatrix &amp;, 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> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix-dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html#l00030">30</a> of file <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matB.<a class="code" href="classdgbmatrix.html#af8a3e6bc4557b151f0049b8224f54580" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) = newmat(j,<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>) = matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>+matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>+1); j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j)-=matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }
  
  matA.<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="a126c457d3d7e16dd0c6af841df491a33"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a126c457d3d7e16dd0c6af841df491a33" args="(const _dsymatrix &amp;, 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> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix-_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html#l00031">31</a> of file <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matB.<a class="code" href="class__dgbmatrix.html#a587a07ecbbfc6bbdb0e05574ea6b088a" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; j++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) = newmat(j,<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>) = matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) -= matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j);
    }
  }

  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a5e102121931a8b18b1e3b082b40ad646"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a5e102121931a8b18b1e3b082b40ad646" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix-dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html#l00025">25</a> of file <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#ae2d6c9f9657d44171a4bfd56276277cb" title="matrix row size">m</a> || matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="classdgsmatrix.html#acefae40b42444c17eebbc23ccee73ae0" title="matrix column size">n</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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">to_dgematrix</a>() );
  <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++){
    newmat(it-&gt;i,it-&gt;j) -= it-&gt;v;
  }
  
  matA.<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="af7068f5e29904707b5260316afa6152b"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="af7068f5e29904707b5260316afa6152b" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix-_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html#l00026">26</a> of file <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#ac45fcfd8d9cddb1fae5b66cd9b7ac98a" title="matrix row size">m</a> || matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>!=matB.<a class="code" href="class__dgsmatrix.html#a8457e7af3a1e2ff7d8c3238260d59936" title="matrix column size">n</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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#acdf4ff1ebc021a28f63351bd6f6a368c">to_dgematrix</a>() );
  <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++){
    newmat(it-&gt;i,it-&gt;j) -= it-&gt;v;
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="ad81c3bfb85f1ebd6f38bed5e746a3c19"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="ad81c3bfb85f1ebd6f38bed5e746a3c19" args="(const _dsymatrix &amp;, const dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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="a4bb60887d2c29b980bf361dc35940800"></a><!-- doxytag: member="_dsymatrix::operator&#45;" ref="a4bb60887d2c29b980bf361dc35940800" args="(const _dsymatrix &amp;, const _dssmatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator- </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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="af51cdd88dbb5453f102e096d05cc73a5"></a><!-- doxytag: member="_dsymatrix::operator*" ref="af51cdd88dbb5453f102e096d05cc73a5" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*dcovector operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dcovector_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-dcovector_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  dsymv_( <span class="charliteral">&#39;l&#39;</span>, mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 1.0, mat.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>,
          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__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a3f3f8a2d27f873d0842cc3599c0cd860"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a3f3f8a2d27f873d0842cc3599c0cd860" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*_dcovector operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dcovector_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-__dcovector_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  dsymv_( <span class="charliteral">&#39;l&#39;</span>, mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 1.0, mat.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>,
          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__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="afda93bc54596d73f3921ca6572790814"></a><!-- doxytag: member="_dsymatrix::operator*" ref="afda93bc54596d73f3921ca6572790814" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*dgematrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html#l00051">51</a> of file <a class="el" href="__dsymatrix-dgematrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a> );
  dsymm_( <span class="charliteral">&#39;l&#39;</span>, <span class="charliteral">&#39;l&#39;</span>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="classdgematrix.html#a60d26e3a499a235310e59a0aebbe0d0a" title="matrix column size">n</a>, 1.0, matA.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 
          matB.<a class="code" href="classdgematrix.html#a006ce8c2f9add380bcb2816fdb13aa01" title="1D array to store matrix data">array</a>, matB.<a class="code" href="classdgematrix.html#a3852aa1d0b876a74c1070e982cf388a2" title="matrix row size">m</a>, 0.0, newmat.array, newmat.m );
  
  matA.<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="a977ef4982892149eb6bb4794e7f72a0f"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a977ef4982892149eb6bb4794e7f72a0f" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dgematrix*dgematrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html#l00049">49</a> of file <a class="el" href="__dsymatrix-__dgematrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a> );
  dsymm_( <span class="charliteral">&#39;l&#39;</span>, <span class="charliteral">&#39;l&#39;</span>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="class__dgematrix.html#ad46b02985eb9537b1bd65937d5b052de" title="matrix column size">n</a>, 1.0, matA.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 
          matB.<a class="code" href="class__dgematrix.html#ad01c168fc4997d06e7c996ba36b12ff6" title="1D array to store matrix data">array</a>, matB.<a class="code" href="class__dgematrix.html#a4c07ace0e1899332aa2d620af184dadf" title="matrix row size">m</a>, 0.0, newmat.array, newmat.m );
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a8d2bf7e9fdf57f365e74ed38313ba8d0"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a8d2bf7e9fdf57f365e74ed38313ba8d0" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html#l00047">47</a> of file <a class="el" href="__dsymatrix-dsymatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>  
  matA.<a class="code" href="class__dsymatrix.html#abfd55ec47814e32e991ddd3c8de6b1f8">complete</a>();
  matB.<a class="code" href="classdsymatrix.html#a09530f2e3c2397751b4d66a6f73ce86b">complete</a>();
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  
  dgemm_( <span class="charliteral">&#39;n&#39;</span>, <span class="charliteral">&#39;n&#39;</span>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 1.0, matA.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>,
          matB.<a class="code" href="classdsymatrix.html#a7efa7c974733dc3111167a9ebb217e1f" title="1D array to store matrix data">array</a>, matB.<a class="code" href="classdsymatrix.html#a4a2ca7ec3f2b7bd62139777fbb6ad1a9" title="matrix column size">n</a>, 0.0, newmat.array, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> );
  
  matA.<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="aa28788252525cda33141d129fd6e6a09"></a><!-- doxytag: member="_dsymatrix::operator*" ref="aa28788252525cda33141d129fd6e6a09" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*_dsymatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00049">49</a> of file <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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>  
  matA.<a class="code" href="class__dsymatrix.html#abfd55ec47814e32e991ddd3c8de6b1f8">complete</a>();
  matB.<a class="code" href="class__dsymatrix.html#abfd55ec47814e32e991ddd3c8de6b1f8">complete</a>();
  <a class="code" href="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat(matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>);
  
  dgemm_( <span class="charliteral">&#39;n&#39;</span>, <span class="charliteral">&#39;n&#39;</span>, matA.<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>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 1.0, matA.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>,
          matB.<a class="code" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a" title="1D array to store matrix data">array</a>, matB.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, 0.0, newmat.array, matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> );
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="afd3c21b75b16dfaed32f26f7332bd591"></a><!-- doxytag: member="_dsymatrix::operator*" ref="afd3c21b75b16dfaed32f26f7332bd591" args="(const _dsymatrix &amp;, 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> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dgematrix*dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html#l00057">57</a> of file <a class="el" href="__dsymatrix-dgbmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="classdgbmatrix.html#a5b6fe4d1cfa71161b09ecf0b59befd0e" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;newmat.<a class="code" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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),j-matB.<a class="code" href="classdgbmatrix.html#a0b4d519e95ddf780bf1151a42018d646" title="upper band width">ku</a>); k&lt;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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) +=matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<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="a1e50366c55b277b6677bdcb87397ee5b"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a1e50366c55b277b6677bdcb87397ee5b" args="(const _dsymatrix &amp;, 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> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix*_dgbmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html#l00059">59</a> of file <a class="el" href="__dsymatrix-__dgbmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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="classdgematrix.html" title="Real Double-precision General Dence Matrix Class.">dgematrix</a> newmat( matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>, matB.<a class="code" href="class__dgbmatrix.html#a95aebc7e2d4ea66b1b7a08e981df6971" title="matrix column size">n</a> );
  newmat.zero();
  
  <span class="keywordflow">for</span>(<span class="keywordtype">long</span> <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;newmat.<a class="code" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f" title="matrix row size">m</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">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),j-matB.<a class="code" href="class__dgbmatrix.html#a116dadc8ed5246ee0b4cd0a9fffc53ea" title="upper band width">ku</a>); k&lt;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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,j) +=matA(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,k)*matB(k,j);
      }
    }
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a577fce177108b8fd791ac8a73813f132"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a577fce177108b8fd791ac8a73813f132" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html#l00047">47</a> of file <a class="el" href="__dsymatrix-dgsmatrix_8hpp_source.html">_dsymatrix-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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,it-&gt;j) += matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,it-&gt;<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>)*it-&gt;v;
    }
  }
  
  matA.<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="ad9c9cafc9a48b0f44b9b39ccf5d149d8"></a><!-- doxytag: member="_dsymatrix::operator*" ref="ad9c9cafc9a48b0f44b9b39ccf5d149d8" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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>_dsymatrix*_dgsmatrix operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html#l00049">49</a> of file <a class="el" href="__dsymatrix-__dgsmatrix_8hpp_source.html">_dsymatrix-_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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a> &lt;&lt; <span class="stringliteral">&quot;x&quot;</span> &lt;&lt; matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" 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__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;matA.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
      newmat(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,it-&gt;j) += matB(<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>,it-&gt;<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>)*it-&gt;v;
    }
  }
  
  matA.<a class="code" href="class__dsymatrix.html#a09e5b4ca8579bef58e01e80c22e799b1">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="a3ea6f20cbc65b317fb261187bf3d8732"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a3ea6f20cbc65b317fb261187bf3d8732" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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="a4b58e624b1d0fd3ceb8976cbbad920a1"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a4b58e624b1d0fd3ceb8976cbbad920a1" args="(const _dsymatrix &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__dsymatrix.html">_dsymatrix</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="aa1c618a4d03afb961c861c64aef29aaf"></a><!-- doxytag: member="_dsymatrix::operator*" ref="aa1c618a4d03afb961c861c64aef29aaf" args="(const _dsymatrix &amp;, const double &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator* </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>_dsymatrix*double operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-double_8hpp_source.html#l00003">3</a> of file <a class="el" href="__dsymatrix-double_8hpp_source.html">_dsymatrix-double.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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      mat.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] *=d;
    }
  }
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a4acbdbaa8ba95f7ef2b461ac116ad4ae"></a><!-- doxytag: member="_dsymatrix::operator/" ref="a4acbdbaa8ba95f7ef2b461ac116ad4ae" args="(const _dsymatrix &amp;, const double &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</a> operator/ </td>
          <td>(</td>
          <td class="paramtype">const <a class="el" href="class__dsymatrix.html">_dsymatrix</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>dsymatrix/double operator </p>

<p>Definition at line <a class="el" href="__dsymatrix-double_8hpp_source.html#l00015">15</a> of file <a class="el" href="__dsymatrix-double_8hpp_source.html">_dsymatrix-double.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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      mat.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] /=d;
    }
  }
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<a class="anchor" id="a53c70446bb321088fa2e2a79eba8b58c"></a><!-- doxytag: member="_dsymatrix::operator*" ref="a53c70446bb321088fa2e2a79eba8b58c" args="(const double &amp;, const _dsymatrix &amp;)" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname"><a class="el" href="class__dsymatrix.html">_dsymatrix</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__dsymatrix.html">_dsymatrix</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*_dsymatrix operator </p>

<p>Definition at line <a class="el" href="double-__dsymatrix_8hpp_source.html#l00003">3</a> of file <a class="el" href="double-__dsymatrix_8hpp_source.html">double-_dsymatrix.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__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>=0; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>&lt;mat.<a class="code" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b" title="matrix column size">n</a>; <a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>++){
    <span class="keywordflow">for</span>(<span class="keywordtype">long</span> j=0; j&lt;=<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>; j++){
      mat.<a class="code" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc" title="array of pointers of column head addresses">darray</a>[j][<a class="code" href="class__dsymatrix.html#a0697d2173dfc389bf1312d7bf6a5bd3b">i</a>] *=d;
    }
  }
  
  <span class="keywordflow">return</span> mat;
}
</pre></div>
</div>
</div>
<hr/><h2>Member Data Documentation</h2>
<a class="anchor" id="a0ae0bb2812e15cb1a0667fd6f6f6b44f"></a><!-- doxytag: member="_dsymatrix::m" ref="a0ae0bb2812e15cb1a0667fd6f6f6b44f" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long const&amp; <a class="el" href="class__dsymatrix.html#a0ae0bb2812e15cb1a0667fd6f6f6b44f">_dsymatrix::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="__dsymatrix_8hpp_source.html#l00009">9</a> of file <a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="__dgsmatrix-__dsymatrix_8hpp_source.html#l00052">operator*()</a>, <a class="el" href="__dgsmatrix-__dsymatrix_8hpp_source.html#l00003">operator+()</a>, and <a class="el" href="__dgsmatrix-__dsymatrix_8hpp_source.html#l00026">operator-()</a>.</p>

</div>
</div>
<a class="anchor" id="a8b56cdf252a15da3e211b558abd46c1b"></a><!-- doxytag: member="_dsymatrix::n" ref="a8b56cdf252a15da3e211b558abd46c1b" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">long <a class="el" href="class__dsymatrix.html#a8b56cdf252a15da3e211b558abd46c1b">_dsymatrix::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="__dsymatrix_8hpp_source.html#l00010">10</a> of file <a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dsymatrix-misc_8hpp_source.html#l00158">_()</a>, <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00003">_dsymatrix()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00022">complete()</a>, <a class="el" href="dsymatrix-constructor_8hpp_source.html#l00029">dsymatrix::dsymatrix()</a>, <a class="el" href="__dsymatrix-calc_8hpp_source.html#l00015">i()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dsymatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00059">operator*()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00049">dgematrix::operator*=()</a>, <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00015">dsymatrix::operator+=()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00003">dgematrix::operator+=()</a>, <a class="el" href="__dgbmatrix-__dsymatrix_8hpp_source.html#l00031">operator-()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00038">dsymatrix::operator-=()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00026">dgematrix::operator-=()</a>, <a class="el" href="__dsymatrix-double_8hpp_source.html#l00015">operator/()</a>, <a class="el" href="__dsymatrix-io_8hpp_source.html#l00028">operator&lt;&lt;()</a>, <a class="el" href="dsymatrix-misc_8hpp_source.html#l00066">dsymatrix::shallow_copy()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00022">to_dgematrix()</a>, <a class="el" href="__dsymatrix-cast_8hpp_source.html#l00003">to_zhematrix()</a>, and <a class="el" href="__dsymatrix-io_8hpp_source.html#l00050">write()</a>.</p>

</div>
</div>
<a class="anchor" id="a656ab414ee51638aedc9bf4960fe428a"></a><!-- doxytag: member="_dsymatrix::array" ref="a656ab414ee51638aedc9bf4960fe428a" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double* <a class="el" href="class__dsymatrix.html#a656ab414ee51638aedc9bf4960fe428a">_dsymatrix::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="__dsymatrix_8hpp_source.html#l00011">11</a> of file <a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dsymatrix-misc_8hpp_source.html#l00158">_()</a>, <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00003">_dsymatrix()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00012">destroy()</a>, <a class="el" href="dsymatrix-constructor_8hpp_source.html#l00029">dsymatrix::dsymatrix()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dgematrix-__dsymatrix_8hpp_source.html#l00049">operator*()</a>, <a class="el" href="dgematrix-__dsymatrix_8hpp_source.html#l00049">dgematrix::operator*=()</a>, <a class="el" href="__dsymatrix-unary_8hpp_source.html#l00010">operator-()</a>, <a class="el" href="dsymatrix-misc_8hpp_source.html#l00066">dsymatrix::shallow_copy()</a>, and <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00029">~_dsymatrix()</a>.</p>

</div>
</div>
<a class="anchor" id="a9db2d341c99952e75adb7e183450c3dc"></a><!-- doxytag: member="_dsymatrix::darray" ref="a9db2d341c99952e75adb7e183450c3dc" args="" -->
<div class="memitem">
<div class="memproto">
      <table class="memname">
        <tr>
          <td class="memname">double** <a class="el" href="class__dsymatrix.html#a9db2d341c99952e75adb7e183450c3dc">_dsymatrix::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="__dsymatrix_8hpp_source.html#l00012">12</a> of file <a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a>.</p>

<p>Referenced by <a class="el" href="dsymatrix-misc_8hpp_source.html#l00158">_()</a>, <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00003">_dsymatrix()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00022">complete()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00012">destroy()</a>, <a class="el" href="dsymatrix-constructor_8hpp_source.html#l00029">dsymatrix::dsymatrix()</a>, <a class="el" href="__dsymatrix-misc_8hpp_source.html#l00003">nullify()</a>, <a class="el" href="__dsymatrix-io_8hpp_source.html#l00003">operator()()</a>, <a class="el" href="__dsymatrix-double_8hpp_source.html#l00003">operator*()</a>, <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00003">operator+()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00015">dsymatrix::operator+=()</a>, <a class="el" href="__dsymatrix-__dsymatrix_8hpp_source.html#l00026">operator-()</a>, <a class="el" href="dsymatrix-__dsymatrix_8hpp_source.html#l00038">dsymatrix::operator-=()</a>, <a class="el" href="__dsymatrix-double_8hpp_source.html#l00015">operator/()</a>, <a class="el" href="dsymatrix-misc_8hpp_source.html#l00066">dsymatrix::shallow_copy()</a>, and <a class="el" href="__dsymatrix-constructor_8hpp_source.html#l00029">~_dsymatrix()</a>.</p>

</div>
</div>
<hr/>The documentation for this class was generated from the following files:<ul>
<li><a class="el" href="__dsymatrix_8hpp_source.html">_dsymatrix.hpp</a></li>
<li><a class="el" href="__dsymatrix-cast_8hpp_source.html">_dsymatrix-cast.hpp</a></li>
<li><a class="el" href="__dsymatrix-constructor_8hpp_source.html">_dsymatrix-constructor.hpp</a></li>
<li><a class="el" href="__dsymatrix-io_8hpp_source.html">_dsymatrix-io.hpp</a></li>
<li><a class="el" href="__dsymatrix-misc_8hpp_source.html">_dsymatrix-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>
