Examples with WrappedException org.matheclipse.core.eval.exception.WrappedException used on opensource projects

Example 1 with WrappedException

use of org.matheclipse.core.eval.exception.WrappedException in project symja_android_library by axkr.

the class NMinimize method simplexSolver.

protected static IExpr simplexSolver(VariablesSet vars, LinearObjectiveFunction f, OptimizationData... optData) {
    try {
        SimplexSolver solver = new SimplexSolver();
        PointValuePair solution = solver.optimize(optData);
        double[] values = solution.getPointRef();
        IAST result = F.List(F.num(f.value(values)), F.List(values));
        return result;
    } catch (MathIllegalStateException oe) {
        throw new WrappedException(oe);
    // if (Config.SHOW_STACKTRACE) {
    // oe.printStackTrace();
    // }
    }
}

Example 2 with WrappedException

use of org.matheclipse.core.eval.exception.WrappedException in project symja_android_library by axkr.

the class Roots method findRoots.

/**
	 * <p>
	 * Given a set of polynomial coefficients, compute the roots of the
	 * polynomial. Depending on the polynomial being considered the roots may
	 * contain complex number. When complex numbers are present they will come
	 * in pairs of complex conjugates.
	 * </p>
	 * 
	 * @param coefficients
	 *            coefficients of the polynomial.
	 * @return the roots of the polynomial
	 */
@Nonnull
public static IAST findRoots(double... coefficients) {
    int N = coefficients.length - 1;
    // Construct the companion matrix
    RealMatrix c = new Array2DRowRealMatrix(N, N);
    double a = coefficients[N];
    for (int i = 0; i < N; i++) {
        c.setEntry(i, N - 1, -coefficients[i] / a);
    }
    for (int i = 1; i < N; i++) {
        c.setEntry(i, i - 1, 1);
    }
    try {
        IAST roots = F.List();
        EigenDecomposition ed = new EigenDecomposition(c);
        double[] realValues = ed.getRealEigenvalues();
        double[] imagValues = ed.getImagEigenvalues();
        for (int i = 0; i < N; i++) {
            roots.append(F.chopExpr(F.complexNum(realValues[i], imagValues[i]), Config.DEFAULT_ROOTS_CHOP_DELTA));
        }
        return roots;
    } catch (Exception ime) {
        throw new WrappedException(ime);
    }
}