Changeset fd0df1 in git


Ignore:
Timestamp:
Jun 10, 2011, 12:44:10 PM (13 years ago)
Author:
Stefan Steidel <steidel@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'c7af8613769b29c741d6c338945669719f1fc4f8')
Children:
3ffd6d2e233ebf7267b3b68dc7d6d7d506305516
Parents:
1035602ff5b8a29d54e29d4d0b863eab2f4ecab4
Message:
Some optimization regarding the number of processes used im modStd, some minor documentation improvements.

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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/modstd.lib

    r103560 rfd0df1  
    368368      for(j = 1; j <= size(I[i]); j++)
    369369      {
    370          if((leadcoef(I[i][j]) mod p) == 0) { return(0); }
     370         if((numerator(leadcoef(I[i][j])) mod p) == 0) { return(0); }
     371         if((denominator(leadcoef(I[i][j])) mod p) == 0) { return(0); }
    371372      }
    372373   }
     
    650651@*       - variant = 5, if I is not homogeneous, 1-block-ordering and
    651652                        exactness = 1,
    652 @*       - variant = 3, if I is not homogeneous, complicated ordering (lp or
     653@*       - variant = 6, if I is not homogeneous, complicated ordering (lp or
    653654                        > 1 block) and exactness = 1.
    654655EXAMPLE: example modStd; shows an example
     
    668669   int index = 1;
    669670   int i,k,c;
    670    int pd = printlevel-voice+2;
    671671   int j = 1;
    672672   int pTest, sizeTest;
     
    682682         int n1 = #[1];
    683683         int exactness = 1;
    684          int n2 = 10;
    685          int n3 = 10;
     684         if(n1 >= 10)
     685         {
     686            int n2 = n1 + 1;
     687            int n3 = n1;
     688         }
     689         else
     690         {
     691            int n2 = 10;
     692            int n3 = 10;
     693         }
    686694      }
    687695      if(size(#) == 2)
     
    689697         int n1 = #[1];
    690698         int exactness = #[2];
    691          int n2 = 10;
    692          int n3 = 10;
     699         if(n1 >= 10)
     700         {
     701            int n2 = n1 + 1;
     702            int n3 = n1;
     703         }
     704         else
     705         {
     706            int n2 = 10;
     707            int n3 = 10;
     708         }
    693709      }
    694710      if(size(#) == 4)
     
    696712         int n1 = #[1];
    697713         int exactness = #[2];
    698          int n2 = #[3];
    699          int n3 = #[4];
     714         if(n1 >= #[3])
     715         {
     716            int n2 = n1 + 1;
     717         }
     718         else
     719         {
     720            int n2 = #[3];
     721         }
     722         if(n1 >= #[4])
     723         {
     724            int n3 = n1;
     725         }
     726         else
     727         {
     728            int n3 = #[4];
     729         }
    700730      }
    701731   }
     
    706736      int n2 = 10;
    707737      int n3 = 10;
     738   }
     739
     740   if(printlevel >= 10)
     741   {
     742      "n1 = "+string(n1)+", n2 = "+string(n2)+", n3 = "+string(n3);
    708743   }
    709744
     
    834869      if(h)
    835870      {
    836          variant = 1; if(printlevel >= 10) { "variant = 1"; }
     871         variant = 1;
     872         if(printlevel >= 10) { "variant = 1"; }
    837873         rl[1] = L[5];
    838874         def @r = ring(rl);
     
    961997   {
    962998      tt = timer; rt = rtimer;
    963 
     999     
     1000      if(printlevel >= 10) { "size(L) = "+string(size(L)); }
     1001         
    9641002      if(n1 > 1)
    9651003      {
    966          if(printlevel >= 10) { "size(L) = "+string(size(L)); }
    9671004         while(j <= size(L) + 1)
    9681005         {
     
    11251162               }
    11261163            }
    1127             if(pd > 2) { "pTest o.k. but result wrong"; }
    1128          }
    1129          if(pd > 2) { "pTest o.k. but result wrong"; }
     1164         }
    11301165      }
    11311166
     
    11831218   size(reduce(J1,J2));
    11841219   size(reduce(J2,J1));
    1185 
    1186    /*
    1187    ring R4 = 0,x(1..4),wp(1,-1,-1,1);
    1188    ideal I = cyclic(4);
    1189    ideal J1 = modStd(I,1,0);
    1190    */
    11911220}
    11921221
Note: See TracChangeset for help on using the changeset viewer.