fold - E. A. Robinson’s subroutine to perform a convolution
call fold(la, a, lb, b, lc, c)
Basic Seismic
Utilities (BSU) subroutine originally published by Robinson
(1967) which computes a convolution between "a"
and "b", returning the output in "c".
Fortran Version.
Arguments
la (int)
Number of
samples in a
a (float)
One Signal
lb (int)
Number of
samples in b
b (float)
The other signal
lc (int)
Number of
samples in c (lc=la+lb-1)
c (float)
The result of convolving a with b.
NOTE:
The dimension, "lc" should be set by the calling
program as indicated above.
cross(3), dot(3), eureka(3), shape(3)
Robinson, E.A.,
1967, "Multichannel Time Series Analysis
with Digital Computer Programs", Holden-Day, 298p.
No known bugs.
Copyright © 1967 by Holden-Day Inc.
Holden-Day Inc., Liquidation Trust, 1259 S.W. 14th Street Boca Raton, FL 33486 Phone: 561.750.9229 Fax: 561.394.6809
This subroutine is included for distribution with Basic Seismic Utilities (BSU) by joint permission of the author, Enders Robinson, and Holden-Day Inc. Liquidation Trust, Frederick H. Murphy, Trustee. Being part of the BSU package, it is governed by the terms of the GPL license, and is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
E. A. Robinson