<div dir="ltr">Just an update: I ran with -info, and saw that the 186x6 matrix uses compressed rows routines, whereas 171x21 does not.<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/10/9 Stefano Zampini <span dir="ltr"><<a href="mailto:s.zampini@cineca.it" target="_blank">s.zampini@cineca.it</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div>I just noticed unexpected behaviour of MatMult_SeqBAIJ_3 in PCBDDC code (on branch stefano_zampini/pcbddc-baijfixes) with rectangular matrices having number of rows >> number of cols and number of cols =2*blocksize. I guess this problem will show up for other values of block sizes.<br>

<br>When some rows of the matrix are empty, the corresponding entry in 
the resulting vector is not zero (as expected, mathematically speaking);
 instead, the values are not overwritten at all. Is it the expected 
behaviour?<br>
<br>Attached are two matrices . The one with sizes 186x6 shows the problem on my laptop (ubuntu, gnu compilers,  --with-debugging, , CFLAGS -g -O0)<br><br></div>I'm not familiar with PETSc low-level implementation of BAIJ, so I cannot dig anymore into the problem apart saying that: <br clear="all">

<br>A) if "converted" to SeqAIJ with a call <br><br>MatConvert(Mymat,MATSEQAIJ,MAT_REUSE_MATRIX,&Mymat)<br><br>the output of MatMult is always correct.<br><br>B) the problem disappear by increasing the number of columns (the matrix with sizes 171x21 does not show the problem).<br>

<br><br></div><div>Hope you can reproduce the problem,<br></div><div><div><div><div>-- <br><div dir="ltr"><pre cols="72"><font><span style="font-family:garamond,serif"><span style="font-family:garamond,serif"></span>Ph. D. Stefano Zampini<br>


CINECA  
SuperComputing Applications and Innovations Department - SCAI
Via dei Tizii, 6 00185 Roma - ITALY
------------------------------------------------------------------------------------------------------------------------
Email: <a href="mailto:s.zampini@cineca.it" target="_blank"><span>s.zampini@cineca.it</span></a>
SkypeID: stefano.zampini
GoogleTalk: <a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>
Tel: <a value="+390644486720">+39 06.44486.707<br></a></span></font><a value="+390644486720"><font><span style="font-family:garamond,serif">------------------------------------------------------------------------------------------------------------------------
</span></font><br></a></pre></div>
</div></div></div></div></div>
</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><pre cols="72"><font><span style="font-family:garamond,serif"><span style="font-family:garamond,serif"></span>Ph. D. Stefano Zampini<br>CINECA  
SuperComputing Applications and Innovations Department - SCAI
Via dei Tizii, 6 00185 Roma - ITALY
------------------------------------------------------------------------------------------------------------------------
Email: <a href="mailto:s.zampini@cineca.it" target="_blank"><span>s.zampini@cineca.it</span></a>
SkypeID: stefano.zampini
GoogleTalk: <a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>
Tel: <a value="+390644486720">+39 06.44486.707<br></a></span></font><a value="+390644486720"><font><span style="font-family:garamond,serif">------------------------------------------------------------------------------------------------------------------------
</span></font><br></a></pre></div>
</div>