Changeset 8a57d4 in git for Singular/LIB/fpaprops.lib


Ignore:
Timestamp:
Mar 15, 2018, 3:51:09 PM (6 years ago)
Author:
Karim Abou Zeid <karim23697@…>
Branches:
(u'spielwiese', 'd0474371d8c5d8068ab70bfb42719c97936b18a6')
Children:
634a329f13a9b02b22a7736c597585dbbe8a2a48
Parents:
7e5164822dfad56a986d4d205e77a60f62611fc3
Message:
Fix lpSubstitute for '1'
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/fpaprops.lib

    r7e51648 r8a57d4  
    10511051    for (int j = 1; j <= size(ivfi); j++) {
    10521052      int varindex = ivfi[j];
    1053       int subindex = lpIndexOf(s1, var(varindex));
    1054       if (subindex > 0) {
    1055         s2[subindex] = lpNF(s2[subindex],G);
    1056         fis = lpMult(fis, s2[subindex]);
    1057       } else {
    1058         fis = lpMult(fis, lpNF(iv2lp(varindex),G));
    1059       }
    1060       /*fis = lpNF(fis,G);*/
    1061       kill varindex; kill subindex;
     1053      if (varindex > 0) {
     1054        int subindex = lpIndexOf(s1, var(varindex));
     1055        if (subindex > 0) {
     1056          s2[subindex] = lpNF(s2[subindex],G);
     1057          fis = lpMult(fis, s2[subindex]);
     1058        } else {
     1059          fis = lpMult(fis, lpNF(iv2lp(varindex),G));
     1060        }
     1061        /*fis = lpNF(fis,G);*/
     1062        kill subindex;
     1063      }
     1064      kill varindex;
    10621065    } kill j;
    10631066    kill ivfi;
     
    11501153    for (int j = 1; j <= size(ivfi); j++) {
    11511154      int varindex = ivfi[j];
    1152       int subindex = lpIndexOf(s1, var(varindex));
    1153       if (subindex > 0) {
    1154         tmpDegBound = tmpDegBound + deg(s2[subindex]);
    1155       } else {
    1156         tmpDegBound = tmpDegBound + 1;
    1157       }
    1158       kill varindex; kill subindex;
     1155      if (varindex > 0) {
     1156        int subindex = lpIndexOf(s1, var(varindex));
     1157        if (subindex > 0) {
     1158          tmpDegBound = tmpDegBound + deg(s2[subindex]);
     1159        } else {
     1160          tmpDegBound = tmpDegBound + 1;
     1161        }
     1162        kill subindex;
     1163      }
     1164      kill varindex;
    11591165    } kill j;
    11601166    if (tmpDegBound > maxDegBound) {
Note: See TracChangeset for help on using the changeset viewer.