トリック

Question

私は巨大な銀河カタログ上の特定の計算を行いコンパイル言語（Fortran 95の）を使用してコードを開発しています。いくつかの変更を実装するたびに、コードをコンパイルして実行し、ディスクから銀河のデータをASCIIファイルで読み取るだけで約3分かかります。これは時間の無駄です。

は、配列データを含む変数が異なるコンパイル間のメモリに保持されるので、それは、異なるだろう、私はIDLまたは、MATLABでこのプロジェクトを始めました。

はしかし、私は何かが偽のRAMパーティションか何かでファイルを持っているように、ディスクからの読み取り狼狽することをスピードアップするために行うことができると思います。代わりに、私はあなたがバイナリのものにASCIIデータベースに切り替える提案するRAMディスクの詳細に入るの

Answer 1

。ここで

program writeArr 
    use,intrinsic :: ISO_Fortran_env, only: REAL64 
    implicit none 
    real(REAL64),allocatable :: tmp(:,:) 
    integer :: uFile, i 

    allocate(tmp(10000,10000)) 

    ! Formatted read 
    open(unit=uFile, file='ASCII.txt',form='formatted', & 
     status='replace',action='write') 
    do i=1,size(tmp,1) 
    write(uFile,*) tmp(:,i) 
    enddo !i 
    close(uFile) 

    ! Unformatted read 
    open(unit=uFile, file='binary.bin',form='unformatted', & 
     status='replace',action='write') 
    write(uFile) tmp 
    close(uFile) 

end program 

サイズの面で結果がされています：ここでは非常に単純化した例... ASCII（ASCII.txt）として格納され、ランダムな数字の配列、およびバイナリ日付など（binary.bin）がある

:> ls -lah ASCII.txt binary.bin 
-rw-rw-r--. 1 elias elias 2.5G Feb 20 20:59 ASCII.txt 
-rw-rw-r--. 1 elias elias 763M Feb 20 20:59 binary.bin 

したがって、保存率で約3.35という係数を保存します。 は今楽しい部分が来る：で...

program readArr 
    use,intrinsic :: ISO_Fortran_env, only: REAL64 
    implicit none 
    real(REAL64),allocatable :: tmp(:,:) 
    integer :: uFile, i 
    integer :: count_rate, iTime1, iTime2 

    allocate(tmp(10000,10000)) 

    ! Get the count rate 
    call system_clock(count_rate=count_rate) 

    ! Formatted write 
    open(unit=uFile, file='ASCII.txt',form='formatted', & 
     status='old',action='read') 

    call system_clock(iTime1) 
    do i=1,size(tmp,1) 
    read(uFile,*) tmp(:,i) 
    enddo !i 
    call system_clock(iTime2) 
    close(uFile) 
    print *,'ASCII read ',real(iTime2-iTime1,REAL64)/real(count_rate,REAL64) 

    ! Unformatted write 
    open(unit=uFile, file='binary.bin',form='unformatted', & 
     status='old',action='read') 
    call system_clock(iTime1) 
    read(uFile) tmp 
    call system_clock(iTime2) 
    close(uFile) 
    print *,'Binary read ',real(iTime2-iTime1,REAL64)/real(count_rate,REAL64) 

end program 

戻ってそれを読んだ結果がそう

ASCII read 37.250999999999998  
Binary read 1.5460000000000000 

、> 24倍です！

ので、代わりの何か他のものの考え方、最初のバイナリファイル形式に切り替えてください。