“…Fourth, if, in addition to p‐box parameters, some parameters are characterized by their precise CDFs, the optimization step is embedded in a Monte Carlo sampling loop ( S0 and S8 ); thereby increasing the number of optimizations by a factor of N (the total number of Monte Carlo samples). To decrease the computational cost, practitioners may opt to use more efficient optimization methods (Deng et al., 2013), fast‐to‐evaluate approximations of the original model or meta‐models (Ellis et al., 2020), parallelization of step S6 by distributing the optimization task across computing units (e.g., central processing unit cores) in combination with using high‐performance computing (Kurgalin & Borzunov, 2019), and a more efficient design‐of‐experiment for steps S3 and S4 (Schöbi & Sudret, 2017). Nevertheless, we expect a higher computational burden since a PBA imposes fewer restrictions (i.e., we do not assume a functional form), leading to a larger region of uncertainty over which a model needs to be evaluated.…”