/*
 * $Id$
 */

package edu.jas.application;


import java.io.IOException;
import java.io.Reader;
import java.io.StringReader;
import java.util.ArrayList;
import java.util.List;

import org.apache.log4j.BasicConfigurator;

import edu.jas.arith.BigDecimal;
import edu.jas.arith.BigInteger;
import edu.jas.arith.BigRational;
import edu.jas.arith.ModInteger;
import edu.jas.arith.ModIntegerRing;
import edu.jas.arith.Product;
import edu.jas.arith.ProductRing;
import edu.jas.gb.GroebnerBase;
import edu.jas.gbufd.GBFactory;
import edu.jas.gbufd.RGroebnerBasePseudoSeq;
import edu.jas.gbufd.RReductionSeq;
import edu.jas.kern.ComputerThreads;
import edu.jas.poly.GenPolynomial;
import edu.jas.poly.GenPolynomialRing;
import edu.jas.poly.GenPolynomialTokenizer;
import edu.jas.poly.PolynomialList;
import edu.jas.poly.TermOrder;


/**
 * Examples for application usage.
 * @author Christoph Zengler.
 * @author Heinz Kredel.
 */

public class Examples {


    /**
     * main.
     */
    public static void main(String[] args) {
        BasicConfigurator.configure();
        if (args.length > 0) {
            example1();
            example2();
            example3();
            example4();
        }
        example5();
        example6();
        ComputerThreads.terminate();
    }


    /**
     * example1. cyclic n-th roots polynomial systems.
     */
    public static void example1() {
        int n = 4;

        BigInteger fac = new BigInteger();
        GenPolynomialRing ring = new GenPolynomialRing(fac, n); //,var);
        System.out.println("ring = " + ring + "\n");

        List> cp = new ArrayList>(n);
        for (int i = 1; i <= n; i++) {
            GenPolynomial p = cyclicPoly(ring, n, i);
            cp.add(p);
            System.out.println("p[" + i + "] = " + p);
            System.out.println();
        }
        System.out.println("cp = " + cp + "\n");

        List> gb;
        //GroebnerBase sgb = new GroebnerBaseSeq();
        GroebnerBase sgb = GBFactory.getImplementation(fac);
        gb = sgb.GB(cp);
        System.out.println("gb = " + gb);

    }


    static GenPolynomial cyclicPoly(GenPolynomialRing ring, int n, int i) {

        List> X = /*(List>)*/ring
                        .univariateList();

        GenPolynomial p = ring.getZERO();
        for (int j = 1; j <= n; j++) {
            GenPolynomial pi = ring.getONE();
            for (int k = j; k < j + i; k++) {
                pi = pi.multiply(X.get(k % n));
            }
            p = p.sum(pi);
            if (i == n) {
                p = p.subtract(ring.getONE());
                break;
            }
        }
        return p;
    }


    /**
     * example2. abtract types:
     * List>>>.
     */
    public static void example2() {
        List>>> L = null;
        L = new ArrayList>>>();

        BigRational bfac = new BigRational(1);
        GenPolynomialRing pfac = null;
        pfac = new GenPolynomialRing(bfac, 3);

        List> F = null;
        F = new ArrayList>();

        GenPolynomial p = null;
        for (int i = 0; i < 2; i++) {
            p = pfac.random(5, 4, 3, 0.4f);
            if (!p.isConstant()) {
                F.add(p);
            }
        }
        //System.out.println("F = " + F);

        Ideal id = new Ideal(pfac, F);
        id.doGB();
        if (id.isONE() || id.isZERO()) {
            System.out.println("id zero or one = " + id);
            return;
        }
        ResidueRing rr = new ResidueRing(id);
        System.out.println("rr = " + rr);

        ProductRing> pr = null;
        pr = new ProductRing>(rr, 3);

        String[] vars = new String[] { "a", "b" };
        GenPolynomialRing>> fac;
        fac = new GenPolynomialRing>>(pr, 2, vars);

        GenPolynomial>> pp;
        for (int i = 0; i < 1; i++) {
            pp = fac.random(2, 4, 4, 0.4f);
            if (!pp.isConstant()) {
                L.add(pp);
            }
        }
        System.out.println("L = " + L);

        //PolynomialList>> Lp = null;
        //Lp = new PolynomialList>>(fac,L);
        //System.out.println("Lp = " + Lp);

        GroebnerBase>> bb = new RGroebnerBasePseudoSeq>>(
                        pr);

        System.out.println("isGB(L) = " + bb.isGB(L));

        List>>> G = null;

        G = bb.GB(L);
        System.out.println("G = " + G);
        System.out.println("isGB(G) = " + bb.isGB(G));
    }


    /**
     * example3. abtract types: GB of List>>.
     */
    public static void example3() {
        List>> L = null;
        L = new ArrayList>>();

        BigRational bfac = new BigRational(1);
        GenPolynomialRing pfac = null;
        pfac = new GenPolynomialRing(bfac, 2);

        List> F = null;
        F = new ArrayList>();

        GenPolynomial p = null;
        for (int i = 0; i < 2; i++) {
            p = pfac.random(5, 5, 5, 0.4f);
            //p = pfac.parse("x0^2 -2" );
            if (!p.isConstant()) {
                F.add(p);
            }
        }
        //System.out.println("F = " + F);

        Ideal id = new Ideal(pfac, F);
        id.doGB();
        if (id.isONE() || id.isZERO()) {
            System.out.println("id zero or one = " + id);
            return;
        }
        ResidueRing rr = new ResidueRing(id);
        System.out.println("rr = " + rr);

        String[] vars = new String[] { "a", "b" };
        GenPolynomialRing> fac;
        fac = new GenPolynomialRing>(rr, 2, vars);

        GenPolynomial> pp;
        for (int i = 0; i < 2; i++) {
            pp = fac.random(2, 4, 6, 0.2f);
            if (!pp.isConstant()) {
                L.add(pp);
            }
        }
        System.out.println("L = " + L);

        GroebnerBase> bb;
        //bb = new GroebnerBasePseudoSeq>(rr);
        bb = GBFactory.getImplementation(rr);

        System.out.println("isGB(L) = " + bb.isGB(L));

        List>> G = null;

        G = bb.GB(L);
        System.out.println("G = " + G);
        System.out.println("isGB(G) = " + bb.isGB(G));
    }


    /**
     * example4. abtract types: comprehensive GB of
     * List>>.
     */
    public static void example4() {
        int kl = 2;
        int ll = 3;
        int el = 3;
        float q = 0.2f; //0.4f
        GenPolynomialRing cfac;
        GenPolynomialRing> fac;

        List>> L;

        ComprehensiveGroebnerBaseSeq bb;

        GenPolynomial> a;
        GenPolynomial> b;
        GenPolynomial> c;

        BigRational coeff = new BigRational(kl);
        String[] cv = { "a", "b" };
        cfac = new GenPolynomialRing(coeff, 2, cv);
        String[] v = { "x", "y" };
        fac = new GenPolynomialRing>(cfac, 2, v);
        bb = new ComprehensiveGroebnerBaseSeq(coeff);

        L = new ArrayList>>();

        a = fac.random(kl, ll, el, q);
        b = fac.random(kl, ll, el, q);
        c = a; //c = fac.random(kl, ll, el, q );

        if (a.isZERO() || b.isZERO() || c.isZERO()) {
            return;
        }

        L.add(a);
        System.out.println("CGB exam L = " + L);
        L = bb.GB(L);
        System.out.println("CGB( L )   = " + L);
        System.out.println("isCGB( L ) = " + bb.isGB(L));

        L.add(b);
        System.out.println("CGB exam L = " + L);
        L = bb.GB(L);
        System.out.println("CGB( L )   = " + L);
        System.out.println("isCGB( L ) = " + bb.isGB(L));

        L.add(c);
        System.out.println("CGB exam L = " + L);
        L = bb.GB(L);
        System.out.println("CGB( L )   = " + L);
        System.out.println("isCGB( L ) = " + bb.isGB(L));
    }


    /**
     * example5. comprehensive GB of
     * List>> and GB for regular ring.
     */
    public static void example5() {
        int kl = 2;
        int ll = 4;
        int el = 3;
        float q = 0.3f; //0.4f
        GenPolynomialRing cfac;
        GenPolynomialRing> fac;

        List>> L;

        ComprehensiveGroebnerBaseSeq bb;

        GenPolynomial> a;
        GenPolynomial> b;
        GenPolynomial> c;

        BigRational coeff = new BigRational(kl);
        String[] cv = { "a", "b" };
        cfac = new GenPolynomialRing(coeff, 2, cv);
        String[] v = { "x", "y" };
        fac = new GenPolynomialRing>(cfac, 2, v);
        bb = new ComprehensiveGroebnerBaseSeq(coeff);

        L = new ArrayList>>();

        a = fac.random(kl, ll, el, q);
        b = fac.random(kl, ll, el, q);
        c = a; //c = fac.random(kl, ll, el, q );

        if (a.isZERO() || b.isZERO() || c.isZERO()) {
            return;
        }

        L.add(a);
        L.add(b);
        L.add(c);
        System.out.println("CGB exam L = " + L);
        GroebnerSystem sys = bb.GBsys(L);
        boolean ig = bb.isGB(sys.getCGB());
        System.out.println("CGB( L )   = " + sys.getCGB());
        System.out.println("isCGB( L ) = " + ig);

        List>>> Lr, bLr;
        RReductionSeq>> res = new RReductionSeq>>();

        Lr = PolyUtilApp. toProductRes(sys.list);
        bLr = res.booleanClosure(Lr);

        System.out.println("booleanClosed(Lr)   = " + bLr);

        if (bLr.size() > 0) {
            GroebnerBase>> rbb = new RGroebnerBasePseudoSeq>>(
                            bLr.get(0).ring.coFac);
            System.out.println("isRegularGB(Lr) = " + rbb.isGB(bLr));
        }
    }


    /**
     * Example GBase and real root.
     */
    @SuppressWarnings("unchecked")
    public static void example6() {
        BigRational coeff = new BigRational();
        GroebnerBase gb = GBFactory.getImplementation(coeff);

        String exam = "(x,y,z) L " + "( " + "( x^2 - 2 ), ( y^2 - 3 ), ( z^2 + x * y )" + ") ";
        Reader source = new StringReader(exam);
        GenPolynomialTokenizer parser = new GenPolynomialTokenizer(source);
        PolynomialList F = null;

        try {
            F = (PolynomialList) parser.nextPolynomialSet();
        } catch (ClassCastException e) {
            e.printStackTrace();
            return;
        } catch (IOException e) {
            e.printStackTrace();
            return;
        }
        System.out.println("F = " + F);

        List> G = gb.GB(F.list);

        PolynomialList Gp = new PolynomialList(F.ring, G);
        System.out.println("G = " + Gp);

        // compute real roots of the ideal
        Ideal I = new Ideal(Gp);
        List> Ir = PolyUtilApp. realAlgebraicRoots(I);
        for (IdealWithRealAlgebraicRoots R : Ir) {
            R.doDecimalApproximation();
            for (List Dr : R.decimalApproximation()) {
                System.out.println(Dr.toString());
            }
            System.out.println();
        }
    }


    /**
     * example7. Coefficients in Boolean residue class ring.
     */
    public static void example7() {
        String[] vars = { "v3", "v2", "v1" };

        ModIntegerRing z2 = new ModIntegerRing(2);
        GenPolynomialRing z2p = new GenPolynomialRing(z2, vars.length, new TermOrder(
                        TermOrder.INVLEX), vars);
        List> fieldPolynomials = new ArrayList>();

        //add v1^2 + v1, v2^2 + v2, v3^2 + v3 to fieldPolynomials
        for (int i = 0; i < vars.length; i++) {
            GenPolynomial var = z2p.univariate(i);
            fieldPolynomials.add(var.multiply(var).sum(var));
        }


        Ideal fieldPolys = new Ideal(z2p, fieldPolynomials);
        ResidueRing ring = new ResidueRing(fieldPolys);
        String[] mvars = { "mv3", "mv2", "mv1" };
        GenPolynomialRing> ringp = new GenPolynomialRing>(ring,
                        mvars.length, mvars);

        List>> polynomials = new ArrayList>>();

        GenPolynomial> v1 = ringp.univariate(0);
        GenPolynomial> v2 = ringp.univariate(1);
        GenPolynomial> v3 = ringp.univariate(2);
        GenPolynomial> notV1 = v1.sum(ringp.ONE);
        GenPolynomial> notV2 = v2.sum(ringp.ONE);
        GenPolynomial> notV3 = v3.sum(ringp.ONE);

        //v1*v2
        GenPolynomial> p1 = v1.multiply(v2);

        //v1*v2 + v1 + v2 + 1
        GenPolynomial> p2 = notV1.multiply(notV2);

        //v1*v3 + v1 + v3 + 1
        GenPolynomial> p3 = notV1.multiply(notV3);

        polynomials.add(p1);
        polynomials.add(p2);
        polynomials.add(p3);

        //GroebnerBase> gb = new GroebnerBasePseudoSeq>(ring);
        GroebnerBase> gb = GBFactory.getImplementation(ring);
        List>> G = gb.GB(polynomials);

        System.out.println(G);
    }


    /**
     * example8. Coefficients in Boolean residue class ring with cuppling of
     * variables.
     */
    public static void example8() {
        String[] vars = { "v3", "v2", "v1" };

        ModIntegerRing z2 = new ModIntegerRing(2);
        GenPolynomialRing z2p = new GenPolynomialRing(z2, vars.length, new TermOrder(
                        TermOrder.INVLEX), vars);
        List> fieldPolynomials = new ArrayList>();

        //add v1^2 + v1, v2^2 + v2, v3^2 + v3 to fieldPolynomials
        for (int i = 0; i < vars.length; i++) {
            GenPolynomial var = z2p.univariate(i);
            fieldPolynomials.add(var.multiply(var).sum(var));
        }


        Ideal fieldPolys = new Ideal(z2p, fieldPolynomials);
        ResidueRing ring = new ResidueRing(fieldPolys);
        String[] mvars = { "mv3", "mv2", "mv1" };
        GenPolynomialRing> ringp = new GenPolynomialRing>(ring,
                        mvars.length, mvars);

        List>> polynomials = new ArrayList>>();

        GenPolynomial> v1 = ringp.univariate(0);
        GenPolynomial> v2 = ringp.univariate(1);
        GenPolynomial> v3 = ringp.univariate(2);
        GenPolynomial> notV1 = v1.sum(ringp.ONE);
        GenPolynomial> notV2 = v2.sum(ringp.ONE);
        GenPolynomial> notV3 = v3.sum(ringp.ONE);

        //v1*v2
        GenPolynomial> p1 = v1.multiply(v2);

        //v1*v2 + v1 + v2 + 1
        GenPolynomial> p2 = notV1.multiply(notV2);

        //v1*v3 + v1 + v3 + 1
        GenPolynomial> p3 = notV1.multiply(notV3);

        polynomials.add(p1);
        polynomials.add(p2);
        polynomials.add(p3);

        List> gens = ring.generators();
        System.out.println("gens = " + gens);
        GenPolynomial> mv3v3 = v3.subtract(gens.get(1));
        GenPolynomial> mv2v2 = v2.subtract(gens.get(2));
        GenPolynomial> mv1v1 = v1.subtract(gens.get(3));

        System.out.println("mv3v3 = " + mv3v3);
        System.out.println("mv2v2 = " + mv2v2);
        System.out.println("mv1v1 = " + mv1v1);

        polynomials.add(mv3v3);
        polynomials.add(mv2v2);
        polynomials.add(mv1v1);

        //GroebnerBase> gb = new GroebnerBasePseudoSeq>(ring);
        GroebnerBase> gb = GBFactory.getImplementation(ring);

        List>> G = gb.GB(polynomials);

        System.out.println(G);
    }


    /**
     * example9. Groebner base and dimension.
     */
    public static void example9() {
        String[] vars = { "d1", "d2", "d3", "p1a", "p1b", "p1c", "p2a", "p2b", "p2c", "p3a", "p3b", "p3c",
                "p4a", "p4b", "p4c", "A", "B", "C", "D" };

        BigRational br = new BigRational();
        GenPolynomialRing pring = new GenPolynomialRing(br, vars);
        //GenPolynomialRing pring = new GenPolynomialRing(br, vars.length, new TermOrder(TermOrder.INVLEX), vars);

        GenPolynomial e1 = pring.parse("A*p1a+B*p1b+C*p1c+D"); // (1)
        GenPolynomial e2 = pring.parse("A*p2a+B*p2b+C*p2c+D"); // (2)
        GenPolynomial e3 = pring.parse("A*p3a+B*p3b+C*p3c+D"); // (3)
        GenPolynomial e4 = pring.parse("A*p4a+B*p4b+C*p4c+D"); // (4)
        GenPolynomial e5 = pring.parse("p2a-p3a"); // (5)
        GenPolynomial e6 = pring.parse("p2b-p3b"); // (6)
        GenPolynomial e7 = pring.parse("p2c-p3c"); // (7)
        GenPolynomial e8 = pring.parse("(p2a-p1a)^2+(p2b-p1b)^2+(p2c-p1c)^2-d1^2"); // (8)
        GenPolynomial e9 = pring.parse("(p4a-p3a)^2+(p4b-p3b)^2+(p4c-p3c)^2-d2^2"); // (9)

        List> cp = new ArrayList>(9);
        cp.add(e1);
        cp.add(e2);
        cp.add(e3);
        cp.add(e4);
        cp.add(e5);
        cp.add(e6);
        cp.add(e7);
        cp.add(e8);
        cp.add(e9);

        GenPolynomial e10 = pring.parse("(p4a-p1a)^2+(p4b-p1b)^2+(p4c-p1c)^2-d3^2"); // (10)
        cp.add(e10);

        List> gb;
        GroebnerBase sgb = GBFactory.getImplementation(br);
        gb = sgb.GB(cp);
        //System.out.println("gb = " + gb);

        PolynomialList pl = new PolynomialList(pring, gb);
        Ideal id = new Ideal(pl, true);
        System.out.println("cp = " + cp);
        System.out.println("id = " + id);

        Dimension dim = id.dimension();
        System.out.println("dim = " + dim);
    }

}