the last block column</span></div>
<div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;        {</div>
<div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;                IT * blocknnz = <span class="keyword">new</span> IT[A.nbr];  <span class="comment">// nnz per block responsible</span></div>
<div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;                <span class="keywordflow">for</span>(IT i=0; i&lt;A.nbr; ++i)</div>
<div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160;                {</div>
<div class="line"><a name="l00422"></a><span class="lineno">  422</span>&#160;                        blocknnz[i] = A.top[i][j] - A.top[i][j-1];</div>
<div class="line"><a name="l00423"></a><span class="lineno">  423</span>&#160;                }</div>
<div class="line"><a name="l00424"></a><span class="lineno">  424</span>&#160;                sum[j-1] = std::accumulate(blocknnz, blocknnz + (A.nbr-1), 0);         <span class="comment">// ignore the last block row</span></div>
<div class="line"><a name="l00425"></a><span class="lineno">  425</span>&#160;                <span class="keyword">delete</span> [] blocknnz;</div>
<div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;        }</div>
<div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;        <span class="keywordtype">float</span> colave = std::accumulate(sum.begin(), sum.end(), 0.0) / <span class="keyword">static_cast&lt;</span><span class="keywordtype">float</span><span class="keyword">&gt;</span>(A.nbc-1);</div>
<div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;        vector&lt;float&gt;::iterator colmax = std::max_element(sum.begin(), sum.end());</div>
<div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;        <span class="keywordflow">return</span> (*colmax) / colave;</div>
<div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;}</div>
<div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;</div>
<div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;</div>
<div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160;<span class="preprocessor"></span></div>
<div class="ttc" id="csbsym_8h_html"><div class="ttname"><a href="csbsym_8h.html">csbsym.h</a></div></div>
<div class="ttc" id="friends_8h_html_ad6f53d26802d2b6d175082e55c083e84"><div class="ttname"><a href="friends_8h.html#ad6f53d26802d2b6d175082e55c083e84">bicsb_gespmv</a></div><div class="ttdeci">void bicsb_gespmv(const BiCsb&lt; NT, IT &gt; &amp;A, const RHS *__restrict x, LHS *__restrict y)</div><div class="ttdef"><b>Definition:</b> <a href="friends_8h_source.html#l00113">friends.h:113</a></div></div>
<div class="ttc" id="class_bm_csb_html_ac23f83142d9655235ef568971050c98c"><div class="ttname"><a href="class_bm_csb.html#ac23f83142d9655235ef568971050c98c">BmCsb::rowsize</a></div><div class="ttdeci">IT rowsize() const </div><div class="ttdef"><b>Definition:</b> <a href="bmcsb_8h_source.html#l00033">bmcsb.h:33</a></div></div>
<div class="ttc" id="class_bm_sym_html_ab4b8f74f857c58cea834c929d56f786c"><div class="ttname"><a href="class_bm_sym.html#ab4b8f74f857c58cea834c929d56f786c">BmSym::isPar</a></div><div class="ttdeci">bool isPar() const </div><div class="ttdef"><b>Definition:</b> <a href="bmsym_8h_source.html#l00065">bmsym.h:65</a></div></div>
<div class="ttc" id="class_bm_csb_html"><div class="ttname"><a href="class_bm_csb.html">BmCsb</a></div><div class="ttdef"><b>Definition:</b> <a href="bmcsb_8h_source.html#l00021">bmcsb.h:21</a></div></div>
<div class="ttc" id="class_bi_csb_html_a8b8b5238b5805fdf28ac3e5ac3cd415d"><div class="ttname"><a href="class_bi_csb.html#a8b8b5238b5805fdf28ac3e5ac3cd415d">BiCsb::rowsize</a></div><div class="ttdeci">IT rowsize() const </div><div class="ttdef"><b>Definition:</b> <a href="bicsb_8h_source.html#l00034">bicsb.h:34</a></div></div>
<div class="ttc" id="friends_8h_html_a8e7ce87d0bb073f78c728ba9f1510034"><div class="ttname"><a href="friends_8h.html#a8e7ce87d0bb073f78c728ba9f1510034">csbsym_gespmv</a></div><div class="ttdeci">void csbsym_gespmv(const CsbSym&lt; NT, IT &gt; &amp;A, const NT *__restrict x, NT *__restrict y)</div><div class="ttdef"><b>Definition:</b> <a href="friends_8h_source.html#l00281">friends.h:281</a></div></div>
<div class="ttc" id="utility_8h_html_acc2a6872d719d865052a48ab5cffa13d"><div class="ttname"><a href="utility_8h.html#acc2a6872d719d865052a48ab5cffa13d">prescan</a></div><div class="ttdeci">unsigned prescan(unsigned *a, MTYPE *const M, int n)</div><div class="ttdef"><b>Definition:</b> <a href="utility_8h_source.html#l00191">utility.h:191</a></div></div>
<div class="ttc" id="friends_8h_html_a29fb3844b386e91a3e7a10e