convolve2D: Fast Column-wise Convolution of a Matrix

In waveslim: Basic Wavelet Routines for One-, Two-, and Three-Dimensional Signal Processing

Fast Column-wise Convolution of a Matrix

Description

Use the Fast Fourier Transform to perform convolutions between a sequence and each column of a matrix.

Usage

convolve2D(x, y, conj = TRUE, type = c("circular", "open"))

Arguments

x	MxN matrix.
y	numeric sequence of length N.
conj	logical; if TRUE, take the complex conjugate before back-transforming (default, and used for usual convolution).
type	character; one of circular, open (beginning of word is ok). For circular, the two sequences are treated as circular, i.e., periodic. For open and filter, the sequences are padded with zeros (from left and right) first; filter returns the middle sub-vector of open, namely, the result of running a weighted mean of x with weights y.

Details

This is a corrupted version of convolve made by replacing fft with mvfft in a few places. It would be nice to submit this to the R Developers for inclusion.

Author(s)

B. Whitcher

See Also

convolve

waveslim documentation built on June 22, 2024, 9:43 a.m.