This subroutine contains a potential hidden memory scaling obstacle, which is the evaluation of the boundary amplitudes in GBTOlib. Even though GBTOlib internally stores just a compact array of orbital boundary amplitudes (dimension Npw × Norb), its interface only allows retrieval of the channel boundary amplitudes (Nchan × Norb). For calculations with large numbers of channels (= target states), this is a significantly larger matrix, which should be distributed. However, this is not possible with the current GBTOlib interface.