Hi, Thanks. Kind regards.
El vie., 28 de feb. de 2020 a la(s) 10:04, hzh...@mcs.anl.gov ( hzh...@mcs.anl.gov) escribió: > 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. >> >>