Overview of our global algorithm to infer agnostic and self-consistent planetary interior models. The difference between $\tilde{\rho }$ used by the compositional algorithm and ρ compatible with the interior pressure and planetary mass decreases until it vanishes using our iterative process. Traditionally, empirical models only use a single iteration, leading to inconsistencies. Figure B 1 contains more details.