------- Comment #6 from dominiq at lps dot ens dot fr 2007-02-27 15:36 ------- I don't know if this reduced test account for all the problem, but it exhibits at least one bug:
module perdida_m implicit none contains subroutine perdida (dt, strain_tensor) real (kind = 8), intent(in) :: dt real (kind = 8), dimension(:,:), intent(in) :: strain_tensor end subroutine perdida end module perdida_m program iztaccihuatl use perdida_m implicit none real (kind = 8) :: dt, time, coefficient real (kind = 8), dimension(3,3,10) :: strain_tensor dt = 5.555555555555556D-006 strain_tensor = 0.0_8 time = dt coefficient = sin(time) print *, coefficient coefficient = cos(time) call perdida (dt, strain_tensor(:,:,1)) end program iztaccihuatl gives 1.947453947074031E-314 If I comment the line 'call perdida (dt, strain_tensor(:,:,1))', I get 5.555555555526978E-006 As far as I can tell, all the ingredients in this test case are necessary to trigger the bug. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=30969