Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Write generational output #552

Open
wants to merge 3 commits into
base: master
Choose a base branch
from
Open

Write generational output #552

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
146a219
changes to allow the moead_gen algorithm to write the population at e…
kirbyherm May 31, 2023
340ff79
minor description issues
kirbyherm May 31, 2023
cde7220
appending instead of overwriting
kirbyherm May 31, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion include/pagmo/algorithms/moead_gen.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -93,7 +93,7 @@ class PAGMO_DLL_PUBLIC moead_gen
moead_gen(unsigned gen = 1u, std::string weight_generation = "grid", std::string decomposition = "tchebycheff",
population::size_type neighbours = 20u, double CR = 1.0, double F = 0.5, double eta_m = 20.,
double realb = 0.9, unsigned limit = 2u, bool preserve_diversity = true,
unsigned seed = pagmo::random_device::next());
unsigned seed = pagmo::random_device::next(), std::string outfile = "");

// Algorithm evolve method
population evolve(population) const;
Expand Down Expand Up @@ -213,6 +213,7 @@ class PAGMO_DLL_PUBLIC moead_gen
bool m_preserve_diversity;
mutable detail::random_engine_type m_e;
unsigned m_seed;
std::string m_outfile;
unsigned m_verbosity;
mutable log_type m_log;
boost::optional<bfe> m_bfe;
Expand Down
30 changes: 26 additions & 4 deletions src/algorithms/moead_gen.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -51,14 +51,17 @@ see https://www.gnu.org/licenses/. */
// the other serialization headers.
#include <boost/serialization/optional.hpp>


#include <cstdio>

namespace pagmo
{

moead_gen::moead_gen(unsigned gen, std::string weight_generation, std::string decomposition, population::size_type neighbours,
double CR, double F, double eta_m, double realb, unsigned limit, bool preserve_diversity, unsigned seed)
double CR, double F, double eta_m, double realb, unsigned limit, bool preserve_diversity, unsigned seed, std::string outfile)
: m_gen(gen), m_weight_generation(weight_generation), m_decomposition(decomposition), m_neighbours(neighbours),
m_CR(CR), m_F(F), m_eta_m(eta_m), m_realb(realb), m_limit(limit), m_preserve_diversity(preserve_diversity),
m_e(seed), m_seed(seed), m_verbosity(0u)
m_e(seed), m_seed(seed), m_outfile(outfile), m_verbosity(0u)
{
// Sanity checks
if (m_weight_generation != "random" && m_weight_generation != "grid" && m_weight_generation != "low discrepancy") {
Expand Down Expand Up @@ -175,8 +178,26 @@ population moead_gen::evolve(population pop) const

// Main Generational MOEA/D loop --------------------------------------------------------------------------------------------
for (decltype(m_gen) gen = 1u; gen <= m_gen; ++gen) {
// 0a - custom logging of entire population at each gen
if (m_outfile != ""){
FILE *ostrm;
ostrm = fopen(m_outfile.c_str(), "a");
for (decltype(pop.size()) i = 0u; i < pop.size(); ++i) {
vector_double pop_gen_x = pop.get_x()[i];
for (decltype(pop_gen_x.size()) j = 0u; j < pop_gen_x.size(); ++j) {
fprintf(ostrm, "%f, ", pop_gen_x[j]);
}
vector_double pop_gen_f = pop.get_f()[i];
for (decltype(pop_gen_f.size()) j = 0u; j < pop_gen_f.size(); ++j) {
if(j < pop.get_f()[i].size()-1)
fprintf(ostrm, "%f, ", pop_gen_f[j]);
else
fprintf(ostrm, "%f\n", pop_gen_f[j]);
}
}
fclose(ostrm);
}
// 0 - Logs and prints (verbosity modes > 1: a line is added every m_verbosity generations)

if (m_verbosity > 0u) {
// Every m_verbosity generations print a log line
if (gen % m_verbosity == 1u || m_verbosity == 1u) {
Expand Down Expand Up @@ -383,6 +404,7 @@ std::string moead_gen::get_extra_info() const
stream(ss, "\n\tDistribution index: ", m_eta_m);
stream(ss, "\n\tChance for diversity preservation: ", m_realb);
stream(ss, "\n\tSeed: ", m_seed);
stream(ss, "\n\tOutfile: ", m_outfile);
stream(ss, "\n\tVerbosity: ", m_verbosity);
return ss.str();
}
Expand All @@ -392,7 +414,7 @@ template <typename Archive>
void moead_gen::serialize(Archive &ar, unsigned)
{
detail::archive(ar, m_gen, m_weight_generation, m_decomposition, m_neighbours, m_CR, m_F, m_eta_m, m_realb, m_limit,
m_preserve_diversity, m_e, m_seed, m_verbosity, m_log, m_bfe);
m_preserve_diversity, m_e, m_seed, m_outfile, m_verbosity, m_log, m_bfe);
}

std::vector<population::size_type>
Expand Down