数组越界错误与数组输出太慢

sqs · 发表于 6 天前

主要问题就是1. ES赋值后进行输出就死暗示数组越界了，注释掉那一句话后就没有这个问题了，2.但是如果不对ES赋值，对对角化得到的本征态进行输出，对整个二维数组输出的时候非常慢。请问问题出在哪里呀

[Fortran] 纯文本查看 复制代码

    !ES = ES_H(:, N:N+1)
    !ES = reshape([ES_H(:, N), ES_H(:, N+1)], [2*N, 2])

    !call ES_TRIPOD(theta, Phi0)

    ! the output file
    open(1, file = 'eb.dat' , status = 'replace')

    !do i = 1, 2*N
        write(1,*)EV_H(1) 
        ! write(*,*)i
    !enddo

    close(1)

[求助] 数组越界错误与数组输出太慢