Optimization and parallelization of the thermal–hydraulic subchannel code CTF for high-fidelity multi-physics applications