Changeset 08e898 in git


Ignore:
Timestamp:
Nov 14, 2022, 9:03:14 PM (3 months ago)
Author:
Johannes Hoffmann <johannes.hoffmann@…>
Branches:
(u'spielwiese', 'ad2543eab51733612ba7d118afc77edca719600e')
Children:
4df27e2988096e614ab54338dbdd2bb990a110e3
Parents:
0a4a20dfba04555d929927f769e18fd086517dd8
Message:
fixed edge case of total localization
File:
1 edited

Legend:

Unmodified
Added
Removed
  • Singular/LIB/olga.lib

    r0a4a20d r08e898  
    404404            locData = normalizeRational(locData);
    405405        }
    406         ideal I = p;
    407         intvec modLocData = intvecComplement(locData, 1..nvars(basering));
    408         I = eliminateNC(I, modLocData);
    409         if (size(I)) {
    410             return(1);
     406        int n = nvars(basering);
     407        if (size(locData) < n) { // there are variables to eliminate
     408            ideal I = p;
     409            intvec modLocData = intvecComplement(locData, 1..nvars(basering));
     410            I = eliminateNC(I, modLocData);
     411            if (size(I)) {
     412                return(1);
     413            }
     414        } else {
     415            return(p != 0);
    411416        }
    412417    }
     
    14711476        ERROR("Weyl rational isInS negative failed");
    14721477    }
     1478    intvec w = 4,2,3,4,1;
     1479    if (!isInS(x*y*Dx*Dy,2,w)) {
     1480        ERROR("Weyl total rational isInS positive failed");
     1481    }
     1482    if (isInS(0, 2, w)) {
     1483        ERROR("Weyl total rational isInS negative failed");
     1484    }
    14731485    print("    isInS OK");
    14741486}
Note: See TracChangeset for help on using the changeset viewer.