How can I access fmincon's function count within a custom anonymous function?

7 次查看（过去 30 天）
James
James 2017-4-20
编辑： Matt J 2017-4-22
I am trying to use fmincon to optimise an aerofoil shape with CFD (computational fluid dynamics). I have defined my own anonymous function to generate the aerofoil shape and run the CFD case on a remote cluster. I want to access fmincon's iteration number and function count within my function so I can name the input files to track the algorithm's progress. This is particularly important when I run fmincon in parallel as any nasty workarounds (global variables or writing files to disc) fail horribly and you end up with files with the same name.
I see that the build in plotting functions like optimplotfval must know the iteration number but I can't work out how to use them in my case.
If you can't access the function count directly, manually indexing the calls to the function in some other robust way would be fine.
I have included an example below, please can anyone help me?
function [x,fval] = optimise_aerofoil()
% Calculate input parameters (x) that yield the best aerofoil design
% measured by the output (fval)
% Define some constants that are needed to run the case
a = 1; b = 2; c = 3;
% Initialise the first design and bounds
x0 = [1 2 3 4 5]; lb = [0 0 0 0 0]; ub = [10 10 10 10 10];
% Define the anonymous function handle to run the cfd and return
% the output value
f = @(x)run_cfd(x,a,b,c);
% Choose optimiser settings
options = optimoptions('fmincon','Display','iter','UseParallel',1);
% Run the optimiser
[x,fval] = fmincon(f,x0,[],[],[],[],lb,ub,[],options);
end
function f = run_cfd(x,a,b,c)
% Generate a design, save it, run the CFD and evaluate the design
% Generate the design from x
xyz = generate_design(x,a,b,c);
%%%%%%%HERE IS WHERE I WOULD LIKE TO KNOW F-COUNT %%%%%%%
iter_count = ?;
% Save the input file to CFD
filename = ['Run_' num2str(iter_count) '.dat'];
dlmwrite(filename,xyz);
% Run the CFD
run_cfd(filename)
% Evaluate the design
f = evaluate_design(filename);
end
  2 个评论
Matt J
Matt J 2017-4-20
I don't really see how the code you've posted can run without error messages. The run_cfd() function is multiply defined - both as your objective function to fmincon and as some external routine
run_cfd(filename)
that takes filename as input.
James
James 2017-4-21
Sorry, I made a mistake in calling both functions run_cfd in my example. I should have called the second one something else.
...
dlmwrite(filename,xyz);
% Run the CFD on a remote machine
run_remote_cfd_case(filename)
% Evaluate the design
...

请先登录，再进行评论。

请先登录，再回答此问题。

回答（3 个）

Steven Lord
Steven Lord 2017-4-20
Use an output function. According to the documentation that lists the fields passed into the OutputFcn, fmincon does pass the field "iteration" into the OutputFcn in the optimValues struct array.
  5 个评论
显示 3更早的评论
Steven Lord
Steven Lord 2017-4-21
I'm not sure if this will work in your parallel use case, but maybe the memoize function introduced in release R2017a can help you in reusing previously computed results.
Matt J
Matt J 2017-4-21
编辑：Matt J 2017-4-22
@James,
It appears that you have two purposes here. One purpose is to generate unique filenames to avoid conflicts between parallel calls to run_remote_cfd_case(). I believe this is already covered by suggestions made by Walter and me.
The second purpose is to save the results of each iteration so that you can resume the optimization if it aborts prematurely. The latter is best done using an Output Function as Steve suggests, here. There is an example of how to do this in the documentation on Output Functions, but basically, it would look like this for you
% Choose optimiser settings
options = optimoptions('fmincon','Display','iter',...
'UseParallel',1,'OutputFcn',@outfun);
function stop = outfun(x,optimValues,state)
if strcmp(state,'iter')
filename = ['Iteration_' num2str(optimValues.iteration)];
save(filename,'x');
end
end

请先登录，再进行评论。

Walter Roberson
Walter Roberson 2017-4-20
Perhaps you could use tempname() to generate unique names? Though I am not certain that parallel workers have JVM
With R2017a onward, you might be able to use the new send() / poll() facilities.
  2 个评论
James
James 2017-4-21
Thank you, this is a nice idea.
Ideally I'd like to have names which are repeatable though. Running the CFD is quite expensive and sometimes tempermental. If I restart fmincon I'd like to to recognise that it has already converged some designs and re-use the output.
Matt J
Matt J 2017-4-21
编辑：Matt J 2017-4-21
If I restart fmincon I'd like to to recognise that it has already converged some designs and re-use the output.
For that purpose, it is a bad idea to save all evaluations of your objective function. Some evaluations result in no progress, or worse. You should be saving the final result of each iteration, using the output function method suggested by Steve, and as demonstrated in this example.

请先登录，再进行评论。

Matt J
Matt J 2017-4-21
编辑：Matt J 2017-4-21
It might make the most sense to generate a temporary filename based on labindex()
filename = ['Run_' num2str(labindex) '.dat'];
The iteration count and funccount will not let you distinguish between parallel workers, as labindex does.
  2 个评论
Walter Roberson
Walter Roberson 2017-4-21
"However, inside a parfor-loop, labindex always returns a value of 1 on all workers in all iterations."
Matt J
Matt J 2017-4-21
编辑：Matt J 2017-4-21
OK, but this will work in parfor
t=getCurrentTask(); t.ID;
filename = ['Run_' num2str(t.ID) '.dat'];
For some reason, though, Edric Ellis does prefer tempname()

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Performance and Memory 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by