Changeset 02af91 in git for factory/ffops.h

Timestamp:

Jun 18, 1996, 8:57:00 AM (28 years ago)

Author:

Rüdiger Stobbe <stobbe@…>

Branches:

(u'spielwiese', '17f1d200f27c5bd38f5dfc6e8a0879242279d1d8')

Children:

cfa45ec3cccf336939c446d85de2dca6283593b1

Parents:

35930b6cf4435468b76d893cb0f2be0e692720cf

Message:

"Now the functionality to handle prime numbers that are bigger than 2^15.
This results in various changes and some new functions.
"


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

File:

• factory/ffops.h (modified) (4 diffs)

factory/ffops.h

@@ -1 +1 @@
 // emacs edit mode for this file is -*- C++ -*-
-// $Id: ffops.h,v 1.0 1996-05-17 10:59:40 stobbe Exp $
+// $Id: ffops.h,v 1.1 1996-06-18 06:57:00 stobbe Exp $
 
 #ifndef INCL_FFOPS_H
@@ -7 +7 @@
 /*
 $Log: not supported by cvs2svn $
+Revision 1.0  1996/05/17 10:59:40  stobbe
+Initial revision
+
 */
 
@@ -15 +18 @@
 extern int ff_halfprime;
 extern short * ff_invtab;
+extern bool ff_big;
 
 int ff_newinv ( const int );
+int ff_biginv ( const int );
 void ff_setprime ( const int );
 
 inline int ff_norm ( const int a )
 {
-  int n = a % ff_prime;
-  if ( n < 0 )
-    return n + ff_prime;
-  else
-    return n;
+    int n = a % ff_prime;
+    if ( n < 0 )
+        return n + ff_prime;
+    else
+        return n;
 }
 
@@ -38 +43 @@
 inline int ff_longnorm ( const long a )
 {
-  int n = (int)(a % (long)ff_prime);
-  if ( n < 0 )
-    return n + ff_prime;
-  else
-    return n;
+    int n = (int)(a % (long)ff_prime);
+    if ( n < 0 )
+        return n + ff_prime;
+    else
+        return n;
+}
+
+inline int ff_bignorm ( const long long int a )
+{
+    int n = (int)(a % (long long int)ff_prime);
+    if ( n < 0 )
+        return n + ff_prime;
+    else
+        return n;
 }
 
 inline int ff_add ( const int a, const int b )
 {
-  return ff_norm( a + b );
+    return ff_norm( a + b );
 }
 
 inline int ff_sub ( const int a, const int b )
 {
-  return ff_norm( a - b );
+    return ff_norm( a - b );
 }
 
 inline int ff_neg ( const int a )
 {
-  return ff_norm( -a );
+    return ff_norm( -a );
 }
 
 inline int ff_mul ( const int a, const int b )
 {
-  return ff_longnorm ( (long)a * (long)b );
+    if ( ff_big )
+        return ff_bignorm( (long long int)a * (long long int)b );
+    else
+        return ff_longnorm ( (long)a * (long)b );
 }
 
 inline int ff_inv ( const int a )
 {
-  register b;
-  if ( (b = (int)(ff_invtab[a])) )
-    return b;
-  else
-    return ff_newinv( a );
+    if ( ff_big )
+        return ff_biginv( a );
+    else {
+        register b;
+        if ( (b = (int)(ff_invtab[a])) )
+            return b;
+        else
+            return ff_newinv( a );
+    }
+    
 }
 
 inline int ff_div ( const int a, const int b )
 {
-  return ff_mul( a, ff_inv( b ) );
+    return ff_mul( a, ff_inv( b ) );
 }
 
 #endif /* INCL_FFOPS_H */
-
-
-
-
-
-