From b8e92e9f6dcb3320123908c845d81ce0db19a2cc Mon Sep 17 00:00:00 2001
From: JanThorbecke <janth@xs4all.nl>
Date: Thu, 6 Feb 2020 11:47:36 +0100
Subject: [PATCH] fconv scaling with dt added
---
utils/fconv.c | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/utils/fconv.c b/utils/fconv.c
index 42686ef..d44c92a 100644
--- a/utils/fconv.c
+++ b/utils/fconv.c
@@ -507,7 +507,7 @@ void corr(float *data1, float *data2, float *cov, int nrec, int nsam, float dt,
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/(float)optn;
+ scl = dt/(float)optn;
crmfft(&ccov[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,cov,nsam);
@@ -578,7 +578,7 @@ void corr3(float *data1, float *data2, float *cov, int nrec, int nsam, float dt)
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/(float)optn;
+ scl = dt/(float)optn;
crmfft(&ccov[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,cov,ntout);
@@ -689,7 +689,7 @@ void deconv(float *data1, float *data2, float *decon, int nrec, int nsam,
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/(float)optn;
+ scl = dt/(float)optn;
crmfft(&cdec[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,decon,nsam);
@@ -777,7 +777,7 @@ void power(float *data1, float *data2, float *cov, int nrec, int nsam, float dt,
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/(float)optn;
+ scl = dt/(float)optn;
crmfft(&acov[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,cov,nsam);
@@ -866,7 +866,7 @@ void convol(float *data1, float *data2, float *con, int nrec, int nsam, float dt
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/((float)(optn));
+ scl = dt/((float)(optn));
crmfft(&ccon[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,con,nsam);
@@ -979,7 +979,7 @@ void cohr(float *data1, float *data2, float *cov, int nrec, int nsam, float dt,
/* inverse frequency-time FFT and scale result */
sign = 1;
- scl = 1.0/(float)optn;
+ scl = dt/(float)optn;
crmfft(&ccov[0], &rdata1[0], optn, nrec, nfreq, optn, sign);
scl_data(rdata1,optn,nrec,scl,cov,nsam);
--
GitLab