SAD部屋 - Ticket-25 Diff

! SOL+MULT errata under BZ closed to 0

:Priority:Normal
:Reporter:Akio Morita
:Status:New
:Assigned to:?
:Version:?
:Milestone:?
:Created:2011-04-08

!! Description

BZ=0近傍での SOL+MULTの動作がおかしいとの報告有り

BZ多項式からなる 無限小間の除算があると推測される。

現時点で判明しているのは
*src/tmulti.fの a12, a14の定義

!! Attachement files
{{attach_map('Ticket-25')}}
{{attach_form('','')}}

!! Changelog
!!!Akio Morita (2011-04-08 (金) 11:58:58)
*src/tmulti.f a14の定義
a22 >= 0にて
a14 = a24^2 / bzp / (1 + a22)
a12 = a24 / bzp
なので、
a14 = a12 * a24 / (1 + a22)
と置換して a14の特異性を a12に押し込めた

a12 = sin(-bzp * const.) / bzp
に関しては修正が必要
!!!Akio Morita (2011-04-11 (月) 14:54:02)
*src/tquase.f(r3581)
babs=sqrt(bxs**2+bys**2+bzs**2)が、前提条件bxs != 0かつ byz != 0を
満足する場合でも underflowする場合がある。
babs=hypot(bzs,hypot(bxs,bys))へ置換

*src/tquase.f(r3582)
alb=1/babsが babs〜0で overflowする
babsに相関して 0近傍の値となる項と albの積が NaNになる

前提条件から babs > 0が保証されるので、foo * albを foo / babsに置換
!!!Akio Morita (2011-04-12 (火) 01:30:40)
*src/tsolque.f(r3584, r3585)
bzp^2 underflowに伴う NaN(bzp^2/bzp^2)の修正
cod(5)及び trans1(5,6)
!!!Akio Morita (2011-04-12 (火) 01:32:33)
*src/tsolqu.f
tsolqueのcod(5)と同様のbzp^2 underflowが存在するはず。要解析

{{its_edit_ticket_form}}