Changeset 364784a in git
 Timestamp:
 Jan 14, 2007, 11:09:42 PM (17 years ago)
 Branches:
 (u'fiekerDuVal', '117eb8c30fc9e991c4decca4832b1d19036c4c65')(u'spielwiese', 'b21a664aa22dc6e196223af8a74ad4885e83547c')
 Children:
 9baa295fa0016bec13c69de439789cba6e7053d2
 Parents:
 7b8cb0e2c52093f9cb60de8aeb4a829b4a0cac0d
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

Singular/LIB/control.lib
r7b8cb0 r364784a 1 version="$Id: control.lib,v 1.3 5 20060720 22:17:11levandov Exp $";1 version="$Id: control.lib,v 1.36 20070114 22:09:42 levandov Exp $"; 2 2 category="System and Control Theory"; 3 3 info=" … … 660 660 RETURN: list (of strings) 661 661 PURPOSE: determine parametric expressions which have been assumed to be nonzero in the process of computing the Groebner basis 662 NOTE: the output list consists of two strings. The first string contains the variables only, whereas the second string contains polynomials.662 NOTE: the output list consists of strings. The first string contains the variables only, whereas each further string contain a single polynomial in parameters. 663 663 $* We strongly recommend to switch on the redSB and redTail options. 664 664 @* The procedure is effective with the lift procedure for modules with parameters … … 817 817 ideal I = ideal(M); 818 818 I = simplify(I,2+4); 819 int s = size(I);819 int s = ncols(I); 820 820 ideal Den; 821 821 poly p; … … 833 833 string newvars = parstr(basering); 834 834 def save = basering; 835 string NewRing = "ring @NR = " +string(char(basering))+",("+newvars+"),Dp;";835 string NewRing = "ring @NR =(" +string(char(basering))+"),("+newvars+"),Dp;"; 836 836 execute(NewRing); 837 837 ideal F; 838 838 ideal Den = imap(save,Den); 839 839 Den = simplify(Den,2); 840 int s1 = size(Den);840 int s1 = ncols(Den); 841 841 for (i=1; i<=s1; i++) 842 842 { 843 if ( Den[i]!=1)843 if (normalize(Den[i]) !=1) 844 844 { 845 845 F= F, factorize(Den[i],1); … … 848 848 F = simplify(F, 2+4+8); 849 849 ideal @L = F; 850 @L = simplify(@L,2); 850 851 list SL; 851 852 int c,j; 852 853 string Mono; 853 854 c = 1; 854 for (j=1; j<= size(@L);j++)855 for (j=1; j<= ncols(@L);j++) 855 856 { 856 857 if (leadcoef(@L[j]) <0) … … 858 859 @L[j] = 1*@L[j]; 859 860 } 860 if ( (@L[j]  lead(@L[j]))==0) //@L[j] is a monomial861 if (((@L[j]  lead(@L[j]))==0 ) && (@L[j]!=0) ) //@L[j] is a monomial 861 862 { 862 863 Mono = Mono + string(@L[j])+ ","; // concatenation … … 870 871 if (Mono!="") 871 872 { 872 Mono = Mono[1..size(Mono)1]; // delete the last semicolon873 Mono = Mono[1..size(Mono)1]; // delete the last colon 873 874 } 874 875 SL[1] = Mono;
Note: See TracChangeset
for help on using the changeset viewer.