Changeset 5fe545 in git

Timestamp:

May 18, 2005, 5:39:19 PM (19 years ago)

Author:

Hans Schönemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

b340969c053d9f828187f38bac97b1feb96e654e

Parents:

7835d2ba9bdd252cdfa90863a2ff862d818f9ac1

Message:

*hannes: ngcPower fixed


git-svn-id: file:///usr/local/Singular/svn/trunk@8218 2c84dea3-7e68-4137-9b89-c4e89433aadc

File:

• kernel/gnumpc.cc (modified) (3 diffs)

kernel/gnumpc.cc

@@ -2 +2 @@
 *  Computer Algebra System SINGULAR     *
 ****************************************/
-/* $Id: gnumpc.cc,v 1.2 2005-01-18 08:53:45 Singular Exp $ */
+/* $Id: gnumpc.cc,v 1.3 2005-05-18 15:39:19 Singular Exp $ */
 /*
 * ABSTRACT: computations with GMP complex floating-point numbers
@@ -290 +290 @@
     return;
   }
-  if ( exp == 1 )
+  else if ( exp == 1 )
   {
     nNew(u);
@@ -306 +306 @@
     return;
   }
-  ngcPower(x,exp-1,u);
-  gmp_complex *n=new gmp_complex();
-  *n=*(gmp_complex*)x;
-  *(gmp_complex*)(*u) *= *(gmp_complex*)n;
-  delete n;
+  else if (exp == 2)
+  {
+    nNew(u);
+    if ( x == NULL )
+    {
+      gmp_complex* n = new gmp_complex();
+      *u=(number)n;
+    }
+    else
+    {
+      gmp_complex* n = new gmp_complex();
+      *n= *(gmp_complex*)x;
+      *u=(number)n;
+      *(gmp_complex*)(*u) *= *(gmp_complex*)n;
+    }
+    return;
+  }
+  if (exp&1==1)
+  {
+    ngcPower(x,exp-1,u);
+    gmp_complex *n=new gmp_complex();
+    *n=*(gmp_complex*)x;
+    *(gmp_complex*)(*u) *= *(gmp_complex*)n;
+    delete n;
+  }
+  else
+  {
+    number w;
+    nNew(&w);
+    ngcPower(x,exp/2,&w);
+    ngcPower(w,2,u);
+    nDelete(&w);
+  }
 }