Changeset 747a1e in git for Singular/LIB
- Timestamp:
- May 28, 2010, 12:16:16 PM (14 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '4bd32dfef92ec9f5ed8dceee82d14318ae147107')
- Children:
- 12cca3ff11b4a22801170a24e54b4b0fc1eda12d
- Parents:
- 8017e90b3d8ba6d6c93419d574d5ec289a1d885a
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/random.lib
r8017e9 r747a1e 126 126 //----------------------------- set defaults ---------------------------------- 127 127 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; }}}} 131 131 //------------------ use sparsemat for creation of sparseid ------------------- 132 132 int ii; matrix i[1][k]; intmat m; … … 137 137 u=1; 138 138 } 139 for ( ii=u; ii<=o; ii =ii+1)139 for ( ii=u; ii<=o; ii++) 140 140 { 141 141 m = sparsemat(size(maxideal(ii)),k,p,b); … … 169 169 int ii; ideal i; intmat m; ideal id; 170 170 171 for ( ii=k; ii>0; ii =ii-1)171 for ( ii=k; ii>0; ii--) 172 172 { 173 173 id = maxideal(random(u, o)); // monomial basis of some degree … … 207 207 if( p<40 ) 208 208 { 209 for( ii=1; ii<=t; ii =ii+1)209 for( ii=1; ii<=t; ii++ ) 210 210 { r=( random(1,100)>p ); v[1,ii]=r*random(1,b); h=h+r; } 211 211 } … … 248 248 int pe=50;int u=0;int pp=75;int b=100; 249 249 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]; }}}} 253 253 //------------------- use sparsemat and sparseid ----------------------------- 254 254 intmat I = sparsemat(n,m,pe,1); 255 for(ii= 1; ii<=n;ii++)255 for(ii=n; ii>0;ii--) 256 256 { 257 257 id = sparseid(m,u,o,pp,b); 258 for(jj= 1; jj<=m; jj++)258 for(jj=m; jj>0; jj--) 259 259 { 260 260 if( I[ii,jj] !=0) … … 286 286 //----------------------------- set defaults ---------------------------------- 287 287 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; }}}} 291 291 int ii; poly f; 292 292 //----------------- 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]; } 294 294 return(f); 295 295 } … … 320 320 if( n<=m ) { min=n-1; M[n,n]=1; } 321 321 else { min=m; } 322 for( ii=1; ii<=min; ii =ii+1)322 for( ii=1; ii<=min; ii++ ) 323 323 { 324 324 l=r+1; r=r+n-ii;
Note: See TracChangeset
for help on using the changeset viewer.