PurityPaperMora.g

#This is the first example shown in "An Axiomatic Setup for Algorithmic Homological Algebra and an Alternative Approach to Localization"
#in the version where Mora's algorithm fails (using ByASmallerPresentation)

LoadPackage( "RingsForHomalg" );;
R := HomalgFieldOfRationalsInDefaultCAS( ) *  "a,b,c,d,e";;
K := HomalgMatrix( "[\
2*a+c+d+e-2,2*a+c+d+e-2,2*a+c+d+e-2,0,\
2*c*d+d^2+2*c*e+2*d*e-3*e^2-2*c-d-6*e+1,\
2*c*d+d^2+2*c*e+2*d*e-3*e^2-2*c-d-6*e+1,\
2*c*d+d^2+2*c*e+2*d*e-3*e^2-2*c-d-6*e+1,0,\
-4*a+2*b-c-d-e+2,-4*a+2*b-c-d-e+2,\
-c+d+e+2,4*a*d-2*b*d+2*d^2+2*d*e,\
c^2-d-1,c^2-d-1,c^2-d-1,0,\
4*d*e^2-d^2+2*c*e+4*d*e-3*e^2-2*c+3*d-6*e+5,\
4*d*e^2-d^2+2*c*e+4*d*e-3*e^2-2*c+3*d-6*e+5,\
4*d*e^2-d^2+2*c*e+4*d*e-3*e^2-2*c+3*d-6*e+5,0,\
0,b^2+a+c+d+e,0,b^2*e+a*e+c*e+d*e+e^2,\
0,b^2*d+a*d+c*d+d^2+d*e,0,0,\
0,a*b^2+a^2+a*c+a*d+a*e,0,0,\
4*b^3*d-4*d^3-12*d^2*e-32*c*e^2+12*e^3+21*d^2\
-42*c*e+40*d*e+27*e^2+8*c-9*d+16*e-17,\
-4*a*b*d-4*b*c*d-4*b*d^2-4*d^3-4*b*d*e-12*d^2*e-32*c*e^2+\
12*e^3+21*d^2-42*c*e+40*d*e+27*e^2+8*c-9*d+16*e-17,\
-12*d^2*e-32*c*e^2+12*e^3+21*d^2-42*c*e+44*d*e+\
27*e^2+8*c-9*d+16*e-17,-4*b^3+4*d^2+4*e\
]", 9 , 4 , R );;

LoadPackage( "LocalizeRingForHomalg" );;
R0Mora:=LocalizePolynomialRingAtZeroWithMora(R);

LoadPackage( "Modules" );
K1:=LeftPresentation(K);

K0:=R0Mora*K1;

ByASmallerPresentation(K0);

Display(K0);

filt0:=PurityFiltration(K0);

Display(filt0);
Display(Source(IsomorphismOfFiltration(filt0)));