benruijl / forcer Goto Github PK

test.frm
__________________________________________________________
#--
table sparse rtop(2);
table sparse top(1);
table sparse c(2);
Auto S x;
S w,d;
CF rat,cc;

PolyRatFun rat;

* Include Transformations 
Fill rtop(1,1) = + cc(1,1)*top(2) 
                 + cc(1,2)*top(3) 
                 + cc(1,3)*top(4);
Fill rtop(2,4) = + cc(2,1)*top(3);
Fill rtop(2,2) = + cc(2,2)*top(5);
Fill rtop(3,6) = top(7);
Fill rtop(3,3) = top(5);

* Include Coefficients 
Fill c(1,1) = rat((-2*d^3+14*d^2-32*d+24),(6*d^3-69*d^2+264*d-336));
Fill c(1,2) = rat(((50*d^4-550*d^3+2248*d^2-4040*d+2688)*w-40*d^4+452*d^3-1892*d^2+3472*d-2352),
                            ((18*d^4-267*d^3+1482*d^2-3648*d+3360)*w-18*d^4+267*d^3-1482*d^2+3648*d-3360));
Fill c(1,3) = rat(((2*d^3-14*d^2+32*d-24)*w),((18*d^3-195*d^2+702*d-840)*w-18*d^3+195*d^2-702*d+840));
Fill c(2,1) = rat((-5*d+14),((d-4)*w));
Fill c(2,2) = rat((5*d-14),(3*d-10));

L F = top(1);

#do i=1,3
    id top(?a) = rtop(`i',?a);
    id cc(?a) = c(?a);
    id rtop(x1?,x2?) = top(x2);
    .sort
#enddo
print +s;
.end