gcc/testsuite/gfortran.dg/reassoc_12.f90

   1 ! { dg-do compile }
   2 ! { dg-options "-O2 -ffast-math" }
   3 ! PR middle-end/57370
   4
   5  SUBROUTINE xb88_lr_adiabatic_lda_calc(e_ndrho_ndrho_ndrho, &
   6                                        grad_deriv,npoints, sx)
   7     IMPLICIT REAL*8 (t)
   8     INTEGER, PARAMETER :: dp=8
   9     REAL(kind=dp), DIMENSION(1:npoints) :: e_ndrho_ndrho_ndrho, &
  10                                            e_ndrho_ndrho_rho
  11       DO ii=1,npoints
  12           IF( grad_deriv >= 2 .OR. grad_deriv == -2 ) THEN
  13             t1425 = t233 * t557
  14             t1429 = beta * t225
  15             t1622 = t327 * t1621
  16             t1626 = t327 * t1625
  17             t1632 = t327 * t1631
  18             t1685 = t105 * t1684
  19             t2057 = t1636 + t8 * (t2635 + t3288)
  20           END IF
  21           IF( grad_deriv >= 3 .OR. grad_deriv == -3 ) THEN
  22             t5469 = t5440 - t5443 - t5446 - t5449 - &
  23                     t5451 - t5454 - t5456 + t5459  - &
  24                     t5462 + t5466 - t5468
  25             t5478 = 0.240e2_dp * t1616 * t973 * t645 * t1425
  26             t5489 = 0.1600000000e2_dp * t1429 * t1658
  27             t5531 = 0.160e2_dp * t112 * t1626
  28             t5533 = 0.160e2_dp * t112 * t1632
  29             t5537 = 0.160e2_dp * t112 * t1622
  30             t5541 = t5472 - t5478 - t5523 + t5525 + &
  31                     t5531 + t5533 + t5535 + t5537 + &
  32                     t5540
  33             t5565 = t112 * t1685
  34             t5575 = t5545 - t5548 + t5551 + t5553 - &
  35                     t5558 + t5560 - t5562 + t5564 - &
  36                     0.80e1_dp * t5565 + t5568 + t5572 + &
  37                     t5574
  38             t5611 = t5579 - t5585 + t5590 - t5595 + &
  39                     t5597 - t5602 + t5604 + t5607 + &
  40                     t5610
  41             t5613 = t5469 + t5541 + t5575 + t5611
  42             t6223 = t6189 - &
  43                     0.3333333336e0_dp  * t83 * t84 * t5613 + &
  44                     t6222
  45             t6227 = - t8 * (t5305 + t6223)
  46             e_ndrho_ndrho_rho(ii) = e_ndrho_ndrho_rho(ii) + &
  47                      t6227 * sx
  48             t6352 = t5440 - t5443 - t5446 - t5449 - &
  49                     t5451 - t5454 + &
  50                     0.40e1_dp * t102  * t327 * t2057 * t557 - &
  51                     t5456 + t5459 - t5462 + t5466 - &
  52                     t5468
  53             t6363 = t5480 - t5489 + &
  54                     0.9600000000e2_dp  * t1054 * t640 * t3679
  55             t6367 = t5472 - t5474 - t5478 - t5523 + &
  56                     t5525 + t5531 + t5533 + t5535 + &
  57                     t5537 - 0.20e1_dp * t102 * t105 * t6363 + &
  58                     t5540
  59             t6370 = t5545 - t5548 + t5551 + t5553 - &
  60                     t5558 + t5560 - t5562 + t5564  - &
  61                     0.40e1_dp * t5565 + &
  62                     t5568 + t5572 + t5574
  63             t6373 = t5579 - t5585 + t5590 - t5595 + &
  64                     t5597 - t5602 + t5604 + t5607 + &
  65                     t5610
  66             t6375 = t6352 + t6367 + t6370 + t6373
  67             t6380 = - 0.3333333336e0_dp * t83 * t84 * t6375 + t5701
  68             t6669 = -t4704 - t8 * (t6344 + t6380 + t6665)
  69             e_ndrho_ndrho_ndrho(ii) = e_ndrho_ndrho_ndrho(ii) + &
  70                      t6669 * sx
  71           END IF
  72       END DO
  73   END SUBROUTINE xb88_lr_adiabatic_lda_calc
  74