view gcc/testsuite/gfortran.fortran-torture/compile/pr36078.f90 @ 158:494b0b89df80 default tip

...
author Shinji KONO <kono@ie.u-ryukyu.ac.jp>
date Mon, 25 May 2020 18:13:55 +0900
parents 04ced10e8804
children
line wrap: on
line source

   subroutine foo(func,p,eval)
      real(kind=kind(1.0d0)), dimension(3,0:4,0:4,0:4) :: p
      logical(kind=kind(.true.)), dimension(5,5,5) :: eval
      interface
         subroutine func(values,pt)
            real(kind=kind(1.0d0)), dimension(:), intent(out) :: values
            real(kind=kind(1.0d0)), dimension(:,:), intent(in) :: pt
         end subroutine
      end interface
      real(kind=kind(1.0d0)), dimension(125,3) :: pt
      integer(kind=kind(1)) :: n_pt

      n_pt = 1
      pt(1:n_pt,:) = &
         reshape( &
            pack( &
               transpose(reshape(p,(/3,125/))), &
               spread(reshape(eval,(/125/)),dim=2,ncopies=3)), &
            (/n_pt,3/))

   end subroutine
   end