Wrong 1x1 triangular matrix-vector product
Submitted by Chen-Pang He
Assigned to Nobody
Link to original bugzilla bug (#839)
Description
Created attachment 473
Demo reproducing the bug
Our internal::triangular_matrix_vector_product is specialized for triangular matrix-vector product. When 1x1 times 1xn, it ends prematurely by computing only one entry.
The below is the output from the program compiled with the attachment.
Matrix:
0.144509
0.0973679
-0.228021
0.190039
-0.105256
-0.131808
0.0489002
0.160188
Triangular:
0.144509
0
0
0
0
0
0
0
?TRMV:
0.144509
0
0
0
0
0
0
0
?TRMM:
0.144509
0.0973679
-0.228021
0.190039
-0.105256
-0.131808
0.0489002
0.160188
Attachment 473, "Demo reproducing the bug":
tri.cpp