Currently, in calc_redshift, if subtract_host is True, the function just takes self.dm_excess - self.dm_hostest as the new input DM. However, one must consider the cosmological effect in the host galaxy DM (eg. http://doi.org/10.1086/380598, http://arxiv.org/abs/2107.10858). The new input DM should be DM{excess}-DM_{host}/(1+z). This 1+z factor needs to be carried to the later solving process.
Currently, in calc_redshift, if subtract_host is True, the function just takes self.dm_excess - self.dm_hostest as the new input DM. However, one must consider the cosmological effect in the host galaxy DM (eg. http://doi.org/10.1086/380598, http://arxiv.org/abs/2107.10858). The new input DM should be DM{excess}-DM_{host}/(1+z). This 1+z factor needs to be carried to the later solving process.