------- Comment #10 from jv244 at cam dot ac dot uk  2007-07-02 10:21 -------

> Joost, could you please test this 4.2 patch with CP2K? It fixes the testcase
> you provided (and PR31966, FWIW).
> 

the testcase passes, the CP2K tests that were previously failing now seem to
pass as well .... many thanks for fixing this issue

> BTW: Could you create the testcase that would be suitable for gfortran
> testsuite?

this aborts before and passes after the patch:

      SUBROUTINE T(nsubcell,sab_max,subcells)
       INTEGER, PARAMETER :: dp=KIND(0.0D0)
       REAL(dp) :: sab_max(3), subcells,nsubcell(3)
       nsubcell(:) = MIN(MAX(1,NINT(0.5_dp*subcells/sab_max(:))),20)
      END SUBROUTINE

      INTEGER, PARAMETER :: dp=KIND(0.0D0)
      REAL(dp) :: sab_max(3), subcells,nsubcell(3)
      subcells=2.0_dp
      sab_max=0.590060749244805_dp
      CALL T(nsubcell,sab_max,subcells)
      IF (ANY(nsubcell.NE.2.0_dp)) CALL ABORT()
     END


-- 


http://gcc.gnu.org/bugzilla/show_bug.cgi?id=32533

Reply via email to