I can make multiple instances to run multiple simulations in parallel, but matrix factorisation within the program itself is single-threaded. The stiffness matrix of these structures needs to be reformulated frequently (every time a structural element yields / gets damaged), so this is the largest time-sink, I think.