[petsc-users] MatMatMult

[hzhang at mcs.anl.gov]

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.
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200228/489fabcf/attachment.html>