Changeset b8ad3a in git for Singular/LIB/puiseuxexpansions.lib
- Timestamp:
- Jan 4, 2024, 4:39:13 PM (4 months ago)
- Branches:
- (u'spielwiese', '98550b669234b32be762076c32b3be2c35188ac4')
- Children:
- 4b0877fd8b06e461f8c2450538a142d631260806
- Parents:
- dff47b7996ab7bb15ce9bd372c7fc531a02308f7
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/puiseuxexpansions.lib
rdff47b rb8ad3a 479 479 @c = puiseuxStep(fTemp, slN, 1); 480 480 fc = factorize(@c); 481 481 482 482 for(j = 2; j <= size(fc[1]); j++) 483 483 { … … 505 505 { 506 506 sizeErg = size(erg); 507 } else 507 } else 508 508 { 509 509 sizeErg = 1; … … 538 538 { 539 539 cs[size(cs) + 1] = list((csT[k][1] + c1)*(var(1)^(csT[k][2]*slN)), csT[k][2]*slD); 540 } else 540 } else 541 541 { 542 542 cs[size(cs) + 1] = list(csT[k][1] + c1, csT[k][2]*slD); … … 547 547 // NEW CODING FOR CONJUGATED EXPANSIONS 548 548 cs[size(cs)][6] = insert(csT[k][6], cod); 549 549 550 550 if(newExt == 1) 551 551 { … … 565 565 } 566 566 } 567 } else 567 } else 568 568 { 569 569 def RT = csT[k]; … … 928 928 929 929 fNew1 = subst(f, var(1), var(1)^slD); 930 930 931 931 for(ii = 2; ii <= size(fJe[1]); ii++) 932 932 { … … 948 948 @c = puiseuxStep(fTemp, slN, 1); 949 949 fc = factorize(@c); 950 950 951 951 // Stopping criterium 952 952 stop = 0; … … 957 957 stop = 1; 958 958 } 959 } else 959 } else 960 960 { 961 961 if(slN >= slD * maxDeg) … … 968 968 { 969 969 //minP = fc[1][j]; 970 970 971 971 // The minpoly of Y^5-1 is Y-1 and the minpoly of Y^5-2 is Y^5-2. 972 972 minPBefore = subst(fJe[1][ii],var(1),1); … … 981 981 { 982 982 sizeErg = size(erg); 983 } else 983 } else 984 984 { 985 985 sizeErg = 1; … … 1014 1014 { 1015 1015 cs[size(cs) + 1] = list((csT[k][1] + c1)*(var(1)^(csT[k][2]*slN)), csT[k][2]*slD); 1016 } else 1016 } else 1017 1017 { 1018 1018 cs[size(cs) + 1] = list(csT[k][1] + c1, csT[k][2]*slD); … … 1023 1023 // NEW CODING FOR CONJUGATED EXPANSIONS 1024 1024 cs[size(cs)][6] = insert(csT[k][6], cod); 1025 1025 1026 1026 if(newExt == 1) 1027 1027 { … … 1039 1039 } 1040 1040 } 1041 } else 1041 } else 1042 1042 { 1043 1043 def RT = csT[k];
Note: See TracChangeset
for help on using the changeset viewer.