# HG changeset patch # User Jordi GutiƩrrez Hermoso # Date 1336406370 14400 # Node ID 34f067bcac1230469205f7f211456fd8f7feaf4c # Parent 28e53daab1f8fda81d8032283a112e9b20e992ea# Parent 6980b0f35df97c065a5917c484862bff0701bb88 maint: Periodic merge of stable to default diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/cconv2.f --- a/libcruft/blas-xtra/cconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/cconv2.f Mon May 07 11:59:30 2012 -0400 @@ -70,7 +70,7 @@ do k = 1,na-nb+1 do j = 1,nb do i = 1,mb - call caxpy(ma-mb+1,b(i,j),a(mb+1-i,k+j-1),1,c(1,k),1) + call caxpy(ma-mb+1,b(i,j),a(mb+1-i,k+nb-j),1,c(1,k),1) end do end do end do diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/csconv2.f --- a/libcruft/blas-xtra/csconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/csconv2.f Mon May 07 11:59:30 2012 -0400 @@ -76,7 +76,7 @@ do j = 1,nb do i = 1,mb btmp = b(i,j) - call caxpy(ma-mb+1,btmp,a(mb+1-i,k+j-1),1,c(1,k),1) + call caxpy(ma-mb+1,btmp,a(mb+1-i,k+nb-j),1,c(1,k),1) end do end do end do diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/dconv2.f --- a/libcruft/blas-xtra/dconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/dconv2.f Mon May 07 11:59:30 2012 -0400 @@ -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 diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/sconv2.f --- a/libcruft/blas-xtra/sconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/sconv2.f Mon May 07 11:59:30 2012 -0400 @@ -70,7 +70,7 @@ do k = 1,na-nb+1 do j = 1,nb do i = 1,mb - call saxpy(ma-mb+1,b(i,j),a(mb+1-i,k+j-1),1,c(1,k),1) + call saxpy(ma-mb+1,b(i,j),a(mb+1-i,k+nb-j),1,c(1,k),1) end do end do end do diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/zconv2.f --- a/libcruft/blas-xtra/zconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/zconv2.f Mon May 07 11:59:30 2012 -0400 @@ -70,7 +70,7 @@ do k = 1,na-nb+1 do j = 1,nb do i = 1,mb - call zaxpy(ma-mb+1,b(i,j),a(mb+1-i,k+j-1),1,c(1,k),1) + call zaxpy(ma-mb+1,b(i,j),a(mb+1-i,k+nb-j),1,c(1,k),1) end do end do end do diff -r 28e53daab1f8 -r 34f067bcac12 libcruft/blas-xtra/zdconv2.f --- a/libcruft/blas-xtra/zdconv2.f Mon May 07 08:52:17 2012 -0400 +++ b/libcruft/blas-xtra/zdconv2.f Mon May 07 11:59:30 2012 -0400 @@ -76,7 +76,7 @@ do j = 1,nb do i = 1,mb btmp = b(i,j) - call zaxpy(ma-mb+1,btmp,a(mb+1-i,k+j-1),1,c(1,k),1) + call zaxpy(ma-mb+1,btmp,a(mb+1-i,k+nb-j),1,c(1,k),1) end do end do end do diff -r 28e53daab1f8 -r 34f067bcac12 src/DLD-FUNCTIONS/conv2.cc --- a/src/DLD-FUNCTIONS/conv2.cc Mon May 07 08:52:17 2012 -0400 +++ b/src/DLD-FUNCTIONS/conv2.cc Mon May 07 11:59:30 2012 -0400 @@ -262,9 +262,18 @@ %% Test cases from Bug #34893 %!assert (conv2 ([1:5;1:5], [1:2], "same"), [4 7 10 13 10; 4 7 10 13 10]) %!assert (conv2 ([1:5;1:5]', [1:2]', "same"), [4 7 10 13 10; 4 7 10 13 10]') -%!#assert (conv2 ([1:5;1:5], [1:2], "valid"), [4 7 10 13; 4 7 10 13]) +%!assert (conv2 ([1:5;1:5], [1:2], "valid"), [4 7 10 13; 4 7 10 13]) %!assert (conv2 ([1:5;1:5]', [1:2]', "valid"), [4 7 10 13; 4 7 10 13]') +%!test +%! rand ("seed", 42); +%! x = rand (100); +%! y = ones (5); +%! A = conv2 (x, y)(5:end-4,5:end-4); +%! B = conv2 (x, y, "valid"); +%! assert (B, A); ## Yes, this test is for *exact* equivalence. + + %% Test input validation %!error conv2 () %!error conv2 (1)