diff libcruft/blas-xtra/zdconv2.f @ 10946:1094868ca10d

fix bugs in inner convolution
author Jaroslav Hajek <highegg@gmail.com>
date Tue, 07 Sep 2010 12:23:01 +0200
parents 5af0b4bb384d
children fd0a3ac60b0e
line wrap: on
line diff
--- a/libcruft/blas-xtra/zdconv2.f
+++ b/libcruft/blas-xtra/zdconv2.f
@@ -76,7 +76,7 @@
         do j = 1,nb
           do i = 1,mb
             btmp = b(i,j)
-            call zaxpy(ma-mb+1,btmp,a(i,k+j-1),1,c(1,k),1)
+            call zaxpy(ma-mb+1,btmp,a(mb+1-i,k+j-1),1,c(1,k),1)
           end do
         end do
       end do