It utilises the MDAnalysis function radius_of_gyration(pbc=True), but allows users to specify the atom positions and their corresponding resids on which to perform this calculation at each time step. Figure shows the comparison between using the MDAnalysis radius of gyration function alone compared to the PySoftK rgyr, which captures the right radius of gyration of the micelle when computed on the whole coordinates from make_micelle_whole.