Skip to content
This repository has been archived by the owner on Aug 27, 2024. It is now read-only.

Simulator property tests #70

Merged
merged 7 commits into from
Jun 26, 2024
Merged

Simulator property tests #70

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
4f5decd
tests created
NyMFB42 Jun 26, 2024
21d4ce6
Update crncc_test.fsproj
NyMFB42 Jun 26, 2024
a418b36
Update Compiler.fs
NyMFB42 Jun 26, 2024
8e604d0
Merge branch 'main' into Simulator-property-test
MikaelFangel Jun 26, 2024
d4d96d2
Update Compiler.fs
NyMFB42 Jun 26, 2024
a694174
Merge branch 'Simulator-property-test' of https://github.com/MikaelFa…
NyMFB42 Jun 26, 2024
85b12b0
Update Compiler.fs
NyMFB42 Jun 26, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
308 changes: 154 additions & 154 deletions scripts/PlotSimulator.fsx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10,214 +10,214 @@ open CRN.AST
//Load
let ldreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")]),
ExprS.Expr(
[Species("A")
Species("B")
(
[("A")]),
(
[("A")
("B")
]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")]
(
[("B")]
),
ExprS.Expr([]), 1.0
([]), 1.0
)]
//Addition
let addreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")]),
ExprS.Expr(
[Species("A")
Species("C")
(
[("A")]),
(
[("A")
("C")
]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")]),
ExprS.Expr(
[Species("B")
Species("C")
(
[("B")]),
(
[("B")
("C")
]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("C")]
(
[("C")]
),
ExprS.Expr([]), 1.0
([]), 1.0
)]

//Subtraction
let subreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")]),
ExprS.Expr(
[Species("A")
Species("C")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")]),
ExprS.Expr(
[Species("B")
Species("H")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("C")]),
ExprS.Expr([]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("C")
Species("H")]),
ExprS.Expr(
(
[("A")]),
(
[("A")
("C")]), 1.0)
ReactionS.Reaction(
(
[("B")]),
(
[("B")
("H")]), 1.0)
ReactionS.Reaction(
(
[("C")]),
([]), 1.0)
ReactionS.Reaction(
(
[("C")
("H")]),
(
[]), 1.0)
]

//multiplication
let mulreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")
Species("B")]),
ExprS.Expr(
[Species("A")
Species("B")
Species("C")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("C")]
(
[("A")
("B")]),
(
[("A")
("B")
("C")]), 1.0)
ReactionS.Reaction(
(
[("C")]
),
ExprS.Expr([]), 1.0
([]), 1.0
)]

//division
let divreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")]),
ExprS.Expr(
[Species("A")
Species("C")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")
Species("C")]
(
[("A")]),
(
[("A")
("C")]), 1.0)
ReactionS.Reaction(
(
[("B")
("C")]
),
ExprS.Expr(
[Species("B")]), 1.0
(
[("B")]), 1.0
)]

//square root
let sqrtreaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")]),
ExprS.Expr(
[Species("A")
Species("B")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")
Species("B")]
(
[("A")]),
(
[("A")
("B")]), 1.0)
ReactionS.Reaction(
(
[("B")
("B")]
),
ExprS.Expr(
(
[]), 0.5
)]

//clock simulator
let clock3reaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("A")
Species("B")]),
ExprS.Expr(
[Species("B")
Species("B")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("B")
Species("C")]
(
[("A")
("B")]),
(
[("B")
("B")]), 1.0)
ReactionS.Reaction(
(
[("B")
("C")]
),
ExprS.Expr(
[Species("C")
Species("C")]), 1.0)
(
[("C")
("C")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("C")
Species("A")]
(
[("C")
("A")]
),
ExprS.Expr(
[Species("A")
Species("A")]), 1.0
(
[("A")
("A")]), 1.0
)]
let clock9reaction =
[ReactionS.Reaction(
ExprS.Expr(
[Species("X1")
Species("X2")]),
ExprS.Expr(
[Species("X2")
Species("X2")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X2")
Species("X3")]
(
[("X1")
("X2")]),
(
[("X2")
("X2")]), 1.0)
ReactionS.Reaction(
(
[("X2")
("X3")]
),
ExprS.Expr(
[Species("X3")
Species("X3")]), 1.0)
(
[("X3")
("X3")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X3")
Species("X4")]
(
[("X3")
("X4")]
),
ExprS.Expr(
[Species("X4")
Species("X4")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X4")
Species("X5")]),
ExprS.Expr(
[Species("X5")
Species("X5")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X5")
Species("X6")]
(
[("X4")
("X4")]), 1.0)
ReactionS.Reaction(
(
[("X4")
("X5")]),
(
[("X5")
("X5")]), 1.0)
ReactionS.Reaction(
(
[("X5")
("X6")]
),
ExprS.Expr(
[Species("X6")
Species("X6")]), 1.0)
(
[("X6")
("X6")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X6")
Species("X7")]
(
[("X6")
("X7")]
),
ExprS.Expr(
[Species("X7")
Species("X7")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X7")
Species("X8")]),
ExprS.Expr(
[Species("X8")
Species("X8")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X8")
Species("X9")]
(
[("X7")
("X7")]), 1.0)
ReactionS.Reaction(
(
[("X7")
("X8")]),
(
[("X8")
("X8")]), 1.0)
ReactionS.Reaction(
(
[("X8")
("X9")]
),
ExprS.Expr(
[Species("X9")
Species("X9")]), 1.0)
(
[("X9")
("X9")]), 1.0)
ReactionS.Reaction(
ExprS.Expr(
[Species("X9")
Species("X1")]
(
[("X9")
("X1")]
),
ExprS.Expr(
[Species("X1")
Species("X1")]), 1.0)
(
[("X1")
("X1")]), 1.0)
]
//let initial = Map.ofList ["X1", 1.0-8.0*(0.1e-11); "X2", 0.1e-11; "X3", 0.1e-11; "X4", 0.1e-11; "X5", 0.1e-11; "X6", 0.1e-11; "X7", 0.1e-11; "X8", 0.1e-11; "X9", 0.1e-11]
(*
Expand Down Expand Up @@ -245,7 +245,7 @@ for a in 0..60 do
printfn ""
*)

let initial = Map["A", 5.0; "B", 5.0; "C", 0.0; "H", 0.0]
let (names, states) = (solveODEFast initial 0.01 subreaction)
let initial = Map["A", 1.304408512; "B", 0.02074312655; "C", 0.0]
let (names, states) = (solveODEFast initial 0.01 divreaction)
Seq.map (CRN.Simulator.ArraytoMap names) states
|> CRN.Visualization.plotState (fun s -> seq [ "C" ] |> Seq.contains s) 2000
|> CRN.Visualization.Library.plotState (fun s -> seq [ "C" ] |> Seq.contains s) 20000
7 changes: 3 additions & 4 deletions test/crncc_test/Compiler.fs
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -42,17 +42,16 @@ let ``Compiler: division`` (a:PositiveInt, b:PositiveInt) =
let a = int a
let b = int b
if b = 0 then Prop.classify true "b is zero" true
else if (a / b) >= 13 then Prop.classify true "Timer ran out" true else
else if (a / b) >= 9 then Prop.classify true "Timer ran out" true else
let result = testParser "division.crn"
let result = result |> Result.bind typecheck
let result = Result.bind (fun x -> Ok (compile (Map.ofList [("a0", a);("b0",b)]) x)) result
match result with
| Error a -> Prop.classify true "Error" (Result.isOk result)
| Ok s -> let simulated = convertOut (fst s) 0.025 (snd s)
| Ok s -> let simulated = convertOut (fst s) 0.01 (snd s)
let quot = float (a / b)
let res = float (a % b)
printf "a: %A b: %A q: %A r: %A" a b quot res
let existsState = (Seq.take 150000 simulated) |> Seq.exists (fun map -> abs (Map.find "q" map - quot) <= 0.5
let existsState = (Seq.take 200000 simulated) |> Seq.exists (fun map -> abs (Map.find "q" map - quot) <= 0.5
&& abs (Map.find "r" map - res) <= 0.5 )
Prop.ofTestable existsState

Expand Down
Loading
Loading