https://gcc.gnu.org/bugzilla/show_bug.cgi?id=121185
--- Comment #12 from Jürgen Reuter <juergen.reuter at desy dot de> --- (In reply to kargls from comment #11) > (In reply to kargls from comment #10) > > (In reply to Jürgen Reuter from comment #8) > > > Created attachment 61945 [details] > > > Reproducer, single file, first part (polarization test) > > > > This is what I get > > > > % ./whizard_test > > Running test: polarization_2Test: polarization_2 > > ... success. > > Success. > > *** Test Summary *** > > Success: > > Total = 1 > > Success = 1 > > Failure = 0 > > *** End of test Summary *** > > > > with > > > > GNU Fortran (GCC) 16.0.0 20250712 (experimental) > > Copyright (C) 2025 Free Software Foundation, Inc. > > > > There is no other output. > > Rebuilding with an up-to-date gfortran, I see > > ./whizard_test --check polarizations > Running test: polarization_2 Trace of matrix square = 3.2500000000000000 > > Not permissible as density matrix > Trace of matrix square = 3.3399999999999999 > Not permissible as density matrix > Trace of matrix square = 3.9705327314630523 > Not permissible as density matrix > ... failure. See: err-output/polarization_2.out > *** Test Summary *** > Failure: > Total = 1 > Success = 0 > Failure = 1 > *** End of test Summary *** > hotrats:kargl[369] gfcx --version > GNU Fortran (GCC) 16.0.0 20250723 (experimental) > > ./whizard_test --check polarizations > Running test: polarization_2 Trace of matrix square = 3.2500000000000000 > > Not permissible as density matrix > Trace of matrix square = 3.3399999999999999 > Not permissible as density matrix > Trace of matrix square = 3.9705327314630523 > Not permissible as density matrix > ... failure. See: err-output/polarization_2.out > *** Test Summary *** > Failure: > Total = 1 > Success = 0 > Failure = 1 > *** End of test Summary *** > > % gfcx --version > GNU Fortran (GCC) 16.0.0 20250723 (experimental) Thanks, Steve, for reproducing this. I still need to understand where this difference comes about, but it seems some data structures get not handled correctly.