Commit 55e55d1d authored by thomas.forbriger's avatar thomas.forbriger
Browse files

croposp [FIX]: fix phase computation

parent 08765297
......@@ -668,12 +668,12 @@ int main(int iargc, char* argv[])
{
if ((k!=i) && (l!=i))
{
psd::TDseries::Tcoc cpsd1, cpsd2;
cpsd1=psd::arg(CPSD_vector[pairs(k,i)].series);
if (pairs.swap(k,i)) { cpsd1 = -1.*cpsd1; }
cpsd2=psd::arg(CPSD_vector[pairs(l,i)].series);
if (pairs.swap(l,i)) { cpsd2 = -1.*cpsd1; }
phase_vector[triples(i,k,l)].series=cpsd1-cpsd2;
psd::TDCseries::Tcoc cpsd1, cpsd2;
cpsd1=CPSD_vector[pairs(k,i)].series;
if (pairs.swap(k,i)) { cpsd1 = psd::conj(cpsd1); }
cpsd2=CPSD_vector[pairs(l,i)].series;
if (pairs.swap(l,i)) { cpsd2 = psd::conj(cpsd2); }
phase_vector[triples(i,k,l)].series=psd::arg(cpsd1/cpsd2);
phase_vector[triples(i,k,l)].label
=opt.prefix_phase+" "+vector_of_metadata[k].label
+"-"+vector_of_metadata[l].label+" (ref: "
......
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