◆ _zhpr2()

Rank 2 操作: A <- αxy^H + conjg(α)yx^H + A (エルミート行列) (圧縮形式) (BLAS 2)

目的: 本ルーチンはエルミート rank 2 操作を行う.
A <- αxy^H + conjg(α)yx^H + A

ここで, αはスカラー, xおよびyはnベクトル, また, Aは圧縮形式のn×nエルミート行列である.

引数

[in]	uplo	行列Aの上三角部分あるいは下三角部分のどちらを配列ap[]に格納するかを指定. = 'U': Aの上三角部分をap[]に格納. = 'L': Aの下三角部分をap[]に格納.
[in]	n	行列Aの行および列数. (n >= 0) (n = 0 の場合, 処理を行わずに戻る)
[in]	alpha	スカラーα.
[in]	x[]	配列 x[lx] (lx >= 1 + (n - 1)*abs(incx)) ベクトル x.
[in]	incx	配列内におけるxの要素間隔. (incx != 0)
[in]	y[]	配列 y[ly] (ly >= 1 + (n - 1)*abs(incy)) ベクトル y.
[in]	incy	配列内におけるyの要素間隔. (incy != 0)
[in,out]	ap[]	配列 ap[lap] (lap >= n(n + 1)/2) [in] 圧縮形式のn×nエルミート行列 A. uploに従い上三角部分あるいは下三角部分を格納する. 対角要素の虚数部は設定不要で, 0と仮定される. [out] 圧縮形式の出力行列 (= αxy^H + conjg(α)yx^H + A). uploに従い上三角部分あるいは下三角部分が格納される. 対角要素の虚数部は0に設定される.