Emmanuel: This is a bug in petsc. I've pushed a fix https://gitlab.com/petsc/petsc/-/commit/fd2a003f2c07165526de5c2fa5ca4f3c85618da7
You can edit it in your petsc library, or add MatAssemblyBegin/End in your application code until petsc-release is patched. Thanks for reporting it and sending us the test! Hong Thanks for the answer. > > Emmanuel: >> You can create a dense C with the required parallel layout without >> calling MatAssemblyBegin() and MatAssemblyEnd(). >> Did you get error without calling these routines? >> > > Yes, the output is (after create the C dense matrix and do not assembly > it, run1 - see attached file -): > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > [0]PETSC ERROR: Object is in wrong state > [0]PETSC ERROR: Not for unassembled matrix > [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html > for trouble shooting. > [0]PETSC ERROR: Petsc Release Version 3.12.4, Feb, 04, 2020 > [0]PETSC ERROR: ./comp on a arch-linux-c-opt-O2 named ayala by ayala Thu > Feb 27 16:47:15 2020 > [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 > [0]PETSC ERROR: #1 MatNorm() line 5123 in > /home/ayala/Documents/PETSc/petsc-3.12.4/src/mat/interface/matrix.c > > 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. >> Hong >> > > I attach my test file, or maybe I'm doing something wrong. I tested this > file on my laptop ubuntu 18 > > Kind regards. > >
