<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="ltr"><div dir="ltr">Emmanuel:<br></div><div>This is a bug in petsc. I've pushed a fix <a href="https://gitlab.com/petsc/petsc/-/commit/fd2a003f2c07165526de5c2fa5ca4f3c85618da7">https://gitlab.com/petsc/petsc/-/commit/fd2a003f2c07165526de5c2fa5ca4f3c85618da7</a></div><div><br></div><div>You can edit it in your petsc library, or add <span class="gmail-n">MatAssemblyBegin/End in your application code until petsc-release is patched.</span></div><div><span class="gmail-n">Thanks for reporting it and sending us the test!</span></div><div><span class="gmail-n">Hong</span></div><div><span class="gmail-n"><br></span></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"><div dir="ltr"><div>Thanks for the answer.</div><div><br></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"><div dir="ltr"><div dir="ltr">Emmanuel:<br></div><div>You can create a dense C with the required parallel layout without calling MatAssemblyBegin() and MatAssemblyEnd().</div><div>Did you get error without calling these routines?</div></div></blockquote><div><br></div><div>Yes, the output is (after create the C dense matrix and do not assembly it, run1 - see attached file -):</div><div><br></div><div>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: Object is in wrong state<br>[0]PETSC ERROR: Not for unassembled matrix<br>[0]PETSC ERROR: See <a href="https://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">https://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.12.4, Feb, 04, 2020 <br>[0]PETSC ERROR: ./comp on a arch-linux-c-opt-O2 named ayala by ayala Thu Feb 27 16:47:15 2020<br>[0]PETSC ERROR: Configure options --with-debugging=0 COPTFLAGS="-O2 -march=native -mtune=native" CXXOPTFLAGS="-O2 -march=native -mtune=native" FOPTFLAGS="-O2 -march=native -mtune=native" --download-mpich=1 --download-fblaslapack=1 --with-cxx-dialect=C++11<br>[0]PETSC ERROR: #1 MatNorm() line 5123 in /home/ayala/Documents/PETSc/petsc-3.12.4/src/mat/interface/matrix.c</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>We only updated the help manu, not internal implementation. In the next release, we'll introduce new set of API to consolidate the API of mat-mat-operations.</div><div>Hong</div></div></blockquote><div><br></div><div>I attach my test file, or maybe I'm doing something wrong. I tested this file on my laptop ubuntu 18</div><div><br></div><div>Kind regards.<br></div><div><br></div></div></div>
</blockquote></div></div>