Giter VIP home page Giter VIP logo

Comments (6)

bluescarni avatar bluescarni commented on June 12, 2024

@mattapattu the problem is not about private members, rather it seems like your Qfunc class does not satisfy all the requirements of a user-defined problem. The requirements are listed here:

https://esa.github.io/pagmo2/docs/cpp/problem.html

Every UDP must implement at least the following two methods: [...]

from pagmo2.

mattapattu avatar mattapattu commented on June 12, 2024

I have defined fitness and get_bounds methods in Qfunc class:

class Qfunc {
private:
  // Declare the data members that are needed for the function
  Rcpp::List& inputWeightArrayList;
  const Eigen::MatrixXd probMatrix;
  const Rcpp::DataFrame binned_spiketrain;
  const std::vector<Eigen::MatrixXd> log_xi;
  const Eigen::MatrixXd log_gamma;
  const int curr_state;
  const int next_state;
  const int n_observations;
  const int n_states;
public:
  // Constructor to initialize the data members
  Qfunc(Rcpp::List& inputWeightArrayList_,
        const Eigen::MatrixXd probMatrix_,
        const Rcpp::DataFrame& binned_spiketrain_,
        const std::vector<Eigen::MatrixXd> log_xi_,
        const Eigen::MatrixXd log_gamma_,
        const int curr_state_,
        const int next_state_,
        const int n_observations_,
        const int n_states_);
  

  pagmo::vector_double fitness(const pagmo::vector_double &) const;
  
  std::pair<pagmo::vector_double, pagmo::vector_double> get_bounds() const;

  Eigen::VectorXd derivateQfunc(const Eigen::VectorXd x);

  double computeQval(const Eigen::VectorXd x);

  Eigen::MatrixXd getProbMatrix();

  double Evaluate(const arma::mat& x);

  double EvaluateWithGradient(const arma::mat& x, arma::mat& g);
  
};

Two mandatory functions are defined like this (as dummy functions):

vector_double Qfunc::fitness(const vector_double &x) const
{
        return {x};
}

std::pair<vector_double, vector_double> Qfunc::get_bounds() const
{
        return {{-10, -10, -10, -10}, {10, 10, 10, 10}};
}

from pagmo2.

bluescarni avatar bluescarni commented on June 12, 2024

@mattapattu user-defined problems must also satisfy a few other requirements, see the is_udp type trait:

https://esa.github.io/pagmo2/docs/cpp/miscellanea/type_traits.html#_CPPv4I0EN5pagmo6is_udpE

Your Qfunc does not appear to have a default constructor.

from pagmo2.

mattapattu avatar mattapattu commented on June 12, 2024

Thanks, adding a default constructor fixed the issue. But now when I run archipelago.evolve I get a seg fault:

problem prob{Qfunc(inputWeightArrayList, eigen_ProbMatrix, binned_spiketrain, eigen_log_xi, eigen_log_gamma, curr_state, next_state, n_observations, n_states)};
 pagmo::algorithm algo{sade(100)};
 archipelago archi{16u, algo, prob, 20u};
 archi.evolve(10);
 archi.wait_check();

 for (const auto &isl : archi) {
      std::cout << "champion:" <<isl.get_population().champion_f()[0] << '\n';
  }

It gives below error:

 *** caught segfault ***
address 0x30, cause 'memory not mapped'
Error during wrapup: getAttrib: invalid type (double) for TAG
Error: no more error handlers available (recursive errors?); invoking 'abort' restart
Browse[2]> Error: bad value
Warning: stack imbalance in 'lazyLoadDBfetch', 87 then 86
Warning: stack imbalance in '<-', 84 then 83
Warning: stack imbalance in '{', 80 then 79
Warning: stack imbalance in 'if', 78 then 77
Error during wrapup: bad value
Error: no more error handlers available (recursive errors?); invoking 'abort' restart

 *** caught segfault ***
address 0x6, cause 'memory not mapped'
Error during wrapup: long vectors not supported yet: ../../src/include/Rinlinedfuns.h:537
Error: no more error handlers available (recursive errors?); invoking 'abort' restart
Browse[2]>
 *** caught segfault ***
address (nil), cause 'memory not mapped'

I tried running the same code with the builtin "problem prob{schwefel(30)};" and it works fine. However, when I use my own UDP, I get a segmentation fault. I am not sure what is causing this error, as I have used DE from Ensmallen package successfully on the same objective function before. Is it possible to debug where the error is happening ?

from pagmo2.

bluescarni avatar bluescarni commented on June 12, 2024

@mattapattu archipelagos use multithreading, which requires your problem to be thread-safe. My guess would be that your problem implementation is using internally some global state that is being modified concurrently by multiple problem instances.

In order to confirm this, you should try to create an archipelago with a single island (instead of 16) and see if this keeps on crashing.

from pagmo2.

mattapattu avatar mattapattu commented on June 12, 2024

Thanks for the clarification, this could be an issue with R which I using to load input files for optimization. R seems to call gc() automatically and therefore is not recommended to use in multithreads.

from pagmo2.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.