mirror of
https://asciireactor.com/otho/psdlag-agn.git
synced 2024-11-22 06:05:06 +00:00
replaced src with new version again because it was unclear that all source files were up-to-date.
This commit is contained in:
parent
2d9d4ebf44
commit
5c3cf1ea3b
47
src/.smhist
47
src/.smhist
@ -1,47 +0,0 @@
|
|||||||
data check.dat
|
|
||||||
read f 1
|
|
||||||
read pow 2
|
|
||||||
limits f pow
|
|
||||||
box
|
|
||||||
ptype 10 3
|
|
||||||
points f pow
|
|
||||||
quit
|
|
||||||
data check.dat
|
|
||||||
read f 1
|
|
||||||
read pow 2
|
|
||||||
limits f pow
|
|
||||||
ptype 10 3
|
|
||||||
points f pow
|
|
||||||
box
|
|
||||||
quit
|
|
||||||
data check.dat
|
|
||||||
read f 1
|
|
||||||
read pow 2
|
|
||||||
limits f pow
|
|
||||||
box
|
|
||||||
ptype 10 3
|
|
||||||
points f pow
|
|
||||||
set lgf = log10(f)
|
|
||||||
set lgf = alog10(f)
|
|
||||||
set lgf = lg(f)
|
|
||||||
erase
|
|
||||||
box
|
|
||||||
erase
|
|
||||||
limits lgf pow
|
|
||||||
points lgf pow
|
|
||||||
box
|
|
||||||
erase
|
|
||||||
limits lgf 0. 3.
|
|
||||||
box
|
|
||||||
points lgf pow
|
|
||||||
quit
|
|
||||||
data check.dat
|
|
||||||
read 1 f
|
|
||||||
read f 1
|
|
||||||
read pow 2
|
|
||||||
set lgf = lg(f)
|
|
||||||
limits lgf pow
|
|
||||||
box
|
|
||||||
ptype 10 3
|
|
||||||
points lgf pow
|
|
||||||
quit
|
|
32
src/README
32
src/README
@ -1,32 +0,0 @@
|
|||||||
|
|
||||||
|
|
||||||
|
|
||||||
To compile, modify the file makefile to so that libdir and incdir point to where
|
|
||||||
alglib library and header files are installed respectively, then run: make
|
|
||||||
|
|
||||||
|
|
||||||
FILES:
|
|
||||||
- README: this file
|
|
||||||
- makefile: a make file for compiling
|
|
||||||
- main.cpp: the main file, that process the input file, print help message if
|
|
||||||
needed and reads the light curves files. check readLC to see the light curve
|
|
||||||
format.
|
|
||||||
- mod.cpp: contains the class mod, which is a base class that uses the likelihood method. This on its own does nothing, but it is meant to be inherited by others, e.g. psd , lag , lag10
|
|
||||||
etc ... It is here where the likelihood is calculated, and maximized.
|
|
||||||
|
|
||||||
- psd.cpp: a child class of mod to calculate psds. psd10 is similar by it
|
|
||||||
calculate the log10 of psd. This is better than calculating the psd directly.
|
|
||||||
|
|
||||||
- lag.cpp: a child class of mod to calculate cross spectra and phase lags.
|
|
||||||
similar to psd. there is also lag10 that calcualte the log10 of the cross
|
|
||||||
spectrum plus phase lag.
|
|
||||||
|
|
||||||
- psdlag.cpp: a child class of mod to calculate psd and cross spec and phase
|
|
||||||
lag at the same time.
|
|
||||||
|
|
||||||
- inc: contains the header files.
|
|
||||||
|
|
||||||
- inc/mod.hpp: the header for mod.cpp. This file has the Mod class also, which
|
|
||||||
is a class that holds multiple mod's used for fitting multiple light curve
|
|
||||||
simultaneusly, each mod fits for one light curve, and the total log-likelihood
|
|
||||||
is the sum of individual likelihood (i.e. product of probabilities)
|
|
@ -1,6 +0,0 @@
|
|||||||
2.5e-5 4.681
|
|
||||||
5.5e-5 3.358
|
|
||||||
8.5e-5 2.911
|
|
||||||
1.5e-4 2.650
|
|
||||||
2.5e-4 1.864
|
|
||||||
3.5e-4 0.605
|
|
@ -1,6 +0,0 @@
|
|||||||
1.022e-01 9.694e-01
|
|
||||||
3.260e-01 9.192e-01
|
|
||||||
5.118e-01 9.776e-01
|
|
||||||
3.254e-01 9.269e-01
|
|
||||||
1.000e+00 9.659e-01
|
|
||||||
1.000e+00 9.891e-01
|
|
@ -2,10 +2,8 @@
|
|||||||
|
|
||||||
#libdir='/eos/azoghbi/soft/usr/lib'
|
#libdir='/eos/azoghbi/soft/usr/lib'
|
||||||
#incdir='/eos/azoghbi/soft/usr/include'
|
#incdir='/eos/azoghbi/soft/usr/include'
|
||||||
libdir='/home/caes/science/psdlag-agn/src'
|
libdir='/home/caes/science/psdlag-agn/src/inc'
|
||||||
incdir='/home/caes/science/psdlag-agn/src/inc'
|
incdir='/home/caes/science/psdlag-agn/src/inc'
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
psdlag:
|
psdlag:
|
||||||
g++ *cpp -o psdlag -O3 -Wall -lalglib -I${incdir} -L${libdir}
|
g++ *cpp -o psdlag -O3 -Wall -lalglib -I${incdir} -L${libdir}
|
||||||
|
@ -1,2 +0,0 @@
|
|||||||
|
|
||||||
|
|
Loading…
Reference in New Issue
Block a user