diff libcruft/blas-xtra/dconv2.f @ 14596:ad1d7a53b1ca stable

perform convolution properly independently of matrix orientation (bug #34893) * cconv2.f (cconv2i), csconv2.f (csconv2i), dconv2.f (dconv2i), sconv2.f (sconv2i), zconv2.f (zconv2i), zdconv2.f (zdconv2i): Fix iteration over the a matrix, so it's properly reversed. * conv2.cc: Uncomment test to check for this, and add a new xtest.
author Sunghyun Cho <sodomau@gmail.com>
date Tue, 17 Apr 2012 12:47:34 +0100
parents 72c96de7a403
children
line wrap: on
line diff
--- a/libcruft/blas-xtra/dconv2.f	Wed May 02 17:29:40 2012 -0400
+++ b/libcruft/blas-xtra/dconv2.f	Tue Apr 17 12:47:34 2012 +0100
@@ -70,7 +70,7 @@
       do k = 1,na-nb+1
         do j = 1,nb
           do i = 1,mb
-            call daxpy(ma-mb+1,b(i,j),a(mb+1-i,k+j-1),1,c(1,k),1)
+            call daxpy(ma-mb+1,b(i,j),a(mb+1-i,k+nb-j),1,c(1,k),1)
           end do
         end do
       end do