Quelle  normaliz.gi   Sprache: Shell

##
InstallMethod( ViewString,
               "for a Normaliz cone",
               [ IsNormalizCone ],
function( r )
    # TODO: May print more information when present
    return "<a Normaliz cone>";
end );

BindGlobal("_NmzRecordOfConeProperties", _NmzConePropertiesNamesRecord());

InstallGlobalFunction("NmzConeProperty",
function(cone, prop)
    local result, t, poly, tmp, denom;
    result := _NmzConeProperty(cone, prop);
    if prop = _NmzRecordOfConeProperties.Grading or prop = "Grading" then
        denom := NmzConeProperty(cone, _NmzRecordOfConeProperties.GradingDenom);;
        return result / denom;
    fi;
    if prop = _NmzRecordOfConeProperties.HilbertSeries or prop = "HilbertSeries" then
        t := Indeterminate(Integers, "t");
        poly := UnivariatePolynomial(Integers, result[1], t);
        poly := poly * t^result[3];
        tmp := Collected(result[2]);
        return [poly, tmp];
    fi;
    if prop = _NmzRecordOfConeProperties.HilbertQuasiPolynomial or prop = "HilbertQuasiPolynomial" then
        t := Indeterminate(Rationals, "t");
        denom := Remove(result);
        poly := List(result, coeffs -> UnivariatePolynomial(Rationals, coeffs, t));
        return poly / denom;
    fi;
    return result;
end);

InstallGlobalFunction("NmzPrintConeProperties",
function(cone)
    local prop, val;
    if not IsNormalizCone(cone) then
        Error("First argument must be a Normaliz cone object");
        return;
    fi;
    for prop in NmzKnownConeProperties(cone) do
        val := NmzConeProperty(cone, prop);
        Print(prop," = ");
        if IsMatrix(val) then
            Print("\n");
        fi;
        Display(val);
    od;
end);

BindGlobal("_NmzPrintSomeConeProperties",
function(cone, excluded)
    local prop, val;
    if not IsNormalizCone(cone) then
        Error("First argument must be a Normaliz cone object");
        return;
    fi;
    excluded := Union(excluded, [ "ExplicitHilbertSeries", "Triangulation" ]);
    # iterate over properties in alphabetical order, to reduce
    # fluctuations between Normaliz versions
    for prop in NmzKnownConeProperties(cone) do
        if prop in excluded then
            continue;
        fi;
        val := NmzConeProperty(cone, prop);
        Print(prop," = ");
        if IsMatrix(val) then
            Print("\n");
        fi;
        Display(val);
    od;
end);

InstallGlobalFunction("NmzBasisChange",
function(cone)
    local result;
    result := NmzConeProperty( cone, "Sublattice" );
    return rec(
        Embedding := result[1],
        Projection := result[2],
        Annihilator := result[3],
        );
end);


#
#
#
InstallGlobalFunction("NmzCone",
function(arg)
    local func, opts_rec, opts_list, cone;
    if Length(arg) = 1 then
        if IsList(arg[1]) then
            opts_list := arg[1];
        elif IsRecord(arg[1]) then
            opts_list := [];
            for opts_rec in RecNames(arg[1]) do
              Add(opts_list, opts_rec);
              Add(opts_list, arg[1].(opts_rec));
            od;
        else
            # TODO: better error message
            Error("Unsupported input");
        fi;
    else
        opts_list := arg;
    fi;

    cone := _NmzCone(opts_list);
    return cone;
end);

# TODO: extend NmzCone to allow this syntax:
##cone := NmzCone(rec(integral_closure := M, grading := [ 0, 0, 1 ]));;



#
#
#
InstallGlobalFunction("NmzCompute",
function(arg)
    local cone, propsToCompute;
    if not Length(arg) in [1,2] then
        Error("Wrong number of arguments, expected 1 or 2");
        return fail;
    fi;
    cone := arg[1];
    if not IsNormalizCone(cone) then
        Error("First argument must be a Normaliz cone object");
        return fail;
    fi;
    if Length(arg) = 1 then
        propsToCompute := [];
        if ValueOption("dual") = true then Add(propsToCompute, "DualMode"); fi;
        if ValueOption("DualMode") = true then Add(propsToCompute, "DualMode"); fi;
        if ValueOption("HilbertBasis") = true then Add(propsToCompute, "HilbertBasis"); fi;
        # TODO: add more option names? or just support arbitrary ones, by using
        # iterating over the (undocumented!) OptionsStack???
        if Length(propsToCompute) = 0 then
            propsToCompute := [ "DefaultMode" ];
        fi;
    else
        if IsString(arg[2]) then
            propsToCompute := [arg[2]];
        else
            propsToCompute := arg[2];
        fi;
    fi;
    return _NmzCompute(cone, propsToCompute);
end);