Examples with DblDbl_Dbl primal.primitive.DblDbl_Dbl used on opensource projects

Example 1 with DblDbl_Dbl

use of primal.primitive.DblDbl_Dbl in project suite by stupidsing.

the class Bfgs method lineSearch.

private double lineSearch(Dbl_Dbl phi, Dbl_Dbl phiGradient, double alphax) {
    var c1 = .0001d;
    var c2 = .1d;
    var alpha0 = 0d;
    var v0 = phi.apply(alpha0);
    var g0 = phiGradient.apply(alpha0);
    // TODO
    DblDbl_Dbl interpolate = (a0, a1) -> (a0 + a1) * .5d;
    // TODO
    DblDbl_Dbl choose = (a0, a1) -> (a0 + a1) * .5d;
    DblDbl_Dbl zoom = (alphaLo, alphaHi) -> {
        for (var iter = 0; iter < 16; iter++) {
            var alpha = interpolate.apply(alphaLo, alphaHi);
            var v = phi.apply(alpha);
            double g;
            if (v0 + c1 * alpha * g0 < v || phi.apply(alphaLo) <= v)
                alphaHi = alpha;
            else if (abs(g = phiGradient.apply(alpha)) <= -c2 * g0)
                return alpha;
            else {
                if (0d <= g * (alphaHi - alphaLo))
                    alphaHi = alphaLo;
                alphaLo = alpha;
            }
        }
        return alphaLo;
    };
    var alphap = alpha0;
    var vp = v0;
    var alpha = choose.apply(alphap, alphax);
    for (var iter = 0; iter < 16; iter++) {
        var v = phi.apply(alpha);
        if (v0 + c1 * alpha * g0 < v || 0 < iter && vp <= v)
            return zoom.apply(alphap, alpha);
        var g = phiGradient.apply(alpha);
        if (abs(g) <= -c2 * g0)
            break;
        else if (0d <= g)
            return zoom.apply(alpha, alphap);
        else {
            alphap = alpha;
            vp = v;
            alpha = choose.apply(alpha, alphax);
        }
    }
    return alpha;
}

Example 2 with DblDbl_Dbl

use of primal.primitive.DblDbl_Dbl in project suite by stupidsing.

the class PlotMain method run.

private boolean run() {
    var size = 1024;
    var scale = 1d / (size + 1);
    DblDbl_Dbl variety = (x, y) -> {
        return y * y - (x + .25f) * (x + .15f) * (x + .05f) * (x - .05f) * (x - .15f) * (x - .25f);
    };
    IntInt_Int fp = (fx, fy) -> {
        var x0 = fx * scale - .5d;
        var y0 = fy * scale - .5d;
        var value = variety.apply(x0, y0);
        if (Double.isNaN(value))
            return -2;
        else if (value < 0)
            return -1;
        else
            return 1;
    };
    return // 
    new Render().renderPixels(size, size, (fx, fy) -> {
        var b0 = fp.apply(fx, fy);
        var b1 = fp.apply(fx + 1, fy);
        var b2 = fp.apply(fx, fy + 1);
        var b3 = fp.apply(fx + 1, fy + 1);
        var c = b0 != b1 || b1 != b2 || b2 != b3 ? 1d : 0d;
        return new R3(c, c, c);
    }).view();
}