|
|
|
|
|
x*******3 发帖数: 8342 | 6
我想说的是小白兔....
为啥说成了5d2
找豆腐去。。。 |
|
|
|
t******d 发帖数: 3 | 9 ~{GW0.5D0"5$#:~}
~{!!!!GkTJPmNRUbQy3F:tDc#,NR@4WTL+F=Qs1K06R;8vC@@v5D9z6H!*!*!*VP9z!#NR5DC{WV=P
EK=p2(#,>M=PNR2(2(:CAK!#~}
~{!!!!<8LlG0NR4SIn[ZD31(IO5CV*AKDc5DTbSv#,DG8v:]PD5DE.HKR*EWFzDc!#?4WE1(V=#,NR
0Y~}
~{8P=;
=g9c4s@:Gr0.:CU_5DQ[Gr#,DcRcH;>vH;5XIaFzAKWT<:5DP!+Iq0!#!
~}
~{[email protected];V9K}R;8v#,K}UbVV>uNr5MOB#,WTK=WT@{5DHK2;V55CDcH%0.#,@k>M@k#,
~{!!!!NR#,2(2(#,>MJGUbQyR;8v:CE.HK!#=qDj~}38~{Kj5DNR? |
|
t******d 发帖数: 3 | 10 ~{GW0.5D0"5$#:~}
~{!!!!GkTJPmNRUbQy3F:tDc#,NR@4WTL+F=Qs1K06R;8vC@@v5D9z6H!*!*!*VP9z!#NR5DC{WV=P
EK=p2(#,>M=PNR2(2(:CAK!#~}
~{!!!!<8LlG0NR4SIn[ZD31(IO5CV*AKDc5DTbSv#,DG8v:]PD5DE.HKR*EWFzDc!#?4WE1(V=#,NR
0Y~}
~{8P=;
=g9c4s@:Gr0.:CU_5DQ[Gr#,DcRcH;>vH;5XIaFzAKWT<:5DP!+Iq0!#!
~}
~{[email protected];V9K}R;8v#,K}UbVV>uNr5MOB#,WTK=WT@{5DHK2;V55CDcH%0.#,@k>M@k#,
~{!!!!NR#,2(2(#,>MJGUbQyR;8v:CE.HK!#=qDj~}38~{Kj5DNR? |
|
|
l****o 发帖数: 121 | 12 ~{E\!#~}
~{PT1pNJLb#,E.5D2;R*!#K}SVE\5@MrJ/H%URbV@o5DAl5<#,R*Gs51Da9C!#~}
~{bV3$=2#,DaSP5%N;5D=iI\PEC4#?NRCGR*8y>]~}
~{=iI\PE?MDQ0lAK#,~}
~{03CG5381Ji
~{0!#??O6(;9R*UY?*H+O54s;a!#0U0U0U#,Ev6(>MHkJ@P^PP@-!#~}
~{SZJGK}>M3#3#TZE
~{:s@4K}P^3IR;P)9{#,>-3#V85
~{K}5DVXR*SoB<:\6`#,HCEYR;5c!#~}
1~{!#A50.JGN*AK7VJV~}
2~{!#3T795DU}H7 |
|
c*******h 发帖数: 1096 | 13 DOUBLE PRECISION AA
TOL = 1.0D-18
AA = 0.5D0/TOL
AA 会是多少?我在gdb里面print出来是0 |
|
b*******h 发帖数: 2585 | 14 0.5D0=0
AA=0/somthing=0.0000000000000000000000000000000000000000000 |
|
j****x 发帖数: 943 | 15 0.5D0=0 ? why? It's perfectly legal. |
|
r****y 发帖数: 1437 | 16
此程序four1.for用于追赶法加FFT解椭圆方程会出问题, 需要使用
realft.for如下
SUBROUTINE realft(data,n,isign)
INTEGER isign,n
REAL data(n)
CU USES four1
INTEGER i,i1,i2,i3,i4,n2p3
REAL c1,c2,h1i,h1r,h2i,h2r,wis,wrs
DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
theta=3.141592653589793d0/dble(n/2)
c1=0.5
if (isign.eq.1) then
c2=-0.5
call four1(data,n/2,+1)
else
c2=0.5
theta=-theta
endif
wpr=-2.0d0*sin(0.5d0*theta) |
|