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 75