Re: [Pdl-general] [Pdl-devel] lu_backsub

2021-06-07 Thread Luis Mochan
Hi Ed, On Tue, Jun 01, 2021 at 12:57:51AM +, Ed . wrote: > That is what I was aiming to achieve (and have captured in hopefully-thorough > tests). My thinking is that the linear algebra stuff is founded on matrices, > with at least 2 dimensions. Therefore, getting something back that has at

Re: [Pdl-general] PDL and FFTW3

2021-06-07 Thread Ingo Schmid
Hi Daniel, I think what you want to do is 1d real valued convolution? Have you looked at PDL::Image2D / PDL::ImageND and the ocnv*  functions? Reshape your q and t with secnod dim 1 and try conv2d. Of course you can do rfft1 / irfft1 directly. Instead of using cmul try a simple *. Slicing is