Commit 59d66e99 authored by thomas.forbriger's avatar thomas.forbriger
Browse files

cross [WP][FIX]: assign correct spectrum

program operates correctly in convolution mode; though results are not yet as
expected in cross-correlation mode
parent 3a035edc
......@@ -243,11 +243,11 @@ FourierProcessor::Tspectrum
FourierProcessor::Tspectrum result(sref.size());
result=FourierProcessor::Tspectrum::Tvalue(0.);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BR(sref);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BD(sref);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BD(ssig);
aff::Iterator<FourierProcessor::Tspectrum> IR(result);
while (BR.valid() && BD.valid() && IR.valid())
{
*IR = (*BR) * (*BD);
(*IR) = ((*BR) * (*BD));
++IR;
++BR;
++BD;
......@@ -268,7 +268,7 @@ FourierProcessor::Tspectrum
FourierProcessor::Tspectrum result(sref.size());
result=FourierProcessor::Tspectrum::Tvalue(0.);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BR(sref);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BD(sref);
aff::Browser<FourierProcessor::Tspectrum::Tcoc> BD(ssig);
aff::Iterator<FourierProcessor::Tspectrum> IR(result);
while (BR.valid() && BD.valid() && IR.valid())
{
......@@ -548,8 +548,7 @@ int main(int iargc, char* argv[])
TFXX_debug(opt.debug,
"main()", "prepare FFTreference "
<< TFXX_value(reference.header.wid2().dt));
Tts::Tdttimeseries dtreference(reference,
reference.header.wid2().dt);
Tts::Tdttimeseries dtreference(reference);
TFXX_debug(opt.debug,
"main()", TFXX_value(dtreference.header.dt));
FFTreference.setseries(dtreference);
......@@ -569,7 +568,7 @@ int main(int iargc, char* argv[])
{
TFXX_debug(opt.debug,
"main()", "prepare FFTdata");
Tts::Tdttimeseries dtdata(data, data.header.wid2().dt);
Tts::Tdttimeseries dtdata(data);
FFTdata.setseries(dtdata);
}
if (opt.convolve)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment