Mercurial > octave-nkf

comparison libcruft/blas-xtra/zconv2.f @ 14596:ad1d7a53b1ca stable

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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
comparison
equal deleted inserted replaced
14590:6250e1232c9c 14596:ad1d7a53b1ca
68 integer i,j,k 68 integer i,j,k
69 external zaxpy 69 external zaxpy
70 do k = 1,na-nb+1 70 do k = 1,na-nb+1
71 do j = 1,nb 71 do j = 1,nb
72 do i = 1,mb 72 do i = 1,mb
73 call zaxpy(ma-mb+1,b(i,j),a(mb+1-i,k+j-1),1,c(1,k),1) 73 call zaxpy(ma-mb+1,b(i,j),a(mb+1-i,k+nb-j),1,c(1,k),1)
74 end do 74 end do
75 end do 75 end do
76 end do 76 end do
77 end subroutine 77 end subroutine