From 31751ccb4d382fd69c9d965c04c4ae318393f54e Mon Sep 17 00:00:00 2001
From: skylarmt <apocalypselabs@gmail.com>
Date: Mon, 13 Apr 2015 12:30:02 -0600
Subject: [PATCH] powermod(a, -1, m) now gives modular inverse.

---
 src/net/apocalypselabs/symat/Functions.java | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/src/net/apocalypselabs/symat/Functions.java b/src/net/apocalypselabs/symat/Functions.java
index 52c7d3b..aa7acb5 100644
--- a/src/net/apocalypselabs/symat/Functions.java
+++ b/src/net/apocalypselabs/symat/Functions.java
@@ -175,13 +175,17 @@ public class Functions {
     }
 
     /**
-     * Take (a^b) mod m.
+     * Take (a^b) mod m.  If b is -1, finds the modular inverse of a mod m.
      * @param a Number
      * @param b Exponent
      * @param m Modulus
      * @return pow(a,b)%m.
      */
     public double powermod(double a, double b, double m) {
+        if (b == -1) {
+            BigInteger in = new BigInteger(String.valueOf((int)a));
+            return in.modInverse(new BigInteger(String.valueOf((int)m))).doubleValue();
+        }
         return (pow(a, b) % m + m) % m;
     }