Add Sqrt

lib/FloatedBigDigit32.cpp

@@ -3106,6 +3106,36 @@ int FloatedBigDigit32::PowerDiv_boost(FloatedBigDigit32* V) {
 
 }
 
+/****************************************************************************/
+
+int FloatedBigDigit32::Sqrt() {
+
+    if(this->isOver()) {
+        this->overflow();
+        return floatedBigDigitERR;
+    }
+
+    if(this->isZero()) {
+        this->clear();
+        return floatedBigDigitOK;
+    }
+
+    if(this->isMinus()) {
+        this->overflow();
+        return floatedBigDigitERR;
+    }
+
+    FloatedBigDigit32* F = new FloatedBigDigit32();
+
+    F->set(2);
+
+    int ret = this->PowerDiv(F);
+
+    delete F;
+
+    return ret;
+
+} 
 
 /****************************************************************************/
 /****************************************************************************/

lib/FloatedBigDigit32.h

@@ -197,6 +197,8 @@ class FloatedBigDigit32 {
 
         int PowerDiv(FloatedBigDigit32* V);
 
+        int Sqrt();
+
     /****************************************************************************/
 
         int SetSin(FloatedBigDigit32* V);

test/sqrt.cpp

@@ -0,0 +1,37 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+
+#include "FloatedBigDigit32.h"
+
+int main(int argc,const char **argv) {
+
+    char buff[10000];
+
+    if(argc!=3) {
+        printf("USAGE: sqrt KETA value \n");
+        return 0;
+    }
+
+    long KETA = atol(argv[1]);
+
+    SetFloatedBigDigit32Keta((int)KETA);
+
+    FloatedBigDigit32* VAL = new FloatedBigDigit32();
+
+    VAL->Set(argv[2]);
+
+    int ret = VAL->Sqrt();
+
+    VAL->toString(buff,10000,false);
+
+    printf("%s R:%d\n",buff,ret);
+
+    delete VAL;
+
+    FreeFloatedBigDigit32();
+
+    return 0;
+
+}
+