Commit 14132770 authored by Andrea taracchini's avatar Andrea taracchini

SEOBNRv3_opt Fix debugging output

Original: a2bdb25616e989edbd3db6a4663b4f5c4911e35d
parent 838f37e4
...@@ -3597,7 +3597,7 @@ int XLALSimIMRSpinEOBWaveformAll( ...@@ -3597,7 +3597,7 @@ int XLALSimIMRSpinEOBWaveformAll(
*hIMRlmJTSHiOutput = hIMRlmJTSHi; *hIMRlmJTSHiOutput = hIMRlmJTSHi;
if (debugPK){ if (debugPK){
out = fopen( "JIMRWavesHi.dat", "w" ); out = fopen( "JIMRWavesHi.dat", "w" );
for ( i = 0; i < retLenHi + retLenRDPatchHi; i++ ) for ( i = 0; i < (INT4) hIMR22JTSHi->data->length; i++ )
{ {
fprintf( out, fprintf( out,
"%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n", "%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n",
...@@ -3719,7 +3719,7 @@ int XLALSimIMRSpinEOBWaveformAll( ...@@ -3719,7 +3719,7 @@ int XLALSimIMRSpinEOBWaveformAll(
if (debugPK){ if (debugPK){
out = fopen( "JIMRWaves.dat", "w" ); out = fopen( "JIMRWaves.dat", "w" );
for ( i = 0; i < retLenLow + retLenRDPatchLow; i++ ) for ( i = 0; i < (INT4) hIMR22JTS->data->length; i++ )
{ {
fprintf( out, fprintf( out,
"%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n", "%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n",
...@@ -3849,23 +3849,23 @@ int XLALSimIMRSpinEOBWaveformAll( ...@@ -3849,23 +3849,23 @@ int XLALSimIMRSpinEOBWaveformAll(
*hIMRoutput = XLALSphHarmTimeSeriesAddMode( *hIMRoutput, hIMR2m1ITS, 2, -1 ); *hIMRoutput = XLALSphHarmTimeSeriesAddMode( *hIMRoutput, hIMR2m1ITS, 2, -1 );
*hIMRoutput = XLALSphHarmTimeSeriesAddMode( *hIMRoutput, hIMR2m2ITS, 2, -2 ); *hIMRoutput = XLALSphHarmTimeSeriesAddMode( *hIMRoutput, hIMR2m2ITS, 2, -2 );
XLALSphHarmTimeSeriesSetTData( *hIMRoutput, tlistRDPatch ); XLALSphHarmTimeSeriesSetTData( *hIMRoutput, tlistRDPatch );
}
if (debugPK){ if (debugPK){
out = fopen( "IWaves.dat", "w" ); out = fopen( "IWaves.dat", "w" );
for ( i = 0; i < retLenLow + retLenRDPatchLow; i++ ) for ( i = 0; i < (INT4) hIMR22ITS->data->length; i++ )
{ {
fprintf( out, fprintf( out,
"%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n", "%.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e %.16e\n",
tlistRDPatch->data[i], tlistRDPatch->data[i],
creal(hIMR22ITS->data->data[i]), cimag(hIMR22ITS->data->data[i]), creal(hIMR22ITS->data->data[i]), cimag(hIMR22ITS->data->data[i]),
creal(hIMR21ITS->data->data[i]), cimag(hIMR21ITS->data->data[i]), creal(hIMR21ITS->data->data[i]), cimag(hIMR21ITS->data->data[i]),
creal(hIMR20ITS->data->data[i]), cimag(hIMR20ITS->data->data[i]), creal(hIMR20ITS->data->data[i]), cimag(hIMR20ITS->data->data[i]),
creal(hIMR2m1ITS->data->data[i]), cimag(hIMR2m1ITS->data->data[i]), creal(hIMR2m1ITS->data->data[i]), cimag(hIMR2m1ITS->data->data[i]),
creal(hIMR2m2ITS->data->data[i]), cimag(hIMR2m2ITS->data->data[i]) ); creal(hIMR2m2ITS->data->data[i]), cimag(hIMR2m2ITS->data->data[i]) );
} }
fclose( out ); fclose( out );
} }
}
/* ********************************************************************************* /* *********************************************************************************
* ********************************************************************************* * *********************************************************************************
......
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