Does the parallel version of the code do load-balancing?
No, the code does not have parallel load-balancing capabilities at the moment. It can do static load-balancing i.e. “optimal” partitioning of the domain so that an initial mesh is divided in roughly equal-sized subdomains while minimising the size of the communication boundaries. I don’t personally use the parallel capability of the code very much. What I usually need is several different sequential computations with a different set of parameters. This is of course the ideal case for “parallelism”. I don’t usually require to run “one shot” very large parallel computations.