Implemented efficiecy computation

Author: bwbush

analysis/markov/Linleios/Evolve.lean

@@ -67,3 +67,12 @@ def ebDistribution : Probabilities → HashMap Nat Probability :=
           $ singleton ⟨ state.ebCount , p ⟩
     )
     ∅
+
+def ebEfficiency (states : Probabilities) : Float :=
+  let rbCount := states.keys.head!.rbCount
+  let ebCount :=
+    HashMap.fold
+      (fun acc state p =>acc + state.ebCount.toFloat * p)
+      0
+      states
+  ebCount / (rbCount.toFloat - 1)

analysis/markov/Main.lean

@@ -4,7 +4,7 @@ import Linleios
 
 def env : Environment := makeEnvironment 0.05 0.95 0.90 600 0.75 1 4 7
 
-def sn := simulate env default 1e-6 10
+def sn := simulate env default 1e-6 25
 
 def main : IO Unit :=
   do
@@ -14,4 +14,6 @@ def main : IO Unit :=
     IO.println ""
     print $ ebDistribution sn
     IO.println ""
+    print $ ebEfficiency sn
+    IO.println ""
     print $ 1 - totalProbability sn