<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hong,<div class="">You should probably cherry pick <a href="https://gitlab.com/petsc/petsc/commit/93d7d1d6d29b0d66b5629a261178b832a925de80?merge_request_iid=2069" class="">https://gitlab.com/petsc/petsc/commit/93d7d1d6d29b0d66b5629a261178b832a925de80?merge_request_iid=2069</a> (and remove the MatNest part).</div><div class="">This fixes a similar issue in MatTransposeMatMult with nontrivial LDAs.</div><div class="">Since this commit is part of a feature MR that is unlikely to be ready for tomorrow, this fix (as of now) is also unlikely to be in master for the release.</div><div class=""><br class=""></div><div class="">Thanks,</div><div class="">Pierre<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 23 Sep 2019, at 6:02 PM, Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov" class="">hzhang@mcs.anl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">Barry:</div>
<div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
As a hack for this release could you have the Numeric portion of the multiply routines check if the symbolic data is there and if not just call the symbolic an attach the needed data? You might need to have a utility function that does all the symbolic
part except the allocation of the matrix and then call this from the numeric part as well as the real symbolic part.<br class="">
</blockquote>
<div class=""> </div>
<div class="">I'm working on this now. I was not aware of MatSeqDenseSetLDA() which changes pattern of data access in seqdense matrix. </div>
<div class=""><span style="font-family:Roboto,RobotoDraft,Helvetica,Arial,sans-serif;font-size:14px" class="">Pierre's patch:</span><br class="">
</div>
<div class="">"change Bm here <a href="https://www.mcs.anl.gov/petsc/petsc-dev/src/mat/impls/aij/mpi/mpimatmatmult.c.html#line549" target="_blank" class="">https://www.mcs.anl.gov/petsc/petsc-dev/src/mat/impls/aij/mpi/mpimatmatmult.c.html#line549</a> to the LDA of B" </div>
<div class="">fix this bug. I'll further test it and submit a pull request. </div>
<div class="">Then, I'll check slepc's bug report.</div>
<div class="">Hong</div>
</div>
</div>
</div>
</div></blockquote></div><br class=""></div></body></html>