[petsc-users] Problem with BLASdot in 3.19.4

[Thuc Bui]

Dear Petsc users and experts,

 

If someone can direct me how to track this bug, I would really appreciate
it.

 

The Petsc DLL library version 3.19.4 was built on Windows 10 with Visual
Studio 2022, and with Microsoft MPI 10.1.2 and Intel MKL 2020.3.279. The
same code works fine with Petsc 3.18.6 using the same versions of MS MPI and
Intel MKL, and built with Visual Studio 2019.

 

When my code calls PetscCall(KSPSolve(...)), it gets to
PetscCall(KSPGMRESCycle(...)), PetscCall(VecNormalize(VEC_VV(0), &res)),
PetscCall(VecNorm(x, NORM_2, &norm)) and then crashes in VecNorm_Seq() at:

 

                PetscCallBLAS("BLASdot", ztmp[type == NORM_1_AND_2] =
PetscSqrtReal(PetscRealPart(BLASdot_(&bn, xx, &one, xx, &one))));

 

I tried to step into BLASdot_, but was unable to. I assume BLASdot_ belong
to Intel MKL, and its library does not contain debug information. After the
code exits, the errors are shown below the hash line.

 

Should I install the latest Intel MKL to work with Petsc 3.19.4?

 

Many thanks in advance for your help,

Thuc Bui

Senior R&D Engineer

Calabazas Creek Research, Inc

(650) 948-5361

 

 

###############################################

 

[0]PETSC ERROR:
------------------------------------------------------------------------

[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,
probably memory access out of range

[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger

[0]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind and
https://petsc.org/release/faq/

[0]PETSC ERROR: ---------------------  Stack Frames
------------------------------------

[0]PETSC ERROR: The line numbers in the error traceback are not always
exact.

[0]PETSC ERROR: #1 BLASdot()

[0]PETSC ERROR: #2 VecNorm_Seq() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\vec\vec\impls\seq\bvec2.c:213

[0]PETSC ERROR: #3 VecNorm() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\vec\vec\interface\rvector.c:23
7

[0]PETSC ERROR: #4 VecNormalize() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\vec\vec\interface\rvector.c:31
8

[0]PETSC ERROR: #5 KSPGMRESCycle() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\ksp\ksp\impls\gmres\gmres.c:11
0

[0]PETSC ERROR: #6 KSPSolve_GMRES() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\ksp\ksp\impls\gmres\gmres.c:22
7

[0]PETSC ERROR: #7 KSPSolve_Private() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\ksp\ksp\interface\itfunc.c:898

[0]PETSC ERROR: #8 KSPSolve() at
D:\Users\bbwannabe\Documents\Petsc\latest\src\ksp\ksp\interface\itfunc.c:107
0

 

job aborted:

[ranks] message

 

[0] application aborted

aborting MPI_COMM_WORLD (comm=0x44000000), error 59, comm rank 0

 

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