Question

I have a simulation written in C++ in which I need to maintain a variable number of agents, and I am having trouble deciding how to implement it well. Every agent looks something similar to:

class Agent{
public:
    Vector2f pos;
    float health;
    float data[DATASIZE];
    vector<Rule> rules;
}

I need to maintain a variable number of agents in my simulation such that:

1. Preferably, there is no upper bound on the number of agents
2. I can easily add an Agent
3. I can easily remove any agent under some condition (say health<0)
4. I can easily iterate all agents and do something (say health–)
5. Preferably, I can parallelize the work using openMP, because many updates are somewhat costly, but completely independent of other agents.
6. (edit) the order of the agents doesn’t matter at all

What kind of container or design principles should I use for the agents? Until now I was using a vector, but I think it pretty hard to erase from this structure: something I need to do quite often, as things die all the time. Are there any alternatives I should look at? I thought of something like List, but I don’t think they can be parallelized because they are implemented as linked lists with iterator objects?

Thank you

Answer 1

Until now I was using a vector, but I think it pretty hard to erase from this structure: something I need to do quite often, as things die all the time.

How many do you actually expect to die per each step of your simulation? What seems like “all the time” to a human could still be considered very infrequent to a computer. For instance, if each step of your simulation processes thousands of agents but on average only 1 agent dies every few steps, then agent death is a minor incident. With those kind of numbers, your program spends far more time processing live agents than it does dealing with dead agents and so worrying about the performance of removing a dead agent may not be worth while at all. If making agent removal more efficient would end up making normal agent iteration and processing less efficient (yet agent removal is relatively rare), then that would probably be a poor trade-off.

On the other hand, if large numbers of agents are born and die every simulation step, then you might want to make sure those events can be handled efficiently. So it really depends on the kind of numbers you expect to be dealing with.

My general advice would be to proceed with using std::vector (as long as it fits the rest of your design) unless you really expect a significant number of agent deaths per step compared to the number of agents in total.