Since it seems that CD_ORIGINDEX is not available for loops, the only choice is to simply use the loop normals already computed by depsgraph after evaluating modifiers. This revealed a bug where the Auto Smooth settings would be lost from the mesh after complex modifiers, or after edit mesh to mesh conversion, so restoring them is needed to get correct results.