Changeset 747a1e in git for Singular/LIB


Ignore:
Timestamp:
May 28, 2010, 12:16:16 PM (14 years ago)
Author:
Hans Schoenemann <hannes@…>
Branches:
(u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '4bd32dfef92ec9f5ed8dceee82d14318ae147107')
Children:
12cca3ff11b4a22801170a24e54b4b0fc1eda12d
Parents:
8017e90b3d8ba6d6c93419d574d5ec289a1d885a
Message:
optim.

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

Legend:

Unmodified
Added
Removed
  • Singular/LIB/random.lib

    r8017e9 r747a1e  
    126126//----------------------------- set defaults ----------------------------------
    127127   if( size(#)>=3 ) { int o=#[1]; int p=#[2]; int b=#[3]; }
    128    if( size(#)==2 ) { int o=#[1]; int p=#[2]; int b=30000; }
    129    if( size(#)==1 ) { int o=#[1]; int p=75; int b=30000; }
    130    if( size(#)==0 ) { int o=u; int p=75; int b=30000; }
     128   else {if( size(#)==2 ) { int o=#[1]; int p=#[2]; int b=30000; }
     129   else {if( size(#)==1 ) { int o=#[1]; int p=75; int b=30000; }
     130   else {if( size(#)==0 ) { int o=u; int p=75; int b=30000; }}}}
    131131//------------------ use sparsemat for creation of sparseid -------------------
    132132   int ii; matrix i[1][k]; intmat m;
     
    137137      u=1;
    138138   }
    139    for ( ii=u; ii<=o; ii=ii+1)
     139   for ( ii=u; ii<=o; ii++)
    140140   {
    141141       m = sparsemat(size(maxideal(ii)),k,p,b);
     
    169169   int ii; ideal i; intmat m; ideal id;
    170170
    171    for ( ii=k; ii>0; ii=ii-1)
     171   for ( ii=k; ii>0; ii--)
    172172   {
    173173       id = maxideal(random(u, o)); // monomial basis of some degree
     
    207207   if( p<40 )
    208208   {
    209       for( ii=1; ii<=t; ii=ii+1 )
     209      for( ii=1; ii<=t; ii++ )
    210210      { r=( random(1,100)>p ); v[1,ii]=r*random(1,b); h=h+r; }
    211211   }
     
    248248   int pe=50;int u=0;int pp=75;int b=100;
    249249   if( size(#)==4 ) { u=#[1]; pe=#[2]; pp=#[3]; b=#[4]; }
    250    if( size(#)==3 ) { u=#[1]; pe=#[2]; pp=#[3]; }
    251    if( size(#)==2 ) { u=#[1]; pe=#[2]; }
    252    if( size(#)==1 ) { u=#[1]; }
     250   else { if( size(#)==3 ) { u=#[1]; pe=#[2]; pp=#[3]; }
     251   else { if( size(#)==2 ) { u=#[1]; pe=#[2]; }
     252   else {if( size(#)==1 ) { u=#[1]; }}}}
    253253//------------------- use sparsemat and sparseid  -----------------------------
    254254   intmat I = sparsemat(n,m,pe,1);
    255    for(ii=1; ii<=n;ii++)
     255   for(ii=n; ii>0;ii--)
    256256   {
    257257     id = sparseid(m,u,o,pp,b);
    258      for(jj=1; jj<=m; jj++)
     258     for(jj=m; jj>0; jj--)
    259259     {
    260260        if( I[ii,jj] !=0)
     
    286286//----------------------------- set defaults ----------------------------------
    287287   if( size(#)>=3 ) { int o=#[1]; int p=#[2]; int b=#[3]; }
    288    if( size(#)==2 ) { int o=#[1]; int p=#[2]; int b=30000; }
    289    if( size(#)==1 ) { int o=#[1]; int p=75; int b=30000; }
    290    if( size(#)==0 ) { int o=u; int p=75; int b=30000; }
     288   else {if( size(#)==2 ) { int o=#[1]; int p=#[2]; int b=30000; }
     289   else {if( size(#)==1 ) { int o=#[1]; int p=75; int b=30000; }
     290   else {if( size(#)==0 ) { int o=u; int p=75; int b=30000; }}}}
    291291   int ii; poly f;
    292292//----------------- use sparseid for creation of sparsepoly -------------------
    293    for( ii=u; ii<=o; ii=ii+1 ) { f=f+sparseid(1,ii,ii,p,b)[1]; }
     293   for( ii=u; ii<=o; ii++ ) { f=f+sparseid(1,ii,ii,p,b)[1]; }
    294294   return(f);
    295295}
     
    320320   if( n<=m ) { min=n-1; M[n,n]=1; }
    321321   else { min=m; }
    322    for( ii=1; ii<=min; ii=ii+1 )
     322   for( ii=1; ii<=min; ii++ )
    323323   {
    324324      l=r+1; r=r+n-ii;
Note: See TracChangeset for help on using the changeset viewer.