-
Notifications
You must be signed in to change notification settings - Fork 0
/
main.jl
42 lines (37 loc) · 1.14 KB
/
main.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
include("SensitivityAlgorithms.jl")
abstolset = 1e-6;
reltolset = 1e-5;
#############
# PKA model #
#############
modelName = "PKA";
include("./models/PKA.jl");
tspan = (0.0,600.0);
nIter = 100;
avgTimeFS_PKA, avgTimePBS_PKA, avgTimeExp_PKA = timeSensitivityAlgorithms(modelName, tspan, nIter);
################
# CaMKII model #
################
modelName = "CaMKIIs";
include("./models/CaMKIIs.jl");
tspan = (0.0,600.0);
nIter = 100;
avgTimeFS_CaMKII, avgTimePBS_CaMKII, avgTimeExp_CaMKII = timeSensitivityAlgorithms(modelName, tspan, nIter);
########################
# Chua's circuit model #
########################
modelName = "ChuaCircuit";
include("./models/ChuaCircuit.jl");
tspan = (0.0,10.0);
nIter = 10;
avgTimeFS_Chua, avgTimePBS_Chua, avgTimeExp_Chua = timeSensitivityAlgorithms(modelName, tspan, nIter);
########################
# Random Linear System #
########################
include("./models/RandomLinearSystem.jl")
startDim = 5;
endDim = 20;
stepSize = 5;
nIter = 10;
tspan = (0.0, 10.0);
avgTimeFS_RandomLinearSystem, avgTimePBS_RandomLinearSystem, avgTimeExp_RandomLinearSystem = runtimePerDimension(startDim, endDim, stepSize, nIter, tspan);