|
|
fnpeaks, mpeaks
Zbigniew Kołaczkowski, Astronomical Institute, University of Wroclaw
We deliver source codes of two programs designed for efficient calculation
of Fourier amplitude spectra. The first one, FNPEAKS, calculates amplitude
spectrum for a single-object time-series data. In order to speed-up the
calculations, the algorithm uses the modification proposed by Kurtz
(1985, MNRAS 213,773) and implements SSE instructions available
in all modern processors.
The output from the program is the amplitude spectrum (option -f)
and the list of highest peaks in this spectrum.
The second program, MPEAKS, is designed for the same task, but for a sample
of objects with the same time vector. A typical example for this kind of
data is differential photometry for objects located in the same field-of-view
of CCD detector. The program performs calculations using LAPACK and ATLAS,
libraries optimized for linear algebra,
(http://www.netlib.org/lapack/, http://math-atlas.sourceforge.net).
In fact, the time vector for analysed objects does not
have to be exactly the same for all objects; it has only to be a subsample
of a predefined time vector. This allows using the program in a case when
some point-rejection procedures are applied in the analysis.
The output (for each object) is the same as for FNPEAKS.
FNPEAKS
Computing of Fourier amplitude spectrum for time series data.
Usage: fnpeaks [-f] <DataFile> <StartFr> <EndFr> <StepFr>
option -f : writing spectrum into file (with extension .trf).
---------------------------------------------------------------------
MPEAKS
Usage: mpeaks [<options>] <time vector> <datacat>
<options>:
-f=<first_fr>,<last_fr>,<fr_step> (default: 0,30,computed)
-f=<first_fr>,<last_fr> (<fr_step> is computed)
-fc=<centr_fr>,<delta_fr>,<fr_step>
-fc=<centr_fr>,<delta_fr> (<fr_step> is computed)
-fc=<centr_fr> (<delta_fr>=0.2, <fr_step> is computed)
-append (data are appended to *.max and *.trf files)
-trf (spectrum is written into *.trf file)
-dt=<max_time_diff> (in days)
-minobs=<min_no_of_obs>
-clean (light-curve cleaning)
-k=<kappa_clean>
Input file format: <time> <signal> [omitted fields]
Most prominent peaks are written into *.max file
Some important parameters are written into log-file "mpeaks.log"
Source codes
Questions/remarks: kolaczkowski at astro.uni.wroc.pl
|
