<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Nevermind, this is the wrong fix.<div class="">The proper fix is in PETSc. It should not error out if the matrix is also symmetric.</div><div class="">Indeed, complex symmetric Hermitian => complex with no imaginary part.</div><div class="">Thus all operations like MatMult, MatMultHermitianTranspose, Cholesky… will work for bs > 1, since all is filled with zeroes.</div><div class="">I will take care of this, I’m c/c’ing petsc-dev so that they don’t have to “reverse engineer” the trivial change to MatSetOption_SeqSBAIJ.</div><div class=""><br class=""></div><div class="">Sorry about the noise.</div><div class=""><br class=""></div><div class="">Thank you,</div><div class="">Pierre</div><div class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 10 Sep 2019, at 8:37 AM, Pierre Jolivet <<a href="mailto:pierre.jolivet@enseeiht.fr" class="">pierre.jolivet@enseeiht.fr</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello,<div class="">Could you consider not setting MAT_HERMITIAN here <a href="http://slepc.upv.es/documentation/current/src/sys/classes/st/interface/stsles.c.html#line276" class="">http://slepc.upv.es/documentation/current/src/sys/classes/st/interface/stsles.c.html#line276</a> when using SBAIJ matrices with bs > 1?</div><div class="">This makes PETSc error out with</div><div class=""><div class="">#<span class="Apple-tab-span" style="white-space:pre"> </span>[1]PETSC ERROR: No support for this operation for this object type</div><div class="">#<span class="Apple-tab-span" style="white-space:pre"> </span>[1]PETSC ERROR: No support for Hermitian with block size greater than 1</div></div><div class=""><br class=""></div><div class="">The change does not bring any regression, since PETSc is always giving an error without it, but on the contrary, it improves the range of applicability of SLEPc, e.g., for complex Hermitian problems with SBAIJ matrices and bs > 1 that _don’t_ require the flag MAT_HERMITIAN set to true.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Pierre</div></div></div></blockquote></div><br class=""></div></body></html>