24531904_j

fma ─ 浮動小数点積和 (Floating-point Multiply Add)
書式 : (qp) fma.pc.sf f 1 = f 3 , f 4 , f 2 F1
説明 : FR f3 と FR f4 と の積が無限の精度で計算され、 次に FR f2 がこの積にやはり無限の精度で
加算される。 次に、 結果の値が、 FPSR.sf.rc で指定される丸めモー ド と pc ( さらに
FPSR.sf.pc および FPSR.sf.wre) で指定される精度にしたが って丸められ、 その結果が FR
f1 に格納される。
FR f3 、 FR f4 、 FR f2 のどれか 1 つでも NaTVal である場合は、 FR f1 は計算結果ではな く
NaTVal に設定される。
f2 が f0 の場合は、 積和演算の代わ り に IEEE の乗算が行われる。 3:77 ページの 「fmpy —
浮動小数点乗算 (Floating-point Multiply)」 を参照のこ と。
オペコ ー ド の pc のニーモニ ッ ク値は 3:50 ページの表 2-22 に示している。 sf のニーモ
ニック値は3:50 ページの表 2-23 に示している。 ステータス ・ フ ィ ール ド の pc、 wre、 rc
のエンコ ーデ ィ ング と意味については、 第 1 巻、 第 1 部の表 5-5 および第 1 巻、 第 1 部
の表 5-6 を参照のこ と。
操作 : if (PR[qp]) {
fp_check_target_register(f 1);
if (tmp_isrcode = fp_reg_disabled(f 1 , f 2 , f 3 , f 4 ))
disabled_fp_register_fault(tmp_isrcode, 0);
}
if (fp_is_natval(FR[f 2 ]) || fp_is_natval(FR[f 3 ]) ||
fp_is_natval(FR[f 4 ])) {
FR[f 1] = NATVAL;
fp_update_psr(f 1 );
} else {
tmp_default_result = fma_exception_fault_check(f 2, f 3, f 4,
pc, sf, &tmp_fp_env);
if (fp_raise_fault(tmp_fp_env))
fp_exception_fault(fp_decode_fault(tmp_fp_env));
}
if (fp_is_nan_or_inf(tmp_default_result)) {
FR[f 1] = tmp_default_result;
} else {
tmp_res = fp_mul(fp_reg_read(FR[f 3 ]), fp_reg_read(FR[f 4 ]));
if (f 2 != 0)
tmp_res = fp_add(tmp_res, fp_reg_read(FR[f 2 ]), tmp_fp_env);
FR[f 1 ] = fp_ieee_round(tmp_res, &tmp_fp_env);
}
fp_update_fpsr(sf, tmp_fp_env);
fp_update_psr(f 1);
if (fp_raise_traps(tmp_fp_env))
fp_exception_trap(fp_decode_trap(tmp_fp_env));
FP 例外 : 無効操作 (Invalid Operation: V) アンダーフロー (Underflow: U)
デノーマル/ アンノ ーマル ・ オペラン ド オーバーフロー (Overflow: O)
(Denormal/Unnormal Operand: D) 不正確結果 (Inexact: I)
ソフ ト ウェア・アシスト ソフ ト ウェア・アシスト
(Software Assist: SWA) フォルト (Software Assist: SWA) ト ラップ
割り込み : 無効操作フ ォ ル ト 浮動小数点例外フ ォ ル ト
浮動小数点レ ジ ス タ無効フ ォ ル ト 浮動小数点例外 ト ラ ッ プ
3:70 第 3 巻 : 命令リファレンス
fma