1 files changed, 17 insertions, 0 deletions
diff --git a/src/zmodsqr.c b/src/zmodsqr.c
new file mode 100644
index 0000000..36e9ed1
--- /dev/null
+++ b/src/zmodsqr.c
@@ -0,0 +1,17 @@
+/* See LICENSE file for copyright and license details. */
+#include "internals"
+
+
+void
+zmodsqr(z_t a, z_t b, z_t c)
+{
+	/* TODO What is the fastest way to do zmodsqr? */
+	if (a == c) {
+		zset(libzahl_tmp_modsqr, c);
+		zsqr(a, b);
+		zmod(a, a, libzahl_tmp_modsqr);
+	} else {
+		zsqr(a, b);
+		zmod(a, a, c);
+	}
+}