COPY_K(m, b, incb, buffer, 1);
}
- for (is = 0; is < m; is += DTB_ENTRIES){
+/*FIXME the GEMV unrolling performed here was found to be broken, see issue 1332 */
+/* Multiplying DTB size by 100 is just a quick-and-dirty hack to disable it for now[B */
- min_i = MIN(m - is, DTB_ENTRIES);
+ for (is = 0; is < m; is += DTB_ENTRIES * 100){
+
+ min_i = MIN(m - is, DTB_ENTRIES * 100);
#ifndef TRANSA
if (is > 0){
+fprintf(stderr,"WARNING unrolling of the trmv_U loop may give wrong results\n");
GEMV_N(is, min_i, 0, dp1,
a + is * lda, lda,
B + is, 1,