Minimizing with differential evolutionMinimizing a function of many coordinatesMinimizing a function with...

Has a sovereign Communist government ever run, and conceded loss, on a fair election?

Why is gluten-free baking possible?

Source permutation

Does an unused member variable take up memory?

What materials can be used to make a humanoid skin warm?

What is the generally accepted pronunciation of “topoi”?

What is this diamond of every day?

What problems would a superhuman have who's skin is constantly hot?

How can I find out information about a service?

School performs periodic password audits. Is my password compromised?

MySQL importing CSV files really slow

Does a difference of tense count as a difference of meaning in a minimal pair?

What is better: yes / no radio, or simple checkbox?

Why couldn't the separatists legally leave the Republic?

Can I negotiate a patent idea for a raise, under French law?

From an axiomatic set theoric approach why can we take uncountable unions?

Vocabulary for giving just numbers, not a full answer

What ability score modifier does a javelin's damage use?

Plausibility of Mushroom Buildings

Getting the || sign while using Kurier

When Schnorr signatures are part of Bitcoin will it be possible validate each block with only one signature validation?

Are all players supposed to be able to see each others' character sheets?

Expressing logarithmic equations without logs

Do I really need to have a scientific explanation for my premise?



Minimizing with differential evolution


Minimizing a function of many coordinatesMinimizing a function with some restrictionsMinimizing Multiple FunctionsProblem in minimizing expressionSolving 4 coupled differential equation and minimizing the solutionProblem when minimizing user-defined function in Mathematica with Minimize[]Minimizing with constraintsMinimizing functions with parametersMinimizing a conditional function with parametersMinimizing a function problem













4












$begingroup$


A differential evolution algorithm is given here. I would like to get this kind of animation. I thought I could use NMinimize, given
DifferentialEvolution as an option, but it turns out that does not work as I espected.



Is it possible to extract intermediate step in DifferentialEvolution, or do I have to implement algorithm myself?



f[x_, y_] := 

  -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - E^(0.5 (Cos[2 π x] + Cos[2 π y])) + E + 20



p1 = 

  Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

  PerformanceGoal -> "Quality",

  ColorFunction -> "WatermelonColors", 

  Mesh -> None, 

  BoxRatios -> {1, 1, 1}];   



p2 = 

   DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

     ColorFunction -> "WatermelonColors", 

     PlotPoints -> 200, 

     PerformanceGoal -> "Quality", 

     Frame -> False, 

     PlotRangePadding -> None];



p3 = Plot3D[0, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], Mesh -> None];



Show[p1, p3, PlotRange -> {0, 15}]


enter image description here



When I use StepMonitor to track iterations as follows, it does not work.



{fit, intermediates} = 

  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

    MaxIterations -> 1000, 

    Method -> {"DifferentialEvolution", "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

    StepMonitor :> Sow[{x, y}]]];



Table[

  ListPlot[Take[intermediates[[1, i ;; i + 10]]], 

    Frame -> True, ImageSize -> 350, AspectRatio -> 1], 

  {i, 10, 1000, 100}]


EDIT
Here is the result when we used @Michael E2 solution. Cool!!



f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

  E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20



p1 = Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   PerformanceGoal -> "Quality", ColorFunction -> "WatermelonColors", 

   Mesh -> None, BoxRatios -> {1, 1, 1}];



p2 = DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   ColorFunction -> "WatermelonColors", PerformanceGoal -> "Quality", 

   Frame -> False, PlotRangePadding -> None];



p3 = Plot3D[-0.5, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], 

   Mesh -> None];



p4 = Show[p1, p3, PlotRange -> {-0.5, 15}]



Block[{f}, 

  f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

    E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;

  {fit, intermediates} = 

   Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

     MaxIterations -> 30, 

     Method -> {"DifferentialEvolution", 

       "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

     StepMonitor :> 

      Sow[{Optimization`NMinimizeDump`vecs, 

        Optimization`NMinimizeDump`vals}]]];] // Quiet



Multicolumn[

 Table[Show[p4, 

   ListPointPlot3D[{Append[#, 0] & /@ intermediates[[1, i, 1]]}, 

    PlotRange -> {{-5, 5}, {-5, 5}, {-5, 5}}, Boxed -> False, 

    PlotStyle -> Directive[AbsolutePointSize[3], Black]]], {i, 1, 30, 

   2}], 5, Appearance -> "Horizontal"]


enter image description here






mathematical-optimization







share|improve this question











$endgroup$












  • $begingroup$
    Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
    $endgroup$
    – Michael E2
    50 mins ago


















4












$begingroup$


A differential evolution algorithm is given here. I would like to get this kind of animation. I thought I could use NMinimize, given
DifferentialEvolution as an option, but it turns out that does not work as I espected.



Is it possible to extract intermediate step in DifferentialEvolution, or do I have to implement algorithm myself?



f[x_, y_] := 

  -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - E^(0.5 (Cos[2 π x] + Cos[2 π y])) + E + 20



p1 = 

  Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

  PerformanceGoal -> "Quality",

  ColorFunction -> "WatermelonColors", 

  Mesh -> None, 

  BoxRatios -> {1, 1, 1}];   



p2 = 

   DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

     ColorFunction -> "WatermelonColors", 

     PlotPoints -> 200, 

     PerformanceGoal -> "Quality", 

     Frame -> False, 

     PlotRangePadding -> None];



p3 = Plot3D[0, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], Mesh -> None];



Show[p1, p3, PlotRange -> {0, 15}]


enter image description here



When I use StepMonitor to track iterations as follows, it does not work.



{fit, intermediates} = 

  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

    MaxIterations -> 1000, 

    Method -> {"DifferentialEvolution", "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

    StepMonitor :> Sow[{x, y}]]];



Table[

  ListPlot[Take[intermediates[[1, i ;; i + 10]]], 

    Frame -> True, ImageSize -> 350, AspectRatio -> 1], 

  {i, 10, 1000, 100}]


EDIT
Here is the result when we used @Michael E2 solution. Cool!!



f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

  E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20



p1 = Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   PerformanceGoal -> "Quality", ColorFunction -> "WatermelonColors", 

   Mesh -> None, BoxRatios -> {1, 1, 1}];



p2 = DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   ColorFunction -> "WatermelonColors", PerformanceGoal -> "Quality", 

   Frame -> False, PlotRangePadding -> None];



p3 = Plot3D[-0.5, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], 

   Mesh -> None];



p4 = Show[p1, p3, PlotRange -> {-0.5, 15}]



Block[{f}, 

  f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

    E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;

  {fit, intermediates} = 

   Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

     MaxIterations -> 30, 

     Method -> {"DifferentialEvolution", 

       "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

     StepMonitor :> 

      Sow[{Optimization`NMinimizeDump`vecs, 

        Optimization`NMinimizeDump`vals}]]];] // Quiet



Multicolumn[

 Table[Show[p4, 

   ListPointPlot3D[{Append[#, 0] & /@ intermediates[[1, i, 1]]}, 

    PlotRange -> {{-5, 5}, {-5, 5}, {-5, 5}}, Boxed -> False, 

    PlotStyle -> Directive[AbsolutePointSize[3], Black]]], {i, 1, 30, 

   2}], 5, Appearance -> "Horizontal"]


enter image description here






mathematical-optimization







share|improve this question











$endgroup$












  • $begingroup$
    Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
    $endgroup$
    – Michael E2
    50 mins ago
















4












4








4


2



$begingroup$


A differential evolution algorithm is given here. I would like to get this kind of animation. I thought I could use NMinimize, given
DifferentialEvolution as an option, but it turns out that does not work as I espected.



Is it possible to extract intermediate step in DifferentialEvolution, or do I have to implement algorithm myself?



f[x_, y_] := 

  -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - E^(0.5 (Cos[2 π x] + Cos[2 π y])) + E + 20



p1 = 

  Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

  PerformanceGoal -> "Quality",

  ColorFunction -> "WatermelonColors", 

  Mesh -> None, 

  BoxRatios -> {1, 1, 1}];   



p2 = 

   DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

     ColorFunction -> "WatermelonColors", 

     PlotPoints -> 200, 

     PerformanceGoal -> "Quality", 

     Frame -> False, 

     PlotRangePadding -> None];



p3 = Plot3D[0, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], Mesh -> None];



Show[p1, p3, PlotRange -> {0, 15}]


enter image description here



When I use StepMonitor to track iterations as follows, it does not work.



{fit, intermediates} = 

  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

    MaxIterations -> 1000, 

    Method -> {"DifferentialEvolution", "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

    StepMonitor :> Sow[{x, y}]]];



Table[

  ListPlot[Take[intermediates[[1, i ;; i + 10]]], 

    Frame -> True, ImageSize -> 350, AspectRatio -> 1], 

  {i, 10, 1000, 100}]


EDIT
Here is the result when we used @Michael E2 solution. Cool!!



f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

  E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20



p1 = Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   PerformanceGoal -> "Quality", ColorFunction -> "WatermelonColors", 

   Mesh -> None, BoxRatios -> {1, 1, 1}];



p2 = DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   ColorFunction -> "WatermelonColors", PerformanceGoal -> "Quality", 

   Frame -> False, PlotRangePadding -> None];



p3 = Plot3D[-0.5, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], 

   Mesh -> None];



p4 = Show[p1, p3, PlotRange -> {-0.5, 15}]



Block[{f}, 

  f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

    E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;

  {fit, intermediates} = 

   Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

     MaxIterations -> 30, 

     Method -> {"DifferentialEvolution", 

       "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

     StepMonitor :> 

      Sow[{Optimization`NMinimizeDump`vecs, 

        Optimization`NMinimizeDump`vals}]]];] // Quiet



Multicolumn[

 Table[Show[p4, 

   ListPointPlot3D[{Append[#, 0] & /@ intermediates[[1, i, 1]]}, 

    PlotRange -> {{-5, 5}, {-5, 5}, {-5, 5}}, Boxed -> False, 

    PlotStyle -> Directive[AbsolutePointSize[3], Black]]], {i, 1, 30, 

   2}], 5, Appearance -> "Horizontal"]


enter image description here






mathematical-optimization







share|improve this question











$endgroup$




A differential evolution algorithm is given here. I would like to get this kind of animation. I thought I could use NMinimize, given
DifferentialEvolution as an option, but it turns out that does not work as I espected.



Is it possible to extract intermediate step in DifferentialEvolution, or do I have to implement algorithm myself?



f[x_, y_] := 

  -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - E^(0.5 (Cos[2 π x] + Cos[2 π y])) + E + 20



p1 = 

  Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

  PerformanceGoal -> "Quality",

  ColorFunction -> "WatermelonColors", 

  Mesh -> None, 

  BoxRatios -> {1, 1, 1}];   



p2 = 

   DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

     ColorFunction -> "WatermelonColors", 

     PlotPoints -> 200, 

     PerformanceGoal -> "Quality", 

     Frame -> False, 

     PlotRangePadding -> None];



p3 = Plot3D[0, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], Mesh -> None];



Show[p1, p3, PlotRange -> {0, 15}]


enter image description here



When I use StepMonitor to track iterations as follows, it does not work.



{fit, intermediates} = 

  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

    MaxIterations -> 1000, 

    Method -> {"DifferentialEvolution", "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

    StepMonitor :> Sow[{x, y}]]];



Table[

  ListPlot[Take[intermediates[[1, i ;; i + 10]]], 

    Frame -> True, ImageSize -> 350, AspectRatio -> 1], 

  {i, 10, 1000, 100}]


EDIT
Here is the result when we used @Michael E2 solution. Cool!!



f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

  E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20



p1 = Plot3D[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   PerformanceGoal -> "Quality", ColorFunction -> "WatermelonColors", 

   Mesh -> None, BoxRatios -> {1, 1, 1}];



p2 = DensityPlot[f[x, y], {x, -5, 5}, {y, -5, 5}, 

   ColorFunction -> "WatermelonColors", PerformanceGoal -> "Quality", 

   Frame -> False, PlotRangePadding -> None];



p3 = Plot3D[-0.5, {x, -5, 5}, {y, -5, 5}, PlotStyle -> Texture[p2], 

   Mesh -> None];



p4 = Show[p1, p3, PlotRange -> {-0.5, 15}]



Block[{f}, 

  f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

    E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;

  {fit, intermediates} = 

   Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

     MaxIterations -> 30, 

     Method -> {"DifferentialEvolution", 

       "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

     StepMonitor :> 

      Sow[{Optimization`NMinimizeDump`vecs, 

        Optimization`NMinimizeDump`vals}]]];] // Quiet



Multicolumn[

 Table[Show[p4, 

   ListPointPlot3D[{Append[#, 0] & /@ intermediates[[1, i, 1]]}, 

    PlotRange -> {{-5, 5}, {-5, 5}, {-5, 5}}, Boxed -> False, 

    PlotStyle -> Directive[AbsolutePointSize[3], Black]]], {i, 1, 30, 

   2}], 5, Appearance -> "Horizontal"]


enter image description here






mathematical-optimization




mathematical-optimization






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited 1 hour ago







Okkes Dulgerci

















asked 3 hours ago









Okkes DulgerciOkkes Dulgerci

5,2691917




5,2691917












  • $begingroup$
    Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
    $endgroup$
    – Michael E2
    50 mins ago




















  • $begingroup$
    Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
    $endgroup$
    – Michael E2
    50 mins ago


















$begingroup$
Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
$endgroup$
– Michael E2
50 mins ago






$begingroup$
Note that blocking f (Block[{f}, ...]) isn't necessary. It was just to prevent f from being defined, which is a habit I have with single-lettter symbols on SE, esp. ones I use like f, x, etc. -- thanks for the accept!
$endgroup$
– Michael E2
50 mins ago












1 Answer
1






active

oldest

votes


















3












$begingroup$

Here's a way:



Block[{f},

 f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 

   E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;

 {fit, intermediates} = 

  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 

    MaxIterations -> 30, 

    Method -> {"DifferentialEvolution", 

      "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 

    StepMonitor :> 

     Sow[{Optimization`NMinimizeDump`vecs, 

       Optimization`NMinimizeDump`vals}]]];

 ]



Manipulate[

 Graphics[{

   PointSize[Medium],

   Point[intermediates[[1, n, 1]], 

    VertexColors -> 

     ColorData["Rainbow"] /@ 

      Rescale[intermediates[[1, n, 2]], 

       MinMax[intermediates[[1, All, 2]]]]]

   },

  PlotRange -> 5, Frame -> True],

 {n, 1, Length@intermediates[[1]], 1}

 ]


enter image description here



You can find out about things like Optimization`NMinimizeDump`vecs by inspecting the code for Optimization`NMinimizeDump`CoreDE.






share|improve this answer









$endgroup$













    Your Answer





    StackExchange.ifUsing("editor", function () {
    return StackExchange.using("mathjaxEditing", function () {
    StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
    StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
    });
    });
    }, "mathjax-editing");

    StackExchange.ready(function() {
    var channelOptions = {
    tags: "".split(" "),
    id: "387"
    };
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function() {
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled) {
    StackExchange.using("snippets", function() {
    createEditor();
    });
    }
    else {
    createEditor();
    }
    });

    function createEditor() {
    StackExchange.prepareEditor({
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: false,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: null,
    bindNavPrevention: true,
    postfix: "",
    imageUploader: {
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    },
    onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    });


    }
    });














    draft saved

    draft discarded


















    StackExchange.ready(
    function () {
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f193009%2fminimizing-with-differential-evolution%23new-answer', 'question_page');
    }
    );

    Post as a guest















    Required, but never shown

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    3












    $begingroup$

    Here's a way:



    Block[{f},
    
 f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 
    
   E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;
    
 {fit, intermediates} = 
    
  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 
    
    MaxIterations -> 30, 
    
    Method -> {"DifferentialEvolution", 
    
      "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 
    
    StepMonitor :> 
    
     Sow[{Optimization`NMinimizeDump`vecs, 
    
       Optimization`NMinimizeDump`vals}]]];
    
 ]
    

    
Manipulate[
    
 Graphics[{
    
   PointSize[Medium],
    
   Point[intermediates[[1, n, 1]], 
    
    VertexColors -> 
    
     ColorData["Rainbow"] /@ 
    
      Rescale[intermediates[[1, n, 2]], 
    
       MinMax[intermediates[[1, All, 2]]]]]
    
   },
    
  PlotRange -> 5, Frame -> True],
    
 {n, 1, Length@intermediates[[1]], 1}
    
 ]


    enter image description here



    You can find out about things like Optimization`NMinimizeDump`vecs by inspecting the code for Optimization`NMinimizeDump`CoreDE.






    share|improve this answer









    $endgroup$


















      3












      $begingroup$

      Here's a way:



      Block[{f},
      
 f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 
      
   E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;
      
 {fit, intermediates} = 
      
  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 
      
    MaxIterations -> 30, 
      
    Method -> {"DifferentialEvolution", 
      
      "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 
      
    StepMonitor :> 
      
     Sow[{Optimization`NMinimizeDump`vecs, 
      
       Optimization`NMinimizeDump`vals}]]];
      
 ]
      

      
Manipulate[
      
 Graphics[{
      
   PointSize[Medium],
      
   Point[intermediates[[1, n, 1]], 
      
    VertexColors -> 
      
     ColorData["Rainbow"] /@ 
      
      Rescale[intermediates[[1, n, 2]], 
      
       MinMax[intermediates[[1, All, 2]]]]]
      
   },
      
  PlotRange -> 5, Frame -> True],
      
 {n, 1, Length@intermediates[[1]], 1}
      
 ]


      enter image description here



      You can find out about things like Optimization`NMinimizeDump`vecs by inspecting the code for Optimization`NMinimizeDump`CoreDE.






      share|improve this answer









      $endgroup$
















        3












        3








        3





        $begingroup$

        Here's a way:



        Block[{f},
        
 f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 
        
   E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;
        
 {fit, intermediates} = 
        
  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 
        
    MaxIterations -> 30, 
        
    Method -> {"DifferentialEvolution", 
        
      "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 
        
    StepMonitor :> 
        
     Sow[{Optimization`NMinimizeDump`vecs, 
        
       Optimization`NMinimizeDump`vals}]]];
        
 ]
        

        
Manipulate[
        
 Graphics[{
        
   PointSize[Medium],
        
   Point[intermediates[[1, n, 1]], 
        
    VertexColors -> 
        
     ColorData["Rainbow"] /@ 
        
      Rescale[intermediates[[1, n, 2]], 
        
       MinMax[intermediates[[1, All, 2]]]]]
        
   },
        
  PlotRange -> 5, Frame -> True],
        
 {n, 1, Length@intermediates[[1]], 1}
        
 ]


        enter image description here



        You can find out about things like Optimization`NMinimizeDump`vecs by inspecting the code for Optimization`NMinimizeDump`CoreDE.






        share|improve this answer









        $endgroup$



        Here's a way:



        Block[{f},
        
 f[x_, y_] := -20 E^(-0.2 Sqrt[0.5 (x^2 + y^2)]) - 
        
   E^(0.5 (Cos[2 [Pi] x] + Cos[2 [Pi] y])) + E + 20;
        
 {fit, intermediates} = 
        
  Reap[NMinimize[{f[x, y], -5 <= x <= 5, -5 <= y <= 5}, {x, y}, 
        
    MaxIterations -> 30, 
        
    Method -> {"DifferentialEvolution", 
        
      "InitialPoints" -> Tuples[Range[-5, 5], 2]}, 
        
    StepMonitor :> 
        
     Sow[{Optimization`NMinimizeDump`vecs, 
        
       Optimization`NMinimizeDump`vals}]]];
        
 ]
        

        
Manipulate[
        
 Graphics[{
        
   PointSize[Medium],
        
   Point[intermediates[[1, n, 1]], 
        
    VertexColors -> 
        
     ColorData["Rainbow"] /@ 
        
      Rescale[intermediates[[1, n, 2]], 
        
       MinMax[intermediates[[1, All, 2]]]]]
        
   },
        
  PlotRange -> 5, Frame -> True],
        
 {n, 1, Length@intermediates[[1]], 1}
        
 ]


        enter image description here



        You can find out about things like Optimization`NMinimizeDump`vecs by inspecting the code for Optimization`NMinimizeDump`CoreDE.







        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered 2 hours ago









        Michael E2Michael E2

        148k12198478




        148k12198478






























            draft saved

            draft discarded




















































            Thanks for contributing an answer to Mathematica Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid …



            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.


            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f193009%2fminimizing-with-differential-evolution%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            -

            Popular posts from this blog

            Why not use the yoke to control yaw, as well as pitch and roll? Announcing the arrival of...

            Image
            Does the Mueller report show a conspiracy between Russia and the Trump Campaign? Relating to the President and obstruction, were Mueller's conclusions preordained? How much damage would a cupful of neutron star matter do to the Earth? Nose gear failure in single prop aircraft: belly landing or nose landing? Should a wizard buy fine inks every time he want to copy spells into his spellbook? Can you force honesty by using the Speak with Dead and Zone of Truth spells together? A trigonometry question from STEP examination Special flights What is the "studentd" process? systemd and copy (/bin/cp): no such file or directory...
            Read more

            Couldn't open a raw socket. Error: Permission denied (13) (nmap)Is it possible to run networking commands...

            Image
            Why didn't Boeing produce its own regional jet? What is the most common color to indicate the input-field is disabled? Theorists sure want true answers to this! What is the opposite of "eschatology"? What do you call someone who asks many questions? Is there a hemisphere-neutral way of specifying a season? Is this draw by repetition? Bug? String.split chopping off empty trailing value Is it inappropriate for a student to attend their mentor's dissertation defense? How dangerous is XSS Why do I get negative height? Is it possible to map ...
            Read more

            VNC viewer RFB protocol error: bad desktop size 0x0I Cannot Type the Key 'd' (lowercase) in VNC Viewer...

            Image
            risk of flooding in petra in november How to get directions in deep space? How could a planet have erratic days? 15% tax on $7.5k earnings. Is that right? Did the UK lift the requirement for registering SIM cards? Is it allowed to activate the ability of multiple planeswalkers in a single turn? Stack Interview Code methods made from class Node and Smart Pointers Non-trope happy ending? What is the highest possible scrabble score for placing a single tile Open a doc from terminal, but not by its name The IT department bottlenecks progress, how should I handle this? ...
            Read more