KAGALI PE code development Hyung Won Lee Inje

KAGALI PE code development Hyung Won Lee, Inje University, 29 August 2017 coll. With Jeongcho Kim(Inje U. ), Chunglee Kim(KASI) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 1

Contents Motivation Parameter Estimation with Bayesian Inference Implementation in LALSuite Strategy for KAGALI Development for KAGALI Preliminary test run Discussions 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 2

Motivation Implement MCMC PE for KAGALI Based on LALSuite Minimum full functionality for MCMC Parallel tempered Improvement as many as possible from LALSuite Modularize as much as possible for easy maintenance and readable Test for i. KAGRA data 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 3

Parameter Estimation with Bayesian inference 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 4

Parameter Estimation with Bayesian inference MCMC Samples Post Process 2017 -08 -29 PTMCMCOutput. #. nn. h 5 PTMCMCOutput. #. h 5 cbc. Bayes. Post. Proc. py KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 5

Waveform Templates There are many waveform templates Concentrate Taylor. F 2 with • Amplitude correction up to 2. 5 PN • Effect of eccentricity up to 3. 0 PN 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 6

Likelihood Calculation 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 7

MCMC Process 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 8

Parallel Tempering 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 9

Implementation in lalsuite LAL MCMC Pipeline Directory structure of lalsuite Basic entities of LALSuite LAL applications Likelihood, proposal and prior functions Simulation related functions and waveforms Post Processing python script 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 10

LAL MCMC Pipeline Start Initialization MPI synchronize PT MCMC Clearance End 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 11

main LALInference. Init. Run. State() LALInference. Inject. Signal() LALInference. Apply. Calibration. Errors() init_ptmcmc() LALInference. Init. CBCPrior() LALInference. Init. Likelihood() LALInference. Draw. Threads() run. State->algorithm 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 12

Init. Run. State() 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 13

LALInference. Inject. Inspiral. Signal() 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 14

LALInference. Apply. Calibration. Errors() 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 15

init_ptmcmc() Set Algorithm function pointer Set parallel swapping function pointer Set CBC model Set number of varying parameters Set temperature ladder Set CBCthreads if necessary Set proposal function and proposal arguments for all threads Set adaptation parameters 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 16

LALInference. Inspiral. Prior() Set various prior functions 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 17

PTMCMCAlgorithm() 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 18

Template LALInference. Template. XLALSim. Inspiral. Choose. Waveform 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 19

Strategy for KAGALI Base on LALSuite for many functions Modularize if possible Use structures • • • Command parameters(KGLProcess. Params. Table) IFO information(KGLInference. IFOData) Waveform parameters(KGLWaveform. Model. Params) MCMC status(KGLInference. Run. State) Simulation parameters Use simple text file for injections Implement lalinference 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 20

Development for KAGALI Use a new branch “mcmc” for development Library and application directory are separated (version 4) • • • root/kagali/inference/src root/kagali/kglcommon/clock/src root/kagali/kglcommon/detector/src root/kagali/kglcommon/std/src root/kagali/waveform/src root/kagaliapps/cbc/mpisrc Cloned on KISTI and local Ubuntu server Build and run on both systems 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 21

Main features Use a unique union structure for waveform parameters (KGLWaveform. Model. Params) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 22

Structures KGLInference. IFOData KGLDetector KGLWaveform. Model. Params. Type KGLWaveform. Model. Params. Value KGLProcess. Params. Table KGLProcess. Params. Pair KGLWaveform. Model. Params. Item KGLInference. Algorithm KGLInference. Run. State KGLWaveform. Model. Params. Index KGLWaveform. Model. Params. Value KGLInference. Algorithm. One. Step 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 23

2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 24

2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 25

Development for KAGALI kagali/inference/src 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 26

Development for KAGALI kagali/kglcommon/clock/src 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 27

Development for KAGALI kagali/kglcommon/detector/src 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 28

Development for KAGALI kagali/kglcommon/std/src 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 29

Development for KAGALI kagali/waveform/src 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 30

Development for KAGALI kagaliapps/cbc/mpisrc 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 31

KGLProcess. Params. Table. h. in 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 32

KGLProcess. Params. Table. h. in KGLProcess. Params. Pair KGLProcess. Params. Table 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 33

KGLCommand. Utils. c 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 34

KGLInference. MCMC. c 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 35

inference. MCMC. c 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 36

inference. MCMC. c 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 37

Makefile. am 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 38

build KAGALI 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 39

build script 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 40

build script(cont’d) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 41

Build results nohup 2017 -08 -29 . /build-branch. sh & KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 42

Sample Run(condor) /data/kagra/home/hwlee/projects/LAL/s cratch/kagali/mcmc. sub condor_submit mcmc. sub 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 43

Sample Run(condor) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 44

Sample Run(condor) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 45

Sample Run(condor) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 46

Sample Run(condor) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 47

Sample run(mpirun) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 48

Sample run(con’d) 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 49

Progress from the last Defines some structures • Waveform Parameters, Datetime, etc. Implement basic placeholder functions Implement few more placeholder for functions Fixing some nil point errors 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 50

Discussions Basic data structures are ready Basic program architecture fixed Need to write source codes and hurry up Hanyang, SNU Seoul Toyama KISTI, NIMS Daejeon Pusan U. Gimhae 2017 -08 -29 Osaka KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 51

Thanks 2017 -08 -29 KAGALI PE code development, The 17 th KAGRA F 2 F meeting, Toyama 52