1 files changed, 58 insertions, 1 deletions

diff --git a/src/lib/libcrypto/bn/bn_sqr.c b/src/lib/libcrypto/bn/bn_sqr.c
index ff254764e3..74d5eded94 100644
--- a/src/lib/libcrypto/bn/bn_sqr.c
+++ b/src/lib/libcrypto/bn/bn_sqr.c
@@ -1,4 +1,4 @@
-/* $OpenBSD: bn_sqr.c,v 1.21 2023/01/21 14:10:46 jsing Exp $ */
+/* $OpenBSD: bn_sqr.c,v 1.22 2023/01/23 12:09:06 jsing Exp $ */
 /* Copyright (C) 1995-1998 Eric Young (eay@cryptsoft.com)
  * All rights reserved.
  *
@@ -56,6 +56,7 @@
  * [including the GNU Public Licence.]
  */
 
+#include <assert.h>
 #include <stdio.h>
 #include <string.h>
 
@@ -178,6 +179,62 @@ bn_sqr_comba8(BN_ULONG *r, const BN_ULONG *a)
 }
 #endif
 
+#ifndef HAVE_BN_SQR_WORDS
+#if defined(BN_LLONG) || defined(BN_UMULT_HIGH)
+void
+bn_sqr_words(BN_ULONG *r, const BN_ULONG *a, int n)
+{
+        assert(n >= 0);
+        if (n <= 0)
+                return;
+
+#ifndef OPENSSL_SMALL_FOOTPRINT
+        while (n & ~3) {
+                sqr(r[0], r[1], a[0]);
+                sqr(r[2], r[3], a[1]);
+                sqr(r[4], r[5], a[2]);
+                sqr(r[6], r[7], a[3]);
+                a += 4;
+                r += 8;
+                n -= 4;
+        }
+#endif
+        while (n) {
+                sqr(r[0], r[1], a[0]);
+                a++;
+                r += 2;
+                n--;
+        }
+}
+#else /* !(defined(BN_LLONG) || defined(BN_UMULT_HIGH)) */
+void
+bn_sqr_words(BN_ULONG *r, const BN_ULONG *a, int n)
+{
+        assert(n >= 0);
+        if (n <= 0)
+                return;
+
+#ifndef OPENSSL_SMALL_FOOTPRINT
+        while (n & ~3) {
+                sqr64(r[0], r[1], a[0]);
+                sqr64(r[2], r[3], a[1]);
+                sqr64(r[4], r[5], a[2]);
+                sqr64(r[6], r[7], a[3]);
+                a += 4;
+                r += 8;
+                n -= 4;
+        }
+#endif
+        while (n) {
+                sqr64(r[0], r[1], a[0]);
+                a++;
+                r += 2;
+                n--;
+        }
+}
+#endif
+#endif
+
 /* tmp must have 2*n words */
 void
 bn_sqr_normal(BN_ULONG *r, const BN_ULONG *a, int n, BN_ULONG *tmp)