Dear Andy,
The relationship file that I referenced is for the 2020 PUMAs (posted on the census website) but they are not in use yet for the current 2021 PUMS data (I use the API). The relationship file in the post has a lot of PUMAs that do not exist in the 2021 PUMS data. So the PUMS data is still using the 2012 (or so) PUMAs. You have to use the GEOCORR rcrosswalk/relationship file to get 2020 tracts in a 2012 PUMA. This is the best that you can do until the 2022 ? PUMS data comes out which I assume will use the 2020 PUMAs I think that there is some "apples to oranges" going on here. The 2017-2021 ACS detail tables use the 2020 tracts and the 2021 PUMs uses the 2012 PUMAs. Very confusing.
I'm using Small Area Estimation (SAE) to take a multiway 2021 PUMS cross tabulation and adjusting it with 2021 ACS tract level detail tables for the marginals.
Dave