3d array assignment error in fortran segfault

Question

My code is currently breaking when I try to pass an array to another array. I'm using gfortran as the compiler with no special flags. I'm also writing the file as a ".f95". Here's the backbone of the code that's causing the problem:

    program foo
        integer, parameter :: pp = 16

    call fooyoo

    contains 

    subroutine fooyoo

        real(kind=pp), allocatable :: f(:,:,:), y_ix(:), r(:)
        integer :: ii, jj, nn


        nn = 32
        r = zoo(nn)
        y_ix = r    ! this is done because I'm actually allocating two other 1D arrays

    do ii = 1, nn
        do jj = 1, nn
            write(*,*) "crashes here"            
                f(ii,jj,:) = y_ix
            write(*,*) "nothing here"
        end do
    end do

    end subroutine

    function zoo(nn) result(r)
        integer :: i, nn
        real(kind=pp) :: r(nn)

    do i = 1, nn
        r(i) = i
    end do

    end function
    end program

High Performance Mark · Accepted Answer · 2014-05-01 18:12:33Z

3

Probably

do ii = 1, nn
    do jj = 1, nn
        write(*,*) "crashes on executing the next line because programmer forgot ", &
                   "to allocate f before trying to set values"            
            f(ii,jj,:) = y_ix
        write(*,*) "nothing here"
    end do
end do

In this line

    y_ix = r

automatic allocation works because the expression on the lhs is an allocatable array. In

    f(ii,jj,:) = y_ix

automatic allocation fails because the expression on the lhs is an array section, not an allocatable array.

edited May 1, 2014 at 18:12

answered May 1, 2014 at 18:04

High Performance Mark

78.6k7 gold badges109 silver badges168 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

Indelible Over a year ago

I did this: IF (.NOT. ALLOCATED(f)) ALLOCATE(f(nn,nn,nn)) seems to have worked, but seems inelegant.

High Performance Mark Over a year ago

I don't understand why you bothered to test whether f was allocated. From the code you have shown there is no possibility that it could be. But whatever floats your boat.

Collectives™ on Stack Overflow

3d array assignment error in fortran segfault

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related