Mailing-List: contact user-help@commons.apache.org; run by ezmlm
Precedence: bulk
Reply-To: "Commons Users List" <user@commons.apache.org>
Received-SPF: pass (athena.apache.org: domain of mmakster@gmail.com designates
 209.85.212.47 as permitted sender)
MIME-Version: 1.0
From: Mister Mak <mmakster@gmail.com>
Date: Thu, 23 May 2013 17:17:16 -0400
Message-ID: 
 <CACOOymsXxmM1ahpinWMPug5sy=VKJSv0S8073dCuXo-q5_OPFw@mail.gmail.com>
Subject: Re: [math] Struggling with opimization
To: user@commons.apache.org
Content-Type: multipart/alternative; boundary=089e0158a8883954c104dd6936e4

--089e0158a8883954c104dd6936e4
Content-Type: text/plain; charset=ISO-8859-1

Thanks to Thomas and Gilles for their tips.  The univariate solution
proposed works (for the record:
     UnivariateFunction f = new UnivariateFunction() {public double
value(double x) {return -(x - 5.0) * (x - 5.0) + 2.25; }};
        BrentOptimizer optimizer = new BrentOptimizer(1e-6, 1e-12);
        UnivariatePointValuePair solution = optimizer.optimize(new
UnivariateObjectiveFunction(f), new MaxEval(100), GoalType.MAXIMIZE, new
SearchInterval(-10, 10));
        System.out.println("Min is: " + solution.getValue() + "\tobtained
at: " + solution.getPoint());
).

I am now trying to maximize a function of two variables, and can't get the
following code running:

        MultivariateFunction g =
            new MultivariateFunction() {
                public double value(double[] x) {
                    return -(x[0] - 5.0) * (x[1] - 5.0) + 2.25;
                }
            };
        SimplexOptimizer optimizerMult = new SimplexOptimizer(1e-6, 1e-12);
        PointValuePair solutionMult = optimizerMult.optimize(new
ObjectiveFunction(g), GoalType.MAXIMIZE, new MultiDirectionalSimplex(100));

The last line causes the following NullPointerException:

Exception in thread "main" java.lang.NullPointerException
    at
org.apache.commons.math3.optim.nonlinear.scalar.noderiv.AbstractSimplex.build(AbstractSimplex.java:217)
    at
org.apache.commons.math3.optim.nonlinear.scalar.noderiv.SimplexOptimizer.doOptimize(SimplexOptimizer.java:151)
    at
org.apache.commons.math3.optim.nonlinear.scalar.noderiv.SimplexOptimizer.doOptimize(SimplexOptimizer.java:89)
    at
org.apache.commons.math3.optim.BaseOptimizer.optimize(BaseOptimizer.java:143)
    at
org.apache.commons.math3.optim.BaseMultivariateOptimizer.optimize(BaseMultivariateOptimizer.java:66)
    at
org.apache.commons.math3.optim.nonlinear.scalar.MultivariateOptimizer.optimize(MultivariateOptimizer.java:64)
    at
org.apache.commons.math3.optim.nonlinear.scalar.noderiv.SimplexOptimizer.optimize(SimplexOptimizer.java:122)
    at ca.queensu.montecarlo.Test_MC.<init>(Unknown Source)
    at ca.queensu.montecarlo.Test_MC.main(Unknown Source)

Any suggestions?

Thanks again,

Philippe

--089e0158a8883954c104dd6936e4--