edu.jas.root.RootUtil Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of ext.bundle.math.symja_jas Show documentation
The newest version!
/*
 * $Id: RootUtil.java 4063 2012-07-27 13:05:07Z kredel $
 */

package edu.jas.root;


import java.util.List;

import org.apache.log4j.Logger;

import edu.jas.arith.Rational;
import edu.jas.poly.Complex;
import edu.jas.structure.RingElem;
import edu.jas.structure.RingFactory;


/**
 * Real root utilities. For example real root count.
 * @author Heinz Kredel
 */
public class RootUtil {


    private static final Logger logger = Logger.getLogger(RootUtil.class);


    private static boolean debug = logger.isDebugEnabled();


    /**
     * Count changes in sign.
     * @param  coefficient type.
     * @param L list of coefficients.
     * @return number of sign changes in L.
     */
    public static > long signVar(List L) {
        long v = 0;
        if (L == null || L.isEmpty()) {
            return v;
        }
        C A = L.get(0);
        for (int i = 1; i < L.size(); i++) {
            C B = L.get(i);
            while (B == null || B.signum() == 0) {
                i++;
                if (i >= L.size()) {
                    return v;
                }
                B = L.get(i);
            }
            if (A.signum() * B.signum() < 0) {
                v++;
            }
            A = B;
        }
        return v;
    }


    /**
     * Parse interval for a real root from String.
     * @param s String, syntax: [left, right] or [mid].
     * @return Interval from s.
     */
    public static  & Rational> Interval parseInterval(RingFactory fac, String s) {
        int r = s.length();
        int el = s.indexOf("[");
        if (el >= 0) {
            int ri = s.indexOf("]");
            if (ri > 0) {
                r = ri;
            }
        } else {
            el = -1;
        }
        //System.out.println("s  = " + s);
        String iv = s.substring(el + 1, r).trim();
        //System.out.println("iv = " + iv);
        int k = iv.indexOf(",");
        if (k < 0) {
            k = s.indexOf(" ");
        }
        if (k < 0) {
            C mid = fac.parse(iv);
            return new Interval(mid);
        }
        //System.out.println("k  = " + k + ", len = " + iv.length());
        String ls = iv.substring(0, k).trim();
        String rs = iv.substring(k + 1, iv.length()).trim();
        //System.out.println("ls = " + ls + ", rs = " + rs);
        C left = fac.parse(ls);
        C right = fac.parse(rs);
        if (debug) {
            logger.debug("Interval: left = " + left + ", right = " + right);
        }
        return new Interval(left, right);
    }


    /**
     * Parse rectangle for a complex root from String.
     * @param s String, syntax: [south-west, north-east] or [mid].
     * @return Interval from s.
     */
    public static  & Rational> Rectangle parseRectangle(RingFactory> fac,
                    String s) {
        int r = s.length();
        int el = s.indexOf("[");
        if (el >= 0) {
            int ri = s.indexOf("]");
            if (ri > 0) {
                r = ri;
            }
        } else {
            el = -1;
        }
        //System.out.println("s  = " + s);
        String iv = s.substring(el + 1, r).trim();
        //System.out.println("iv = " + iv);
        int k = iv.indexOf(",");
        if (k < 0) {
            k = s.indexOf(" ");
        }
        if (k < 0) {
            Complex mid = fac.parse(iv);
            return new Rectangle(mid);
        }
        //System.out.println("k  = " + k + ", len = " + iv.length());
        String ls = iv.substring(0, k).trim();
        String rs = iv.substring(k + 1, iv.length()).trim();
        //System.out.println("ls = " + ls + ", rs = " + rs);
        Complex sw = fac.parse(ls);
        Complex ne = fac.parse(rs);
        if (debug) {
            logger.debug("Rectangle: sw = " + sw + ", ne = " + ne);
        }
        return new Rectangle(sw, ne);
    }

}