! { dg-do run }
! { dg-require-effective-target fortran_large_real }
! Test that the internal pack and unpack routines work OK
! for our large real type.

program main
  implicit none
  integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1)
  real(kind=k), dimension(3) :: rk
  complex(kind=k), dimension(3) :: ck

  rk = (/ -1.0_k, 1.0_k, -3.0_k /)
  call sub_rk(rk(1:3:2))
  if (any(rk /= (/ 3.0_k, 1.0_k, 2.0_k/))) STOP 1

  ck = (/ (-1.0_k, 0._k), (1.0_k, 0._k), (-3.0_k, 0._k) /)
  call sub_ck(ck(1:3:2))
  if (any(real(ck) /= (/ 3.0_k, 1.0_k, 2.0_k/))) STOP 2
  if (any(aimag(ck) /= 0._k)) STOP 3

end program main

subroutine sub_rk(r)
  implicit none
  integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1)
  real(kind=k), dimension(2) :: r
  if (r(1) /= -1._k) STOP 4
  if (r(2) /= -3._k) STOP 5
  r(1) = 3._k
  r(2) = 2._k
end subroutine sub_rk

subroutine sub_ck(r)
  implicit none
  integer,parameter :: k = selected_real_kind (precision (0.0_8) + 1)
  complex(kind=k), dimension(2) :: r
  if (r(1) /= (-1._k,0._k)) STOP 6
  if (r(2) /= (-3._k,0._k)) STOP 7
  r(1) = 3._k
  r(2) = 2._k
end subroutine sub_ck