<p>我希望有一个真正的Fortran程序员出现,但是在没有更好的建议的情况下,我只指定<code>x(:)</code>的形状而不是大小,使用临时数组<code>temp(size(x))</code>,并使输出为y<code>allocatable</code>。然后在第一次传递之后,<code>allocate(y(j))</code>并从临时数组中复制值。但是我不能强调我不是一个Fortran程序员,所以我不能说这个语言是否有一个可增长的数组,或者是否有一个用于后者的库。</p>
<pre><code>program test
implicit none
integer:: x(10) = (/1,0,2,0,3,0,4,0,5,0/)
print "(10I2.1)", select(x)
contains
function select(x) result(y)
implicit none
integer, intent(in):: x(:)
integer:: i, j, temp(size(x))
integer, allocatable:: y(:)
j = 0
do i = 1, size(x)
if (x(i) /= 0) then
j = j + 1
temp(j) = x(i)
endif
enddo
allocate(y(j))
y = temp(:j)
end function select
end program test
</code></pre>
<hr/>
<p>编辑:</p>
<p>根据M.S.B.的回答,这里有一个函数的修订版本,它在过度分配的情况下增长了<strike>temp</strike><code>y</code>。<strike>就像它在末尾将结果复制到y一样。</strike>事实证明,我不需要显式地以最终大小分配新数组。相反,它可以通过赋值自动完成。</p>
<pre><code> function select(x) result(y)
implicit none
integer, intent(in):: x(:)
integer:: i, j, dsize
integer, allocatable:: temp(:), y(:)
dsize = 0; allocate(y(0))
j = 0
do i = 1, size(x)
if (x(i) /= 0) then
j = j + 1
if (j >= dsize) then !grow y using temp
dsize = j + j / 8 + 8
allocate(temp(dsize))
temp(:size(y)) = y
call move_alloc(temp, y) !temp gets deallocated
endif
y(j) = x(i)
endif
enddo
y = y(:j)
end function select
</code></pre>