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.

Reply via email to