Problem Extension
Problem extension can be achieved easily via PRTYPE which is of type Val, or another subtype of the superclass can be created.
QuantizedSystemSolver.NLODEProblem — TypeNLODEProblem{PRTYPE,T,Z,Y,CS} This is superclass for all NLODE problems. It is parametric on these:
- The problem type PRTYPE.
+Problem Extension · Quantized System Solver Problem Extension
Problem extension can be achieved easily via PRTYPE which is of type Val, or another subtype of the superclass can be created.
QuantizedSystemSolver.NLODEProblem — TypeNLODEProblem{PRTYPE,T,Z,Y,CS} This is superclass for all NLODE problems. It is parametric on these:
- The problem type PRTYPE.
- The number of continuous variables T
diff --git a/docs/build/algorithmDev/index.html b/docs/build/algorithmDev/index.html
index 2882f29..9658ca2 100644
--- a/docs/build/algorithmDev/index.html
+++ b/docs/build/algorithmDev/index.html
@@ -1,4 +1,5 @@
Algorithm Extension · Quantized System Solver Algorithm Extension
Currently only QSS1,2,3 ; LiQSS1,2,3 ; and mLiQSS1,2,3 exist. Any new algorithm can be added via the name N which is of type Val, with O is the order which also of type Val. For example qss1 is created by:
qss1()=QSSAlgorithm(Val(:qss),Val(1))
QuantizedSystemSolver.ALGORITHM — TypeALGORITHM{N,O} This is superclass for all QSS algorithms. It is parametric on these:
- The name of the algorithm N
+Algorithm Extension · Quantized System Solver Algorithm Extension
Currently only QSS1,2,3 ; LiQSS1,2,3 ; and mLiQSS1,2,3 exist. Any new algorithm can be added via the name N which is of type Val, with O is the order which also of type Val. For example qss1 is created by:
qss1()=QSSAlgorithm(Val(:qss),Val(1))
QuantizedSystemSolver.ALGORITHM — TypeALGORITHM{N,O} This is superclass for all QSS algorithms. It is parametric on these:
- The name of the algorithm N
- The order of the algorithm O
sourceSettings
This document was generated with Documenter.jl version 1.5.0 on Wednesday 23 October 2024. Using Julia version 1.9.2.
Examples · Quantized System Solver Examples
Systems of 2 Linear Time Invariant Differential equations
odeprob = NLodeProblem(quote
Examples · Quantized System Solver Examples
Systems of 2 Linear Time Invariant Differential equations
odeprob = NLodeProblem(quote
name=(sysb53,)
u = [-1.0, -2.0]
@@ -6,12 +7,14 @@
du[2] =1.24*u[1]-0.01*u[2]+0.2
end)
tspan=(0.0,1.0)
This is a great example that shows when we need to use the explicit qss, the implicit liqss, or the modified implicit nmliqss. This is a stiff problem so we need to use the implicit methods, but it also contains larger entries outside the main diagonal of the Jacobian. Therefore, nmliqss should the most appropriate algorithm to use.
sol=solve(odeprob,qss1(),tspan)
+tspan=(0.0,1.0)
This is a great example that shows when we need to use the explicit qss, the implicit liqss, or the modified implicit nmliqss. This is a stiff problem so we need to use the implicit methods, but it also contains larger entries outside the main diagonal of the Jacobian. Therefore, nmliqss should the most appropriate algorithm to use.
sol=solve(odeprob,qss1(),tspan)
save_Sol(sol)
sol=solve(odeprob,qss2(),tspan)
save_Sol(sol)
sol=solve(odeprob,liqss1(),tspan)
save_Sol(sol)
sol=solve(odeprob,liqss2(),tspan)
save_Sol(sol)
sol=solve(odeprob,nmliqss1(),tspan)
save_Sol(sol)
sol=solve(odeprob,nmliqss2(),tspan)
save_Sol(sol)
The nmliqss plot does not finish at the final time because it terminated when it reached the equilibrium in which the values are the same as the values at the final time.
The Tyson Model
function test(solvr,absTol,relTol)
+save_Sol(sol)
The nmliqss plot does not finish at the final time because it terminated when it reached the equilibrium in which the values are the same as the values at the final time.
The Tyson Model
function test(solvr,absTol,relTol)
odeprob = NLodeProblem(quote
name=(tyson,)
u = [0.0,0.75,0.25,0.0,0.0,0.0]
@@ -40,6 +43,25 @@
du[1] = u[2]
du[2] = (1.0-u[1]*u[1])*u[2]-u[1]
+end)
+tspan=(0.0,10.0)
+sol=solve(odeprob,nmliqss2(),tspan)
+save_Sol(sol)
odeprob = NLodeProblem(quote
+ name=(oregonator,)
+ u = [1.0,1.0,0.0]
+ du[1] = 100.8*(9.523809523809524e-5*u[2]-u[1]*u[2]+u[1]*(1.0-u[1]))
+ du[2] =40320.0*(-9.523809523809524e-5*u[2]-u[1]*u[2]+u[3])
+ du[3] = u[1] -u[3]
+end)
+tspan=(0.0,10.0)
+sol=solve(odeprob,nmliqss2(),tspan)
+save_Sol(sol)
Bouncing Ball
odeprob = NLodeProblem(quote
+end
This model also is stiff and it needs a stiff method, but also the normal liqss will produce unnecessary cycles. Hence, the nmliqss is again the most appropriate. 
Oregonator; Vanderpl
odeprob = NLodeProblem(quote
+ name=(vanderpol,)
+ u = [0.0,1.7]
+ du[1] = u[2]
+ du[2] = (1.0-u[1]*u[1])*u[2]-u[1]
+
end)
tspan=(0.0,10.0)
sol=solve(odeprob,nmliqss2(),tspan)
diff --git a/docs/build/userGuide/index.html b/docs/build/userGuide/index.html
index 7aa96ae..d2d3fc6 100644
--- a/docs/build/userGuide/index.html
+++ b/docs/build/userGuide/index.html
@@ -53,17 +53,24 @@
end)
Explanation:
Parameters: These are constants used in the model. C, L, R, U, T, DC, ROn, ROff are physical constants related to the buck converter circuit.
Discrete and continuous variables: These represent the initial states of discrete and continuous variables in the model.
Renaming variables: For convenience, the elements of discrete and u are renamed to more descriptive variable names (rd, rs, nextT, lastT, diodeon, il, uc).
Helper equations: These are intermediate expressions needed for the differential equations. id is the current through the diode.
Differential equations: These represent the system's dynamics. du[1] is the rate of change of the inductor current. du[2] is the rate of change of the capacitor voltage.
Events: These are conditions that modify the continous and discrete variables. They handle the switching behavior of the buck converter and the diode's state transitions.
Step 2: Solve the ODE Problem
Next, we solve the ODE problem using the solve function. We need to specify the problem, the algorithm, and the simulation settings.
# Define the time span for the simulation
tspan = (0.0, 0.001) # Start at 0 seconds, end at 0.001 seconds
# Solve the ODE problem with the chosen algorithm and settings
-sol = solve(odeprob, nmliqss2(), tspan, abstol=1e-4, reltol=1e-3)
Explanation:
Time span: tspan defines the interval over which the solution is computed, from 0.0 to 0.001 seconds.
Solver function: solve is used to compute the solution.
odeprob is the ODE problem we defined.
nmliqss2() specifies the algorithm used to solve the problem (e.g., qss2,nmliqss1... might be other algorithms).
abstol and reltol are the absolute and relative tolerances for the solver, controlling the accuracy of the solution.
Step 3: Query the Solution
After solving the problem, we can query the solution to extract useful information such as variable values at specific times, the number of steps, events, and more.
# Get the value of variable 2 at time 0.0005
+sol = solve(odeprob, nmliqss2(), tspan, abstol=1e-4, reltol=1e-3)
Explanation:
Time span: tspan defines the interval over which the solution is computed, from 0.0 to 0.001 seconds.
Solver function: solve is used to compute the solution.
odeprob is the ODE problem we defined.
nmliqss2() specifies the algorithm used to solve the problem (e.g., qss2,nmliqss1... might be other algorithms).
abstol and reltol are the absolute and relative tolerances for the solver, controlling the accuracy of the solution.
Step 3: Query the Solution
After solving the problem, we can query the solution to extract useful information such as variable values at specific times, the number of steps, events, and more.
+# returns a vector of the values of all variables at time 0.0005 # feature not available in v1.0.1
+sol(0.0005)
+
+# Get the value of variable 2 at time 0.0005
value_at_time = sol(2, 0.0005)
+# Get the value of variable 2 at time 0.0005
+value_at_time = sol(0.0005,idxs=2)# feature not available in v1.0.1
+
# Get the total number of steps to end the simulation
-total_steps = sol.totalSteps
+total_steps = sol.totalSteps # feature available only in v1.0.1 ->sol.stats in later versions
# Get the number of simultaneous steps during the simulation
-simul_step_count = sol.simulStepCount
+simul_step_count = sol.simulStepCount # feature available only in v1.0.1 ->sol.stats in later versions
# Get the total number of events during the simulation
-event_count = sol.evCount
+event_count = sol.evCount # feature available only in v1.0.1 ->sol.stats in later versions
# Get the saved times and variables
saved_times = sol.savedTimes
@@ -79,6 +86,15 @@
plot_Sol(sol,1)
# Display the plot
display(plot_obj)
+# Generate a plot object for all variables of the solution. It requires "using Plots". feature not available in v1.0.1.
+plot(sol)
+# Generate a plot object for variable 1 of the solution
+p1=plot(sol,idxs=(0,2)) # feature not available in v1.0.1
+# Generate a plot object for variable 1 in function of variable 1
+p1=plot(sol,idxs=(1,2)) # feature not available in v1.0.1
+# Generate a 3D plot object for variables 1 and 2 in function of time
+p1=plot(sol,idxs=(0,1,2)) # feature not available in v1.0.1
+
#plot variables 1 and 2 of the solution
diff --git a/docs/src/interface.md b/docs/src/interface.md
index e493440..50236f5 100644
--- a/docs/src/interface.md
+++ b/docs/src/interface.md
@@ -5,11 +5,12 @@
NLodeProblem(odeExprs)
```
```@docs
-solve(prob::NLODEProblem{PRTYPE,T,Z,D,CS},al::QSSAlgorithm{SolverType, OrderType},tspan::Tuple{Float64, Float64};sparsity::Val{Sparsity}=Val(false),saveat=1e-9::Float64,abstol=1e-4::Float64,reltol=1e-3::Float64,maxErr=Inf::Float64,maxStepsAllowed=10000000) where{PRTYPE,T,Z,D,CS,SolverType,OrderType,Sparsity}
+solve(prob::NLODEProblem{PRTYPE,T,Z,D,CS},al::QSSAlgorithm{SolverType, OrderType},tspan::Tuple{Float64, Float64};sparsity::Val{Sparsity}=Val(false),saveat=1e-9::Float64,abstol=1e-4::Float64,reltol=1e-3::Float64,maxErr=Inf::Float64,maxiters=10000000) where{PRTYPE,T,Z,D,CS,SolverType,OrderType,Sparsity}
```
```@docs
-plot_Sol(sol::Sol{T,O},xvars::Int...;note=" "::String,xlims=(0.0,0.0)::Tuple{Float64, Float64},ylims=(0.0,0.0)::Tuple{Float64, Float64},legend=:true::Bool,marker=:circle::Symbol) where{T,O}
+plot_Sol(sol::Sol{T,O},xvars::Int...;note=" "::String,xlims=(0.0,0.0)::Tuple{Float64, Float64},ylims=(0.0,0.0)::Tuple{Float64, Float64},legend=:true::Bool,marker=:circle::Symbol,title="") where{T,O}
```
+
```@docs
save_Sol(sol::Sol{T,O},xvars::Int...;note=" "::String,xlims=(0.0,0.0)::Tuple{Float64, Float64},ylims=(0.0,0.0)::Tuple{Float64, Float64},legend=:true::Bool) where{T,O}
```
@@ -19,4 +20,3 @@ getError(sol::Sol{T,O},index::Int,f::Function) where{T,O}
```@docs
getAverageErrorByRefs(sol::Sol{T,O},solRef::Vector{Any}) where{T,O}
```
-
diff --git a/docs/src/userGuide.md b/docs/src/userGuide.md
index 8d432ec..3406901 100644
--- a/docs/src/userGuide.md
+++ b/docs/src/userGuide.md
@@ -95,17 +95,24 @@ abstol and reltol are the absolute and relative tolerances for the solver, contr
After solving the problem, we can query the solution to extract useful information such as variable values at specific times, the number of steps, events, and more.
```julia
+
+# returns a vector of the values of all variables at time 0.0005 # feature not available in v1.0.1
+sol(0.0005)
+
# Get the value of variable 2 at time 0.0005
value_at_time = sol(2, 0.0005)
+# Get the value of variable 2 at time 0.0005
+value_at_time = sol(0.0005,idxs=2)# feature not available in v1.0.1
+
# Get the total number of steps to end the simulation
-total_steps = sol.totalSteps
+total_steps = sol.totalSteps # feature available only in v1.0.1 ->sol.stats in later versions
# Get the number of simultaneous steps during the simulation
-simul_step_count = sol.simulStepCount
+simul_step_count = sol.simulStepCount # feature available only in v1.0.1 ->sol.stats in later versions
# Get the total number of events during the simulation
-event_count = sol.evCount
+event_count = sol.evCount # feature available only in v1.0.1 ->sol.stats in later versions
# Get the saved times and variables
saved_times = sol.savedTimes
@@ -134,6 +141,15 @@ plot_obj = plot_Sol(sol)
plot_Sol(sol,1)
# Display the plot
display(plot_obj)
+# Generate a plot object for all variables of the solution. It requires "using Plots". feature not available in v1.0.1.
+plot(sol)
+# Generate a plot object for variable 1 of the solution
+p1=plot(sol,idxs=(0,2)) # feature not available in v1.0.1
+# Generate a plot object for variable 1 in function of variable 1
+p1=plot(sol,idxs=(1,2)) # feature not available in v1.0.1
+# Generate a 3D plot object for variables 1 and 2 in function of time
+p1=plot(sol,idxs=(0,1,2)) # feature not available in v1.0.1
+
#plot variables 1 and 2 of the solution
diff --git a/lcov.info b/lcov.info
index 79f15e5..718ca41 100644
--- a/lcov.info
+++ b/lcov.info
@@ -1,230 +1,230 @@
SF:src\Common\Helper_QSSNLDiscreteProblem.jl
-DA:2,24
-DA:9,4
-DA:10,4
-DA:12,4
-DA:13,84
-DA:14,7
-DA:15,7
-DA:16,70
+DA:2,36
+DA:9,5
+DA:10,5
+DA:12,5
+DA:13,114
+DA:14,12
+DA:15,12
+DA:16,90
DA:18,4
DA:19,4
DA:20,2
DA:22,4
DA:23,4
DA:24,4
-DA:26,73
-DA:29,4
-DA:30,8
-DA:31,7
-DA:32,7
-DA:33,14
-DA:34,7
-DA:35,7
-DA:36,14
-DA:37,8
-DA:42,4
-DA:43,4
-DA:44,4
-DA:45,109
-DA:46,9
-DA:47,9
-DA:48,91
-DA:49,6
-DA:50,2
-DA:51,2
-DA:52,0
-DA:54,2
-DA:55,2
-DA:57,6
-DA:59,94
-DA:61,4
-DA:62,8
-DA:63,18
-DA:64,9
-DA:65,9
-DA:66,9
-DA:67,9
-DA:68,9
-DA:69,18
-DA:70,4
-DA:73,12
-DA:74,12
-DA:76,12
-DA:77,147
-DA:78,5
-DA:79,5
-DA:80,5
-DA:81,1
-DA:83,5
-DA:84,5
-DA:85,137
-DA:87,12
-DA:88,12
-DA:89,6
-DA:91,12
-DA:92,12
-DA:94,130
-DA:96,12
-DA:100,6
-DA:101,6
-DA:102,6
-DA:103,6
-DA:104,6
-DA:105,12
-DA:106,24
-DA:107,24
-DA:108,42
-DA:109,12
-DA:110,24
-DA:111,24
-DA:112,24
-DA:113,38
-DA:114,14
-DA:115,20
-DA:116,20
-DA:117,28
-DA:119,14
-DA:120,14
-DA:121,18
-DA:123,28
-DA:124,9
-DA:125,18
-DA:126,24
-DA:127,24
-DA:128,24
-DA:129,42
-DA:130,6
-DA:133,6
-DA:134,6
-DA:135,6
-DA:136,6
-DA:137,12
-DA:138,24
-DA:139,24
-DA:140,42
-DA:141,12
-DA:142,24
-DA:143,24
-DA:144,24
-DA:145,42
-DA:146,22
-DA:147,39
-DA:148,61
-DA:149,22
-DA:150,22
+DA:26,98
+DA:29,5
+DA:30,10
+DA:31,12
+DA:32,12
+DA:33,24
+DA:34,12
+DA:35,12
+DA:36,24
+DA:37,10
+DA:42,7
+DA:43,7
+DA:44,7
+DA:45,279
+DA:46,21
+DA:47,21
+DA:48,237
+DA:49,17
+DA:50,13
+DA:51,13
+DA:52,9
+DA:54,13
+DA:55,13
+DA:57,17
+DA:59,241
+DA:61,7
+DA:62,14
+DA:63,37
+DA:64,19
+DA:65,19
+DA:66,19
+DA:67,19
+DA:68,19
+DA:69,38
+DA:70,7
+DA:73,18
+DA:74,18
+DA:76,18
+DA:77,540
+DA:78,31
+DA:79,31
+DA:80,31
+DA:81,13
+DA:83,31
+DA:84,31
+DA:85,478
+DA:87,50
+DA:88,50
+DA:89,39
+DA:91,50
+DA:92,50
+DA:94,459
+DA:96,18
+DA:100,7
+DA:101,7
+DA:102,7
+DA:103,7
+DA:104,7
+DA:105,14
+DA:106,36
+DA:107,36
+DA:108,65
+DA:109,14
+DA:110,36
+DA:111,36
+DA:112,36
+DA:113,52
+DA:114,34
+DA:115,50
+DA:116,140
+DA:117,158
+DA:119,34
+DA:120,34
+DA:121,58
+DA:123,68
+DA:124,83
+DA:125,166
+DA:126,54
+DA:127,36
+DA:128,36
+DA:129,65
+DA:130,7
+DA:133,7
+DA:134,7
+DA:135,7
+DA:136,7
+DA:137,14
+DA:138,36
+DA:139,36
+DA:140,65
+DA:141,14
+DA:142,36
+DA:143,36
+DA:144,36
+DA:145,65
+DA:146,21
+DA:147,36
+DA:148,57
+DA:149,21
+DA:150,21
DA:151,8
-DA:153,44
+DA:153,42
DA:154,4
DA:155,8
DA:156,28
-DA:157,24
-DA:158,24
-DA:159,42
-DA:160,6
-DA:164,12
-DA:165,12
-DA:166,12
-DA:167,24
-DA:168,48
-DA:169,84
-DA:170,24
-DA:171,48
-DA:172,83
-DA:173,27
-DA:174,40
-DA:175,87
-DA:176,24
-DA:177,33
-DA:178,48
-DA:179,84
-DA:180,12
-LH:139
+DA:157,36
+DA:158,36
+DA:159,65
+DA:160,7
+DA:164,14
+DA:165,14
+DA:166,14
+DA:167,28
+DA:168,72
+DA:169,130
+DA:170,28
+DA:171,72
+DA:172,111
+DA:173,189
+DA:174,222
+DA:175,134
+DA:176,25
+DA:177,35
+DA:178,72
+DA:179,130
+DA:180,14
+LH:140
LF:140
end_of_record
SF:src\Common\Helper_QSSNLProblem.jl
-DA:1,112
-DA:2,112
-DA:3,634
-DA:4,110
-DA:5,110
-DA:6,110
-DA:7,110
-DA:9,497
-DA:11,112
-DA:13,97
-DA:14,97
-DA:15,24
-DA:16,3
-DA:17,18
-DA:18,7
+DA:1,149
+DA:2,149
+DA:3,872
+DA:4,109
+DA:5,109
+DA:6,109
+DA:7,109
+DA:9,688
+DA:11,149
+DA:13,125
+DA:14,125
+DA:15,40
+DA:16,7
+DA:17,26
+DA:18,11
DA:19,4
DA:20,2
DA:21,0
-DA:23,12
-DA:26,85
-DA:28,97
-DA:30,151
-DA:31,151
-DA:32,34
-DA:33,4
-DA:34,26
-DA:35,11
+DA:23,20
+DA:26,105
+DA:28,125
+DA:30,180
+DA:31,180
+DA:32,58
+DA:33,10
+DA:34,38
+DA:35,17
DA:36,4
-DA:37,19
-DA:43,134
-DA:45,151
-DA:47,11
-DA:48,11
+DA:37,31
+DA:43,151
+DA:45,180
+DA:47,16
+DA:48,16
DA:49,6
DA:50,1
DA:51,4
DA:52,1
DA:53,2
DA:54,4
-DA:60,8
-DA:62,11
-DA:64,66
-DA:65,66
-DA:66,557
-DA:67,99
-DA:68,99
-DA:69,327
+DA:60,13
+DA:62,16
+DA:64,81
+DA:65,81
+DA:66,698
+DA:67,125
+DA:68,125
+DA:69,412
DA:70,0
-DA:72,426
-DA:74,66
-DA:76,43
-DA:77,43
-DA:78,677
-DA:79,411
-DA:80,40
-DA:81,267
-DA:82,85
-DA:83,182
-DA:84,11
-DA:85,171
-DA:86,324
-DA:90,480
-DA:91,85
-DA:92,170
-DA:93,310
-DA:94,11
-DA:95,11
-DA:97,677
-DA:99,43
-DA:101,25
-DA:102,25
-DA:103,260
-DA:104,122
-DA:105,41
-DA:106,61
-DA:107,25
-DA:108,36
+DA:72,537
+DA:74,81
+DA:76,65
+DA:77,65
+DA:78,855
+DA:79,615
+DA:80,100
+DA:81,336
+DA:82,99
+DA:83,237
+DA:84,16
+DA:85,221
+DA:86,429
+DA:90,578
+DA:91,99
+DA:92,198
+DA:93,380
+DA:94,16
+DA:95,16
+DA:97,855
+DA:99,65
+DA:101,42
+DA:102,42
+DA:103,582
+DA:104,452
+DA:105,174
+DA:106,154
+DA:107,50
+DA:108,104
DA:109,7
-DA:110,29
+DA:110,97
DA:111,0
-DA:114,260
-DA:116,25
+DA:114,582
+DA:116,42
DA:120,26
DA:121,26
DA:122,26
@@ -366,70 +366,70 @@ DA:206,12
DA:207,72
DA:208,12
DA:209,12
-DA:212,6
-DA:213,6
-DA:214,12
-DA:215,16
-DA:216,7
-DA:217,9
-DA:218,7
-DA:219,9
-DA:220,9
-DA:221,18
-DA:222,9
-DA:224,26
-DA:225,6
-DA:226,6
-DA:227,6
-DA:228,6
-DA:229,6
-DA:230,6
-DA:231,36
-DA:232,6
-DA:233,6
-DA:236,18
-DA:237,18
-DA:238,18
-DA:239,1051
-DA:240,2084
-DA:241,36
-DA:242,34
-DA:243,29
-DA:244,5
-DA:245,5
-DA:246,1014
-DA:247,2028
-DA:248,4040
-DA:249,15150
-DA:250,3030
-DA:252,0
-DA:254,6060
-DA:255,1014
-DA:256,1014
-DA:258,50
-DA:259,18
-DA:261,18
-DA:262,18
-DA:263,18
-DA:264,1051
-DA:265,2084
-DA:267,36
-DA:268,34
-DA:269,58
-DA:270,54
-DA:271,108
-DA:272,5
-DA:273,5
-DA:274,2028
-DA:275,4040
-DA:276,15150
-DA:277,3030
-DA:279,0
-DA:281,6060
-DA:283,1014
-DA:285,50
-DA:286,18
-LH:163
+DA:212,7
+DA:213,7
+DA:214,14
+DA:215,31
+DA:216,12
+DA:217,14
+DA:218,12
+DA:219,19
+DA:220,19
+DA:221,37
+DA:222,19
+DA:224,55
+DA:225,7
+DA:226,7
+DA:227,7
+DA:228,7
+DA:229,7
+DA:230,7
+DA:231,42
+DA:232,7
+DA:233,7
+DA:236,19
+DA:237,19
+DA:238,19
+DA:239,1064
+DA:240,2109
+DA:241,38
+DA:242,38
+DA:243,30
+DA:244,8
+DA:245,8
+DA:246,1026
+DA:247,2052
+DA:248,4096
+DA:249,15330
+DA:250,3066
+DA:252,4
+DA:254,6140
+DA:255,1026
+DA:256,1026
+DA:258,57
+DA:259,19
+DA:261,19
+DA:262,19
+DA:263,19
+DA:264,1064
+DA:265,2109
+DA:267,38
+DA:268,38
+DA:269,60
+DA:270,59
+DA:271,118
+DA:272,8
+DA:273,8
+DA:274,2052
+DA:275,4096
+DA:276,15330
+DA:277,3066
+DA:279,4
+DA:281,6140
+DA:283,1026
+DA:285,57
+DA:286,19
+LH:165
LF:167
end_of_record
SF:src\Common\QSSNL_AbstractTypes.jl
@@ -437,293 +437,323 @@ LH:0
LF:0
end_of_record
SF:src\Common\QSSNLdiscrProblem.jl
-DA:4,12
-DA:22,53
-DA:23,53
-DA:29,6
-DA:30,6
-DA:32,6
-DA:33,6
-DA:34,6
-DA:37,6
-DA:38,6
-DA:39,6
-DA:40,6
-DA:41,6
-DA:42,6
-DA:45,6
-DA:46,6
-DA:47,6
-DA:48,6
-DA:49,6
-DA:50,96
-DA:51,4
-DA:53,60
-DA:54,16
-DA:55,8
-DA:56,16
+DA:4,14
+DA:22,83
+DA:23,83
+DA:29,7
+DA:30,7
+DA:32,7
+DA:33,7
+DA:34,7
+DA:37,7
+DA:38,7
+DA:39,7
+DA:40,7
+DA:41,7
+DA:42,7
+DA:45,7
+DA:46,7
+DA:47,7
+DA:48,7
+DA:49,7
+DA:50,173
+DA:51,5
+DA:53,104
+DA:54,18
+DA:55,9
+DA:56,18
DA:57,4
-DA:58,4
+DA:58,5
DA:59,0
DA:61,0
-DA:63,4
-DA:64,4
-DA:65,4
-DA:66,46
-DA:67,44
-DA:68,60
+DA:63,5
+DA:64,5
+DA:65,5
+DA:66,87
+DA:67,85
+DA:68,106
DA:70,0
-DA:71,4
-DA:73,4
-DA:74,4
-DA:75,4
-DA:76,3
-DA:78,4
-DA:79,40
-DA:80,12
-DA:81,12
-DA:82,12
-DA:83,24
+DA:71,7
+DA:73,7
+DA:74,7
+DA:75,7
+DA:76,4
+DA:78,7
+DA:79,78
+DA:80,18
+DA:81,18
+DA:82,18
+DA:83,36
DA:84,0
-DA:87,12
-DA:88,12
-DA:89,1
-DA:91,12
-DA:97,12
-DA:98,11
-DA:99,11
-DA:100,11
-DA:103,12
-DA:104,12
-DA:106,12
-DA:107,1
-DA:108,1
-DA:110,11
-DA:113,12
-DA:114,12
-DA:116,12
-DA:117,12
-DA:119,12
-DA:120,12
-DA:121,12
-DA:122,12
-DA:124,68
-DA:125,68
-DA:127,46
-DA:128,68
-DA:129,22
-DA:131,12
-DA:134,34
-DA:135,43
+DA:87,18
+DA:88,18
+DA:89,2
+DA:91,18
+DA:97,18
+DA:98,12
+DA:99,12
+DA:100,12
+DA:103,18
+DA:104,18
+DA:106,18
+DA:107,6
+DA:108,6
+DA:110,12
+DA:113,18
+DA:114,18
+DA:116,18
+DA:117,18
+DA:119,18
+DA:120,18
+DA:121,18
+DA:122,18
+DA:124,88
+DA:125,88
+DA:127,68
+DA:128,88
+DA:129,20
+DA:131,24
+DA:134,44
+DA:135,62
DA:136,0
-DA:138,41
-DA:142,34
-DA:144,12
-DA:145,12
-DA:146,12
-DA:147,24
-DA:148,1
-DA:150,4
-DA:152,4
-DA:153,4
-DA:154,0
-DA:156,2
-DA:158,2
-DA:159,8
-DA:160,0
-DA:162,6
-DA:165,2
-DA:167,12
-DA:168,12
-DA:169,12
-DA:170,12
-DA:173,62
-DA:174,6
-DA:176,6
-DA:177,12
-DA:178,16
-DA:179,24
-DA:180,12
-DA:181,16
-DA:183,12
-DA:184,4
-DA:186,18
-DA:187,6
-DA:188,6
-DA:189,6
-DA:191,6
-DA:194,6
-DA:195,10
-DA:196,6
-DA:197,8
-DA:199,6
-DA:200,6
-DA:201,6
-DA:204,6
-DA:205,6
-DA:206,12
-DA:207,18
-DA:208,30
-DA:210,6
-DA:211,6
-DA:212,6
-DA:214,6
-DA:216,12
-DA:217,6
-DA:220,6
-DA:221,6
-DA:222,6
-DA:223,6
-DA:224,42
-DA:225,6
-DA:227,6
-DA:228,6
-DA:229,6
-DA:230,6
-DA:231,6
-DA:232,6
-DA:239,6
-DA:244,6
-DA:246,6
-DA:247,6
-DA:250,6
-DA:251,6
-DA:252,6
-DA:253,6
-DA:256,6
-DA:257,6
-DA:258,6
-DA:259,13
-DA:260,22
-DA:261,9
-DA:262,4
-DA:263,8
-DA:264,4
+DA:138,58
+DA:142,44
+DA:144,18
+DA:145,18
+DA:146,18
+DA:147,36
+DA:148,6
+DA:150,22
+DA:152,21
+DA:153,22
+DA:154,1
+DA:156,10
+DA:158,11
+DA:159,44
+DA:160,1
+DA:162,33
+DA:165,11
+DA:167,18
+DA:168,18
+DA:169,18
+DA:170,18
+DA:173,106
+DA:174,7
+DA:176,7
+DA:177,14
+DA:178,23
+DA:179,32
+DA:180,16
+DA:181,18
+DA:183,16
+DA:184,7
+DA:186,25
+DA:187,7
+DA:188,7
+DA:189,7
+DA:191,7
+DA:194,7
+DA:195,12
+DA:196,11
+DA:197,17
+DA:199,7
+DA:200,7
+DA:201,7
+DA:204,7
+DA:205,7
+DA:206,14
+DA:207,29
+DA:208,51
+DA:210,7
+DA:211,7
+DA:212,7
+DA:214,7
+DA:216,14
+DA:217,7
+DA:220,7
+DA:221,7
+DA:222,7
+DA:223,7
+DA:224,49
+DA:225,7
+DA:227,7
+DA:228,7
+DA:229,7
+DA:230,7
+DA:231,7
+DA:232,7
+DA:239,7
+DA:244,7
+DA:246,7
+DA:247,7
+DA:250,7
+DA:251,7
+DA:252,7
+DA:253,7
+DA:256,7
+DA:257,7
+DA:258,7
+DA:259,18
+DA:260,31
+DA:261,11
+DA:262,6
+DA:263,12
+DA:264,10
DA:265,2
-DA:266,2
-DA:267,2
-DA:268,8
-DA:269,8
-DA:271,8
-DA:272,4
-DA:273,5
-DA:274,6
-DA:277,6
-DA:278,6
-DA:279,6
-DA:280,25
-DA:281,44
-DA:282,9
-DA:283,4
-DA:284,5
-DA:285,6
-LH:173
+DA:266,8
+DA:267,8
+DA:268,32
+DA:269,32
+DA:271,20
+DA:272,6
+DA:273,8
+DA:274,7
+DA:277,7
+DA:278,7
+DA:279,7
+DA:280,38
+DA:281,69
+DA:282,11
+DA:283,18
+DA:284,32
+DA:285,7
+LH:175
LF:180
end_of_record
SF:src\Common\QSS_Algorithm.jl
-DA:4,96
+DA:4,104
DA:11,6
-DA:16,11
+DA:16,12
DA:25,7
-DA:30,11
+DA:30,12
DA:48,6
-DA:52,6
+DA:52,8
LH:7
LF:7
end_of_record
SF:src\Common\QSS_data.jl
-DA:4,96
-DA:31,62
+DA:4,104
+DA:31,68
LH:2
LF:2
end_of_record
SF:src\Common\Scheduler.jl
-DA:2,1156820
-DA:3,1156820
-DA:4,1156820
-DA:5,1156820
-DA:6,1156820
-DA:7,1156820
-DA:8,1156820
-DA:10,1156820
-DA:11,1156820
-DA:12,582033679
-DA:13,7197296
-DA:14,7197296
-DA:16,582033679
-DA:17,108412
-DA:18,108412
-DA:20,116291053
-DA:21,1156820
-DA:22,880802
-DA:23,460662
-DA:24,460662
-DA:26,1442129
-DA:28,1156820
-DA:29,1928
-DA:30,1757
-DA:32,2099
-DA:35,1154892
-DA:36,22350
-DA:38,1132542
-DA:41,1156820
+DA:2,1270205
+DA:3,1270205
+DA:4,1270205
+DA:5,1270205
+DA:6,1270205
+DA:7,1270205
+DA:8,1270205
+DA:10,1270205
+DA:11,1270205
+DA:12,583373620
+DA:13,7523932
+DA:14,7523932
+DA:16,583373620
+DA:17,125614
+DA:18,125614
+DA:20,116547703
+DA:21,1270205
+DA:22,1495420
+DA:23,698187
+DA:24,698187
+DA:26,2557980
+DA:28,1270205
+DA:29,44357
+DA:30,7
+DA:32,88707
+DA:35,1225848
+DA:36,24249
+DA:38,1201599
+DA:41,1270205
DA:42,0
-DA:45,1156820
+DA:45,1270205
LH:30
LF:31
end_of_record
SF:src\Common\Solution.jl
-DA:17,110
-DA:31,55
-DA:33,55
-DA:35,241341233
-DA:36,241341233
-DA:37,235856305
-DA:38,5484928
-DA:39,5484928
-DA:41,0
-DA:45,2002008
-DA:46,2002008
-DA:48,2002008
-DA:50,3872016
-DA:51,116056932
-DA:52,6961732
-DA:53,1740433
-DA:54,1740433
-DA:55,1740433
-DA:57,1740433
-DA:58,114316499
-DA:59,0
-DA:61,0
-DA:63,228503423
-DA:65,261575
-DA:67,7
-DA:68,7
-DA:69,7
-DA:70,7
-DA:71,7
-DA:72,301507
-DA:73,301500
-DA:74,301500
-DA:75,301500
-DA:76,7
-DA:77,7
-DA:78,7
-DA:79,7
-DA:80,7
+DA:2,52
+DA:24,132
+DA:36,66
+DA:38,66
+DA:40,241627180
+DA:41,241627180
+DA:42,236140046
+DA:45,5487134
+DA:51,2002743
+DA:52,2002743
+DA:53,2002743
+DA:54,1
+DA:55,1
+DA:56,2
+DA:57,3
+DA:58,1
+DA:61,2002742
+DA:63,3873484
+DA:64,116198073
+DA:65,6964628
+DA:66,1741157
+DA:67,1741157
+DA:68,1741157
+DA:70,1741157
+DA:71,114456916
+DA:72,3
+DA:74,3
+DA:76,228784244
+DA:78,261582
+DA:81,7
DA:82,7
-DA:83,1027
-DA:84,1027
-DA:85,2054
-DA:86,116301490
-DA:87,4002973
-DA:88,1027
-DA:89,1027
-DA:90,2047
+DA:83,7
+DA:84,7
+DA:85,7
+DA:86,301507
+DA:87,301500
+DA:88,301500
+DA:89,301500
+DA:90,7
DA:91,7
-DA:93,8
-LH:46
-LF:49
+DA:92,7
+DA:94,7
+DA:96,7
+DA:97,1027
+DA:98,1027
+DA:99,2054
+DA:100,2002000
+DA:101,4002973
+DA:102,1027
+DA:103,1027
+DA:104,2047
+DA:105,7
+DA:107,7
+DA:108,7
+DA:109,7
+DA:110,7
+DA:111,7
+DA:112,707
+DA:113,700
+DA:114,700
+DA:115,700
+DA:116,7
+DA:117,7
+DA:120,7
+DA:123,7
+DA:124,7
+DA:125,14
+DA:126,700
+DA:127,1393
+DA:128,7
+DA:129,7
+DA:131,7
+DA:133,39
+DA:134,4
+DA:138,0
+DA:139,0
+DA:140,0
+DA:141,0
+DA:142,0
+LH:74
+LF:79
end_of_record
SF:src\Common\SolutionError.jl
DA:6,2
@@ -754,8 +784,6 @@ DA:34,1
DA:35,1
DA:37,0
DA:39,1
-DA:42,0
-DA:43,0
DA:45,501000
DA:46,501000
DA:54,1
@@ -778,77 +806,121 @@ DA:73,1000
DA:74,1999
DA:75,1
LH:47
-LF:51
+LF:49
end_of_record
SF:src\Common\SolutionPlot.jl
-DA:5,14
-DA:6,7
-DA:7,7
-DA:8,6
-DA:9,14
-DA:10,6
-DA:11,6
-DA:12,8
-DA:13,5
-DA:14,5
-DA:15,3
-DA:16,1
-DA:17,1
-DA:18,2
-DA:19,1
-DA:20,1
-DA:22,1
-DA:23,1
-DA:25,14
-DA:26,19
+DA:13,12
+DA:14,6
+DA:15,6
+DA:17,6
+DA:18,6
+DA:19,5
+DA:20,10
+DA:21,5
+DA:22,5
+DA:23,5
+DA:24,4
+DA:25,4
+DA:27,1
DA:28,1
-DA:29,2
-DA:30,1
-DA:32,1
-DA:34,2
-DA:35,3
-DA:37,9
-DA:38,1
-DA:39,7
-DA:40,1
-DA:41,10
-DA:42,1
-DA:44,4
-DA:46,7
-DA:59,2
-DA:60,1
-DA:61,1
-DA:62,1
-DA:63,1
-DA:65,2
-DA:66,1
-DA:67,1
-DA:68,1
-DA:69,1
-DA:74,10
-DA:75,5
-DA:76,5
-DA:77,5
-DA:78,5
-DA:79,5
-DA:80,5
-DA:81,10
-DA:82,600020
-DA:83,300010
-DA:84,300010
-DA:85,600015
-DA:86,5
+DA:30,10
+DA:31,14
+DA:33,1
+DA:39,2
+DA:40,3
+DA:42,8
+DA:43,1
+DA:44,6
+DA:45,1
+DA:46,8
+DA:47,1
+DA:49,3
+DA:51,6
+DA:59,12
+DA:60,6
+DA:61,6
+DA:64,6
+DA:65,6
+DA:66,2
+DA:67,2
+DA:68,8
+DA:69,2
+DA:70,2
+DA:71,6
+DA:72,2
+DA:73,2
+DA:74,4
+DA:75,2
+DA:76,2
+DA:77,2
+DA:78,1
+DA:79,1
+DA:81,1
+DA:82,1
+DA:84,8
+DA:85,10
+DA:87,0
DA:88,0
-DA:90,7
-DA:91,1
-DA:92,5
-DA:93,1
-DA:94,6
+DA:89,0
+DA:92,4
+DA:93,2
+DA:94,1
DA:95,1
-DA:97,2
-DA:99,5
-LH:65
-LF:66
+DA:97,1
+DA:98,1
+DA:99,3
+DA:101,2
+DA:102,2
+DA:104,1
+DA:105,1
+DA:106,1
+DA:108,1
+DA:109,1
+DA:110,1
+DA:111,3
+DA:114,0
+DA:116,6
+DA:117,0
+DA:118,6
+DA:119,0
+DA:120,12
+DA:121,0
+DA:123,6
+DA:125,6
+DA:138,2
+DA:139,1
+DA:140,1
+DA:141,1
+DA:142,1
+DA:144,2
+DA:145,1
+DA:146,1
+DA:147,1
+DA:148,1
+DA:153,10
+DA:154,5
+DA:155,5
+DA:156,5
+DA:157,5
+DA:158,5
+DA:159,5
+DA:160,10
+DA:161,600020
+DA:162,300010
+DA:163,300010
+DA:164,600015
+DA:165,5
+DA:167,0
+DA:169,7
+DA:170,1
+DA:171,5
+DA:172,1
+DA:173,6
+DA:174,1
+DA:176,2
+DA:178,5
+LH:102
+LF:110
end_of_record
SF:src\Common\TaylorEquationConstruction.jl
DA:5,11
@@ -857,20 +929,20 @@ DA:7,11
DA:8,11
DA:9,11
DA:10,11
-DA:13,42
-DA:14,42
-DA:15,42
-DA:17,1369
-DA:18,48
-DA:19,2
-DA:20,2
-DA:21,2
-DA:22,2
-DA:23,2
-DA:24,2
-DA:25,2
-DA:26,2
-DA:27,44
+DA:13,52
+DA:14,52
+DA:15,52
+DA:17,2166
+DA:18,84
+DA:19,4
+DA:20,4
+DA:21,4
+DA:22,4
+DA:23,4
+DA:24,4
+DA:25,4
+DA:26,4
+DA:27,76
DA:28,3
DA:29,3
DA:30,3
@@ -879,63 +951,63 @@ DA:32,3
DA:33,3
DA:34,3
DA:36,3
-DA:37,38
-DA:38,7
-DA:39,7
-DA:40,7
-DA:41,7
-DA:42,7
-DA:43,7
-DA:44,7
-DA:45,7
-DA:48,23
-DA:49,23
-DA:50,23
-DA:51,23
-DA:52,58
-DA:54,1299
-DA:55,6
-DA:56,6
-DA:57,6
-DA:58,6
-DA:59,6
-DA:61,1287
-DA:62,32
-DA:63,10
-DA:64,10
-DA:65,10
-DA:66,10
-DA:67,10
-DA:68,10
-DA:69,10
-DA:71,10
-DA:72,12
-DA:73,2
-DA:74,2
-DA:75,2
-DA:76,2
-DA:77,2
-DA:78,2
-DA:79,2
-DA:80,2
-DA:82,8
-DA:83,8
-DA:84,8
-DA:85,8
-DA:86,28
-DA:88,1247
-DA:89,5
-DA:90,5
-DA:91,5
-DA:92,5
-DA:93,5
-DA:96,1237
-DA:97,32
-DA:98,32
-DA:99,32
-DA:100,32
-DA:101,32
-DA:102,1173
+DA:37,70
+DA:38,18
+DA:39,18
+DA:40,18
+DA:41,18
+DA:42,18
+DA:43,18
+DA:44,18
+DA:45,18
+DA:48,30
+DA:49,30
+DA:50,30
+DA:51,30
+DA:52,85
+DA:54,2056
+DA:55,7
+DA:56,7
+DA:57,7
+DA:58,7
+DA:59,7
+DA:61,2042
+DA:62,64
+DA:63,11
+DA:64,11
+DA:65,11
+DA:66,11
+DA:67,11
+DA:68,11
+DA:69,11
+DA:71,11
+DA:72,42
+DA:73,6
+DA:74,6
+DA:75,6
+DA:76,6
+DA:77,6
+DA:78,6
+DA:79,6
+DA:80,6
+DA:82,30
+DA:83,30
+DA:84,30
+DA:85,30
+DA:86,77
+DA:88,1948
+DA:89,6
+DA:90,6
+DA:91,6
+DA:92,6
+DA:93,6
+DA:96,1936
+DA:97,43
+DA:98,43
+DA:99,43
+DA:100,43
+DA:101,43
+DA:102,1850
DA:103,12
DA:104,12
DA:105,12
@@ -945,24 +1017,24 @@ DA:108,12
DA:109,12
DA:110,12
DA:111,12
-DA:112,1149
-DA:113,8
-DA:114,8
-DA:115,8
-DA:116,8
-DA:117,8
-DA:118,1132
+DA:112,1826
+DA:113,24
+DA:114,24
+DA:115,24
+DA:116,24
+DA:117,24
+DA:118,1777
DA:119,2
DA:120,2
DA:121,2
DA:122,2
-DA:124,1129
+DA:124,1774
DA:125,3
DA:126,3
DA:127,3
DA:128,3
DA:129,3
-DA:130,1123
+DA:130,1768
DA:131,2
DA:132,2
DA:133,2
@@ -971,7 +1043,7 @@ DA:136,2
DA:137,2
DA:138,2
DA:139,2
-DA:141,1118
+DA:141,1763
DA:143,1
DA:144,1
DA:145,1
@@ -984,179 +1056,176 @@ DA:153,1
DA:154,1
DA:155,1
DA:156,1
-DA:157,1102
+DA:157,1739
DA:159,0
-DA:160,1228
-DA:164,42
-DA:166,42
+DA:160,1941
+DA:164,52
+DA:166,52
LH:137
LF:138
end_of_record
SF:src\Interface\QSS_Solve.jl
DA:2,2
DA:3,1
-DA:18,96
-DA:19,48
-DA:22,48
-DA:24,48
-DA:25,48
-DA:26,48
-DA:28,48
-DA:29,17
-DA:31,2081
-DA:32,31
-DA:33,31
-DA:43,96
-DA:44,5588890
-DA:45,5588794
-DA:47,96
-DA:51,48
-DA:52,48
-DA:53,2121
-DA:54,48
-DA:55,48
-DA:56,48
-DA:57,48
-DA:58,2121
-DA:59,2121
-DA:60,2121
-DA:61,2121
-DA:62,48
-DA:63,125
-DA:64,48
-DA:65,48
-DA:66,125
-DA:67,65
-DA:68,48
-DA:69,25
-DA:70,42
-DA:71,48
-DA:72,2121
-DA:73,2121
-DA:74,5321
-DA:75,2121
-DA:76,5321
-DA:77,2121
-DA:78,2121
-DA:79,2121
-DA:80,2121
-DA:81,4194
-DA:82,48
-DA:83,48
-DA:84,463
-DA:85,306
-DA:86,48
-DA:87,48
-DA:90,31
-DA:91,31
-DA:92,31
-DA:93,31
-DA:94,2081
-DA:95,2081
-DA:96,4131
-DA:97,31
-DA:98,31
+DA:18,104
+DA:19,52
+DA:22,52
+DA:24,52
+DA:25,52
+DA:26,52
+DA:28,52
+DA:29,18
+DA:31,2105
+DA:32,34
+DA:42,104
+DA:43,11690655
+DA:44,11690551
+DA:46,104
+DA:50,52
+DA:51,2151
+DA:52,52
+DA:53,52
+DA:54,52
+DA:55,52
+DA:56,2151
+DA:57,2151
+DA:58,2151
+DA:59,2151
+DA:60,52
+DA:61,137
+DA:62,52
+DA:63,52
+DA:64,137
+DA:65,83
+DA:66,52
+DA:67,43
+DA:68,74
+DA:71,52
+DA:72,2151
+DA:73,2151
+DA:74,5411
+DA:75,2151
+DA:76,5411
+DA:77,2151
+DA:78,2151
+DA:79,2151
+DA:80,2151
+DA:81,4250
+DA:82,52
+DA:83,52
+DA:84,529
+DA:85,346
+DA:86,52
+DA:89,34
+DA:90,34
+DA:91,34
+DA:92,34
+DA:93,2105
+DA:94,2105
+DA:95,4176
+DA:96,34
+DA:97,34
+DA:100,2068
DA:101,2068
-DA:102,2068
-LH:65
-LF:65
+LH:62
+LF:62
end_of_record
SF:src\Interface\indexMacro.jl
-DA:6,18
-DA:7,18
-DA:8,18
-DA:9,18
-DA:10,18
-DA:11,18
-DA:12,18
-DA:13,18
-DA:14,18
-DA:15,18
-DA:16,18
-DA:17,1002
-DA:18,2
-DA:19,4
-DA:20,4
+DA:6,19
+DA:7,19
+DA:8,19
+DA:9,19
+DA:10,19
+DA:11,19
+DA:12,19
+DA:13,19
+DA:14,19
+DA:15,19
+DA:16,19
+DA:17,1015
+DA:18,3
+DA:19,6
+DA:20,5
DA:21,2
-DA:23,2
-DA:24,1000
-DA:25,1000
-DA:27,6
-DA:30,18
-DA:33,18
-DA:41,18
-DA:42,18
-DA:43,18
-DA:44,18
-DA:45,18
-DA:46,18
-DA:47,18
-DA:48,18
-DA:49,18
-DA:50,18
-DA:51,18
-DA:52,115
-DA:53,196
-DA:54,98
-DA:55,12
-DA:56,131
-DA:57,8
-DA:58,8
-DA:59,78
-DA:60,5
-DA:61,4
-DA:62,6
-DA:63,6
-DA:65,2
-DA:66,19
-DA:67,19
-DA:68,22
+DA:23,3
+DA:24,1013
+DA:25,1013
+DA:27,7
+DA:30,19
+DA:33,19
+DA:41,19
+DA:42,19
+DA:43,19
+DA:44,19
+DA:45,19
+DA:46,19
+DA:47,19
+DA:48,19
+DA:49,19
+DA:50,19
+DA:51,19
+DA:52,158
+DA:53,264
+DA:54,132
+DA:55,24
+DA:56,170
+DA:57,26
+DA:58,26
+DA:59,82
+DA:60,6
+DA:61,5
+DA:62,8
+DA:63,8
+DA:65,3
+DA:66,29
+DA:67,29
+DA:68,34
DA:69,2
DA:70,2
-DA:71,7
-DA:74,73
-DA:75,20
+DA:71,8
+DA:74,76
+DA:75,21
DA:76,16
DA:77,16
DA:78,16
DA:79,16
-DA:81,20
-DA:83,53
-DA:84,30
-DA:85,23
+DA:81,21
+DA:83,55
+DA:84,31
+DA:85,24
DA:86,5
-DA:87,98
+DA:87,132
DA:90,0
-DA:91,5
-DA:92,12
-DA:93,12
-DA:94,12
-DA:95,12
-DA:97,12
-DA:99,12
-DA:100,11
-DA:101,1
-DA:102,1
-DA:103,1
-DA:106,133
-DA:107,18
+DA:91,8
+DA:92,18
+DA:93,18
+DA:94,18
+DA:95,18
+DA:97,18
+DA:99,18
+DA:100,12
+DA:101,6
+DA:102,6
+DA:103,6
+DA:106,177
+DA:107,19
LH:77
LF:78
end_of_record
SF:src\QuantizedSystemSolver.jl
-DA:1,6095
+DA:1,6132
LH:1
LF:1
end_of_record
SF:src\Utils\rootfinders\SimUtils.jl
-DA:2,373096
-DA:3,373096
-DA:4,373096
+DA:2,379920
+DA:3,379920
+DA:4,379920
DA:5,4778
-DA:7,368318
-DA:9,373096
-DA:10,184159
-DA:13,373096
+DA:7,375142
+DA:9,379920
+DA:10,187571
+DA:13,379920
DA:15,72844
DA:16,72844
DA:17,72844
@@ -1165,75 +1234,75 @@ DA:20,25
DA:22,72844
DA:23,0
DA:25,72844
-DA:27,37874
-DA:28,37874
-DA:29,37874
-DA:30,37874
-DA:31,37874
-DA:32,75748
+DA:27,77933
+DA:28,77933
+DA:29,77933
+DA:30,77933
+DA:31,77933
+DA:32,155866
DA:33,0
DA:34,0
DA:36,0
DA:38,0
DA:39,0
-DA:43,37874
-DA:44,37874
+DA:43,77933
+DA:44,77933
DA:45,0
-DA:48,37874
-DA:49,37874
-DA:50,37874
-DA:51,21968
-DA:53,15906
-DA:55,37874
-DA:56,37874
-DA:57,15906
-DA:61,37874
-DA:63,1163883
-DA:64,1163883
-DA:65,1163883
-DA:66,1163883
-DA:67,1163883
-DA:68,2327282
-DA:69,492
-DA:70,278
-DA:72,214
-DA:74,492
-DA:75,103
-DA:79,1163391
-DA:80,1163391
-DA:81,517762
-DA:84,645629
-DA:85,645629
-DA:86,645629
-DA:87,331870
-DA:89,313759
-DA:91,645629
-DA:92,645629
-DA:93,321714
-DA:97,1163883
-DA:99,165333
-DA:100,165333
-DA:101,165333
-DA:102,165333
-DA:103,165333
-DA:104,165337
-DA:105,165329
-DA:106,1
-DA:108,165328
-DA:110,165329
-DA:111,82583
-DA:115,4
-DA:116,4
-DA:117,2
-DA:120,2
-DA:121,2
-DA:122,2
-DA:123,2
+DA:48,77933
+DA:49,77933
+DA:50,77933
+DA:51,30005
+DA:53,47928
+DA:55,77933
+DA:56,77933
+DA:57,47928
+DA:61,77933
+DA:63,2637054
+DA:64,2637054
+DA:65,2637054
+DA:66,2637054
+DA:67,2637054
+DA:68,5273617
+DA:69,499
+DA:70,264
+DA:72,235
+DA:74,499
+DA:75,115
+DA:79,2636555
+DA:80,2636555
+DA:81,906407
+DA:84,1730148
+DA:85,1730148
+DA:86,1730148
+DA:87,982503
+DA:89,747645
+DA:91,1730148
+DA:92,1730148
+DA:93,735593
+DA:97,2637054
+DA:99,404931
+DA:100,404931
+DA:101,404931
+DA:102,404931
+DA:103,404931
+DA:104,404941
+DA:105,404921
+DA:106,0
+DA:108,404921
+DA:110,404921
+DA:111,219182
+DA:115,10
+DA:116,10
+DA:117,5
+DA:120,5
+DA:121,5
+DA:122,5
+DA:123,5
DA:125,0
-DA:127,2
-DA:128,2
+DA:127,5
+DA:128,5
DA:129,0
-DA:133,165333
+DA:133,404931
DA:135,94032
DA:136,94032
DA:137,94032
@@ -1409,7 +1478,7 @@ DA:349,14
DA:350,3
DA:352,11
DA:354,47043
-LH:238
+LH:237
LF:260
end_of_record
SF:src\dense\NL_integrators\NL_LiQSS_Integrator.jl
@@ -1540,329 +1609,328 @@ DA:148,0
DA:149,0
DA:152,4865
DA:153,4865
-DA:154,4865
-DA:155,6027
-DA:156,4865
-DA:158,4865
-DA:159,4865
-DA:160,6027
-DA:161,4865
-DA:163,4865
-DA:164,9730
-DA:165,9730
-DA:166,0
-DA:167,0
-DA:169,14595
-DA:170,14595
-DA:171,7063
-DA:172,4865
-DA:173,4865
-DA:174,4865
-DA:176,111757
-DA:177,111757
-DA:178,111757
-DA:179,9
-LH:140
-LF:149
-end_of_record
-SF:src\dense\NL_integrators\NL_LiQSS_discreteIntegrator.jl
-DA:1,3
-DA:2,3
-DA:3,0
-DA:5,3
-DA:6,3
-DA:7,3
-DA:8,3
-DA:9,3
-DA:10,3
-DA:11,3
-DA:12,3
-DA:13,3
-DA:14,3
-DA:15,3
-DA:16,3
-DA:17,3
-DA:18,3
-DA:19,3
-DA:20,3
-DA:21,3
-DA:22,3
-DA:23,3
-DA:24,3
-DA:25,3
-DA:26,3
-DA:27,3
-DA:28,3
-DA:30,3
-DA:31,3
-DA:32,3
-DA:33,3
-DA:34,3
-DA:35,3
-DA:37,3
-DA:38,3
-DA:40,3
-DA:41,3
-DA:42,28
-DA:43,28
-DA:45,3
-DA:46,34
-DA:47,34
-DA:48,31
-DA:49,3
-DA:50,3
-DA:51,3
-DA:53,8
-DA:54,3
-DA:56,3
-DA:57,3
-DA:58,5
-DA:59,5
-DA:60,19
-DA:61,33
-DA:62,5
-DA:63,32
-DA:64,19
-DA:66,19
-DA:67,19
-DA:68,33
-DA:69,7
-DA:70,3
-DA:71,3
-DA:72,3
-DA:73,12
-DA:74,12
-DA:75,12
-DA:76,12
-DA:77,19
-DA:78,24
-DA:79,12
-DA:80,3
-DA:81,3
-DA:82,2
-DA:83,1
-DA:84,1
-DA:85,1
-DA:86,1
-DA:90,9
-DA:92,21
-DA:93,3
-DA:94,8
-DA:95,4
-DA:96,4
-DA:97,4
-DA:98,4
-DA:99,5
-DA:105,3
-DA:106,3
-DA:107,3
-DA:108,3
-DA:109,3
-DA:110,3
-DA:111,3
-DA:112,3
-DA:113,3
-DA:114,0
-DA:116,28003
-DA:117,28003
-DA:118,0
-DA:120,28003
-DA:121,28003
-DA:122,28003
-DA:123,28003
-DA:124,28003
-DA:126,3
-DA:128,28000
-DA:129,28000
-DA:131,28000
-DA:132,23453
-DA:133,23453
-DA:134,23453
-DA:135,23453
-DA:136,24776
-DA:137,23453
-DA:138,23453
-DA:139,23453
-DA:140,43419
-DA:141,46906
-DA:142,23453
-DA:143,23453
-DA:144,23453
-DA:145,40333
-DA:146,40333
-DA:147,20509
-DA:148,19632
-DA:150,40333
-DA:151,40333
-DA:152,20509
-DA:153,19632
-DA:155,79031
-DA:156,40333
-DA:157,40333
-DA:158,40333
-DA:159,63786
-DA:160,46906
+DA:154,4865
+DA:155,6027
+DA:156,4865
+DA:158,4865
+DA:159,4865
+DA:160,6027
+DA:161,4865
+DA:163,4865
+DA:164,9730
+DA:165,9730
+DA:166,0
+DA:167,0
+DA:169,14595
+DA:170,14595
+DA:171,7063
+DA:172,4865
+DA:173,4865
+DA:174,4865
+DA:176,111757
+DA:177,111757
+DA:178,111757
+DA:179,9
+DA:180,9
+LH:141
+LF:150
+end_of_record
+SF:src\dense\NL_integrators\NL_LiQSS_discreteIntegrator.jl
+DA:1,5
+DA:2,5
+DA:3,5
+DA:4,5
+DA:5,5
+DA:6,5
+DA:7,5
+DA:8,5
+DA:9,5
+DA:10,5
+DA:11,5
+DA:12,5
+DA:13,5
+DA:14,5
+DA:15,5
+DA:16,5
+DA:17,5
+DA:18,5
+DA:19,5
+DA:20,5
+DA:21,5
+DA:22,5
+DA:23,5
+DA:24,5
+DA:25,5
+DA:26,5
+DA:28,5
+DA:29,5
+DA:30,5
+DA:31,5
+DA:32,5
+DA:33,5
+DA:35,5
+DA:36,5
+DA:38,5
+DA:39,5
+DA:40,52
+DA:41,52
+DA:44,5
+DA:45,62
+DA:46,62
+DA:47,57
+DA:48,5
+DA:49,5
+DA:50,5
+DA:52,18
+DA:53,5
+DA:55,5
+DA:57,5
+DA:58,9
+DA:59,9
+DA:60,41
+DA:61,73
+DA:62,9
+DA:63,100
+DA:64,41
+DA:66,41
+DA:67,41
+DA:68,73
+DA:69,13
+DA:71,5
+DA:72,5
+DA:73,5
+DA:74,23
+DA:75,23
+DA:76,23
+DA:77,23
+DA:78,38
+DA:79,46
+DA:80,23
+DA:81,5
+DA:82,5
+DA:83,6
+DA:84,2
+DA:85,2
+DA:86,2
+DA:87,2
+DA:91,18
+DA:93,41
+DA:94,5
+DA:95,24
+DA:96,9
+DA:97,9
+DA:98,9
+DA:99,9
+DA:100,13
+DA:106,5
+DA:107,5
+DA:108,5
+DA:109,5
+DA:110,5
+DA:111,5
+DA:112,5
+DA:113,5
+DA:114,5
+DA:115,43460
+DA:116,43460
+DA:117,43460
+DA:118,43460
+DA:119,43460
+DA:120,43460
+DA:121,43460
+DA:123,5
+DA:125,43455
+DA:126,43455
+DA:128,43455
+DA:129,38744
+DA:130,38744
+DA:131,38744
+DA:132,38744
+DA:133,43062
+DA:134,38744
+DA:135,38744
+DA:136,38744
+DA:137,73061
+DA:138,77488
+DA:139,38744
+DA:140,38744
+DA:141,38744
+DA:142,69051
+DA:143,69051
+DA:144,37976
+DA:145,33342
+DA:147,69051
+DA:148,69051
+DA:149,37976
+DA:150,33342
+DA:152,142320
+DA:153,69051
+DA:154,69051
+DA:155,69051
+DA:156,107795
+DA:157,76471
+DA:158,1017
+DA:159,1958
+DA:160,1958
DA:161,0
DA:162,0
-DA:163,0
-DA:164,0
-DA:165,0
-DA:167,0
-DA:168,0
-DA:169,0
-DA:170,0
-DA:171,0
-DA:173,4547
-DA:174,4543
-DA:175,4543
-DA:176,4702
-DA:177,4543
-DA:178,4543
-DA:179,9085
-DA:180,9086
-DA:181,4543
-DA:182,4702
-DA:183,4861
-DA:184,4543
-DA:189,9089
-DA:190,4543
-DA:191,4543
-DA:192,4543
-DA:193,4543
-DA:194,30
-DA:195,15
-DA:196,15
-DA:197,15
-DA:198,15
-DA:201,4543
-DA:202,4543
-DA:203,4543
-DA:204,4702
-DA:205,4543
-DA:207,4543
-DA:208,4543
-DA:209,4702
-DA:210,4543
-DA:213,4543
-DA:214,9086
-DA:215,9086
-DA:216,13629
-DA:217,9089
-DA:218,4543
-DA:219,4543
-DA:220,4543
-DA:221,9086
-DA:222,4543
-DA:223,4543
-DA:224,4543
-DA:225,4543
-DA:226,4543
-DA:228,4543
-DA:229,9086
-DA:230,9086
-DA:231,0
-DA:232,0
-DA:234,13629
-DA:235,9089
-DA:236,4543
-DA:237,4543
-DA:238,4543
-DA:239,9086
-DA:240,9086
+DA:164,2975
+DA:165,3916
+DA:166,1017
+DA:167,1017
+DA:168,2034
+DA:170,4711
+DA:171,4688
+DA:172,4688
+DA:173,4992
+DA:174,4688
+DA:175,4688
+DA:176,9375
+DA:177,9376
+DA:178,4688
+DA:179,4992
+DA:180,5296
+DA:181,4688
+DA:186,9395
+DA:187,4688
+DA:188,4688
+DA:189,4688
+DA:190,4688
+DA:191,162
+DA:192,73
+DA:193,73
+DA:194,73
+DA:195,58
+DA:198,4688
+DA:199,4688
+DA:200,4688
+DA:201,4992
+DA:202,4688
+DA:204,4688
+DA:205,4688
+DA:206,4992
+DA:207,4688
+DA:210,4688
+DA:211,9384
+DA:212,9408
+DA:213,14072
+DA:214,9395
+DA:215,4688
+DA:216,4688
+DA:217,4688
+DA:218,9376
+DA:219,4688
+DA:220,4688
+DA:221,4688
+DA:222,4688
+DA:223,4688
+DA:225,4688
+DA:226,9384
+DA:227,9384
+DA:228,0
+DA:229,0
+DA:231,14072
+DA:232,9395
+DA:233,4688
+DA:234,4688
+DA:235,4688
+DA:236,9376
+DA:237,9239
+DA:238,137
+DA:239,137
+DA:240,137
DA:241,0
DA:242,0
-DA:243,0
-DA:244,0
-DA:245,0
-DA:247,0
-DA:248,0
-DA:249,0
-DA:250,0
-DA:251,0
-DA:255,8
-DA:256,0
-DA:257,0
-DA:258,0
-DA:259,0
-DA:261,0
-DA:263,8
-DA:264,4
-DA:265,4
-DA:266,4
-DA:267,0
-DA:268,0
-DA:271,4
-DA:272,0
-DA:273,0
-DA:275,4
-DA:276,4
+DA:244,274
+DA:245,274
+DA:246,137
+DA:247,137
+DA:248,137
+DA:252,29
+DA:253,17
+DA:254,17
+DA:255,34
+DA:256,17
+DA:258,34
+DA:260,48
+DA:261,23
+DA:262,23
+DA:263,13
+DA:264,0
+DA:265,0
+DA:268,23
+DA:269,0
+DA:270,0
+DA:272,23
+DA:273,23
+DA:274,0
+DA:275,0
DA:277,0
DA:278,0
-DA:280,0
-DA:281,0
-DA:283,4
-DA:284,4
-DA:285,4
-DA:286,8
+DA:280,23
+DA:281,23
+DA:282,23
+DA:283,46
+DA:284,0
+DA:285,0
+DA:286,0
DA:287,0
-DA:288,0
DA:289,0
-DA:290,0
-DA:292,0
-DA:293,4
-DA:294,5
-DA:295,12
-DA:296,18
-DA:297,24
-DA:298,12
-DA:299,24
-DA:300,12
-DA:301,15
-DA:302,4
-DA:303,4
-DA:304,13
-DA:305,13
-DA:306,4
-DA:307,4
-DA:309,13
-DA:310,13
-DA:311,2
-DA:312,2
-DA:314,13
-DA:315,147
-DA:316,26
-DA:317,26
-DA:318,3
-DA:319,3
-DA:322,111
-DA:323,23
-DA:324,13
-DA:325,13
-DA:326,13
+DA:290,23
+DA:291,26
+DA:292,96
+DA:293,153
+DA:294,192
+DA:295,96
+DA:296,192
+DA:297,96
+DA:298,116
+DA:299,23
+DA:300,24
+DA:301,85
+DA:302,85
+DA:303,14
+DA:304,9
+DA:306,85
+DA:307,85
+DA:308,8
+DA:309,5
+DA:311,85
+DA:312,1011
+DA:313,170
+DA:314,170
+DA:315,9
+DA:316,6
+DA:319,759
+DA:320,167
+DA:321,85
+DA:322,85
+DA:323,85
+DA:324,107
+DA:325,29
+DA:326,17
DA:327,17
-DA:328,8
+DA:328,17
DA:329,0
DA:330,0
-DA:331,0
-DA:332,0
-DA:333,0
-DA:335,0
-DA:336,0
-DA:337,0
-DA:338,0
-DA:339,0
-DA:341,28000
-DA:342,28000
-DA:343,27996
-DA:344,27996
-DA:346,4
-DA:347,13
-DA:348,13
-DA:349,13
-DA:351,28000
-DA:352,3
-LH:242
-LF:295
+DA:332,34
+DA:333,34
+DA:334,17
+DA:335,17
+DA:336,17
+DA:338,43455
+DA:339,43455
+DA:340,43432
+DA:341,43432
+DA:343,24
+DA:344,85
+DA:345,85
+DA:346,85
+DA:348,43455
+DA:349,5
+DA:350,5
+LH:272
+LF:293
end_of_record
SF:src\dense\NL_integrators\NL_QSS_Integrator.jl
DA:1,13
@@ -2004,190 +2072,191 @@ DA:152,64023
DA:153,64023
DA:154,64023
DA:155,13
-LH:133
-LF:139
+DA:156,13
+LH:134
+LF:140
end_of_record
SF:src\dense\NL_integrators\NL_QSS_discreteIntegrator.jl
-DA:1,4
-DA:2,4
-DA:3,4
-DA:4,4
-DA:5,4
-DA:6,4
-DA:7,4
-DA:8,4
-DA:9,4
-DA:10,4
-DA:11,4
-DA:12,4
-DA:13,4
-DA:14,4
-DA:15,4
-DA:16,4
-DA:17,4
-DA:18,4
-DA:19,4
-DA:20,4
-DA:21,4
-DA:22,4
-DA:23,4
-DA:24,4
-DA:25,4
-DA:26,4
-DA:27,14
-DA:28,4
-DA:29,4
-DA:30,4
-DA:31,7
-DA:32,7
-DA:33,23
-DA:34,39
-DA:35,7
-DA:36,72
-DA:37,23
-DA:38,23
-DA:39,23
-DA:40,39
-DA:41,10
-DA:42,4
-DA:43,4
-DA:44,4
-DA:45,14
-DA:46,14
-DA:47,14
-DA:48,14
-DA:49,21
-DA:50,28
-DA:51,14
-DA:52,3
-DA:53,3
-DA:54,2
-DA:55,1
-DA:56,1
-DA:57,1
-DA:58,0
-DA:62,11
-DA:64,24
-DA:65,4
-DA:66,38
-DA:67,7
-DA:68,7
-DA:69,7
-DA:70,7
-DA:71,10
-DA:72,4
-DA:73,4
-DA:74,4
-DA:75,4
-DA:76,4
-DA:77,4
-DA:78,208164
-DA:79,208164
-DA:80,208164
-DA:81,208164
-DA:82,208164
-DA:83,208164
-DA:84,208164
-DA:85,4
-DA:87,208160
-DA:88,208160
-DA:89,208160
-DA:90,203553
-DA:91,203553
-DA:92,203553
-DA:93,369296
-DA:94,203553
-DA:95,203553
-DA:96,238853
-DA:97,407106
-DA:101,203553
-DA:102,369296
-DA:103,535039
-DA:104,203553
-DA:105,203553
-DA:106,203553
-DA:107,398244
-DA:108,398244
-DA:109,362789
-DA:110,197443
-DA:112,398244
-DA:113,398244
-DA:114,362393
-DA:115,197245
-DA:117,398244
-DA:118,796488
-DA:119,796488
-DA:120,37821
-DA:121,36471
-DA:123,1194732
-DA:124,3427120
-DA:125,398244
-DA:126,398244
-DA:127,398244
-DA:128,601797
-DA:129,242790
-DA:130,164316
-DA:131,164316
-DA:132,164316
-DA:133,328632
-DA:134,1643160
-DA:135,164316
-DA:136,164316
-DA:137,328632
-DA:138,4607
-DA:139,4552
-DA:140,4711
-DA:141,4552
-DA:142,4552
-DA:143,9103
-DA:144,9104
-DA:145,4552
-DA:146,4711
-DA:147,4870
-DA:148,4552
-DA:153,9107
-DA:154,4552
-DA:155,4552
-DA:156,4552
-DA:157,4552
-DA:158,44
-DA:159,22
-DA:160,22
-DA:161,22
-DA:162,22
-DA:165,4552
-DA:166,4552
-DA:167,4552
-DA:168,4711
-DA:169,4552
-DA:171,4552
-DA:172,4552
-DA:173,4711
-DA:174,4552
-DA:176,4552
-DA:177,9104
-DA:178,9104
-DA:179,13656
-DA:180,9107
-DA:181,4552
-DA:182,4552
-DA:183,4552
-DA:184,9104
-DA:185,4552
-DA:186,4552
-DA:187,4552
-DA:188,4552
-DA:189,4552
-DA:190,4552
-DA:191,9104
-DA:192,9104
-DA:193,13656
-DA:194,9107
-DA:195,4552
-DA:196,4552
-DA:197,4552
-DA:198,9104
-DA:199,9104
+DA:1,5
+DA:2,5
+DA:3,5
+DA:4,5
+DA:5,5
+DA:6,5
+DA:7,5
+DA:8,5
+DA:9,5
+DA:10,5
+DA:11,5
+DA:12,5
+DA:13,5
+DA:14,5
+DA:15,5
+DA:16,5
+DA:17,5
+DA:18,5
+DA:19,5
+DA:20,5
+DA:21,5
+DA:22,5
+DA:23,5
+DA:24,5
+DA:25,5
+DA:26,5
+DA:27,20
+DA:28,5
+DA:29,5
+DA:30,5
+DA:31,9
+DA:32,9
+DA:33,35
+DA:34,61
+DA:35,9
+DA:36,120
+DA:37,35
+DA:38,35
+DA:39,35
+DA:40,61
+DA:41,13
+DA:42,5
+DA:43,5
+DA:44,5
+DA:45,20
+DA:46,20
+DA:47,20
+DA:48,20
+DA:49,32
+DA:50,40
+DA:51,20
+DA:52,4
+DA:53,4
+DA:54,6
+DA:55,2
+DA:56,2
+DA:57,2
+DA:58,1
+DA:62,16
+DA:64,35
+DA:65,5
+DA:66,50
+DA:67,10
+DA:68,10
+DA:69,10
+DA:70,10
+DA:71,15
+DA:72,5
+DA:73,5
+DA:74,5
+DA:75,5
+DA:76,5
+DA:77,5
+DA:78,209801
+DA:79,209801
+DA:80,209801
+DA:81,209801
+DA:82,209801
+DA:83,209801
+DA:84,209801
+DA:85,5
+DA:87,209796
+DA:88,209796
+DA:89,209796
+DA:90,205183
+DA:91,205183
+DA:92,205183
+DA:93,372556
+DA:94,205183
+DA:95,205183
+DA:96,242086
+DA:97,410366
+DA:101,205183
+DA:102,372556
+DA:103,539929
+DA:104,205183
+DA:105,205183
+DA:106,205183
+DA:107,402317
+DA:108,402317
+DA:109,368295
+DA:110,200196
+DA:112,402317
+DA:113,402317
+DA:114,365117
+DA:115,198607
+DA:117,402317
+DA:118,808107
+DA:119,808107
+DA:120,45967
+DA:121,40544
+DA:123,1210424
+DA:124,3443476
+DA:125,402317
+DA:126,402317
+DA:127,402317
+DA:128,607500
+DA:129,244684
+DA:130,165682
+DA:131,167048
+DA:132,167048
+DA:133,332730
+DA:134,1648624
+DA:135,165682
+DA:136,165682
+DA:137,331364
+DA:138,4613
+DA:139,4557
+DA:140,4721
+DA:141,4557
+DA:142,4557
+DA:143,9113
+DA:144,9114
+DA:145,4557
+DA:146,4721
+DA:147,4885
+DA:148,4557
+DA:153,9127
+DA:154,4557
+DA:155,4557
+DA:156,4557
+DA:157,4557
+DA:158,64
+DA:159,27
+DA:160,27
+DA:161,27
+DA:162,27
+DA:165,4557
+DA:166,4557
+DA:167,4557
+DA:168,4721
+DA:169,4557
+DA:171,4557
+DA:172,4557
+DA:173,4721
+DA:174,4557
+DA:176,4557
+DA:177,9119
+DA:178,9134
+DA:179,13676
+DA:180,9127
+DA:181,4557
+DA:182,4557
+DA:183,4557
+DA:184,9114
+DA:185,4557
+DA:186,4557
+DA:187,4557
+DA:188,4557
+DA:189,4557
+DA:190,4557
+DA:191,9119
+DA:192,9119
+DA:193,13676
+DA:194,9127
+DA:195,4557
+DA:196,4557
+DA:197,4557
+DA:198,9114
+DA:199,9114
DA:200,0
DA:201,0
DA:202,0
@@ -2198,54 +2267,54 @@ DA:207,0
DA:208,0
DA:209,0
DA:210,0
-DA:212,78
+DA:212,80
DA:213,32
DA:214,32
DA:215,26
DA:216,13
DA:218,64
-DA:219,518
-DA:220,55
-DA:221,55
-DA:222,32
+DA:219,522
+DA:220,56
+DA:221,56
+DA:222,33
DA:223,0
DA:224,0
-DA:227,55
+DA:227,56
DA:228,0
DA:229,0
-DA:231,55
-DA:232,39
+DA:231,56
+DA:232,40
DA:233,16
DA:234,16
DA:236,0
DA:237,0
-DA:239,55
-DA:240,55
-DA:241,55
-DA:242,110
+DA:239,56
+DA:240,56
+DA:241,56
+DA:242,112
DA:243,0
DA:244,0
DA:245,0
DA:246,0
DA:248,0
-DA:249,55
-DA:250,107
-DA:251,12
-DA:252,18
-DA:253,24
-DA:254,57
-DA:255,12
-DA:256,12
-DA:257,12
-DA:258,15
-DA:259,55
-DA:260,55
+DA:249,56
+DA:250,109
+DA:251,11
+DA:252,17
+DA:253,22
+DA:254,51
+DA:255,11
+DA:256,11
+DA:257,11
+DA:258,14
+DA:259,56
+DA:260,57
DA:261,64
-DA:262,168
+DA:262,171
DA:263,64
DA:264,64
-DA:265,40
-DA:266,21
+DA:265,42
+DA:266,22
DA:268,64
DA:269,300
DA:270,128
@@ -2258,7 +2327,7 @@ DA:278,64
DA:279,64
DA:280,64
DA:281,119
-DA:282,59
+DA:282,61
DA:283,87
DA:284,51
DA:285,51
@@ -2269,17 +2338,18 @@ DA:290,690
DA:291,69
DA:293,69
DA:294,69
-DA:296,208160
-DA:297,208105
-DA:298,208105
-DA:300,55
+DA:296,209796
+DA:297,209740
+DA:298,209740
+DA:300,57
DA:301,64
DA:302,64
DA:303,64
-DA:306,208160
-DA:307,4
-LH:246
-LF:270
+DA:306,209796
+DA:307,5
+DA:308,5
+LH:248
+LF:271
end_of_record
SF:src\dense\NL_integrators\NL_nmLiQSS_Integrator.jl
DA:1,12
@@ -2478,80 +2548,81 @@ DA:231,661544
DA:232,661544
DA:233,661544
DA:234,12
-LH:195
-LF:196
+DA:235,12
+LH:196
+LF:197
end_of_record
SF:src\dense\NL_integrators\NL_nmLiQSS_discreteIntegrator.jl
-DA:1,7
-DA:2,7
-DA:3,7
-DA:4,7
-DA:5,7
-DA:6,7
-DA:7,7
-DA:8,7
-DA:9,7
-DA:10,7
-DA:11,7
-DA:12,7
-DA:13,7
-DA:14,7
-DA:15,7
-DA:16,7
-DA:17,7
-DA:18,7
-DA:19,7
-DA:20,7
-DA:21,7
-DA:22,7
-DA:23,7
-DA:24,7
-DA:27,7
-DA:28,7
-DA:29,7
-DA:30,7
-DA:31,7
-DA:32,7
-DA:33,7
-DA:34,7
-DA:36,7
-DA:37,7
-DA:38,7
-DA:39,7
-DA:40,76
-DA:41,76
-DA:42,7
-DA:43,90
-DA:44,90
-DA:45,83
-DA:46,7
-DA:47,7
-DA:48,7
-DA:50,30
-DA:51,7
-DA:53,7
-DA:54,7
-DA:55,13
-DA:56,13
-DA:57,52
-DA:58,91
-DA:59,13
-DA:60,170
-DA:61,52
-DA:62,52
-DA:63,52
-DA:64,91
-DA:65,19
-DA:66,7
-DA:67,7
-DA:68,7
-DA:69,27
-DA:70,27
-DA:71,27
-DA:72,27
-DA:73,42
-DA:74,54
-DA:75,27
+DA:1,8
+DA:2,8
+DA:3,8
+DA:4,8
+DA:5,8
+DA:6,8
+DA:7,8
+DA:8,8
+DA:9,8
+DA:10,8
+DA:11,8
+DA:12,8
+DA:13,8
+DA:14,8
+DA:15,8
+DA:16,8
+DA:17,8
+DA:18,8
+DA:19,8
+DA:20,8
+DA:21,8
+DA:22,8
+DA:23,8
+DA:24,8
+DA:27,8
+DA:28,8
+DA:29,8
+DA:30,8
+DA:31,8
+DA:32,8
+DA:33,8
+DA:34,8
+DA:36,8
+DA:37,8
+DA:38,8
+DA:39,8
+DA:40,88
+DA:41,88
+DA:42,8
+DA:43,104
+DA:44,104
+DA:45,96
+DA:46,8
+DA:47,8
+DA:48,8
+DA:50,42
+DA:51,8
+DA:53,8
+DA:54,8
+DA:55,15
+DA:56,15
+DA:57,78
+DA:58,141
+DA:59,15
+DA:60,482
+DA:61,78
+DA:62,78
+DA:63,78
+DA:64,141
+DA:65,22
+DA:66,8
+DA:67,8
+DA:68,8
+DA:69,40
+DA:70,40
+DA:71,40
+DA:72,40
+DA:73,55
+DA:74,80
+DA:75,40
DA:76,5
DA:77,5
DA:78,4
@@ -2559,128 +2630,128 @@ DA:79,1
DA:80,1
DA:81,1
DA:82,1
-DA:86,22
-DA:88,47
-DA:89,7
-DA:90,84
-DA:91,15
-DA:92,15
-DA:93,15
-DA:94,15
-DA:95,23
-DA:101,7
-DA:102,7
-DA:103,7
-DA:104,7
-DA:105,7
-DA:106,7
-DA:107,7
-DA:108,7
-DA:109,7
-DA:110,7
-DA:111,83308
-DA:112,83308
-DA:113,83308
-DA:114,83308
-DA:115,83308
-DA:116,83308
-DA:117,83308
+DA:86,35
+DA:88,72
+DA:89,8
+DA:90,156
+DA:91,21
+DA:92,21
+DA:93,21
+DA:94,21
+DA:95,34
+DA:101,8
+DA:102,8
+DA:103,8
+DA:104,8
+DA:105,8
+DA:106,8
+DA:107,8
+DA:108,8
+DA:109,8
+DA:110,8
+DA:111,179600
+DA:112,179599
+DA:113,179599
+DA:114,179599
+DA:115,179599
+DA:116,179599
+DA:117,179599
DA:118,7
-DA:120,83301
-DA:121,83301
-DA:123,83301
-DA:124,82810
-DA:125,82810
-DA:126,82810
-DA:127,82810
-DA:128,85910
-DA:129,82810
-DA:130,82810
-DA:131,82810
-DA:132,147703
-DA:133,165620
-DA:134,82993
-DA:135,165746
-DA:136,165746
-DA:137,172072
-DA:138,165746
-DA:140,248373
-DA:141,82810
-DA:142,82810
-DA:144,82810
-DA:145,82810
-DA:146,165414
-DA:147,332555
-DA:148,332555
-DA:149,5446
-DA:150,2723
-DA:152,497969
-DA:153,165414
-DA:154,168599
-DA:155,165414
-DA:156,165414
-DA:157,168599
-DA:158,165414
-DA:159,165414
-DA:160,81438
-DA:161,108
-DA:162,1080
-DA:163,108
-DA:164,108
-DA:165,108
-DA:166,216
-DA:167,0
-DA:168,0
-DA:169,0
-DA:170,0
-DA:171,0
+DA:120,179592
+DA:121,179592
+DA:123,179592
+DA:124,134942
+DA:125,134942
+DA:126,134942
+DA:127,134942
+DA:128,190174
+DA:129,134942
+DA:130,134942
+DA:131,134942
+DA:132,251801
+DA:133,269884
+DA:134,135125
+DA:135,345207
+DA:136,345207
+DA:137,530612
+DA:138,345016
+DA:140,479966
+DA:141,134942
+DA:142,134942
+DA:144,134942
+DA:145,134942
+DA:146,344884
+DA:147,972681
+DA:148,972681
+DA:149,223178
+DA:150,111589
+DA:152,1317565
+DA:153,344884
+DA:154,348294
+DA:155,344884
+DA:156,344884
+DA:157,348294
+DA:158,344884
+DA:159,344884
+DA:160,208659
+DA:161,47802
+DA:162,573408
+DA:163,47802
+DA:164,47802
+DA:165,47802
+DA:166,174881
+DA:167,79277
+DA:168,158554
+DA:169,79277
+DA:170,79277
+DA:171,79277
DA:172,0
DA:173,0
-DA:175,0
-DA:176,0
-DA:177,0
-DA:178,0
-DA:179,0
-DA:181,0
-DA:182,0
-DA:183,0
-DA:184,0
-DA:185,0
-DA:187,324
-DA:188,186
-DA:189,30
-DA:190,30
-DA:191,30
-DA:192,60
-DA:193,300
-DA:194,30
-DA:195,30
-DA:196,60
-DA:199,248224
-DA:200,82810
-DA:201,108
-DA:207,82810
-DA:208,165414
-DA:209,165414
-DA:210,3072
-DA:211,1536
-DA:213,165414
-DA:214,165414
-DA:215,1615400
-DA:216,165414
-DA:217,165414
-DA:218,165414
-DA:219,248224
-DA:220,92199
-DA:221,73421
-DA:222,73777
-DA:223,73777
-DA:224,147198
-DA:225,731466
-DA:226,73421
-DA:227,73421
-DA:228,146842
-DA:230,491
+DA:175,79277
+DA:176,298050
+DA:177,298050
+DA:178,174250
+DA:179,87125
+DA:181,377327
+DA:182,951324
+DA:183,79277
+DA:184,79277
+DA:185,79277
+DA:187,222683
+DA:188,53014
+DA:189,42590
+DA:190,127710
+DA:191,127710
+DA:192,170300
+DA:193,511020
+DA:194,42590
+DA:195,42590
+DA:196,85180
+DA:199,479826
+DA:200,134942
+DA:201,34971
+DA:207,134942
+DA:208,344884
+DA:209,344884
+DA:210,3306
+DA:211,1653
+DA:213,344884
+DA:214,344884
+DA:215,3766753
+DA:216,344884
+DA:217,344884
+DA:218,344884
+DA:219,479826
+DA:220,147654
+DA:221,122230
+DA:222,219848
+DA:223,219848
+DA:224,342078
+DA:225,1314326
+DA:226,122230
+DA:227,122230
+DA:228,244460
+DA:230,44650
DA:231,382
DA:232,764
DA:233,382
@@ -2744,112 +2815,113 @@ DA:301,274
DA:302,137
DA:303,137
DA:304,137
-DA:307,153
-DA:308,65
-DA:309,65
-DA:310,58
-DA:311,29
-DA:313,130
-DA:314,908
-DA:315,109
-DA:316,109
-DA:317,55
-DA:318,0
-DA:319,0
-DA:322,109
-DA:323,0
-DA:324,0
-DA:326,109
-DA:327,85
-DA:328,24
-DA:329,24
+DA:307,44322
+DA:308,132512
+DA:309,132512
+DA:310,153838
+DA:311,76919
+DA:313,176726
+DA:314,530816
+DA:315,44268
+DA:316,44268
+DA:317,8629
+DA:318,5448
+DA:319,5448
+DA:322,38820
+DA:323,2705
+DA:324,2705
+DA:326,36115
+DA:327,20764
+DA:328,15351
+DA:329,15351
DA:331,0
DA:332,0
-DA:334,109
-DA:335,109
-DA:336,109
-DA:337,218
+DA:334,36115
+DA:335,36115
+DA:336,36115
+DA:337,72230
DA:338,0
DA:339,0
DA:340,0
DA:341,0
DA:343,0
-DA:344,109
-DA:345,198
-DA:346,97
-DA:347,154
-DA:348,194
-DA:349,97
-DA:350,97
-DA:351,97
-DA:352,97
-DA:353,117
-DA:354,109
-DA:355,110
-DA:356,171
-DA:357,171
-DA:358,97
-DA:359,58
-DA:361,171
-DA:362,171
-DA:363,65
-DA:364,42
-DA:366,171
-DA:367,1269
-DA:368,342
-DA:369,342
-DA:370,93
-DA:371,56
-DA:374,1017
-DA:375,1027
-DA:376,171
-DA:377,171
-DA:378,171
-DA:379,279
-DA:380,115
-DA:381,175
-DA:382,103
-DA:383,103
+DA:344,36115
+DA:345,72210
+DA:346,96
+DA:347,153
+DA:348,192
+DA:349,96
+DA:350,96
+DA:351,96
+DA:352,96
+DA:353,116
+DA:354,36115
+DA:355,36116
+DA:356,432243
+DA:357,432243
+DA:358,730609
+DA:359,365314
+DA:361,432243
+DA:362,432243
+DA:363,194375
+DA:364,97197
+DA:366,432243
+DA:367,21172797
+DA:368,1440582
+DA:369,1440582
+DA:370,339753
+DA:371,169886
+DA:374,10802817
+DA:375,5185891
+DA:376,432243
+DA:377,432243
+DA:378,432243
+DA:379,468357
+DA:380,36121
+DA:381,144219
+DA:382,432175
+DA:383,432175
DA:384,0
DA:385,0
-DA:387,206
-DA:388,1254
-DA:389,139
-DA:390,139
-DA:391,139
-DA:393,83301
-DA:394,83301
-DA:395,83192
-DA:396,83192
-DA:398,110
-DA:399,171
-DA:400,171
-DA:401,171
-DA:403,83301
-DA:404,7
-LH:315
-LF:345
+DA:387,576302
+DA:388,1729662
+DA:389,144173
+DA:390,144173
+DA:391,144173
+DA:393,171439
+DA:394,171439
+DA:395,135324
+DA:396,135324
+DA:398,36116
+DA:399,432243
+DA:400,432243
+DA:401,432243
+DA:403,179592
+DA:404,8
+DA:405,8
+LH:335
+LF:346
end_of_record
SF:src\dense\Quantizers\LiQSS_quantizer1.jl
-DA:1,818921
-DA:2,818921
-DA:3,869991
-DA:4,818921
-DA:5,818921
-DA:6,818921
-DA:7,818921
-DA:8,818921
-DA:9,818921
-DA:10,818921
-DA:11,818921
-DA:12,818921
-DA:13,818921
-DA:14,797098
-DA:15,797098
-DA:16,797098
-DA:17,797098
-DA:18,797098
-DA:19,797098
+DA:1,831217
+DA:2,831217
+DA:3,894583
+DA:4,831217
+DA:5,831217
+DA:6,831217
+DA:7,831217
+DA:8,831217
+DA:9,831217
+DA:10,831217
+DA:11,831217
+DA:12,831217
+DA:13,831217
+DA:14,811104
+DA:15,811104
+DA:16,811104
+DA:17,811104
+DA:18,811104
+DA:19,811104
DA:20,778689
DA:21,671914
DA:22,671914
@@ -2858,127 +2930,127 @@ DA:24,90237
DA:25,90237
DA:27,16538
DA:28,795227
-DA:31,18409
-DA:32,11351
-DA:33,11351
-DA:34,7058
-DA:35,7056
-DA:36,7056
-DA:38,2
-DA:39,2
-DA:40,2
+DA:31,32415
+DA:32,20780
+DA:33,20780
+DA:34,11635
+DA:35,11635
+DA:36,11635
+DA:38,0
+DA:39,0
+DA:40,0
DA:42,0
-DA:43,797098
-DA:48,21823
+DA:43,811104
+DA:48,20113
DA:49,19497
-DA:51,2326
-DA:53,21823
+DA:51,616
+DA:53,20113
DA:54,20112
-DA:56,1711
-DA:59,818921
-DA:60,818921
-DA:61,818921
-DA:63,2180126
-DA:64,2180126
-DA:65,2180126
-DA:66,2180126
-DA:67,2180126
-DA:68,2180126
+DA:56,1
+DA:59,831217
+DA:60,831217
+DA:61,831217
+DA:63,2202375
+DA:64,2202375
+DA:65,2202375
+DA:66,2202375
+DA:67,2202375
+DA:68,2202375
DA:69,1
-DA:71,2180125
-DA:72,2176017
-DA:73,2176017
-DA:74,2096734
-DA:75,79283
-DA:76,78882
-DA:77,46381
-DA:79,2208518
-DA:83,4108
-DA:86,2180126
+DA:71,2202374
+DA:72,2201682
+DA:73,2201682
+DA:74,2122387
+DA:75,79295
+DA:76,78894
+DA:77,46391
+DA:79,2234185
+DA:83,692
+DA:86,2202375
DA:87,0
-LH:63
+LH:60
LF:65
end_of_record
SF:src\dense\Quantizers\LiQSS_quantizer2.jl
-DA:2,53086
-DA:3,53086
-DA:4,56572
-DA:5,53086
-DA:7,53086
-DA:8,53086
-DA:9,53086
-DA:10,53086
-DA:11,53086
-DA:12,53086
-DA:13,53086
-DA:14,53086
-DA:15,53086
-DA:16,53086
-DA:17,53086
-DA:18,53086
-DA:19,53086
-DA:20,53086
-DA:21,53086
-DA:22,52564
-DA:23,24123
-DA:24,8217
-DA:25,8217
-DA:27,15906
-DA:28,15906
-DA:29,40029
-DA:32,28441
-DA:33,28441
-DA:34,6473
-DA:35,6473
-DA:37,21968
-DA:38,21968
-DA:39,50409
+DA:2,108318
+DA:3,108318
+DA:4,113807
+DA:5,108318
+DA:7,108318
+DA:8,108318
+DA:9,108318
+DA:10,108318
+DA:11,108318
+DA:12,108318
+DA:13,108318
+DA:14,108318
+DA:15,108318
+DA:16,108318
+DA:17,108318
+DA:18,108318
+DA:19,108318
+DA:20,108318
+DA:21,108318
+DA:22,107493
+DA:23,67860
+DA:24,19932
+DA:25,19932
+DA:27,47928
+DA:28,47928
+DA:29,115788
+DA:32,39633
+DA:33,39633
+DA:34,9628
+DA:35,9628
+DA:37,30005
+DA:38,30005
+DA:39,69638
DA:42,0
DA:43,0
-DA:45,52564
-DA:46,37874
-DA:47,37874
-DA:49,37874
-DA:52,67254
-DA:56,522
-DA:57,509
-DA:58,509
-DA:59,509
-DA:62,13
+DA:45,107493
+DA:46,77933
+DA:47,77933
+DA:49,77933
+DA:52,137053
+DA:56,825
+DA:57,811
+DA:58,811
+DA:59,811
+DA:62,14
DA:64,0
DA:65,0
DA:66,0
-DA:68,13
-DA:69,13
-DA:70,13
-DA:74,53086
-DA:75,53086
-DA:76,53086
-DA:77,53086
-DA:79,242683
-DA:80,242683
-DA:81,242683
-DA:82,242683
-DA:83,242683
-DA:84,242683
-DA:85,242683
-DA:86,242683
-DA:87,242683
-DA:88,3
-DA:90,242680
-DA:91,242680
-DA:92,242680
-DA:93,122862
-DA:94,122862
-DA:95,122862
-DA:96,77516
-DA:98,119818
-DA:99,119197
-DA:100,119197
-DA:101,119197
-DA:102,79525
-DA:105,242680
-DA:106,15
+DA:68,14
+DA:69,14
+DA:70,14
+DA:74,108318
+DA:75,108318
+DA:76,108318
+DA:77,108318
+DA:79,974989
+DA:80,974989
+DA:81,974989
+DA:82,974989
+DA:83,974989
+DA:84,974989
+DA:85,974989
+DA:86,974989
+DA:87,974989
+DA:88,2
+DA:90,974987
+DA:91,974987
+DA:92,974987
+DA:93,662867
+DA:94,662867
+DA:95,662867
+DA:96,387010
+DA:98,312120
+DA:99,311358
+DA:100,311358
+DA:101,311358
+DA:102,202953
+DA:105,974987
+DA:106,13
LH:74
LF:79
end_of_record
@@ -2993,45 +3065,45 @@ DA:9,0
DA:10,114180
DA:14,6
DA:16,114186
-DA:18,172666
-DA:19,172666
-DA:20,172666
-DA:21,168467
-DA:22,168467
-DA:23,168467
+DA:18,174453
+DA:19,174453
+DA:20,174453
+DA:21,170238
+DA:22,170238
+DA:23,170238
DA:25,0
-DA:26,168467
-DA:30,4199
+DA:26,170238
+DA:30,4215
DA:31,4190
DA:32,4190
DA:33,4190
DA:35,0
DA:36,4190
-DA:40,9
-DA:43,172666
+DA:40,25
+DA:43,174453
DA:62,190698
DA:63,190698
DA:64,190698
-DA:65,4150
-DA:67,186548
-DA:68,370707
-DA:69,186548
-DA:70,370707
-DA:71,186548
-DA:72,186548
-DA:73,186548
-DA:76,339628
-DA:77,339628
-DA:78,339628
-DA:79,56
-DA:81,339572
-DA:82,339572
-DA:83,339572
-DA:84,339572
-DA:85,339572
-DA:86,339572
-DA:87,339572
-DA:89,339628
+DA:65,738
+DA:67,189960
+DA:68,377531
+DA:69,189960
+DA:70,377531
+DA:71,189960
+DA:72,189960
+DA:73,189960
+DA:76,344001
+DA:77,344001
+DA:78,344001
+DA:79,80
+DA:81,343921
+DA:82,343921
+DA:83,343921
+DA:84,343921
+DA:85,343921
+DA:86,343921
+DA:87,343921
+DA:89,344001
DA:93,43
DA:94,43
DA:95,43
@@ -3043,56 +3115,58 @@ DA:102,43
DA:103,43
DA:105,0
DA:107,43
-DA:110,59
-DA:111,118
-DA:112,59
-DA:113,59
-DA:114,59
-DA:115,59
+DA:110,124
+DA:111,248
+DA:112,124
+DA:113,124
+DA:114,124
+DA:115,124
DA:117,0
-DA:120,59
+DA:120,124
DA:121,0
-DA:124,59
-DA:125,59
-DA:126,59
-DA:127,59
+DA:124,124
+DA:125,124
+DA:126,124
+DA:127,124
DA:129,0
-DA:131,59
-DA:132,50
-DA:134,9
+DA:131,124
+DA:132,100
+DA:134,24
DA:135,0
-DA:137,9
-DA:141,59
-DA:168,238306
-DA:169,238306
-DA:170,55
-DA:171,55
-DA:172,238251
-DA:173,74
-DA:175,238202
-DA:176,238177
-DA:177,17
-DA:179,238177
-DA:183,238177
-DA:184,238177
-LH:83
-LF:92
+DA:137,24
+DA:141,124
+DA:168,517734
+DA:169,517734
+DA:170,39866
+DA:171,39866
+DA:172,477868
+DA:173,93
+DA:175,477800
+DA:176,477775
+DA:177,11705
+DA:179,477775
+DA:180,477775
+DA:181,0
+DA:183,477775
+DA:184,477775
+LH:84
+LF:94
end_of_record
SF:src\dense\Quantizers\Quantizer_Common.jl
-DA:1,3288252
-DA:3,3288252
-DA:5,1384937
-DA:6,1837819
-DA:8,224562
-DA:9,449124
-DA:10,224562
-DA:21,2400114
-DA:22,2400114
-DA:23,2400114
-DA:25,605347
-DA:26,605347
-DA:27,605347
-DA:28,605347
+DA:1,3298205
+DA:3,3298205
+DA:5,2125748
+DA:6,3307145
+DA:8,281469
+DA:9,562938
+DA:10,281469
+DA:21,2422363
+DA:22,2422363
+DA:23,2422363
+DA:25,1354924
+DA:26,1354924
+DA:27,1354924
+DA:28,1354924
LH:14
LF:14
end_of_record
@@ -3234,133 +3308,133 @@ LH:123
LF:133
end_of_record
SF:src\dense\Quantizers\mLiQSS_quantizer2.jl
-DA:1,92794
-DA:2,185621
-DA:3,92794
-DA:4,185621
-DA:5,92794
-DA:6,92794
-DA:7,92794
-DA:8,927940
-DA:9,185588
-DA:10,185588
-DA:11,92794
-DA:12,371176
-DA:13,371176
-DA:14,92794
-DA:15,92794
-DA:16,92794
-DA:17,92794
-DA:18,92794
-DA:19,92794
-DA:20,92794
-DA:21,92794
-DA:22,92794
-DA:23,92794
-DA:24,92794
-DA:25,92794
-DA:26,92794
-DA:27,92794
-DA:28,92794
-DA:29,92794
-DA:30,92794
-DA:31,92794
-DA:32,92794
-DA:33,92794
-DA:34,92794
-DA:35,92794
-DA:36,92794
-DA:37,92794
-DA:38,92794
-DA:39,92794
-DA:40,92794
-DA:41,92794
-DA:43,161576
-DA:44,70729
-DA:45,45439
-DA:69,184948
-DA:70,713
-DA:71,713
-DA:72,156
-DA:132,92794
-DA:133,45314
-DA:134,45314
-DA:135,45314
-DA:136,45314
-DA:137,45314
-DA:138,45314
-DA:139,45314
-DA:140,45314
-DA:141,45314
-DA:142,45314
-DA:143,45314
-DA:144,45314
-DA:145,45314
-DA:146,45314
-DA:147,45314
-DA:148,45314
-DA:149,45314
-DA:150,45314
-DA:151,45314
-DA:152,45314
-DA:153,45314
-DA:154,45314
-DA:155,45314
-DA:156,226570
-DA:157,45314
-DA:158,45314
-DA:159,45314
-DA:160,45314
-DA:163,83991
-DA:164,28058
-DA:165,14029
-DA:166,14029
-DA:167,14029
-DA:168,70145
-DA:169,14029
-DA:170,14029
-DA:171,14029
-DA:172,14029
-DA:176,45314
-DA:177,90628
-DA:178,0
-DA:179,0
-DA:180,0
-DA:181,0
-DA:182,0
-DA:183,0
-DA:184,0
-DA:185,0
-DA:186,0
-DA:187,0
-DA:188,0
-DA:191,0
-DA:192,45314
+DA:1,220015
+DA:2,440063
+DA:3,220015
+DA:4,440063
+DA:5,220015
+DA:6,220015
+DA:7,220015
+DA:8,2200150
+DA:9,440030
+DA:10,440030
+DA:11,220015
+DA:12,880060
+DA:13,880060
+DA:14,220015
+DA:15,220015
+DA:16,220015
+DA:17,220015
+DA:18,220015
+DA:19,220015
+DA:20,220015
+DA:21,220015
+DA:22,220021
+DA:23,220015
+DA:24,220015
+DA:25,220015
+DA:26,220015
+DA:27,220015
+DA:28,220015
+DA:29,220015
+DA:30,220015
+DA:31,220015
+DA:32,220015
+DA:33,220015
+DA:34,220015
+DA:35,220656
+DA:36,220015
+DA:37,220015
+DA:38,220015
+DA:39,220015
+DA:40,220015
+DA:41,220015
+DA:43,382897
+DA:44,154071
+DA:45,101927
+DA:69,438416
+DA:70,21138
+DA:71,21138
+DA:72,13462
+DA:132,220015
+DA:133,93008
+DA:134,93008
+DA:135,93008
+DA:136,93008
+DA:137,93008
+DA:138,93008
+DA:139,93008
+DA:140,93008
+DA:141,93008
+DA:142,93008
+DA:143,93008
+DA:144,93008
+DA:145,93008
+DA:146,93008
+DA:147,93008
+DA:148,93008
+DA:149,93008
+DA:150,93008
+DA:151,93008
+DA:152,93008
+DA:153,93008
+DA:154,93008
+DA:155,93008
+DA:156,465040
+DA:157,93008
+DA:158,93008
+DA:159,93008
+DA:160,93008
+DA:163,143817
+DA:164,123176
+DA:165,61588
+DA:166,61588
+DA:167,61588
+DA:168,307940
+DA:169,61588
+DA:170,61588
+DA:171,61588
+DA:172,61588
+DA:176,93008
+DA:177,186024
+DA:178,4
+DA:179,4
+DA:180,4
+DA:181,4
+DA:182,4
+DA:183,4
+DA:184,20
+DA:185,4
+DA:186,4
+DA:187,4
+DA:188,4
+DA:191,4
+DA:192,93008
DA:193,0
-DA:195,45314
+DA:195,93008
DA:196,0
DA:197,0
-DA:199,45314
+DA:199,93008
DA:200,0
DA:201,0
-DA:203,45314
+DA:203,93008
DA:204,0
DA:205,0
-DA:207,45314
-DA:208,45314
-DA:209,45314
-DA:210,45314
-DA:211,45314
-DA:212,45314
-DA:213,45314
-DA:215,92794
-LH:101
+DA:207,93008
+DA:208,93008
+DA:209,93008
+DA:210,93008
+DA:211,93008
+DA:212,93008
+DA:213,93008
+DA:215,220015
+LH:113
LF:120
end_of_record
SF:src\ownTaylor\arithmetic.jl
-DA:7,749768
-DA:8,749768
-DA:10,749750
+DA:7,2005557
+DA:8,2005557
+DA:10,2005539
DA:11,1
DA:12,3
DA:13,1
@@ -3380,10 +3454,10 @@ DA:30,1
DA:31,1
DA:32,2
DA:33,1
-DA:35,106
-DA:36,106
-DA:37,106
-DA:38,106
+DA:35,171
+DA:36,171
+DA:37,171
+DA:38,171
DA:40,1
DA:41,1
DA:43,1
@@ -3432,15 +3506,15 @@ DA:94,2
DA:95,5
DA:96,5
DA:97,1
-DA:99,749746
-DA:100,1499492
-DA:101,1499492
-DA:102,749746
-DA:103,749746
-DA:104,749746
-DA:105,749746
-DA:106,749746
-DA:109,749746
+DA:99,2005509
+DA:100,4011018
+DA:101,4011018
+DA:102,2005509
+DA:103,2005509
+DA:104,2005509
+DA:105,2005509
+DA:106,2005509
+DA:109,2005509
DA:111,3
DA:112,3
DA:113,3
@@ -3460,11 +3534,11 @@ LF:97
end_of_record
SF:src\ownTaylor\arithmeticT.jl
DA:1,0
-DA:2,5763
+DA:2,25413
DA:3,0
-DA:4,29822
-DA:5,29822
-DA:7,29822
+DA:4,70741
+DA:5,70741
+DA:7,70741
DA:8,6
DA:9,6
DA:10,5
@@ -3486,9 +3560,9 @@ DA:27,0
DA:28,0
DA:29,0
DA:30,0
-DA:32,29822
-DA:33,29822
-DA:34,29822
+DA:32,70741
+DA:33,70741
+DA:34,70741
DA:37,5610042
DA:38,5610042
DA:39,5610042
@@ -3504,26 +3578,26 @@ DA:48,1
DA:54,5610042
DA:55,5610042
DA:56,5610042
-DA:60,18944
-DA:61,18944
-DA:62,18944
+DA:60,19134
+DA:61,19134
+DA:62,19134
DA:68,6605
DA:69,6605
DA:70,6605
DA:75,1880944
DA:76,3761888
DA:77,1880944
-DA:79,21996
-DA:80,21996
-DA:81,21996
-DA:82,43992
-DA:83,21996
-DA:88,407403
-DA:89,407403
-DA:90,407403
-DA:91,814806
-DA:92,407403
-DA:94,407403
+DA:79,23129
+DA:80,23129
+DA:81,23129
+DA:82,46258
+DA:83,23129
+DA:88,408720
+DA:89,408720
+DA:90,408720
+DA:91,817440
+DA:92,408720
+DA:94,407185
DA:95,1
DA:96,2
DA:97,1
@@ -3536,26 +3610,26 @@ DA:105,1
DA:106,2
DA:107,2
DA:108,2
-DA:114,252674
-DA:115,505348
-DA:116,252674
+DA:114,501225
+DA:115,1002450
+DA:116,501225
DA:118,1
DA:119,1
DA:120,1
-DA:126,1921
-DA:127,3842
-DA:128,1921
-DA:130,61
-DA:131,61
-DA:132,61
-DA:133,61
-DA:134,61
-DA:136,61
-DA:137,1
-DA:138,2
-DA:139,1
-DA:140,2
-DA:141,1
+DA:126,8471
+DA:127,16942
+DA:128,8471
+DA:130,72
+DA:131,72
+DA:132,72
+DA:133,72
+DA:134,72
+DA:136,72
+DA:137,221197
+DA:138,442394
+DA:139,221197
+DA:140,442394
+DA:141,221197
DA:143,1
DA:144,1
DA:145,1
@@ -3568,36 +3642,36 @@ DA:153,1
DA:154,1
DA:155,1
DA:156,1
-DA:159,2281507
-DA:160,2281507
-DA:166,2517203
-DA:167,5034406
-DA:168,2517203
-DA:170,958166
-DA:171,958166
-DA:172,958166
-DA:173,958166
+DA:159,2282824
+DA:160,2282824
+DA:166,4316564
+DA:167,8633128
+DA:168,4316564
+DA:170,958431
+DA:171,958431
+DA:172,958431
+DA:173,958431
DA:175,1865769
DA:176,3731538
DA:177,1865769
DA:178,1865769
-DA:180,237557
-DA:182,237557
-DA:183,237557
-DA:196,269937
-DA:197,539874
-DA:198,269937
+DA:180,516746
+DA:182,516746
+DA:183,516746
+DA:196,1852991
+DA:197,3705982
+DA:198,1852991
DA:200,29826
DA:201,29826
DA:202,29826
DA:203,29826
DA:205,29816
-DA:206,721849
-DA:208,721849
-DA:209,721849
-DA:212,1889859
-DA:213,3779718
-DA:214,1889859
+DA:206,1984175
+DA:208,1984175
+DA:209,1984175
+DA:212,1930778
+DA:213,3861556
+DA:214,1930778
DA:216,7
DA:217,14
DA:218,7
@@ -3613,22 +3687,22 @@ DA:229,1
DA:230,1
DA:231,1
DA:232,1
-DA:235,4369981
-DA:236,10026764
-DA:237,4369981
-DA:238,4369981
-DA:240,3165757
-DA:244,267373
-DA:245,534746
-DA:246,729279
-DA:247,1191185
-DA:248,656439
-DA:249,850972
-DA:250,1191185
-DA:251,267373
-DA:253,65
-DA:257,65
-DA:258,65
+DA:235,6693014
+DA:236,16973614
+DA:237,6693014
+DA:238,6693014
+DA:240,3726078
+DA:244,546563
+DA:245,1093126
+DA:246,1566849
+DA:247,2587135
+DA:248,1494009
+DA:249,1967732
+DA:250,2587135
+DA:251,546563
+DA:253,76
+DA:257,76
+DA:258,76
DA:261,3750006
DA:262,7904174
DA:263,7500012
@@ -3646,96 +3720,96 @@ DA:275,7470082
DA:277,3720080
DA:281,3720080
DA:282,3720080
-DA:298,239479
-DA:299,239479
-DA:305,1216500
-DA:306,1216500
+DA:298,525219
+DA:299,525219
+DA:305,2472276
+DA:306,2472276
DA:308,1
DA:309,1
DA:310,1
-DA:312,734081
-DA:313,1963070
-DA:317,734081
-DA:318,734081
+DA:312,1513315
+DA:313,4300772
+DA:317,1513315
+DA:318,1513315
DA:320,1
DA:321,1
DA:322,1
DA:323,1
-DA:329,749744
-DA:330,749744
-DA:331,749742
-DA:332,749742
-DA:333,1499484
-DA:334,1274053
-DA:335,1274053
-DA:336,1798364
-DA:337,524311
-DA:338,524311
-DA:339,1274053
-DA:340,1274053
+DA:329,2005520
+DA:330,2005520
+DA:331,2005505
+DA:332,2005505
+DA:333,4011010
+DA:334,3785579
+DA:335,3785579
+DA:336,5565653
+DA:337,1780074
+DA:338,1780074
+DA:339,3785579
+DA:340,3785579
DA:342,0
-DA:344,1798364
-DA:345,749742
-DA:361,772983
-DA:362,772983
-DA:363,7653572
-DA:364,7653572
-DA:365,7653572
-DA:366,14534161
-DA:368,2474086
-DA:369,2474086
-DA:370,22292672
-DA:371,22292672
-DA:372,42111258
+DA:344,5565653
+DA:345,2005505
+DA:361,1804818
+DA:362,1804818
+DA:363,19912026
+DA:364,19912026
+DA:365,19912026
+DA:366,38019234
+DA:368,2496335
+DA:369,2496335
+DA:370,22337170
+DA:371,22337170
+DA:372,42178005
LH:216
LF:228
end_of_record
SF:src\ownTaylor\constructors.jl
-DA:15,754491
+DA:15,2010435
DA:18,1
-DA:19,749750
+DA:19,2005539
DA:20,44
DA:21,132
DA:22,44
DA:23,44
DA:25,2
-DA:26,118457401
-DA:27,109380092
+DA:26,155082965
+DA:27,157901189
DA:29,1
-DA:30,7737630
-DA:31,9919975
-DA:32,7737630
+DA:30,8016820
+DA:31,11585373
+DA:32,8016820
DA:33,1
DA:34,1
-DA:37,1499507
-DA:38,4498525
-DA:39,1499507
-DA:40,1499507
+DA:37,4011033
+DA:38,12033143
+DA:39,4011033
+DA:40,4011033
DA:43,6
DA:44,18
DA:45,6
DA:46,6
DA:48,36977
-DA:49,2182345
-DA:50,2182345
-DA:51,4364690
-DA:52,3189556
-DA:53,4196767
-DA:54,2182345
-DA:56,3189556
-DA:63,1079
-DA:64,1079
-DA:65,2158
-DA:66,3237
-DA:67,1079
-DA:68,1079
-DA:70,3200
-DA:71,3200
+DA:49,3568553
+DA:50,3568553
+DA:51,7137106
+DA:52,5939723
+DA:53,8310893
+DA:54,3568553
+DA:56,5939723
+DA:63,1109
+DA:64,1109
+DA:65,2218
+DA:66,3327
+DA:67,1109
+DA:68,1109
+DA:70,3260
+DA:71,3260
DA:72,0
-DA:73,3200
-DA:74,2121
-DA:76,1079
-DA:77,1079
+DA:73,3260
+DA:74,2151
+DA:76,1109
+DA:77,1109
DA:78,0
DA:79,0
DA:84,1
@@ -3883,188 +3957,179 @@ DA:158,2
DA:159,5
DA:160,5
DA:161,1
-DA:180,0
-DA:181,0
-DA:182,0
-DA:183,0
-DA:184,0
-DA:185,0
-DA:187,0
-DA:191,0
-DA:193,0
-DA:194,165
-DA:195,165
-DA:196,55
-DA:197,55
-DA:199,110
-DA:200,165
-DA:201,55
-DA:202,55
-DA:203,110
-DA:204,110
-DA:206,6
-DA:207,6
-DA:208,2
-DA:209,2
-DA:210,4
-DA:211,2
-DA:212,2
+DA:166,165
+DA:167,165
+DA:168,55
+DA:169,55
+DA:171,110
+DA:172,165
+DA:173,55
+DA:174,55
+DA:175,110
+DA:176,110
+DA:178,6
+DA:179,6
+DA:180,2
+DA:181,2
+DA:182,4
+DA:183,2
+DA:184,2
+DA:186,2
+DA:187,2
+DA:188,0
+DA:189,0
+DA:190,2
+DA:191,2
+DA:193,195
+DA:194,195
+DA:195,65
+DA:196,65
+DA:197,65
+DA:199,130
+DA:200,130
+DA:201,130
+DA:202,195
+DA:203,65
+DA:204,65
+DA:205,65
+DA:206,65
+DA:207,130
+DA:208,130
+DA:209,130
+DA:211,6
+DA:212,6
+DA:213,2
DA:214,2
DA:215,2
-DA:216,0
-DA:217,0
-DA:218,2
-DA:219,2
-DA:221,195
-DA:222,195
-DA:223,65
-DA:224,65
-DA:225,65
-DA:227,130
-DA:228,130
-DA:229,130
-DA:230,195
-DA:231,65
-DA:232,65
-DA:233,65
-DA:234,65
-DA:235,130
-DA:236,130
-DA:237,130
-DA:239,6
-DA:240,6
-DA:241,2
-DA:242,2
-DA:243,2
-DA:244,2
-DA:246,4
-DA:247,6
+DA:216,2
+DA:218,4
+DA:219,6
+DA:220,2
+DA:221,2
+DA:222,4
+DA:223,8
+DA:224,4
+DA:226,3
+DA:227,3
+DA:228,1
+DA:229,1
+DA:230,1
+DA:231,1
+DA:233,2
+DA:234,2
+DA:235,0
+DA:236,0
+DA:237,2
+DA:238,2
+DA:239,2
+DA:241,3
+DA:242,3
+DA:243,1
+DA:244,1
+DA:245,1
+DA:246,1
DA:248,2
DA:249,2
-DA:250,4
-DA:251,8
-DA:252,4
-DA:254,3
-DA:255,3
-DA:256,1
-DA:257,1
-DA:258,1
-DA:259,1
+DA:250,0
+DA:251,0
+DA:252,2
+DA:253,2
+DA:254,2
+DA:256,6
+DA:257,6
+DA:258,2
+DA:259,2
+DA:260,2
DA:261,2
-DA:262,2
-DA:263,0
-DA:264,0
-DA:265,2
-DA:266,2
-DA:267,2
-DA:269,3
-DA:270,3
-DA:271,1
-DA:272,1
-DA:273,1
-DA:274,1
-DA:276,2
-DA:277,2
-DA:278,0
-DA:279,0
-DA:280,2
+DA:263,4
+DA:264,4
+DA:265,0
+DA:266,0
+DA:267,8
+DA:268,4
+DA:269,4
+DA:271,6
+DA:272,6
+DA:273,2
+DA:274,2
+DA:275,2
+DA:277,4
+DA:278,4
+DA:279,4
+DA:280,6
DA:281,2
DA:282,2
-DA:284,6
-DA:285,6
-DA:286,2
-DA:287,2
-DA:288,2
-DA:289,2
-DA:291,4
-DA:292,4
-DA:293,0
-DA:294,0
-DA:295,8
-DA:296,4
-DA:297,4
-DA:299,6
-DA:300,6
-DA:301,2
+DA:283,2
+DA:284,2
+DA:285,4
+DA:286,4
+DA:287,4
+DA:289,3
+DA:290,3
+DA:291,1
+DA:292,1
+DA:293,1
+DA:294,1
+DA:296,2
+DA:297,3
+DA:298,1
+DA:299,1
+DA:300,2
+DA:301,4
DA:302,2
-DA:303,2
-DA:305,4
-DA:306,4
-DA:307,4
-DA:308,6
-DA:309,2
-DA:310,2
+DA:304,3
+DA:305,3
+DA:306,1
+DA:307,1
+DA:308,1
+DA:309,1
DA:311,2
DA:312,2
-DA:313,4
-DA:314,4
-DA:315,4
-DA:317,3
-DA:318,3
-DA:319,1
-DA:320,1
+DA:313,0
+DA:314,0
+DA:315,2
+DA:316,2
+DA:317,2
+DA:319,3
+DA:320,3
DA:321,1
DA:322,1
-DA:324,2
-DA:325,3
-DA:326,1
-DA:327,1
-DA:328,2
-DA:329,4
+DA:323,1
+DA:324,1
+DA:326,2
+DA:327,2
+DA:328,0
+DA:329,0
DA:330,2
-DA:332,3
-DA:333,3
-DA:334,1
-DA:335,1
+DA:331,2
+DA:332,2
+DA:334,3
+DA:335,3
DA:336,1
DA:337,1
-DA:339,2
-DA:340,2
-DA:341,0
-DA:342,0
-DA:343,2
-DA:344,2
-DA:345,2
-DA:347,3
-DA:348,3
-DA:349,1
-DA:350,1
+DA:338,1
+DA:339,1
+DA:341,2
+DA:342,2
+DA:343,0
+DA:344,0
+DA:345,4
+DA:346,2
+DA:347,2
+DA:349,2
+DA:350,2
DA:351,1
DA:352,1
-DA:354,2
-DA:355,2
-DA:356,0
-DA:357,0
-DA:358,2
-DA:359,2
-DA:360,2
-DA:362,3
-DA:363,3
-DA:364,1
-DA:365,1
-DA:366,1
-DA:367,1
-DA:369,2
-DA:370,2
-DA:371,0
-DA:372,0
-DA:373,4
-DA:374,2
-DA:375,2
-DA:377,2
-DA:378,2
-DA:379,1
-DA:380,1
-DA:381,1
-DA:383,0
-DA:388,0
-DA:390,0
-DA:391,0
-DA:392,0
-DA:393,0
-DA:395,0
+DA:353,1
+DA:355,0
+DA:360,0
+DA:362,0
+DA:363,0
+DA:364,0
+DA:365,0
+DA:367,0
LH:285
-LF:315
+LF:306
end_of_record
SF:src\ownTaylor\functionsT.jl
DA:9,54
@@ -4168,10 +4233,10 @@ DA:128,1
DA:130,4
DA:131,2
DA:132,2
-DA:135,0
-DA:136,0
-DA:137,0
-LH:97
+DA:135,1
+DA:136,1
+DA:137,1
+LH:100
LF:104
end_of_record
SF:src\ownTaylor\power.jl
@@ -4220,78 +4285,78 @@ DA:55,6
DA:56,5
DA:57,1
DA:59,12
-DA:71,18
-DA:72,9
+DA:60,18
+DA:61,9
+DA:62,0
+DA:63,0
+DA:65,9
+DA:66,9
+DA:67,18
+DA:68,0
+DA:69,0
+DA:71,15
+DA:72,0
DA:73,0
-DA:74,0
-DA:76,9
-DA:77,9
-DA:78,18
-DA:79,0
-DA:80,0
-DA:82,15
+DA:75,9
+DA:76,3
+DA:77,3
+DA:80,6
+DA:81,6
DA:83,0
-DA:84,0
-DA:86,9
+DA:85,9
+DA:86,6
DA:87,3
DA:88,3
+DA:89,3
+DA:90,6
DA:91,6
-DA:92,6
-DA:94,0
-DA:96,9
-DA:97,6
-DA:98,3
-DA:99,3
-DA:100,3
-DA:101,6
-DA:102,6
-DA:104,11
-DA:105,33
-DA:106,22
-DA:107,44
-DA:108,33
-DA:109,11
+DA:93,11
+DA:94,33
+DA:95,22
+DA:96,44
+DA:97,33
+DA:98,11
+DA:100,60756
+DA:101,60756
+DA:102,0
+DA:103,0
+DA:105,60756
+DA:106,60756
+DA:107,121512
+DA:108,60756
+DA:109,60756
+DA:110,60756
DA:111,60756
-DA:112,60756
-DA:113,0
-DA:114,0
-DA:116,60756
-DA:117,60756
-DA:118,121512
-DA:119,60756
-DA:120,60756
-DA:121,60756
-DA:122,60756
-DA:123,30378
-DA:124,30378
+DA:112,30378
+DA:113,30378
+DA:116,2
+DA:117,4
+DA:118,2
+DA:119,2
+DA:120,0
+DA:122,6
+DA:123,2
+DA:124,4
+DA:125,4
+DA:126,6
DA:127,2
-DA:128,4
-DA:129,2
-DA:130,2
-DA:131,0
-DA:133,6
-DA:134,2
-DA:135,4
-DA:136,4
-DA:137,6
-DA:138,2
-DA:140,174
-DA:141,174
-DA:142,90
-DA:143,90
-DA:144,90
+DA:129,174
+DA:130,174
+DA:131,90
+DA:132,90
+DA:133,90
+DA:134,90
+DA:135,0
+DA:137,90
+DA:138,0
+DA:139,0
+DA:140,90
+DA:141,90
+DA:143,0
DA:145,90
-DA:146,0
+DA:146,45
DA:148,90
-DA:149,0
-DA:150,0
-DA:151,90
-DA:152,90
-DA:154,0
-DA:156,90
-DA:157,45
-DA:159,90
-DA:160,90
+DA:149,90
LH:99
LF:117
end_of_record
diff --git a/paper/diagram.png b/paper/diagram.png
new file mode 100644
index 0000000..df2b2bd
Binary files /dev/null and b/paper/diagram.png differ
diff --git a/paper/paper.bib b/paper/paper.bib
new file mode 100644
index 0000000..e70eea4
--- /dev/null
+++ b/paper/paper.bib
@@ -0,0 +1,37 @@
+@misc{TaylorSeries,
+ author = {Benet, Luis and Sanders, David P.},
+ title = {A {Julia} package for {Taylor} expansions in one or
+ more independent variables},
+ howpublished = {\url{https://github.com/JuliaDiff/TaylorSeries.jl}},
+ doi = {10.5281/zenodo.2557003},
+ year = {2014},
+}
+
+@article{julia,
+ author = {Jeff Bezanson and Alan Edelman and Stefan Karpinski and Viral B. Shah},
+ title = {Julia: A Fresh Approach to Numerical Computing},
+ journal = {SIAM Review},
+ volume = {59},
+ number = {1},
+ pages = {65-98},
+ year = {2017},
+ doi = {10.1137/141000671},
+}
+@misc{MacroTools,
+ author = {Mike Innes},
+ title = {MacroTools.jl},
+ url = {https://github.com/FluxML/MacroTools.jl},
+ year = {2015},
+}
+@misc{SymEngine,
+ author = {Development Team},
+ title = {SymEngine.jl},
+ url = {https://github.com/symengine/SymEngine.jl},
+ year = {2015},
+}
+@misc{QSS,
+ author = {Pietro, F. and Migoni, G. and Kofman, E.},
+ title = {Improving linearly implicit quantized state system methods},
+ journal = {Simulation: Transactions of the Society for Modeling and Simulation International},
+ year = {2019},
+}
diff --git a/paper/paper.md b/paper/paper.md
new file mode 100644
index 0000000..77c89ca
--- /dev/null
+++ b/paper/paper.md
@@ -0,0 +1,74 @@
+---
+title: 'QuantizedSystemSolver: A discontinuous ODE system solver for Julia.'
+tags:
+ - Quantised State system
+ - Discontinuities
+ - Stiff Ordinary Differential Equations
+authors:
+ - name: Elmongi Elbellili
+ orcid: 0000-0003-1230-5488
+ affiliation: "1, 2"
+ - name: Daan Huybrechs
+ orcid: 0000-0002-0536-2647
+ affiliation: 2
+ - name: Ben Lauwens
+ orcid: 0000-0003-0761-6265
+ affiliation: 1
+affiliations:
+ - name: Royal Military Academy, Brussels, Belgium
+ index: 1
+ - name: Kuleuven, Leuven, Belgium
+ index: 2
+date: August 2024
+bibliography: paper.bib
+---
+
+# Summary
+The growing intricacy of contemporary engineering systems, typically reduced to stiff differential equations with events, poses a difficulty in digitally simulating them using traditional numerical integration techniques. The Quantized State System (QSS) and the Linearly Implicit Quantized State System (LIQSS) are different methods for tackling such problems. It is an approach that builds the solution by updating the system variables independently as opposed to classic integration methods that update all the system variables every step. [@QSS].
+
+# Statement of need
+Traditional solvers are challenged by frequent discontinuities where the state of the system abruptly changes
+at specific points or intervals. They struggle to accurately capture the dynamics around discontinuities. They either undergo expensive iterations to pinpoint exact discontinuity instances or resort to interpolating their locations, resulting in unreliable outcomes.
+Written in the easy-to-learn Julia language [Julia programming language](https://julialang.org) [@julia],
+and taking advantage of its features such as multiple dispatch and metaprogramming, the QuantizedSystemSolver.jl is a solver that aims to efficiently solve a set of Ordinary differential Equations with a set of events via implementing the QSS and LIQSS methods. It is the first such tool to be published in the Julia ecosystem.
+
+# Package description
+While the package is optimized to be fast, extensibility is not compromised. It is divided into 3 entities that can be extended separately: Problem, Algorithm, and Solution. The rest of the code is to create these entities and glue them together as shown in the Figure. The API was designed to provide an easier way to handle events than the approach provided by the classic integration solvers. Inside an NLodeProblem function, the user may introduce any parameters, variables, equations, and events:
+
+```
+odeprob = NLodeProblem(quote
+ name
+ parameters
+ discrete and continous variables
+ helper expressions
+ differential equations
+ if-statments for events)
+```
+
+The output of this function is an object of type Problem, and it is passed to the solve function along any other configuration arguments such as the algorithm type, the time span and the tolerance. The solve function dispatches on the given algorithm and start the numerical integration.
+
+```
+tspan = (0.0, 0.001)
+sol= solve(odeprob,nmliqss2(),tspan,abstol=1e-4,reltol=1e-3)
+```
+
+A the end, a solution object is produced that can be queried and plotted.
+```
+sol(0.0005,idxs=2)
+sol.stats
+plot(sol)
+```
+
+
+
+In addition, the package contains several shared helper functions used during the integration process by the algorithm such as the scheduler that organizes which variable of the system to update at any specific time of the simulation.
+The solver uses other packages such as [`MacroTools.jl `]( https://github.com/FluxML/MacroTools.jl)[@MacroTools] for user-code parsing, [`SymEngine.jl `]( https://github.com/symengine/SymEngine.jl)[@SymEngine] for Jacobian computation and dependencies extraction, and a modified [`TaylorSeries.jl`](https://github.com/JuliaDiff/TaylorSeries.jl/)[@TaylorSeries] that uses caching to obtain free Taylor variable operations as the current version of TaylorSeries creates a heap allocated object for every operation. The approximation through Taylor variables transforms any complicated equations to polynomials, which makes root finding cheaper, which the QSS methods relies heavily on it.
+
+In conclusion, the package offers extensive functionality to facilitate practical research tasks, all while being lightweight and well documented to be easily used by researchers and students to efficiently model various dynamical systems with discontinuities, as well as further study and improve the newly developed QSS methods. In Fact, Anyone that has a different problem type that can be handled by a specific QSS algorithm that outputs a different solution, they can define their types as subclasses of the abstract types (Problem, Algorithm, and Solution), and use the solver as usual.
+
+
+
+# Acknowledgements
+This research has received no external funding.
+
+# References
diff --git a/src/Common/QSS_data.jl b/src/Common/QSS_data.jl
index 6cd46b2..8eefbb6 100644
--- a/src/Common/QSS_data.jl
+++ b/src/Common/QSS_data.jl
@@ -19,7 +19,7 @@ struct CommonQSS_data{Z}
dQmin ::Float64
dQrel ::Float64
maxErr ::Float64
- maxStepsAllowed ::Int
+ maxiters ::Int
savedTimes :: Vector{Vector{Float64}}
savedVars:: Vector{Vector{Float64}}
diff --git a/src/Common/Solution.jl b/src/Common/Solution.jl
index 8d3b052..02531b8 100644
--- a/src/Common/Solution.jl
+++ b/src/Common/Solution.jl
@@ -1,3 +1,10 @@
+struct Stats
+ totalSteps::Int
+ simulStepCount::Int
+ evCount::Int
+ numSteps ::Vector{Int}
+end
+
"""LightSol{T,O}
A struct that holds the solution of a system of ODEs. It has the following fields:\n
- size: The number of continuous variables T\n
@@ -22,28 +29,34 @@ struct LightSol{T,O}<:Sol{T,O}
algName::String
sysName::String
absQ::Float64
- totalSteps::Int
- simulStepCount::Int
- evCount::Int
- numSteps ::Vector{Int}
+ stats::Stats
ft::Float64
end
-@inline function createSol(::Val{T},::Val{O}, savedTimes:: Vector{Vector{Float64}},savedVars :: Vector{Vector{Float64}},solver::String,nameof_F::String,absQ::Float64,totalSteps::Int#= ,stepsaftersimul::Int =#,simulStepCount::Int,evCount::Int,numSteps ::Vector{Int},ft::Float64#= ,simulStepsVals :: Vector{Vector{Float64}}, simulStepsDers :: Vector{Vector{Float64}} ,simulStepsTimes :: Vector{Vector{Float64}} =#)where {T,O}
+
+@inline function createSol(::Val{T},::Val{O}, savedTimes:: Vector{Vector{Float64}},savedVars :: Vector{Vector{Float64}},solver::String,nameof_F::String,absQ::Float64,stats::Stats,ft::Float64)where {T,O}
# println("light")
- sol=LightSol(Val(T),Val(O),savedTimes, savedVars,solver,nameof_F,absQ,totalSteps#= ,stepsaftersimul =#,simulStepCount,evCount,numSteps,ft#= ,simulStepsVals,simulStepsDers,simulStepsTimes =#)
+ sol=LightSol(Val(T),Val(O),savedTimes, savedVars,solver,nameof_F,absQ,stats,ft#= ,simulStepsVals,simulStepsDers,simulStepsTimes =#)
end
-function getindex(s::Sol, i::Int64)
+function getindex(s::Sol, i::Int64)#helper for calling savedTimes & savedVars
if i==1
return s.savedTimes
- elseif i==2
- return s.savedVars
+ #elseif i==2
else
- error("sol has 2 attributes: time and states")
+ return s.savedVars
+ #= else
+ error("sol has 2 attributes: time and states") =#
end
end
@inline function evaluateSol(sol::Sol{T,O},index::Int,t::Float64)where {T,O}
(t>sol.ft) && error("given point is outside the solution range! Verify where you want to evaluate the solution")
+ if index==0
+ eval_All=Vector{Float64}(undef, T)
+ for index=1:T
+ eval_All[index] =evaluateSol(sol,index,t)
+ end
+ eval_All
+ else
x=sol[2][index][end]
#integratorCache=Taylor0(zeros(O+1),O)
@@ -64,6 +77,7 @@ end
# println("3rd case")
return x #if var never changed then return init cond or if t>lastSavedTime for this var then return last value
end
+end
function solInterpolated(sol::Sol{T,O},step::Float64)where {T,O}
interpTimes=Float64[]
allInterpTimes=Vector{Vector{Float64}}(undef, T)
@@ -76,11 +90,11 @@ function solInterpolated(sol::Sol{T,O},step::Float64)where {T,O}
push!(interpTimes,sol.ft)
numInterpPoints=length(interpTimes)
interpValues=nothing
- if sol isa LightSol
+ # if sol isa LightSol
interpValues=Vector{Vector{Float64}}(undef, T)
- end
+ #end
for index=1:T
- interpValues[index]=[]
+ interpValues[index]=Float64[]
push!(interpValues[index],sol[2][index][1]) #1st element is the init cond (true value)
for i=2:numInterpPoints-1
push!(interpValues[index],evaluateSol(sol,index,interpTimes[i]))
@@ -88,6 +102,43 @@ function solInterpolated(sol::Sol{T,O},step::Float64)where {T,O}
push!(interpValues[index],sol[2][index][end]) #last pt @ft
allInterpTimes[index]=interpTimes
end
- createSol(Val(T),Val(O),allInterpTimes,interpValues,sol.algName,sol.sysName,sol.absQ,sol.totalSteps#= ,sol.stepsaftersimul =#,sol.simulStepCount,sol.evCount,sol.numSteps,sol.ft#= ,sol.simulStepsVals,sol.simulStepsDers,sol.simulStepsVals =#)
+ createSol(Val(T),Val(O),allInterpTimes,interpValues,sol.algName,sol.sysName,sol.absQ,sol.stats,sol.ft#= ,sol.simulStepsVals,sol.simulStepsDers,sol.simulStepsVals =#)
+end
+function solInterpolated(sol::Sol{T,O},index::Int,step::Float64)where {T,O}
+ interpTimes=Float64[]
+ allInterpTimes=Vector{Vector{Float64}}(undef, 1)
+ t=0.0 #later can change to init_time which could be diff than zero
+ push!(interpTimes,t)
+ while t+step3*abs(ẋj) || (dxj*ẋj)<0.0
+ if abs(dxi)>3*abs(ẋi) || abs(dxi)*3(abs(dxj+ẋj)/2)
+ if abs(dxi-dxithrow)>(abs(dxi+dxithrow)/2)
+ iscycle=true
+ end
+ end =#
+ ########condition:Union i union
+ if (abs(dxj)*33*abs(ẋj) || (dxj*ẋj)<0.0)
+ cancelCriteria=1e-6*quani
+ if abs(dxi)>3*abs(dxithrow) || abs(dxi)*310*abs(ẋi) || abs(dxi)*10 2*quani || abs(qj - xj) > 2*quanj) #checking qi-xi is not needed since firstguess just made it less than delta
+ h1 = (abs(quani / ẋi));h2 = (abs(quanj / ẋj));
+ h=min(h1,h2)
+ h_two=h
+ Δ=(1-h*aii)*(1-h*ajj)-h*h*aij*aji
+ #= if Δ==0
+ Δ=1e-12
+ end =#
+ qi = ((1-h*ajj)*(xi+h*uij)+h*aij*(xj+h*uji))/Δ
+ qj = ((1-h*aii)*(xj+h*uji)+h*aji*(xi+h*uij))/Δ
+ end
+ maxIter=1000
+ while (abs(qi - xi) > 2*quani || abs(qj - xj) > 2*quanj) && (maxIter>0)
+ maxIter-=1
+ h1 = h * 2*(0.99*quani / abs(qi - xi));
+ #= Δtemp=(1-h1*aii)*(1-h1*ajj)-h1*h1*aij*aji
+ qitemp = ((1-h1*ajj)*(xi+h1*uij)+h1*aij*(xj+h1*uji))/Δtemp =#
+ h2 = h * 2*(0.99*quanj / abs(qj - xj));
+ #= Δtemp=(1-h2*aii)*(1-h2*ajj)-h2*h2*aij*aji
+ qjtemp = ((1-h2*ajj)*(xi+h2*uij)+h2*aij*(xj+h2*uji))/Δtemp
+ if abs(qitemp - xi) > 1*quani || abs(qjtemp - xj) > 1*quanj
+ println("regle de croix did not work")
+ end =#
+ h=min(h1,h2)
+ h_three=h
+ Δ=(1-h*aii)*(1-h*ajj)-h*h*aij*aji
+ #= if Δ==0
+ Δ=1e-12
+ println("delta liqss1 simulupdate==0")
+ end =#
+ if Δ!=0
+ qi = ((1-h*ajj)*(xi+h*uij)+h*aij*(xj+h*uji))/Δ
+ qj = ((1-h*aii)*(xj+h*uji)+h*aji*(xi+h*uij))/Δ
+ end
+
+ if maxIter < 1 println("maxiter of updateQ = ",maxIter) end
+ end
+ if maxIter < 1
+ return false
+ end
+
+ q[index][0]=qi# store back helper vars
+ trackSimul[1]+=1 # do not have to recomputeNext if qi never changed
+ q[j][0]=qj
+ #= push!(simulqxiVals,abs(qi-xi))
+ push!(simulqxjVals,abs(qj-xj))
+ push!(simuldeltaiVals,quani)
+ push!(simuldeltajVals,quanj) =#
+ tq[j]=simt
+ end #end second dependecy check
+
+ return iscycle
+end =#
+
+
+#iters from 0===>1e-9
+function nmisCycle_and_simulUpdate(cacheRootsi::Vector{Float64},cacheRootsj::Vector{Float64},acceptedi::Vector{Vector{Float64}},acceptedj::Vector{Vector{Float64}},aij::Float64,aji::Float64,respp::Ptr{Float64}, pp::Ptr{NTuple{2,Float64}},trackSimul,::Val{1},index::Int,j::Int,dirI::Float64,dti::Float64, x::Vector{Taylor0},q::Vector{Taylor0}, quantum::Vector{Float64},exactA::Function,d::Vector{Float64},cacheA::MVector{1,Float64},dxaux::Vector{MVector{1,Float64}},qaux::Vector{MVector{1,Float64}},tx::Vector{Float64},tq::Vector{Float64},simt::Float64,ft::Float64)
+
+ cacheA[1]=0.0;exactA(q,d,cacheA,index,index,simt)
+ aii=cacheA[1]
+ cacheA[1]=0.0;exactA(q,d,cacheA,j,j,simt)
+ ajj=cacheA[1]
+ #= exactA(q,cacheA,index,j)
+ aij=cacheA[1]
+ exactA(q,cacheA,j,index)
+ aji=cacheA[1] =#
+
+
+ xi=x[index][0];xj=x[j][0];ẋi=x[index][1];ẋj=x[j][1]
+ qi=q[index][0];qj=q[j][0]
+ quanj=quantum[j];quani=quantum[index]
+ #qaux[j][1]=qj;
+ elapsed = simt - tx[j];x[j][0]= xj+elapsed*ẋj;
+
+ xj=x[j][0]
+ tx[j]=simt
+
+ qiminus=qaux[index][1]
+ #ujj=ẋj-ajj*qj
+ uji=ẋj-ajj*qj-aji*qiminus
+ #uii=dxaux[index][1]-aii*qaux[index][1]
+ uij=dxaux[index][1]-aii*qiminus-aij*qj
+ iscycle=false
+ dxj=aji*qi+ajj*qj+uji #only future qi #emulate fj
+
+
+
+ dxithrow=aii*qi+aij*qj+uij #only future qi
+
+ qjplus=xj+sign(dxj)*quanj #emulate updateQ(j)...
+
+ dxi=aii*qi+aij*qjplus+uij #both future qi & qj #emulate fi
+ #dxj2=ajj*qjplus+aji*qi+uji
+
+
+ #= if abs(dxithrow)<1e-15 && dxithrow!=0.0
+ dxithrow=0.0
+ end =#
+
+ ########condition:Union
+ #= if abs(dxj)*33*abs(ẋj) || (dxj*ẋj)<0.0
+ if abs(dxi)>3*abs(ẋi) || abs(dxi)*3(abs(dxj+ẋj)/2)
+ if abs(dxi-dxithrow)>(abs(dxi+dxithrow)/2)
+ iscycle=true
+ end
+ end =#
+ ########condition:Union i union
+ if (abs(dxj)*33*abs(ẋj) || (dxj*ẋj)<0.0)
+ cancelCriteria=1e-6*quani
+ if abs(dxi)>3*abs(dxithrow) || abs(dxi)*310*abs(ẋi) || abs(dxi)*10 1*quani || abs(qj - xj) > 1*quanj) #checking qi-xi is not needed since firstguess just made it less than delta
+ h1 = (abs(quani / ẋi));h2 = (abs(quanj / ẋj));
+ h=min(h1,h2)
+
+ Δ=(1-h*aii)*(1-h*ajj)-h*h*aij*aji
+ #= if Δ==0
+ Δ=1e-12
+ end =#
+ qi = ((1-h*ajj)*(xi+h*uij)+h*aij*(xj+h*uji))/Δ
+ qj = ((1-h*aii)*(xj+h*uji)+h*aji*(xi+h*uij))/Δ
+ end
+ maxIter=1000
+ while (abs(qi - xi) < 1*quani || abs(qj - xj) < 1*quanj) && (maxIter>0) # maxiters here is better than maxiter in normal iters because here h is acceptable when we reach maxiter
+ maxIter-=1
+
+ h1 = h * (1.0*quani / abs(qi - xi));
+ #= Δtemp=(1-h1*aii)*(1-h1*ajj)-h1*h1*aij*aji
+ qitemp = ((1-h1*ajj)*(xi+h1*uij)+h1*aij*(xj+h1*uji))/Δtemp =#
+ h2 = h * (1.0*quanj / abs(qj - xj));
+ #= Δtemp=(1-h2*aii)*(1-h2*ajj)-h2*h2*aij*aji
+ qjtemp = ((1-h2*ajj)*(xi+h2*uij)+h2*aij*(xj+h2*uji))/Δtemp
+ if abs(qitemp - xi) > 1*quani || abs(qjtemp - xj) > 1*quanj
+ println("regle de croix did not work")
+ end =#
+ h=max(h1,h2)
+ # h_three=h
+ Δ=(1-h*aii)*(1-h*ajj)-h*h*aij*aji
+ #= if Δ==0
+ Δ=1e-12
+ println("delta liqss1 simulupdate==0")
+ end =#
+ if Δ!=0
+ qi = ((1-h*ajj)*(xi+h*uij)+h*aij*(xj+h*uji))/Δ
+ qj = ((1-h*aii)*(xj+h*uji)+h*aji*(xi+h*uij))/Δ
+ end
+
+ if maxIter < 1 println("maxiter of updateQ = ",maxIter) end
+ end
+
+ #h=h_ # go back to last h if exit the loop
+ qi = ((1-h*ajj)*(xi+h*uij)+h*aij*(xj+h*uji))/Δ
+ qj = ((1-h*aii)*(xj+h*uji)+h*aji*(xi+h*uij))/Δ
+
+ if maxIter < 1 # this is not needed in iters from 0
+ return false
+ end
+
+ q[index][0]=qi# store back helper vars
+ trackSimul[1]+=1 # do not have to recomputeNext if qi never changed
+ q[j][0]=qj
+ #= push!(simulqxiVals,abs(qi-xi))
+ push!(simulqxjVals,abs(qj-xj))
+ push!(simuldeltaiVals,quani)
+ push!(simuldeltajVals,quanj) =#
+ tq[j]=simt
+ end #end second dependecy check
+
+ return iscycle
+end
diff --git a/test/Get_Code_Coverage.jl b/test/Get_Code_Coverage.jl
index 8eff685..6a3621a 100644
--- a/test/Get_Code_Coverage.jl
+++ b/test/Get_Code_Coverage.jl
@@ -3,11 +3,11 @@
Pkg.test("QuantizedSystemSolver"; coverage=true) # run this then comment and run next code
=#
# get lcov from .cov files (summary)
- using Coverage
+#= using Coverage
coverage = process_folder()
open("lcov.info", "w") do io
LCOV.write(io, coverage)
-end;
+end; =#
#clean up the folder when not needed
#= using Coverage
diff --git a/test/Tyson.jl b/test/Tyson.jl
new file mode 100644
index 0000000..f14c12b
--- /dev/null
+++ b/test/Tyson.jl
@@ -0,0 +1,48 @@
+
+using QuantizedSystemSolver
+#using XLSX
+using BenchmarkTools
+#using BSON
+#using TimerOutputs
+#using Plots
+function test(case,solvr)
+ absTol=1e-5
+ relTol=1e-2
+
+
+
+
+ # BSON.@load "formalA2/ref_bson/solVect_Tyson_Rodas5Pe-12.bson" solRodas5PVectorTyson
+ odeprob = NLodeProblem(quote
+ name=(tyson,)
+ u = [0.0,0.75,0.25,0.0,0.0,0.0]
+ du[1] = u[4]-1e6*u[1]+1e3*u[2]
+ du[2] =-200.0*u[2]*u[5]+1e6*u[1]-1e3*u[2]
+ du[3] = 200.0*u[2]*u[5]-u[3]*(0.018+180.0*(u[4]/(u[1]+u[2]+u[3]+u[4]))^2)
+ du[4] =u[3]*(0.018+180.0*(u[4]/(u[1]+u[2]+u[3]+u[4]))^2)-u[4]
+ du[5] = 0.015-200.0*u[2]*u[5]
+ du[6] =u[4]-0.6*u[6]
+ end )
+ println("start tyson solving")
+ timenmliqss=0.0
+ tspan=(0.0,25.0)
+ solnmliqss=solve(odeprob,solvr,abstol=absTol,saveat=0.01,reltol=relTol,tspan,maxiters=10000#= ,maxErr=100*relTol =#)
+ save_Sol(solnmliqss)
+ # @show solnmliqss.totalSteps
+ #save_Sol(solnmliqss,note="cancelcriteria-3_xi10dxi"#= xlims=(10.3778695,14.5789) =#)
+
+
+ # save_Sol(solnmliqss,1,note="x1 intrval13 ",xlims=(4.0,4.38),ylims=(0.0007,0.000723))
+ # solnmliqssInterp=solInterpolated(solnmliqss,0.01)
+ err3=0.0
+ # err3=getAverageErrorByRefs(solRodas5PVectorTyson,solnmliqssInterp)
+
+ #timenmliqss=@belapsed solve($odeprob,$solvr,abstol=$absTol,saveat=0.01,reltol=$relTol,$tspan#= ,maxErr=1000*$relTol =#)
+ resnmliqss11E_2= ("$(solnmliqss.algName)",relTol,err3,solnmliqss.stats.totalSteps,solnmliqss.stats.simulStepCount,timenmliqss)
+ @show resnmliqss11E_2
+
+end
+case="order1_"
+
+println("compareBounds")
+test(case,nmliqss1()) #compareBounds
diff --git a/test/bbal.jl b/test/bbal.jl
new file mode 100644
index 0000000..7ef024c
--- /dev/null
+++ b/test/bbal.jl
@@ -0,0 +1,23 @@
+using QuantizedSystemSolver
+function test()
+ odeprob = NLodeProblem(quote
+ name=(sysd0,)
+ u = [3.0,20.0]
+ discrete=[0.0]
+ du[1] = u[2]
+ du[2] = -9.8+discrete[1]*u[1]
+ if -u[1]>0.0
+ u[2]=-0.9*u[2]
+ end
+ end)
+ tspan=(0.0,20.0)
+ # sol=solve(odeprob,qss2(),tspan)
+ sol= solve(odeprob,nmliqss1(),tspan,abstol=1e-3,reltol=1e-2)
+ @show sol.stats.totalSteps
+ save_Sol(sol)
+
+ #getAverageErrorByRefs(solRef::Vector{Any},solmliqss::Sol{T,O})
+end
+test()
+
+
diff --git a/test/diffEQrailgunElectrical_AllCircuit_12.jl b/test/diffEQrailgunElectrical_AllCircuit_12.jl
new file mode 100644
index 0000000..b5c7a9f
--- /dev/null
+++ b/test/diffEQrailgunElectrical_AllCircuit_12.jl
@@ -0,0 +1,259 @@
+using QuantizedSystemSolver
+using BenchmarkTools
+using TimerOutputs
+reset_timer!()
+function test()
+
+ odeprob = @NLodeProblem begin
+ name=(RGElectrical12,)
+ ROn = 1e-5;ROff = 1e1;
+ # Lpr = 4.2e-9#0.48*1e-6#0.45 * 1e-6
+L1 = 0.6e-4 #28.0*1e-6#0.6*1e-6
+#L2 = 4.0e-6;L3 = 1.1*1e-6
+#L23=5.1e-6
+
+R1= 4.0e-3; #R2 = 0.28e-3;R3 = 3.6e-3
+
+C = 3.08e-3#3.08*1e-3
+
+m=0.12
+#= γ = 50.0e6; w = 7.5#25.0*1e-3#15.0*1e-3
+μ = 4.0*3.14*1e-7 =#
+#rr=manualIntg*sqrt(μ/(3.14*γ))/w=manualIntg*coef =manualIntg* 1.1925695879998878e-8
+Rpr0 = 2.5e-6
+FNmec = 680.0; #α = 0.15
+ t0=1.0
+ discrete = [1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,0.0,1e-3]
+ rd1=discrete[1]; rs1=discrete[2]; operate1=discrete[3]; charge1=discrete[4]
+ rd2=discrete[5]; rs2=discrete[6]; operate2=discrete[7]; charge2=discrete[8];
+ rd3=discrete[9]; rs3=discrete[10];operate3=discrete[11]; charge3=discrete[12];
+ rd4=discrete[13];rs4=discrete[14];operate4=discrete[15]; charge4=discrete[16];
+ rd5=discrete[17];rs5=discrete[18];operate5=discrete[19]; charge5=discrete[20];
+ rd6=discrete[21];rs6=discrete[22];operate6=discrete[23]; charge6=discrete[24];
+ rd7=discrete[25];rs7=discrete[26];operate7=discrete[27]; charge7=discrete[28];
+ rd8=discrete[29];rs8=discrete[30];operate8=discrete[31]; charge8=discrete[32];
+ rd9=discrete[33];rs9=discrete[34];operate9=discrete[35]; charge9=discrete[36];
+ rd10=discrete[37];rs10=discrete[38];operate10=discrete[39]; charge10=discrete[40];
+ rd11=discrete[41];rs11=discrete[42];operate11=discrete[43]; charge11=discrete[44];
+ rd12=discrete[45];rs12=discrete[46];operate12=discrete[47]; charge12=discrete[48];
+ manualIntg=discrete[49]; nextT=discrete[50];
+ u = [0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,0.0]
+ is1=u[1] ;uc1=u[2]; il1=u[3] ;is2=u[4] ;uc2=u[5]; il2=u[6] ;is3=u[7] ;uc3=u[8]; il3=u[9] ;is4=u[10] ;uc4=u[11]; il4=u[12];is5=u[13] ;uc5=u[14]; il5=u[15] ;is6=u[16] ;uc6=u[17]; il6=u[18];is7=u[19] ;uc7=u[20]; il7=u[21] ;is8=u[22] ;uc8=u[23]; il8=u[24] ;is9=u[25] ;uc9=u[26]; il9=u[27] ;is10=u[28] ;uc10=u[29]; il10=u[30];is11=u[31] ;uc11=u[32]; il11=u[33];is12=u[34] ;uc12=u[35]; il12=u[36] ;x=u[37]; v=u[38]
+ #α=LR+Lpr;
+ # RR=4.0e-5
+ rr=manualIntg*1.1925695879998878e-8
+ rpr=Rpr0*(sqrt(t0/(t+1e-4))+(t)^16)/(1.0+(t)^16)
+
+ rrpp=(rpr + rr + 4.53e-7v)
+
+ #= rrpp=rrpp
+ x1=(132.97872599999997 + 35.34758999999999x)
+ x2=(-0.10924199999999998 - 11.782529999999998x)
+ x3=678.7486367999996 + 240.3636119999999x - 8.881784197001252e-16(x^2)
+ rd_11=rd_11;rd22=rd22;rd33=rd33;rd44=rd44 =#
+ Il=il1+il2-0.0+il3+il4+il5+il6-0.0+il7+il8+il9+il10-0.0+il11+il12
+ uf=0.1+0.2*exp(-v/100)
+ F=0.5*0.453e-6*Il*Il*(1-uf*0.124)-uf*FNmec
+ x4=(5.1(5.150399999999999 + 5.436000000000001x))
+ x5=(0.0042 + 0.453x)
+
+ rd_11=(-0.00388 - rd1);rd22=(-0.00388 - rd2);rd33=(0.00388 + rd3);rd44=(-0.00388 - rd4); rd55=(-0.00388 - rd5);rd66=(0.00388 + rd6);rd77=(-0.00388 - rd7);rd88=(-0.00388 - rd8); rd99=(0.00388 + rd9);rd100=(-0.00388 - rd10);rd110=(-0.00388 - rd11);rd120=(-0.00388 - rd12)
+ isid=il1*rd_11 +is1*rd1+il2*rd22+is2*rd2-il3*rd33+is3*rd3+il4*rd44+is4*rd4+il5*rd55 + is5*rd5-il6*rd66+is6*rd6+il7*rd77+is7*rd7+il8*rd88+is8*rd8-il9*rd99+is9*rd9+il10*rd100+is10*rd10+il11*rd110+is11*rd11+il12*rd120+is12*rd12
+
+ du[1] =((-(R1+rs1+rd1)*is1+rd1*il1+uc1)/L1)*operate1
+ du[2]=(-is1/C)*charge1*operate1
+ du[3]=operate1*1e6* (-(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il1*rd_11+is1*rd1)*(0.19607843137254904 ))
+
+
+ du[4] =((-(R1+rs2+rd2)*is2+rd2*il2+uc2)/L1)*operate2
+ du[5]=(-is2/C)*charge2*operate2
+ du[6]=operate2*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il2*rd22+is2*rd2)*(0.19607843137254904 ))
+
+ du[7] =((-(R1+rs3+rd3)*is3+rd3*il3+uc3)/L1)*operate3#*charge3
+ du[8]=((-is3/C)*charge3)*operate3
+ du[9]=operate3*1e6* (-(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il3*rd33+is3*rd3)*(0.19607843137254904 ))
+
+ du[10] =((-(R1+rs4+rd4)*is4+rd4*il4+uc4)/L1)*operate4#*charge3
+ du[11]=((-is4/C)*charge4)*operate4
+ du[12]=operate4*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il4*rd44+is4*rd4)*(0.19607843137254904 ))
+
+ du[13] =((-(R1+rs5+rd5)*is5+rd5*il5+uc5)/L1)*operate5#*charge3
+ du[14]=((-is5/C)*charge5)*operate5
+ du[15]=operate5*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il5*rd55+is5*rd5)*(0.19607843137254904 ))
+
+ du[16] =((-(R1+rs6+rd6)*is6+rd6*il6+uc6)/L1)*operate6#*charge3
+ du[17]=((-is6/C)*charge6)*operate6
+ du[18]=operate6*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il6*rd66+is6*rd6)*(0.19607843137254904 ))
+
+ du[19] =((-(R1+rs7+rd7)*is7+rd7*il7+uc7)/L1)*operate7
+ du[20]=(-is7/C)*charge7*operate7
+ du[21]=operate7*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il7*rd77+is7*rd7)*(0.19607843137254904 ))
+
+ du[22] =((-(R1+rs8+rd8)*is8+rd8*il8+uc8)/L1)*operate8
+ du[23]=(-is8/C)*charge8*operate8
+ du[24]=operate8*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il8*rd88+is8*rd8)*(0.19607843137254904 ))
+
+ du[25] =((-(R1+rs9+rd9)*is9+rd9*il9+uc9)/L1)*operate9
+ du[26]=(-is9/C)*charge9*operate9
+ du[27]=operate9*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il9*rd99+is9*rd9)*(0.19607843137254904 ))
+
+ du[28] =((-(R1+rs10+rd10)*is10+rd10*il10+uc10)/L1)*operate10
+ du[29]=(-is10/C)*charge10*operate10
+ du[30]=operate10*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il10*rd100+is10*rd10)*(0.19607843137254904 ))
+
+ du[31] =((-(R1+rs11+rd11)*is11+rd11*il11+uc11)/L1)*operate11
+ du[32]=(-is11/C)*charge11*operate11
+ du[33]=operate11*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il11*rd110+is11*rd11)*(0.19607843137254904 ))
+
+ du[34] =((-(R1+rs12+rd12)*is12+rd12*il12+uc12)/L1)*operate12
+ du[35]=(-is12/C)*charge12*operate12
+ du[36]=operate12*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il12*rd120+is12*rd12)*(0.19607843137254904 ))
+
+
+ du[37]=v # v=u[26]
+ du[38]=F/m
+
+ if t-0.00019>0.0
+ operate4=1.0
+ end
+
+ if t-0.00038>0.0
+ operate5=1.0
+ end
+ if t-0.00057>0.0
+ operate6=1.0
+ end
+ if t-0.00095>0.0
+ operate7=1.0
+ end
+ if t-0.00115>0.0
+ operate8=1.0
+ end
+ if t-0.00142>0.0
+ operate9=1.0
+ end
+ if t-0.00175>0.0
+ operate10=1.0
+ end
+ if t-0.00195>0.0
+ operate11=1.0
+ end
+ if t-0.00210>0.0
+ operate12=1.0
+ end
+
+ if -(uc1)>0.0
+ charge1=0.0 # rs off not needed since charge=0
+ rs1=ROff;
+ rd1=ROn;
+ uc1=0.0
+ #= uc1=0.0;
+ is1=0.0 =#
+ end
+ if -(uc2)>0.0
+ charge2=0.0
+ rs2=ROff;
+ rd2=ROn;
+ #@show rd2
+ uc2=0.0
+
+ end
+
+ if -(uc3)>0.0
+ charge3=0.0
+ rs3=ROff;
+ rd3=ROn;
+ uc3=0.0
+
+ end
+ if -(uc4)>0.0
+ charge4=0.0
+ rs4=ROff;
+ rd4=ROn;
+ uc4=0.0
+
+ end
+
+ if -(uc5)>0.0
+ charge5=0.0
+ rs5=ROff;
+ rd5=ROn;
+ uc5=0.0
+
+ end
+ if -(uc6)>0.0
+ charge6=0.0
+ rs6=ROff;
+ rd6=ROn;
+ uc6=0.0
+
+ end
+ if -(uc7)>0.0
+ charge7=0.0
+ rs7=ROff;
+ rd7=ROn;
+ uc7=0.0
+
+ end
+ if -(uc8)>0.0
+ charge8=0.0
+ rs8=ROff;
+ rd8=ROn;
+ uc8=0.0
+
+ end
+ if -(uc9)>0.0
+ charge9=0.0
+ rs9=ROff;
+ rd9=ROn;
+ uc9=0.0
+ end
+ if -(uc10)>0.0
+ charge10=0.0
+ rs10=ROff;
+ rd10=ROn;
+ uc10=0.0
+ end
+ if -(uc11)>0.0
+ charge11=0.0
+ rs11=ROff;
+ rd11=ROn;
+ uc11=0.0
+ end
+ if -(uc12)>0.0
+ charge12=0.0
+ rs12=ROff;
+ rd12=ROn;
+ uc12=0.0
+ end
+
+ if t-nextT>0
+ manualIntg=manualIntg+v/sqrt(t[0])
+ nextT=nextT+1e-3
+
+ end
+
+ end
+
+
+ # @show odeprob.eqs
+ # println("start solving")
+ tspan = (0.0, 9.0e-3)
+ sol= solve(odeprob,nmliqss2(),tspan,abstol=1e-3,reltol=1e-2)
+ # sol= solve(odeprob,nmliqss1(),tspan,abstol=1e-3,reltol=1e-2)
+
+
+ # save_Sol(sol,9,note="z1",xlims=(0.0,0.2e-8),ylims=(-0.005,0.005))
+ save_SolSum(sol,3,6,9,12,15,18,21,24,27,30,33,36,interp=0.00001) #add interpl preference
+
+
+
+ save_Sol(sol,37)
+ save_Sol(sol,38)
+
+end
+#@time
+#test()
+#@btime
+test()
+print_timer()
diff --git a/test/diffEqsysB53.jl b/test/diffEqsysB53.jl
index e9de14d..9ac912e 100644
--- a/test/diffEqsysB53.jl
+++ b/test/diffEqsysB53.jl
@@ -13,14 +13,16 @@ function odeDiffEquPackage()
#= du[1] = t
du[2] =1.24*u[1]-0.01*u[2]+0.2 =#
- du[1] = t
+ #= du[1] = t
for k in 2:5
du[k]=(u[k]-u[k-1]) ;
- end
+ end =#
+ du[1] = -u[2]
+ du[2]=u[1]
end
- tspan = (0.0,5.0)
- u0= [1.0, 0.0,1.0, 0.0,1.0]
+ tspan = (0.0,6.0)
+ u0= [1.0, 0.0]
prob = ODEProblem(funcName,u0,tspan)
@@ -29,8 +31,8 @@ function odeDiffEquPackage()
solRosenbrock23 = solve(prob,Rosenbrock23(),saveat=0.01,abstol = absTol, reltol = relTol) #1.235 ms (1598 allocations: 235.42 KiB)
- p1=plot!(solRosenbrock23,marker=(:circle),markersize=2)
- savefig(p1, "plot_solRosenbrock23_N8.png")
+ p1=plot!(solRosenbrock23,idxs=(0,1,2))
+ savefig(p1, "plot_solRosenbrock23_vars0_1_2.png")
end
diff --git a/test/diffeqBUCK copy.jl b/test/diffeqBUCK copy.jl
new file mode 100644
index 0000000..602b62f
--- /dev/null
+++ b/test/diffeqBUCK copy.jl
@@ -0,0 +1,45 @@
+using DifferentialEquations
+using Plots
+function odeDiffEquPackage()
+ function f(du, u, p, t)
+ C = 1e-4; L = 1e-4; R = 10;U = 24.0; T = 1e-4; DC = 0.5; ROn = 1e-5;ROff = 1e5;
+ rd=p[1];rs=p[2];
+ il=u[1] ;uc=u[2]
+ id=(il*rs-U)/(rd+rs) # diode's current
+ du[1] =(-id*rd-uc)/L # inductor's current
+ du[2]=(il-uc/R)/C # capacitor's voltage
+ end
+ function condition(out, u, t, integrator) # Event when condition(out,u,t,integrator) == 0
+ out[1] = (t-p[3])
+ out[2] = (t-p[4]-0.5*1e-4)
+ out[3] = (p[5]*((u[1]*p[2]-24.0)/(p[1]+p[2]))+(1.0-p[5])*((u[1]*p[2]-24.0)*p[1]/(p[1]+p[2])))
+ out[4] = -(p[5]*((u[1]*p[2]-24.0)/(p[1]+p[2]))+(1.0-p[5])*((u[1]*p[2]-24.0)*p[1]/(p[1]+p[2])))
+ end
+
+ function affect!(integrator, idx)
+ if idx == 1
+ p[4]=p[3]
+ p[3]=p[3]+1e-4
+ p[2]=1e-5
+ elseif idx == 2
+ p[2]=1e5
+ elseif idx == 3
+ p[1]=1e-5
+ p[5]=1.0
+ elseif idx == 4
+ p[1]=1e5
+ p[5]=0.0
+ end
+ end
+
+
+ cbs = VectorContinuousCallback(condition, affect!, 4)
+ u0 = [0.0, 0.0]
+ tspan = (0.0, 0.0002)
+ p = [1e5,1e-5,1e-4,0.0,0.0,0.0]
+ prob = ODEProblem(f, u0, tspan, p)
+ sol = solve(prob, Rodas5P(), callback = cbs, reltol=1e-6,abstol=1e-9,dtmax=1e-12,maxiters=Int(1e12)#= dt = 1e-3, adaptive = false =#)
+ p1=plot!(sol);
+ savefig(p1, "Rodas5P()_vector_buck_ft0002_")
+ end
+ odeDiffEquPackage()
\ No newline at end of file
diff --git a/test/diffeqBUCK.jl b/test/diffeqBUCK.jl
new file mode 100644
index 0000000..b1ee6f1
--- /dev/null
+++ b/test/diffeqBUCK.jl
@@ -0,0 +1,49 @@
+using DifferentialEquations
+using Plots
+function odeDiffEquPackage()
+ function f(du, u, p, t)
+ C = 1e-4; L = 1e-4; R = 10;U = 24.0; T = 1e-4; DC = 0.5; ROn = 1e-5;ROff = 1e5;
+ rd=p[1];rs=p[2];
+ il=u[1] ;uc=u[2]
+ id=(il*rs-U)/(rd+rs) # diode's current
+ du[1] =(-id*rd-uc)/L # inductor's current
+ du[2]=(il-uc/R)/C # capacitor's voltage
+ end
+ function condition1( u, t, integrator)
+ (t-p[3])
+ end
+ function condition2( u, t, integrator)
+ (t-p[4]-0.5*1e-4)
+ end
+ function condition3( u, t, integrator)
+ (p[5]*((u[1]*p[2]-24.0)/(p[1]+p[2]))+(1.0-p[5])*((u[1]*p[2]-24.0)*p[1]/(p[1]+p[2])))
+ end
+ function affect1!(integrator)
+ p[4]=p[3]
+ p[3]=p[3]+1e-4
+ p[2]=1e-5
+ end
+ function affect2!(integrator)
+ p[2]=1e5
+ end
+ function affect3!(integrator)
+ p[1]=1e-5
+ p[5]=1.0
+ end
+ function affect33!(integrator)
+ p[1]=1e5
+ p[5]=0.0
+ end
+ cb1 = ContinuousCallback(condition1, affect1!,nothing; )
+ cb2 = ContinuousCallback(condition2, affect2!,nothing; )
+ cb3 = ContinuousCallback(condition3, affect3!,affect33!; )
+ cbs = CallbackSet(cb1, cb2,cb3)
+ u0 = [0.0, 0.0]
+ tspan = (0.0, 0.0002)
+ p = [1e5,1e-5,1e-4,0.0,0.0,0.0]
+ prob = ODEProblem(f, u0, tspan, p)
+ sol = solve(prob, ImplicitEuler(), callback = cbs, reltol=1e-2,abstol=1e-3#= ,dtmax=1e-12 =#,maxiters=Int(1e12)#= dt = 1e-3, adaptive = false =#)
+ p1=plot!(sol);
+ savefig(p1, "Rodas5P()_34_buck_ft0002_")
+ end
+ odeDiffEquPackage()
\ No newline at end of file
diff --git a/test/diffeqRG12.jl b/test/diffeqRG12.jl
new file mode 100644
index 0000000..9498b70
--- /dev/null
+++ b/test/diffeqRG12.jl
@@ -0,0 +1,239 @@
+using DifferentialEquations
+using Plots
+function odeDiffEquPackage()
+ function f(du, u, p, t)
+ ROn = 1e-5;ROff = 1e1;
+ # Lpr = 4.2e-9#0.48*1e-6#0.45 * 1e-6
+ L1 = 0.6e-4 #28.0*1e-6#0.6*1e-6
+ #L2 = 4.0e-6;L3 = 1.1*1e-6
+ #L23=5.1e-6
+
+ R1= 4.0e-3; #R2 = 0.28e-3;R3 = 3.6e-3
+
+ C = 3.08e-3#3.08*1e-3
+
+ m=0.12
+ #= γ = 50.0e6; w = 7.5#25.0*1e-3#15.0*1e-3
+ μ = 4.0*3.14*1e-7 =#
+ #rr=manualIntg*sqrt(μ/(3.14*γ))/w=manualIntg*coef =manualIntg* 1.1925695879998878e-8
+ Rpr0 = 2.5e-6
+ FNmec = 680.0; #α = 0.15
+ t0=1.0
+ p = [1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,0.0,1e-3]
+ rd1=p[1]; rs1=p[2]; operate1=p[3]; charge1=p[4]
+ rd2=p[5]; rs2=p[6]; operate2=p[7]; charge2=p[8];
+ rd3=p[9]; rs3=p[10];operate3=p[11]; charge3=p[12];
+ rd4=p[13];rs4=p[14];operate4=p[15]; charge4=p[16];
+ rd5=p[17];rs5=p[18];operate5=p[19]; charge5=p[20];
+ rd6=p[21];rs6=p[22];operate6=p[23]; charge6=p[24];
+ rd7=p[25];rs7=p[26];operate7=p[27]; charge7=p[28];
+ rd8=p[29];rs8=p[30];operate8=p[31]; charge8=p[32];
+ rd9=p[33];rs9=p[34];operate9=p[35]; charge9=p[36];
+ rd10=p[37];rs10=p[38];operate10=p[39]; charge10=p[40];
+ rd11=p[41];rs11=p[42];operate11=p[43]; charge11=p[44];
+ rd12=p[45];rs12=p[46];operate12=p[47]; charge12=p[48];
+ manualIntg=p[49]; nextT=p[50];
+ u = [0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,0.0]
+ is1=u[1] ;uc1=u[2]; il1=u[3] ;is2=u[4] ;uc2=u[5]; il2=u[6] ;is3=u[7] ;uc3=u[8]; il3=u[9] ;is4=u[10] ;uc4=u[11]; il4=u[12];is5=u[13] ;uc5=u[14]; il5=u[15] ;is6=u[16] ;uc6=u[17]; il6=u[18];is7=u[19] ;uc7=u[20]; il7=u[21] ;is8=u[22] ;uc8=u[23]; il8=u[24] ;is9=u[25] ;uc9=u[26]; il9=u[27] ;is10=u[28] ;uc10=u[29]; il10=u[30];is11=u[31] ;uc11=u[32]; il11=u[33];is12=u[34] ;uc12=u[35]; il12=u[36] ;x=u[37]; v=u[38]
+ #α=LR+Lpr;
+ # RR=4.0e-5
+ rr=manualIntg*1.1925695879998878e-8
+ rpr=Rpr0*(sqrt(t0/(t+1e-4))+(t)^16)/(1.0+(t)^16)
+
+ rrpp=(rpr + rr + 4.53e-7v)
+
+ #= rrpp=rrpp
+ x1=(132.97872599999997 + 35.34758999999999x)
+ x2=(-0.10924199999999998 - 11.782529999999998x)
+ x3=678.7486367999996 + 240.3636119999999x - 8.881784197001252e-16(x^2)
+ rd_11=rd_11;rd22=rd22;rd33=rd33;rd44=rd44 =#
+ Il=il1+il2-0.0+il3+il4+il5+il6-0.0+il7+il8+il9+il10-0.0+il11+il12
+ uf=0.1+0.2*exp(-v/100)
+ F=0.5*0.453e-6*Il*Il*(1-uf*0.124)-uf*FNmec
+ x4=(5.1(5.150399999999999 + 5.436000000000001x))
+ x5=(0.0042 + 0.453x)
+
+ rd_11=(-0.00388 - rd1);rd22=(-0.00388 - rd2);rd33=(0.00388 + rd3);rd44=(-0.00388 - rd4); rd55=(-0.00388 - rd5);rd66=(0.00388 + rd6);rd77=(-0.00388 - rd7);rd88=(-0.00388 - rd8); rd99=(0.00388 + rd9);rd100=(-0.00388 - rd10);rd110=(-0.00388 - rd11);rd120=(-0.00388 - rd12)
+ isid=il1*rd_11 +is1*rd1+il2*rd22+is2*rd2-il3*rd33+is3*rd3+il4*rd44+is4*rd4+il5*rd55 + is5*rd5-il6*rd66+is6*rd6+il7*rd77+is7*rd7+il8*rd88+is8*rd8-il9*rd99+is9*rd9+il10*rd100+is10*rd10+il11*rd110+is11*rd11+il12*rd120+is12*rd12
+
+ du[1] =((-(R1+rs1+rd1)*is1+rd1*il1+uc1)/L1)*operate1
+ du[2]=(-is1/C)*charge1*operate1
+ du[3]=operate1*1e6* (-(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il1*rd_11+is1*rd1)*(0.19607843137254904 ))
+
+
+ du[4] =((-(R1+rs2+rd2)*is2+rd2*il2+uc2)/L1)*operate2
+ du[5]=(-is2/C)*charge2*operate2
+ du[6]=operate2*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il2*rd22+is2*rd2)*(0.19607843137254904 ))
+
+ du[7] =((-(R1+rs3+rd3)*is3+rd3*il3+uc3)/L1)*operate3#*charge3
+ du[8]=((-is3/C)*charge3)*operate3
+ du[9]=operate3*1e6* (-(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il3*rd33+is3*rd3)*(0.19607843137254904 ))
+
+ du[10] =((-(R1+rs4+rd4)*is4+rd4*il4+uc4)/L1)*operate4#*charge3
+ du[11]=((-is4/C)*charge4)*operate4
+ du[12]=operate4*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il4*rd44+is4*rd4)*(0.19607843137254904 ))
+
+ du[13] =((-(R1+rs5+rd5)*is5+rd5*il5+uc5)/L1)*operate5#*charge3
+ du[14]=((-is5/C)*charge5)*operate5
+ du[15]=operate5*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il5*rd55+is5*rd5)*(0.19607843137254904 ))
+
+ du[16] =((-(R1+rs6+rd6)*is6+rd6*il6+uc6)/L1)*operate6#*charge3
+ du[17]=((-is6/C)*charge6)*operate6
+ du[18]=operate6*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il6*rd66+is6*rd6)*(0.19607843137254904 ))
+
+ du[19] =((-(R1+rs7+rd7)*is7+rd7*il7+uc7)/L1)*operate7
+ du[20]=(-is7/C)*charge7*operate7
+ du[21]=operate7*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il7*rd77+is7*rd7)*(0.19607843137254904 ))
+
+ du[22] =((-(R1+rs8+rd8)*is8+rd8*il8+uc8)/L1)*operate8
+ du[23]=(-is8/C)*charge8*operate8
+ du[24]=operate8*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il8*rd88+is8*rd8)*(0.19607843137254904 ))
+
+ du[25] =((-(R1+rs9+rd9)*is9+rd9*il9+uc9)/L1)*operate9
+ du[26]=(-is9/C)*charge9*operate9
+ du[27]=operate9*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(-il9*rd99+is9*rd9)*(0.19607843137254904 ))
+
+ du[28] =((-(R1+rs10+rd10)*is10+rd10*il10+uc10)/L1)*operate10
+ du[29]=(-is10/C)*charge10*operate10
+ du[30]=operate10*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il10*rd100+is10*rd10)*(0.19607843137254904 ))
+
+ du[31] =((-(R1+rs11+rd11)*is11+rd11*il11+uc11)/L1)*operate11
+ du[32]=(-is11/C)*charge11*operate11
+ du[33]=operate11*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il11*rd110+is11*rd11)*(0.19607843137254904 ))
+
+ du[34] =((-(R1+rs12+rd12)*is12+rd12*il12+uc12)/L1)*operate12
+ du[35]=(-is12/C)*charge12*operate12
+ du[36]=operate12*1e6* ( -(-12.0*Il*rrpp*(1.0-0.016339869281045753*x4/x5) + (isid))*( x5 / x4) +(il12*rd120+is12*rd12)*(0.19607843137254904 ))
+
+
+ du[37]=v # v=u[26]
+ du[38]=F/m
+ end
+ function condition(out, u, t, integrator) # Event when condition(out,u,t,integrator) == 0
+ out[1] = t-0.00019
+ out[2] = t-0.00038
+ out[3] = t-0.00057
+ out[4] = t-0.00095
+ out[5] = t-0.00115
+ out[6] = t-0.00142
+ out[7] = t-0.00175
+ out[8] = t-0.00195
+ out[9] = t-0.00210
+
+ out[10] = -u[2]
+ out[11] = -u[5]
+ out[12] = -u[8]
+ out[13] = -u[11]
+ out[14] = -u[14]
+ out[15] = -u[17]
+ out[16] = -u[20]
+ out[17] = -u[23]
+ out[18] = -u[26]
+ out[19] = -u[29]
+ out[20] = -u[32]
+ out[21] = -u[35]
+ out[22] = t-p[50]
+
+ end
+
+ function affect!(integrator, idx)
+ if idx == 1
+ p[15]=1.0
+ elseif idx == 2
+ p[19]=1.0
+ elseif idx == 3
+ p[23]=1.0
+ elseif idx == 4
+ p[27]=1.0
+ elseif idx == 5
+ p[31]=1.0
+ elseif idx == 6
+ p[35]=1.0
+ elseif idx == 7
+ p[39]=1.0
+ elseif idx == 8
+ p[43]=1.0
+ elseif idx == 9
+ p[47]=1.0
+ elseif idx == 10
+ p[4]=0.0
+ p[2]=10.0
+ p[1]=1e-5
+ integrator.u[2]=0.0
+ elseif idx == 11
+ p[8]=0.0
+ p[6]=10.0
+ p[5]=1e-5
+ integrator.u[5]=0.0
+ elseif idx == 12
+ p[12]=0.0
+ p[10]=10.0
+ p[9]=1e-5
+ integrator.u[8]=0.0
+ elseif idx == 13
+ p[16]=0.0
+ p[14]=10.0
+ p[13]=1e-5
+ integrator.u[11]=0.0
+ elseif idx == 14
+ p[20]=0.0
+ p[18]=10.0
+ p[17]=1e-5
+ integrator.u[14]=0.0
+ elseif idx == 15
+ p[24]=0.0
+ p[22]=10.0
+ p[21]=1e-5
+ integrator.u[17]=0.0
+ elseif idx == 16
+ p[28]=0.0
+ p[26]=10.0
+ p[25]=1e-5
+ integrator.u[20]=0.0
+ elseif idx == 17
+ p[32]=0.0
+ p[30]=10.0
+ p[29]=1e-5
+ integrator.u[23]=0.0
+ elseif idx == 18
+ p[36]=0.0
+ p[34]=10.0
+ p[33]=1e-5
+ integrator.u[26]=0.0
+ elseif idx == 19
+ p[40]=0.0
+ p[38]=10.0
+ p[37]=1e-5
+ integrator.u[29]=0.0
+ elseif idx == 20
+ p[44]=0.0
+ p[42]=10.0
+ p[41]=1e-5
+ integrator.u[32]=0.0
+ elseif idx == 21
+ p[48]=0.0
+ p[46]=10.0
+ p[45]=1e-5
+ integrator.u[35]=0.0
+ elseif idx == 22
+ p[49]=p[49]+integrator.u[38]/sqrt(integrator.t)
+ p[50]=p[50]+1e-3
+
+ end
+ end
+
+
+ cbs = VectorContinuousCallback(condition, affect!, 22)
+ u0 = [0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,10750.0,0.0,0.0,0.0]
+ tspan = (0.0, 9.0e-3)
+ p = [1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,1.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,1e5,1e-5,0.0,1.0,0.0,1e-3]
+ prob = ODEProblem(f, u0, tspan, p)
+ sol = solve(prob, Rodas5P(), callback = cbs, reltol=1e-3,abstol=1e-4#= ,dtmax=1e-12 =#,maxiters=Int(1e12)#= dt = 1e-3, adaptive = false =#)
+ p1=plot(sol,idxs=1);
+ savefig(p1, "Rodas5P()_vector_rg12_1")
+ p2=plot(sol,idxs=2);
+ savefig(p2, "Rodas5P()_vector_rg12_2")
+ p3=plot(sol,idxs=3);
+ savefig(p3, "Rodas5P()_vector_rg12_3")
+ end
+ odeDiffEquPackage()
\ No newline at end of file
diff --git a/test/exampleTest.jl b/test/exampleTest.jl
index 0356204..6f65a18 100644
--- a/test/exampleTest.jl
+++ b/test/exampleTest.jl
@@ -11,8 +11,11 @@ sol=solve(odeprob,liqss1(),tspan)
sol=solve(odeprob,liqss2(),tspan)
sol=solve(odeprob,nmliqss1(),tspan)
sol=solve(odeprob,nmliqss2(),tspan)
-@show -2.60.0
+ operate4=1.0
+ end
+
+ if t-0.00038>0.0
+ operate5=1.0
+ end
+ if t-0.00057>0.0
+ operate6=1.0
+ end
+ if t-0.00095>0.0
+ operate7=1.0
+ end
+ if t-0.00115>0.0
+ operate8=1.0
+ end
+ if t-0.00142>0.0
+ operate9=1.0
+ end
+ if t-0.00175>0.0
+ operate10=1.0
+ end
+ if t-0.00195>0.0
+ operate11=1.0
+ end
+ if t-0.00210>0.0
+ operate12=1.0
+ end
+
+ if -(uc1)>0.0
+ charge1=0.0 # rs off not needed since charge=0
+ rs1=ROff;
+ rd1=ROn;
+ uc1=0.0
+ #= uc1=0.0;
+ is1=0.0 =#
+ end
+ if -(uc2)>0.0
+ charge2=0.0
+ rs2=ROff;
+ rd2=ROn;
+ #@show rd2
+ uc2=0.0
+
+ end
+
+ if -(uc3)>0.0
+ charge3=0.0
+ rs3=ROff;
+ rd3=ROn;
+ uc3=0.0
+
+ end
+ if -(uc4)>0.0
+ charge4=0.0
+ rs4=ROff;
+ rd4=ROn;
+ uc4=0.0
+
+ end
+
+ if -(uc5)>0.0
+ charge5=0.0
+ rs5=ROff;
+ rd5=ROn;
+ uc5=0.0
+
+ end
+ if -(uc6)>0.0
+ charge6=0.0
+ rs6=ROff;
+ rd6=ROn;
+ uc6=0.0
+
+ end
+ if -(uc7)>0.0
+ charge7=0.0
+ rs7=ROff;
+ rd7=ROn;
+ uc7=0.0
+
+ end
+ if -(uc8)>0.0
+ charge8=0.0
+ rs8=ROff;
+ rd8=ROn;
+ uc8=0.0
+
+ end
+ if -(uc9)>0.0
+ charge9=0.0
+ rs9=ROff;
+ rd9=ROn;
+ uc9=0.0
+ end
+ if -(uc10)>0.0
+ charge10=0.0
+ rs10=ROff;
+ rd10=ROn;
+ uc10=0.0
+ end
+ if -(uc11)>0.0
+ charge11=0.0
+ rs11=ROff;
+ rd11=ROn;
+ uc11=0.0
+ end
+ if -(uc12)>0.0
+ charge12=0.0
+ rs12=ROff;
+ rd12=ROn;
+ uc12=0.0
+ end
+
+ if t-nextT>0
+ manualIntg=manualIntg+v/sqrt(t[0])
+ nextT=nextT+1e-3
+
+ end
+
+ end
+
+
+ # @show odeprob.eqs
+ # println("start solving")
+ tspan = (0.0, 9.0e-3)
+ sol= solve(odeprob,nmliqss2(),tspan,abstol=1e-3,reltol=1e-2)
+ # sol= solve(odeprob,nmliqss1(),tspan,abstol=1e-3,reltol=1e-2)
+
+
+ # save_Sol(sol,9,note="z1",xlims=(0.0,0.2e-8),ylims=(-0.005,0.005))
+ save_SolSum(sol,3,6,9,12,15,18,21,24,27,30,33,36,interp=0.00001) #add interpl preference
+
+
+
+ save_Sol(sol,37)
+ save_Sol(sol,38)
+
+end
+#@time
+#test()
+#@btime
+test()
+print_timer()
diff --git a/test/runtests.jl b/test/runtests.jl
index 8077516..62fd0ed 100644
--- a/test/runtests.jl
+++ b/test/runtests.jl
@@ -1,5 +1,6 @@
using QuantizedSystemSolver
using Test
+using Plots
using BSON
@testset "QuantizedSystemSolver.jl" begin
# Write your tests here.
@@ -45,6 +46,7 @@ using BSON
sol=solve(odeprob,nmliqss1(),tspan)
@test sol.algName == "nmliqss1"
@test 18.80.0
- discrete[1]=0.0
- end
- if t-3.0>0.0
- u[1] = 1.0
- u[2] = 0.0
- u[3] = 1.0
- u[4] = 0.0
- u[5] = 1.0
- discrete[1]=1.0
-
- end
+
end)
tspan=(0.0,6.0)
-sol=solve(odeprob,qss1(),tspan)
-sol=solve(odeprob,liqss1(),tspan)
+
sol=solve(odeprob,nmliqss2(),tspan)
-@test 1.1