001/* 002 * $Id: FactorAlgebraicTest.java 3789 2011-10-01 18:54:43Z kredel $ 003 */ 004 005package edu.jas.ufd; 006 007 008import java.util.SortedMap; 009 010import junit.framework.Test; 011import junit.framework.TestCase; 012import junit.framework.TestSuite; 013 014import edu.jas.arith.BigRational; 015import edu.jas.kern.ComputerThreads; 016import edu.jas.poly.AlgebraicNumber; 017import edu.jas.poly.AlgebraicNumberRing; 018import edu.jas.poly.GenPolynomial; 019import edu.jas.poly.GenPolynomialRing; 020import edu.jas.poly.TermOrder; 021 022 023/** 024 * Factor algebraic tests with JUnit. 025 * @author Heinz Kredel. 026 */ 027 028public class FactorAlgebraicTest extends TestCase { 029 030 031 /** 032 * main. 033 */ 034 public static void main(String[] args) { 035 //BasicConfigurator.configure(); 036 junit.textui.TestRunner.run(suite()); 037 } 038 039 040 /** 041 * Constructs a <CODE>FactorAlgebraicTest</CODE> object. 042 * @param name String. 043 */ 044 public FactorAlgebraicTest(String name) { 045 super(name); 046 } 047 048 049 /** 050 */ 051 public static Test suite() { 052 TestSuite suite = new TestSuite(FactorAlgebraicTest.class); 053 return suite; 054 } 055 056 057 int rl = 3; 058 059 060 int kl = 5; 061 062 063 int ll = 5; 064 065 066 int el = 3; 067 068 069 float q = 0.3f; 070 071 072 @Override 073 protected void setUp() { 074 } 075 076 077 @Override 078 protected void tearDown() { 079 ComputerThreads.terminate(); 080 } 081 082 083 /** 084 * Test dummy for Junit. 085 * 086 */ 087 public void testDummy() { 088 } 089 090 091 /** 092 * Test algebraic factorization. 093 * 094 */ 095 public void testAlgebraicFactorization() { 096 097 TermOrder to = new TermOrder(TermOrder.INVLEX); 098 BigRational cfac = new BigRational(1); 099 String[] alpha = new String[] { "alpha" }; 100 String[] vars = new String[] { "z" }; 101 GenPolynomialRing<BigRational> pfac = new GenPolynomialRing<BigRational>(cfac, 1, to, alpha); 102 GenPolynomial<BigRational> agen = pfac.univariate(0, 2); 103 agen = agen.sum(pfac.getONE()); // x^2 + 1 104 AlgebraicNumberRing<BigRational> afac = new AlgebraicNumberRing<BigRational>(agen, true); 105 GenPolynomialRing<AlgebraicNumber<BigRational>> apfac = new GenPolynomialRing<AlgebraicNumber<BigRational>>( 106 afac, 1, to, vars); // univariate 107 108 //System.out.println("agen = " + agen); 109 //System.out.println("afac = " + afac); 110 //System.out.println("apfac = " + apfac); 111 112 FactorAlgebraic<BigRational> fac = new FactorAlgebraic<BigRational>(afac); 113 114 for (int i = 1; i < 2; i++) { 115 int facs = 0; 116 GenPolynomial<AlgebraicNumber<BigRational>> a; 117 GenPolynomial<AlgebraicNumber<BigRational>> c = apfac.random(2, ll + i, el + i, q); 118 //a = a.monic(); 119 GenPolynomial<AlgebraicNumber<BigRational>> b = apfac.random(2, ll + i, el + i, q); 120 if (b.degree() == 0) { 121 b = b.multiply(apfac.univariate(0)); 122 } 123 //b = b.monic(); 124 //if ( false && ! a.leadingBaseCoefficient().isONE() ) { 125 //continue; 126 //ExpVector e = a.leadingExpVector(); 127 //a.doPutToMap(e,cfac.getONE()); 128 //} 129 if (c.degree() > 0) { 130 facs++; 131 } 132 if (b.degree() > 0) { 133 facs++; 134 } 135 //a = apfac.univariate(0,2).sum( apfac.getONE() ); // x^2 + 1 136 //a = a.multiply(a); 137 //a = a.multiply( apfac.univariate(0,2).subtract( apfac.getONE() ) ); // x^2 - 1 138 //a = apfac.univariate(0,3).subtract( apfac.getONE() ); // x^3 - 1 139 //a = apfac.univariate(0,3).sum( apfac.getONE() ); // x^3 + 1 140 a = c.multiply(b); 141 //a = a.monic(); 142 //System.out.println("\na = " + a); 143 //System.out.println("b = " + b.monic()); 144 //System.out.println("c = " + c.monic()); 145 146 SortedMap<GenPolynomial<AlgebraicNumber<BigRational>>, Long> sm = fac.baseFactors(a); 147 //System.out.println("\na = " + a); 148 //System.out.println("sm = " + sm); 149 if (sm.size() >= facs) { 150 assertTrue("#facs < " + facs, sm.size() >= facs); 151 } else { 152 System.out.println("sm.size() < facs = " + facs); 153 } 154 boolean t = fac.isFactorization(a, sm); 155 //System.out.println("t = " + t); 156 assertTrue("prod(factor(a)) = a", t); 157 } 158 } 159 160}