-
Notifications
You must be signed in to change notification settings - Fork 0
/
run.jl
55 lines (50 loc) · 1.38 KB
/
run.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
using LinearAlgebra
using Distributions
using SpecialFunctions
using Printf
using StatsBase
using Base.Threads
using LsqFit
using SIMD
using Glob
#include("shiftcheck_cgp_main.jl")
include("cgp_main.jl")
T=true;F=false
ktype = logMat52 ## Mat52/logMat52/Mat32/logMat32/RBF/logRBF
Kernel= string(ktype)
multihw=false
const sigRfac = 0.1
const numN = 20000
const mstep = 3000
Monotonic=T
Convex= T
#inps = glob("./data/*.txt")
inps = glob("sample_ncsmdata.dat")
#inps = glob("./data/ncsm_6Li_hw16_lam202_N3LO_6-14.txt")
#inps = glob("./data/Shin_6Li_ncsm_JISP16_2-18.txt")
for paramean in [F,T]
if paramean==T
qT = 2.e-1; qY = 1.e-2; qYfac = 5.e-2
else
qT = 2.e-1; qY = 1.e-2; qYfac = 1.e-1
end
for inpname in inps
if occursin("N3LO",inpname); inttype="N3LO"
elseif occursin("JISP16",inpname); inttype="JISP16"
elseif occursin("NNLOopt",inpname); inttype="NNLOopt"
else; inttype="unknownint";end
println("inttype $inttype parametric mean:$paramean ")
print(" numN:$numN step:$mstep Kernel:$Kernel")
print(" Monotonic:$Monotonic Convex:$Convex ")
if inttype=="N3LO"
xpMax = 38
else
xpMax = 46
end
@time Em,Ev=main(
mstep,numN,sigRfac,ktype,
inpname,inttype,
xpMax,Monotonic,Convex,paramean,
qT,qY,qYfac)
end
end