| 查看: 365 | 回复: 0 | ||
[求助]
放在循环里的数组越界了，求帮忙
|
|
数组放在左边没有越界，放在等式右边就出现了越界。上传的图片中前面用来测试的y(k)好好的，循环里就有问题了，求大神帮忙 subroutine wfft(y,forback) use fft_mod implicit none real*8,parameter :: pi=3.141592654_DP complex*16,intent(inout) :: y(  integer,intent(in) :: forback real*8 :: e integer :: N2,N,k,ktem,kN complex*16 W,temp integer,external :: reverse N=size(y) N2=N kN=N/2 k=0 e=forback*2.0_DP*pi/N W=exp(cmplx(0.0_DP,e)) do k=1,kN,1 y(k)=cmplx(real(y(k)),0.0_DP) write(*,*) k,real(Y(k)) enddo do N2=N2/2 if(N2<1)exit do k=1,kN write(*,*) k,kN,k+kN,y(k) temp=W**(reverse(k))*y(k+KN) y(k+kN)=y(k)-temp y(k)=y(k)+temp enddo enddo endsubroutine  _7P4OEAHC(UFOUR8$Z%~R`4.png  2L1Y{@@TDQ%Z056K0N%[IOC.png |
回复此楼» 猜你喜欢
材料考研调剂
已经有3人回复
-
材料调剂
已经有12人回复
-
英一数一408，总分284，二战真诚求调剂
已经有14人回复
-
085410 一志愿211 22408分数359求调剂
已经有4人回复
-
271求调剂
已经有19人回复
-
385分 生物学（071000）求调剂
已经有3人回复
-
一志愿安徽大学计算机科学与技术学硕，331分求调剂
已经有3人回复
-
318求调剂，计算材料方向
已经有8人回复
-
291求调剂
已经有25人回复
-
一志愿北京科技大学085601材料工程英一数二初试总分335求调剂
已经有6人回复
