WROCŁAWSKA STRONA HELAS-a

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

Webmaster: