This is a Java translation of the Fortran version of Grant Foster's WWZ algorithm. The original (C) notice from the Fortran code is included below. As per the documentation accompanying the program, written (email) permission was sought, and granted, from the AAVSO Director to use the Fortran code in this way. Note that the comment about maximum data points below does not apply here, since we dynamically allocate the arrays based upon the number of observations.
WEIGHTED WAVELET Z-TRANSFORM This is a fortran version of Grant Foster's WWZ1.1.BAS BASIC program, stripped down for speed and readability improvements. It is modestly less flexible than the basic code (esp regarding varying input file formats), but will be easy to modify to suit your needs. In the event you have more than 100000 data points, you can resize all arrays to meet your needs -- the only limitation is your free memory. (Note that all variables are double precision, though.) A description of the mathematics can be found in G. Foster, "Wavelets for Period Analysis of Unevenly Sampled Time Series", Astronomical Journal 112, 1709 (Oct 1996).
-Matthew Templeton, August 15, 2002
(C) Copyright 1996, 2002 by the American Association of Variable Star Observers; all rights reserved.
|
|
|
|