@@ -491,6 +491,13 @@ NumberTMNE List := D -> (
491491-- METHODS FOR DEPENDENCY EQUILIBRIA --
492492-- ***************************************--
493493
494+ -- -----------------------------------------------
495+ -- probabilityRing (List)
496+ --
497+ -- Given a list Di this constructs a ring of joint
498+ -- probabilities for a game of format Di
499+ -- -----------------------------------------------
500+
494501
495502probabilityRing = method (Options => { CoefficientRing => QQ , ProbabilityVariableName => " p" })
496503probabilityRing List := Ring => opts -> Di -> (
@@ -507,12 +514,27 @@ probabilityRing List := Ring => opts -> Di -> (
507514 R)
508515
509516
517+ -- ----------------------------------------------
518+ -- randomGame (List)
519+ --
520+ -- Given a list Di this constructs a random game
521+ -- of format Di
522+ -- ----------------------------------------------
523+
510524randomGame = method (Options => {CoefficientRing => QQ })
511525randomGame List := List => opts -> Di -> (
512526 K := opts.CoefficientRing ;
513527 apply (length Di, i -> randomTensor(K, Di)))
514528
515529
530+ -- -------------------------------------------------------
531+ -- spohnMatrices (Ring, List)
532+ --
533+ -- Given the underlying probability ring R and a game X
534+ -- this constructs the Spohn matrices for X
535+ -- -------------------------------------------------------
536+
537+
516538spohnMatrices = method ()
517539spohnMatrices (Ring , List ) := List => (PR, X) -> (
518540 p := PR#" probabilityVariable" ;
@@ -522,12 +544,25 @@ spohnMatrices (Ring, List) := List => (PR, X) -> (
522544 apply (n, i -> matrix apply (d_i, k -> {sum (select (J, j -> j_i==k), j -> p#j),
523545 sum (select (J, j -> j_i==k), j -> (X_i)#j * p#j) })))
524546
547+ -- ------------------------------------------------------
548+ -- spohnIdeal (Ring, List)
549+ --
550+ -- Given the underlying probability ring R and a game X
551+ -- this constructs the Spohn ideal for X
552+ -- ------------------------------------------------------
553+
525554
526555spohnIdeal = method ()
527556spohnIdeal (Ring , List ) := List => (PR, X) -> (
528557 M := spohnMatrices(PR, X);
529558 sum (M, m -> minors (2, m)))
530559
560+ -- ------------------------------------------------------
561+ -- konstanzMatrix (Ring, List)
562+ --
563+ -- Given the underlying probability ring R and a game X
564+ -- this constructs the Konstanz Matrix for X
565+ -- ------------------------------------------------------
531566
532567konstanzMatrix = method (Options =>{ KonstanzVariableName => " k" })
533568konstanzMatrix (Ring , List ) := Matrix => opts -> (PR, X) -> (
@@ -542,8 +577,7 @@ konstanzMatrix (Ring, List) := Matrix => opts -> (PR, X) -> (
542577 P := vector (apply (J, j -> p#j));
543578 fold ((M0, M1) -> M0 || M1,
544579 apply (n, i -> transpose matrix apply (Di_i,
545- j -> diff (P, (LinearForms_i)_(j, 0 )))))
546- )
580+ j -> diff (P, (LinearForms_i)_(j, 0 ))))))
547581
548582
549583
0 commit comments