Changeset c623f27 in git
- Timestamp:
- Oct 6, 2010, 8:10:20 PM (14 years ago)
- Branches:
- (u'fieker-DuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', '38077648e7239f98078663eb941c3c979511150a')
- Children:
- ce136a3deb909101f53d2844ba9fdd9bdc1c7842
- Parents:
- 3561bd1ffefd8ae194f9d9ec57e0d2d9d54be51c
- Location:
- Singular/LIB
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/dmodapp.lib
r3561bd rc623f27 606 606 if (sl>2) 607 607 { 608 609 608 if (typeof(l[3])<>"string"){wronginput = 1;} 609 if (sl>3){wronginput = 1;} 610 610 } 611 611 } … … 1512 1512 if (size(#)>2) 1513 1513 { 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1514 if (typeof(#[3])=="intvec") 1515 { 1516 if (size(#[3])==2*n && allPositive(#[3])==1) 1517 { 1518 homogweights = #[3]; 1519 } 1520 else 1521 { 1522 print("// Homogenization weight vector must consist of positive entries and be"); 1523 print("// of size " + string(n) + ". Using weight (1,...,1)."); 1524 } 1525 } 1526 1526 } 1527 1527 } … … 1569 1569 } 1570 1570 def Dh = nc_algebra(1,@relD); 1571 setring Dh; kill @Dh; 1571 setring Dh; kill @Dh; 1572 1572 dbprint(ppl-1,"// computing in ring",Dh); 1573 1573 // 2. Compute the initial ideal … … 1886 1886 if (typeof(#[2])=="int" || typeof(#[2])=="number") 1887 1887 { 1888 1889 1888 l0 = int(#[2]); 1889 l0set = 1; 1890 1890 } 1891 1891 if (size(#)>2) 1892 1892 { 1893 1893 if (typeof(#[3])=="ideal") 1894 1894 { 1895 1896 1897 1895 G = #[3]; 1896 Gset = 1; 1897 } 1898 1898 } 1899 1899 } … … 2014 2014 for (k=1; k<=ncols(B); k++) 2015 2015 { 2016 2017 2018 2016 p = B[k]*G[i]; 2017 p = f(p); 2018 M[size(M)+1] = p; 2019 2019 } 2020 2020 } … … 2192 2192 if (M[1,i]<>0) 2193 2193 { 2194 2195 2196 2197 2198 2194 v = M[2..r,i]; 2195 if (v == zm) 2196 { 2197 J[size(J+1)] = M[1,i]; 2198 } 2199 2199 } 2200 2200 } … … 2302 2302 if (v[i] == v[i+1]) 2303 2303 { 2304 2304 ERROR("No double entries allowed in intvec"); 2305 2305 } 2306 2306 } … … 2371 2371 if (v[i] == v[i+1]) 2372 2372 { 2373 2373 ERROR("No double entries allowed in intvec"); 2374 2374 } 2375 2375 } … … 2453 2453 if (typeof(#[2])=="int" || typeof(#[2])=="number") 2454 2454 { 2455 2456 2455 l0 = int(#[2]); 2456 l0set = 1; 2457 2457 } 2458 2458 if (size(#)>2) 2459 2459 { 2460 2460 if (typeof(#[3])=="ideal") 2461 2461 { 2462 2463 2464 2462 G = #[3]; 2463 Gset = 1; 2464 } 2465 2465 } 2466 2466 } … … 2476 2476 if (v == intvec(0)) 2477 2477 { 2478 2478 v[1] = i; 2479 2479 } 2480 2480 else 2481 2481 { 2482 2482 v[size(v)+1] = i; 2483 2483 } 2484 2484 } … … 2616 2616 if (allPositive(#[1])==1) 2617 2617 { 2618 2618 w = #[1]; 2619 2619 } 2620 2620 else 2621 2621 { 2622 2623 2622 print("// Entries of intvec must be strictly positive"); 2623 print("// Using weight " + string(w)); 2624 2624 } 2625 2625 if (size(#)>1) 2626 2626 { 2627 2627 if (typeof(#[2])=="int" || typeof(#[2])=="number") 2628 2628 { 2629 2630 2631 2629 whichengine = int(#[2]); 2630 } 2631 if (size(#)>2) 2632 2632 { 2633 2633 if (typeof(#[3])=="int" || typeof(#[3])=="number") 2634 2634 { 2635 2636 2637 2638 2635 l0 = int(#[3]); 2636 l0set = 1; 2637 } 2638 if (size(#)>3) 2639 2639 { 2640 2640 if (typeof(#[4])=="ideal") 2641 2641 { 2642 2643 2644 2645 2646 2642 G = #[4]; 2643 Gset = 1; 2644 } 2645 } 2646 } 2647 2647 } 2648 2648 } … … 2658 2658 if (v == intvec(0)) 2659 2659 { 2660 2660 v[1] = i; 2661 2661 } 2662 2662 else 2663 2663 { 2664 2664 v[size(v)+1] = i; 2665 2665 } 2666 2666 } … … 2711 2711 while (p<N[j]) 2712 2712 { 2713 2713 j++; 2714 2714 } 2715 2715 N = insertGenerator(N,p,j+1); … … 2762 2762 if(typeof(#[2])=="int" || typeof(#[2])=="number") 2763 2763 { 2764 2765 2764 l0 = int(#[2]); 2765 l0given = 1; 2766 2766 } 2767 2767 } -
Singular/LIB/ncfactor.lib
r3561bd rc623f27 1268 1268 w = D; D=X; X=w; 1269 1269 } 1270 // DONE with assumptions 1270 // DONE with assumptions 1271 1271 //Input successfully checked 1272 1272 intvec lexpofX = leadexp(X); -
Singular/LIB/normal.lib
r3561bd rc623f27 4363 4363 // I = I1 \cap I2 4364 4364 printlevel = printlevel + 1; 4365 ideal JDefault = 0; 4365 ideal JDefault = 0; // Now it uses the default J; 4366 4366 list nor1 = normalM(Id1, decomp, withDelta, denomOption, JDefault, JDefault)[1]; 4367 4367 list nor2 = normalM(Id2, decomp, withDelta, denomOption, JDefault, JDefault)[1]; -
Singular/LIB/primitiv.lib
r3561bd rc623f27 226 226 countz++; 227 227 if (countz>=grad2) 228 228 { ERROR("No primitive element found!! This should NEVER happen!"); } 229 229 } 230 230 g = g +x^countx *z^countz; … … 234 234 mat=coeffs(mat[countz+1,1],x); 235 235 if (size(mat)>countx) 236 236 { 237 237 if (mat[countx+1,1] != 0) 238 238 { 239 239 found=1; // d.h. hier: neues g gefunden 240 240 }}}
Note: See TracChangeset
for help on using the changeset viewer.