[Octave-bug-tracker] [bug #66882] Convolution code path improvements

[Dmitri A. Sergatskov]

Follow-up Comment #70, bug #66882 (group octave):

I see (both with gcc and clang):

octave:1> A=ones(1e5,1);
octave:2> B=A';
octave:3> tic; conv2(B,B); toc
Elapsed time is 0.955019 seconds.
octave:4> tic; conv2(complex(B),complex(B)); toc
Elapsed time is 8.06255 seconds.
octave:5> tic; conv2(complex(A),complex(A)); toc
Elapsed time is 7.91182 seconds.

I also tried to remove fall-back (loops) specialization and I see:

../liboctave/numeric/oct-convn.cc:77:1: note: candidate function not viable:
no known conversion from 'std::complex<double>' to 'const F77_DBLE_CMPLX' (aka
'const _Complex double') for 2nd argument
   77 | blas_axpy (const F77_INT& n, const F77_DBLE_CMPLX& alpha,
      | ^                            ~~~~~~~~~~~~~~~~~~~~~~~~~~~

(and many others like that).

Dmitri.
-- 



    _______________________________________________________

Reply to this item at:

  <https://savannah.gnu.org/bugs/?66882>

_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/

signature.asc
Description: PGP signature