내가 효율적으로 배열을 전치하고

Question

다음과 같은 고려 다른 유형의 배열에 할당 할 수있는 방법 :

program main 
integer, parameter :: n=10, m=20 
integer ints(n,m) 
real floats(m,n) 

!... initialize ints 
! ... 

floats=transpose(ints) 
!... do stuff with floats 

end 

이 gfortran에 대한 문서를보고, transpose(ints) 다음 캐스트됩니다 정수 배열을 반환 것으로 보인다 재판에. 이 작업에서 컴파일러 (gfortran)는 transposed 배열에 대해 임시 배열을 만듭니다.이 배열은 낭비처럼 보입니다 ( gfortran -O3 -Warray-temporaries -o test test.f90으로 컴파일). 또한 배열 "floats"를 integer 배열로 변경하면 경고 메시지가 사라집니다.

임시 배열을 생성하지 않고 임의의 유형에 대해이를 수행 할 수있는 방법이 있습니까? (나는 어딘가에 그것이 중요하다고 읽었 기 때문에 floats(:,:)=transpose(ints)도 시도했다.). 다른 컴파일러와 함께 이런 식으로 동작합니까?

Answer 1

당신은

floats = transpose(real(ints)) 

을 시도 할 수 있지만, gfortran (또는 다른 컴파일러)이을 구현하기 위해 임시 배열을 생성하면 정말 놀라지 않을 것입니다. 그렇지 않다면 더 놀랄 것입니다.

또한 컴파일러가이를 구현하기 위해 임시 배열을 만든 경우 나는 놀라지 않을 것이다, 다시

forall (J=1:N, K=1:M) floats(K, J) = real(ints(J, K)) 

을 시도 할 수 있습니다. 이 서브 루틴이 아닌 함수가 될해야하지만

Answer 2

do i = 1, n 
do j = 1, m 
    floats(j,i) = real(ints(i,j)) 
enddo 
enddo 

당신은 다른 데이터 유형을 처리하기위한 자신의 transpose 인터페이스를 만들 수 있습니다.

interface transpose_ 
    module procedure transpose_ints_to_reals 
end interface 

subroutine transpose_ints_to_reals(ints_in, reals_out) 
... 
end subroutine 

call transpose_(ints,floats)