Attachment 357501 Details for Bug 517578 – demo program with required interface

demo program with required interface

array_display_repaired.f90 (text/plain), 2.02 KB, created by Jeff Bastian on 2009-08-14 19:08:18 UTC

(hide)

Description:

Filename:

MIME Type:

Creator: Jeff Bastian

Created: 2009-08-14 19:08:18 UTC

Size: 2.02 KB

patch

obsolete

>program array_display
>  implicit none
>  integer, parameter:: dp=kind(0.d0)
>  real(dp), pointer       :: p(:)
>  real(dp)                :: x(1000)
>  interface
>    subroutine subr1(y,p)
>      integer, parameter:: dp=kind(0.d0)
>      real(dp), intent(in)    :: y(:)
>      real(dp), intent(in)    :: p(:)
>    end subroutine subr1
>  end interface
>  allocate(p(1000))
>  p(:) = 0._dp
>  x(:) = 0._dp
>  call subr1(x,p)
>  stop
>end program array_display
>subroutine subr1(y,p)
>  integer, parameter:: dp=kind(0.d0)
>  real(dp), intent(in)    :: y(:)
>  real(dp), intent(in)    :: p(:)
>  real(dp)                :: r,s
>  real(dp)                :: t(3),u(3)
>  real(dp)                :: v
>  real(dp)                :: w
>  integer                 :: i
>  real(dp), allocatable   :: ar1(:)
>  real(dp), allocatable   :: ar2(:,:)
>  integer                 :: j
>  interface
>    subroutine subr2(r,t,s,u,v,w,ar1,ar2,j,y,p)
>      integer, parameter:: dp=kind(0.d0)
>      real(dp), intent(in)    :: r
>      real(dp), intent(in)    :: t(3)
>      real(dp), intent(in)    :: s
>      real(dp), intent(out)   :: u(3)
>      real(dp), intent(in)    :: v
>      real(dp), intent(in)    :: w
>      real(dp), intent(inout) :: ar1(:)
>      real(dp), intent(inout) :: ar2(:,:)
>      integer,  intent(inout) :: j
>      real(dp), intent(in)    :: y(:)
>      real(dp), intent(in)    :: p(:)
>    end subroutine subr2
>  end interface
>  i = 4*size(p)
>  allocate(ar1(i),ar2(3,i))
>  call subr2(r,t,s,u,v,w,ar1,ar2,j,y,p)
>  return
>end subroutine subr1
>subroutine subr2(r,t,s,u,v,w,ar1,ar2,j,y,p)
>  integer, parameter:: dp=kind(0.d0)
>  real(dp), intent(in)    :: r
>  real(dp), intent(in)    :: t(3)
>  real(dp), intent(in)    :: s
>  real(dp), intent(out)   :: u(3)
>  real(dp), intent(in)    :: v
>  real(dp), intent(in)    :: w
>  real(dp), intent(inout) :: ar1(:)
>  real(dp), intent(inout) :: ar2(:,:)
>  integer,  intent(inout) :: j
>  real(dp), intent(in)    :: y(:)
>  real(dp), intent(in)    :: p(:)
>  u(:) = 0._dp
>!Stop on the next line and look at ar1, ar2, y, and the contents of p.
>  return 
>end subroutine subr2

Actions: View

Attachments on bug 517578: 357499 | 357501