- This topic has 4개 답변, 3명 참여, and was last updated 9 years, 11 months 전에 by
.
-
글
-
-
% Program caused arithmetic error: Floating illegal operand 등 에러가 나며
가령 M2=0.408557
M4=0.0367701
M6=0.00330931 일때eta=M4^2/(M2*M6) 을 계산하면 0.99999888…… 이지만 print를 하면 1로 인식을 합니다.
그래서 이 뒤에 mu=2*(eta-1)
aa=eta/mu
라고 계산하게되면 aa는 -inf 로 뜨게 됩니다.mu에서 eta-1이 0이 되서 에러가 나는거 같은데
혹시 eta를 1이 아닌 0.999처럼 원래 값이 찍히게 할수는 없는지 하여 글 올립니다. -
-
idl 8.3버전을 사용하고 있으며, 위에서 쉽게 설명하기위해 M2, M4,M6값을 지정해서 계산하는 방식으로 설명하였지만 실제로는 아래와 같이 M2, M4, M6가 계산됩니다.
delta_d=0.2 ; [mm]
M_2=fltarr(orb, ncount)
M2=fltarr(orb)
M_4=fltarr(orb, ncount)
M4=fltarr(orb)
M_6=fltarr(orb, ncount)
M6=fltarr(orb)for aa=0, orb-1 do begin
for i=0, ncount-1 do begin
M_2[aa,i]=(d[aa,i]^2)*nd[aa,i]*delta_d
M2[aa]=M2[aa]+M_2[aa,i]M_4[aa,i]=(d[aa,i]^4)*nd[aa,i]*delta_d
M4[aa]=M4[aa]+M_4[aa,i]M_6[aa,i]=(d[aa,i]^6)*nd[aa,i]*delta_d
M6[aa]=M6[aa]+M_6[aa,i]endfor ; for i
endfor ; for aa위에서 계산된 M2,M4,M6 를 가지고 아래 eta, mu 계산에서 문제가 발생하여 알려주신대로 double을 썼을때와
쓰지 않았을때의 값을 print 하였을때;=== M246 moments ===
; eta=M4^2/(M2*M6)
; mu=((7-11*eta)-(eta^2+14*eta+1)^0.5)/(2*(eta-1))
; lamda=[(M2/M4)*(mu+3)*(mu+4)]^0.5
; No=(M2*lamda^(mu+3))/gamma(mu+3)eta=(M4^2)/(M2*M6) => 기존포맷
eta=double(M4)^2/(double(M2)*double(M6)) => 알려주신 포맷두 각각의 포맷을 print 하였을때 둘다 1값을 프린트하며
mu=((7-11*eta)-(eta^2+14*eta+1)^0.5)/(2*(eta-1))
mu 값에서는 두 값이 서로 미세하지만 차이를 보였습니다.
두 값모두 이론적으로 계산하였을때의 값과 차이를 보여 1)알려주신 방법을 사용하는 것이 맞는것인지 아님
또다른 방법이 있는지 해서 다시 질문드립니다.
2) 또한 위에서 구해진 mu를 통해 lamda, No가 계산되어지는데 아래 부분에서도 double을 써서 계산해야하는것이
맞는지 궁금하니다. 알려주신 double을 쓰면 ‘program caused….’와 같은 에러명령어는 사라지고
잘 수행되는데 값에 대한 신뢰성에 대해 판단이 잘 들지 않아서요 ㅠㅠlamda=((M2/M4)*(mu+3)*(mu+4))^0.5
No=(M2*lamda^(mu+3))/gamma(mu+3) -
-
-
