作者poiuy1203 (po)
看板Fortran
标题Re: [问题] 程式发生错误 找很久了 高手能帮帮忙吗?
时间Mon May 10 14:39:06 2010
※ 引述《poiuy1203 (po)》之铭言：
: ※ 引述《poiuy1203 (po)》之铭言：
: : 程式跑到挂掉
: : 但小弟功力太弱找不到原因
: : 恳请高手帮看指点
: : 画面显示
: : run-time error M6203: MATH
: : -csin:OVERFLOW error
: : 跑到L=2,M=3,N=4时会挂掉
: : 下面附上程式码:
: : SUBROUTINE IILMN(IL,ML1,NN)
: : COMPLEX*16 II(6,5,80),K(6,5,80),Q(6,5,80),AX(6,5),AZ(6,5),DZZ,CI
: : REAL*8 DX(6),Z(6,5)
: : DIMENSION ML1(6),NN(6)
: : COMMON / DAT2 / DX,Z
: : COMMON / DAT4 / AX,AZ
: : COMMON / DAT5 / K,Q
: : COMMON / DAT22 / II
: : CI=CMPLX(0.0,1.0)
: : DO 1,L=1,IL
: : DO 2,N=1,NN(L)
: : II(L,2,N)=1.0
: : 2 CONTINUE
: : 1 CONTINUE
: : DO 3,L=1,IL
: : ML1L=ML1(L)
: : IF(ML1L.EQ.2) GO TO 3
: : DO 4,M=3,ML1L
: : DZZ=(Z(L,M-1)-Z(L,M))*AX(L,M)/AZ(L,M)
: : WRITE(*,*)'ML1L',ML1L
: : DO 5,N=1,NN(L)
: : DZZ=DZZ*K(L,M,N)
: : II(L,M,N)=II(L,M-1,N)*(COS(DZZ)+CI*Q(L,M,N)*SIN(DZZ))
: : WRITE(*,*)L,M,N
: : 5 CONTINUE
: : 4 CONTINUE
: : 3 CONTINUE
: : RETURN
: 目前所知道的是DZZ的值为复数形式
: DZZ=(0.625 , 1.734918577506039E-018)
: 带入SIN(DZZ)和COS(DZZ)时会爆掉
将SIN和COS改成DSIN和DCOS後产生以下问题:
error LNK2001: unresolved external symbol _DCOS@4
error LNK2001: unresolved external symbol _DSIN@4
Debug/teach.exe : fatal error LNK1120: 2 unresolved externals
Error executing link.exe.
因为我其他地方也有一样的问题
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.117.98.63