[petsc-users] Petsc 3.19.2 to 3.18.0 error possibility.

[Barry Smith]

  Could you send us the exact error output that occurs? Cut and paste the run command and the entire error message.

   Also send the configure options you used. Have you tried configuring the later PETSc versions with all optimization turned off; use --with-debugging=1 --with-cflags='-g -O0'   --with-cxxflags='-g -O0' --with-fflags='-g -O0'  Does the same error occur?

  Barry




> On Jun 23, 2023, at 5:30 PM, Paweł Stebliński via petsc-users <petsc-users at mcs.anl.gov> wrote:
> 
> I am micromagnetic (MAGPAR) software developer. Old Magpar version has been using petsc 3.1-p8. I have decided to upgrade to petsc 3.19.2 with avx512 support. Unfortunately there appeared an error during software testing. Error appeares in ranning code after proper compiling and linking. The bug is in a code part which initializes matrix in petsc library versions: 3.19.2, 3.19.1, 3.19.0 and also from 3.18.5 to 3.18.0. If we use petsc version 3.17.5 the error doesn't appear. With this version (3.17.5) all is ok and simulation is running without any errors or throwing exceptions. My guess is linked to avx512 implementation which is good up to 3.17.5 version and buggy in upper mentioned versions with higher numbers. Avx512 is buggy according to tested SeqAij matrices.
> 
> The exception is not thrown if we comment code fragment below.
> 
> ierr = MatCreateSeqAIJ(
>     PETSC_COMM_SELF,
>     nvert,nvert,
>     0,ia,
>     &mat
>   );CHKERRQ(ierr);
>   ierr = MatSetFromOptions(mat);CHKERRQ(ierr);
> 
> ia - is number of nonzeros array which is obtained according parmetis partitioning. There were the same version of parmetis  (3.1.1) in the all considered cases.
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230623/40c2c6b0/attachment.html>