[petsc-users] MatMatMult

[Emmanuel Ayala]

Hi, Thanks.

Kind regards.

El vie., 28 de feb. de 2020 a la(s) 10:04, hzhang at mcs.anl.gov (
hzhang at 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.
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200228/2f4c2f7c/attachment-0001.html>