From e0305643a78a9df7e7e76f28fcb5c9fa8ac3fb2d Mon Sep 17 00:00:00 2001 From: Dave Parker Date: Wed, 16 Nov 2011 21:45:54 +0000 Subject: [PATCH] Add (commented out) code for backwards Gauss-Seidel on DTMCs. git-svn-id: https://www.prismmodelchecker.org/svn/prism/prism/trunk@4168 bbc10eb1-c90d-0410-af57-cb519fbb1720 --- prism/src/explicit/DTMCExplicit.java | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/prism/src/explicit/DTMCExplicit.java b/prism/src/explicit/DTMCExplicit.java index b1838c0e..4d23b012 100644 --- a/prism/src/explicit/DTMCExplicit.java +++ b/prism/src/explicit/DTMCExplicit.java @@ -206,6 +206,15 @@ public abstract class DTMCExplicit extends ModelExplicit implements DTMC maxDiff = diff > maxDiff ? diff : maxDiff; vect[s] = d; } + // Use this code instead for backwards Gauss-Seidel + /*for (s = numStates - 1; s >= 0; s--) { + if (subset.get(s)) { + d = mvMultJacSingle(s, vect); + diff = absolute ? (Math.abs(d - vect[s])) : (Math.abs(d - vect[s]) / d); + maxDiff = diff > maxDiff ? diff : maxDiff; + vect[s] = d; + } + }*/ } return maxDiff; }