<div dir="ltr">Hello,<div><br></div><div>If we look at lines 2330-2331 in file baij2.c, it looks like there are some </div><div>mistakes in assigning the `sum..` variables to the z array, causing</div><div>the function MatMultAdd_SeqBAIJ_11() to not produce the correct</div><div>answer.</div><div><br></div><div>I don't have a good example program to demonstrate this yet - it's</div><div>currently causing problems in a dev branch of pflotan_ogs that</div><div>can produce blocksize 11 matrices. When in parallel, a standard</div><div>matrix-vector multiplication calls MatMultAdd for the off-proc</div><div>contributions, and the result is wrong when this is redirected</div><div>to MatMultAdd_SeqBAIJ_11. Seems to be the root cause of</div><div>several solvers failing such as fgmres.</div><div><br></div><div>Can anyone confirm that these two lines seem incorrect?</div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Daniel</div></div>