Changeset 5600a4 in git for Singular/LIB/modular.lib

Timestamp:

Sep 7, 2016, 6:16:24 PM (8 years ago)

Author:

Hans Schoenemann <hannes@…>

Branches:

(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')

Children:

e4ac9d386d815e4411fb4fcea3826b2c8d70fd37

Parents:

78321f7f1ef646d795a532545c8159f388de4136a88bbe6011f0f3b0eec150334d2c6243aa1f57ad

git-author:

Hans Schoenemann <hannes@mathematik.uni-kl.de>2016-09-07 18:16:24+02:00

git-committer:

GitHub <noreply@github.com>2016-09-07 18:16:24+02:00

Message:

Merge pull request #784 from steenpass/modular

chg: update modular.lib

File:

• Singular/LIB/modular.lib (modified) (6 diffs)

Singular/LIB/modular.lib

@@ -40 +40 @@
 ";
 
+LIB "general.lib";
 LIB "resources.lib";
 LIB "tasks.lib";
@@ -157 +158 @@
         // delete unlucky primes
         indices = deleteUnluckyPrimes(modresults);
+        indices = sort(indices)[1];
+        for (i = size(indices); i > 1; i--) {
+            if (indices[i] == indices[i-1]) {
+                indices = delete(indices, i);
+            }
+        }
         for (i = size(indices); i > 0; i--) {
             modresults = delete(modresults, indices[i]);
@@ -164 +171 @@
         // lift result
         if (N == 1) {
-            result_lift = chinrem(modresults, primes);
+            result_lift = chinrem_recursive(modresults, primes);
         }
         else {
-            result_lift = chinrem(list(result_lift)+modresults,
+            result_lift = chinrem_recursive(list(result_lift)+modresults,
                 list(N)+primes);
         }
@@ -211 +218 @@
     ideal I = x9y2+x10, x2y7-y8;
     modular("std", list(I));
+}
+
+static proc chinrem_recursive(list modresults, list moduli)
+{
+    if (typeof(modresults[1]) != "list") {
+        return(chinrem(modresults, moduli));
+    }
+    // if typeof(modresults[1]) == "list", then we assume
+    // size(modresults[1]) == size(modresults[k])
+    // for k = 2, ..., size(modresults)
+    int n_sublists = size(modresults[1]);
+    int i, j;
+    list modresults_sublists;
+    for (j = n_sublists; j > 0; j--) {
+        modresults_sublists[j] = list();
+        for (i = size(modresults); i > 0; i--) {
+            modresults_sublists[j][i] = modresults[i][j];
+        }
+        task t(j) = "Modular::chinrem_recursive",
+            list(modresults_sublists[j], moduli);
+    }
+    startTasks(t(1..n_sublists));
+    waitAllTasks(t(1..n_sublists));
+    list result;
+    for (j = n_sublists; j > 0; j--) {
+        result[j] = getResult(t(j));
+        killTask(t(j));
+    }
+    return(result);
 }
 
@@ -290 +326 @@
 {
     arg_type = typeof(farey_arg);
-    if (arg_type != "bigint" && arg_type != "ideal"
-        && arg_type != "module" && arg_type != "matrix") {
+    if (arg_type != "bigint" && arg_type != "ideal" && arg_type != "module"
+        && arg_type != "matrix" && arg_type != "list") {
         ERROR("wrong input type");
+    }
+    int i;
+    if (arg_type == "list") {
+        int size_farey_arg = size(farey_arg);
+        for (i = size_farey_arg; i > 0; i--) {
+            task t(i) = "Modular::farey_parallel", list(farey_arg[i], farey_N);
+        }
+        startTasks(t(1..size_farey_arg));
+        waitAllTasks(t(1..size_farey_arg));
+        list result;
+        for (i = size_farey_arg; i > 0; i--) {
+            result[i] = getResult(t(i));
+            killTask(t(i));
+        }
+        return(result);
     }
     if (arg_type == "bigint" || arg_type == "matrix") {
@@ -302 +353 @@
     int chunks = par_range(size_arg);
     intvec range;
-    int i;
     for (i = chunks; i > 0; i--) {
         range = par_range(size_arg, i);