>
|
side := "Nag_RightSide":
uplo := "Nag_Lower":
transt := "Nag_ConjTrans":
diag := "Nag_NonUnitDiag":
m := 9:
n := 9:
alpha := 0.7 -0.9*I:
t := Matrix([[1.146853146853147 +0.04295704295704296*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.02297702297702298 -0.3486513486513487*I , 0.9470529470529471 -0.2927072927072927*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.4925074925074925 -0.4925074925074925*I , 0.3066933066933067 +0.2427572427572428*I , 0.7472527472527473 +0.04295704295704296*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.2627372627372627 -0.3006993006993007*I , -0.4625374625374625 +0.05094905094905095*I , -0.4225774225774226 -0.3486513486513487*I , 0.987012987012987 -0.1568431568431568*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.3326673326673327 +0.4425574425574426*I , 0.4665334665334666 +0.3066933066933067*I , 0.1068931068931069 -0.4925074925074925*I , 0 +0*I , 1.026973026973027 -0.3566433566433567*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.3726273726273726 +0.1068931068931069*I , 0.0969030969030969 +0.2987012987012987*I , 0.2667332667332667 +0.4425574425574426*I , 0 +0*I , 0.05694305694305694 -0.1488511488511489*I , 0.7072927072927073 -0.4925074925074925*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.3426573426573427 -0.1008991008991009*I , 0.4265734265734266 -0.3566433566433567*I , -0.1028971028971029 +0.4505494505494506*I , 0 +0*I , -0.2127872127872128 -0.09290709290709291*I , -0.06293706293706294 -0.3006993006993007*I , 0.8271728271728271 +0.1068931068931069*I , -10000000000 +10000000000*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [0.3866133866133866 -0.1568431568431568*I , 0.1868131868131868 -0.09290709290709291*I , 0.2267732267732268 +0.1068931068931069*I , 0 +0*I , 0.2167832167832168 +0.4985014985014985*I , -0.1328671328671329 +0.4425574425574426*I , -0.4125874125874126 -0.1568431568431568*I , 1.146853146853147 -0.2927072927072927*I , -10000000000 +10000000000*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.1828171828171828 +0.2507492507492508*I , -0.3826173826173826 +0.4985014985014985*I , 0.2567432567432568 -0.1008991008991009*I , 0 +0*I , 0.3466533466533466 +0.04295704295704296*I , 0.4965034965034965 +0.4505494505494506*I , 0.01698301698301698 +0.2507492507492508*I , -0.02297702297702298 +0.0989010989010989*I , 0.5074925074925075 +0.2427572427572428*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ]], datatype=complex[8], order='C_order'):
b := Matrix([[-0.3326673326673327 +0.3066933066933067*I , 0.4665334665334666 +0.4425574425574426*I , -0.1028971028971029 +0.2987012987012987*I , 0 +0*I , -0.1428571428571428 -0.1488511488511489*I , 0.3866133866133866 -0.1568431568431568*I , 0.1868131868131868 -0.09290709290709291*I , 0.4165834165834166 +0.2507492507492508*I , 0.3466533466533466 +0.04295704295704296*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [0.2967032967032967 +0.2987012987012987*I , 0.4265734265734266 +0.1068931068931069*I , 0.2267732267732268 -0.3566433566433567*I , 0 +0*I , -0.4125874125874126 -0.09290709290709291*I , -0.1828171828171828 +0.2507492507492508*I , -0.2527472527472527 +0.04295704295704296*I , -0.4525474525474525 -0.2927072927072927*I , 0.1768231768231768 -0.3486513486513487*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.3726273726273726 -0.3566433566433567*I , 0.4565434565434565 -0.1008991008991009*I , -0.01298701298701299 -0.09290709290709291*I , 0 +0*I , 0.01698301698301698 +0.4985014985014985*I , -0.05294705294705294 -0.2927072927072927*I , -0.4225774225774226 -0.3486513486513487*I , -0.09290709290709291 +0.2427572427572428*I , -0.2927072927072927 -0.4925074925074925*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.3426573426573427 -0.1488511488511489*I , 0.1868131868131868 -0.1568431568431568*I , 0.4165834165834166 +0.4985014985014985*I , 0.3466533466533466 -0.2927072927072927*I , 0.1468531468531468 +0.04295704295704296*I , -0.2227772227772228 +0.0989010989010989*I , 0.1068931068931069 -0.4925074925074925*I , 0.1368631368631369 +0.05094905094905095*I , -0.1328671328671329 +0.4425574425574426*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [0.3866133866133866 -0.09290709290709291*I , -0.3826173826173826 +0.2507492507492508*I , -0.4525474525474525 +0.04295704295704296*I , 0 +0*I , -0.4925074925074925 -0.4925074925074925*I , 0.3066933066933067 +0.2427572427572428*I , 0.3366633366633366 -0.3006993006993007*I , 0.06693306693306693 +0.3066933066933067*I , 0.4965034965034965 +0.4505494505494506*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.05294705294705294 +0.04295704295704296*I , -0.2527472527472527 -0.2927072927072927*I , 0.3766233766233766 -0.3486513486513487*I , 0 +0*I , -0.2627372627372627 -0.3006993006993007*I , 0.4665334665334666 +0.3066933066933067*I , 0.2667332667332667 +0.4425574425574426*I , -0.3026973026973027 +0.2987012987012987*I , -0.1728271728271728 +0.1068931068931069*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.2227772227772228 -0.3486513486513487*I , 0.1068931068931069 +0.2427572427572428*I , -0.09290709290709291 -0.4925074925074925*I , 0 +0*I , -0.3326673326673327 +0.4425574425574426*I , 0.0969030969030969 +0.2987012987012987*I , 0.2267732267732268 +0.1068931068931069*I , 0.02697302697302697 -0.3566433566433567*I , -0.1428571428571428 -0.1008991008991009*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [0.3066933066933067 -0.4925074925074925*I , 0.3366633366633366 +0.05094905094905095*I , 0.06693306693306693 +0.4425574425574426*I , 0 +0*I , 0.2967032967032967 +0.4505494505494506*I , 0.4265734265734266 -0.3566433566433567*I , 0.2567432567432568 -0.1008991008991009*I , -0.2127872127872128 -0.09290709290709291*I , -0.4125874125874126 -0.1568431568431568*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ], [-0.4625374625374625 -0.3006993006993007*I , 0.2667332667332667 +0.3066933066933067*I , -0.3026973026973027 +0.4505494505494506*I , 0 +0*I , -0.3726273726273726 +0.1068931068931069*I , 0.4565434565434565 -0.1488511488511489*I , -0.01298701298701299 -0.1568431568431568*I , 0.2167832167832168 +0.4985014985014985*I , 0.1468531468531468 -0.2927072927072927*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I , 0 +0*I ]], datatype=complex[8], order='C_order'):
NAG:-f16zjc(side, uplo, transt, diag, alpha, t, b, 'm' = m, 'n' = n):
|