1 files changed, 29 insertions, 0 deletions
diff --git a/src/benchmarks/cfrac/psqrt.c b/src/benchmarks/cfrac/psqrt.c
new file mode 100644
index 0000000..00531a6
--- /dev/null
+++ b/src/benchmarks/cfrac/psqrt.c
@@ -0,0 +1,29 @@
+#include "precision.h"
+
+/*
+ *  Square root
+ */
+precision psqrt(y)
+   precision y;
+{
+   int i;
+   precision x = pUndef, lastx = pUndef;
+
+   i = pcmpz(pparm(y));
+   if (i == 0) {				/* if y == 0 */
+      pset(&lastx, pzero);
+   } else if (i < 0) {				/* if y negative */
+      pset(&x, errorp(PDOMAIN, "psqrt", "negative argument"));
+   } else {
+      pset(&x, y);
+      do {
+	 pset(&lastx, x);
+	 pset(&x, phalf(padd(x, pdiv(y, x))));
+      } while (plt(x, lastx));
+   }
+
+   pdestroy(x);
+
+   pdestroy(y);
+   return presult(lastx);
+}