[petsc-users] Bug report VecNorm
Pierre Jolivet
pierre at joliv.et
Sun Dec 10 12:47:43 CST 2023
> On 10 Dec 2023, at 8:40 AM, Stephan Köhler <stephan.koehler at math.tu-freiberg.de> wrote:
>
> Dear PETSc/Tao team,
>
> there is a bug in the voector interface: In the function
> VecNorm, see, eg. https://petsc.org/release/src/vec/vec/interface/rvector.c.html#VecNorm line 197 the check for consistency in line 214 is done on the wrong communicator. The communicator should be PETSC_COMM_SELF.
> Otherwise the program may hang when PetscCheck is executed.
I think the communicator should not be changed, but instead, the check/conditional should be changed, à la PetscValidLogicalCollectiveBool().
Thanks,
Pierre
> Please find a minimal example attached.
>
>
> Kind regards,
> Stephan Köhler
> --
> Stephan Köhler
> TU Bergakademie Freiberg
> Institut für numerische Mathematik und Optimierung
>
> Akademiestraße 6
> 09599 Freiberg
> Gebäudeteil Mittelbau, Zimmer 2.07
>
> Telefon: +49 (0)3731 39-3188 (Büro)
> <minimal_ex_vec_norm.cpp><OpenPGP_0xC9BF2C20DFE9F713.asc>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20231210/190b0f40/attachment.html>
More information about the petsc-users
mailing list