>
|
uplo := "Nag_Lower":
trans := "Nag_ConjTrans":
diag := "Nag_NonUnitDiag":
n := 5:
incx := -2:
t := Matrix([[1.416583416583417, -10000000000, -10000000000, -10000000000, -10000000000, 1.346653346653347, -0.09290709290709291, 0.05694305694305694, 0.1768231768231768, 0.9370629370629371, 0.06693306693306693, -0.2127872127872128, -0.2927072927072927, -0.1328671328671329, 0.8271728271728271, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000], [-0.4525474525474525, 0.6973026973026972, -10000000000, -10000000000, -10000000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000], [0.3766233766233766, 0.02697302697302697, 1.346653346653347, -10000000000, -10000000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000], [-0.09290709290709291, 0.05694305694305694, 0.1768231768231768, 0.9370629370629371, -10000000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000], [0.06693306693306693, -0.2127872127872128, -0.2927072927072927, -0.1328671328671329, 0.8271728271728271, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000, 11200000000]], datatype=float[8], order='C_order'):
x := Vector([0.5969030969030968, -10000000000, 0, -10000000000, 0.9565434565434565, -10000000000, 0.6868131868131868, -10000000000, 0.2472527472527473], datatype=float[8]):
alpha := 1:
NAG:-f16pjc(uplo, trans, diag, n, alpha, t, x, incx):
|