In the paper it emphasizes:

"When the rate k is adjusted, the offset parameter m must also be adjusted to connect the endpoints of the segments. The correct adjustment at changepoint j is easily computed as:"

$$\gamma_j = \left(s_j - m - \sum_{l

My understanding is, adjustment makes sure function is still continuous.

I have 2 questions regarding this model that are not explained in the paper:

How do we come up with the result function is not continuous without the adjustment factor, and how do we compute the adjustment factor?

The answer has been accepted, but I just want to point out a more straightforward way to arrive there. The key to make sense of the formula is to recognize that both final growth rates and offsets are cumulative sums of respective adjustments, plus the initial rate and offset ($$k$$ and $$m$$). It's instantly clear from the code.
$$g(t) = \frac{C}{1+exp(-k(t-b))}$$
Denote $$k_i$$, $$b_i$$ and $$s_i$$ values of growth rates, offsets and starting times for segment $$i$$. We need this relation to hold:
$$k_i(s_i - b_i) = k_{i+1}(s_i - b_{i+1})$$ where $$s_i$$ is a changepoint.
Assuming we already know offsets up to the current point, we need to find next offset $$\gamma_{i+1}$$ such that
$$k_i(s_i - b_i) = (k_i + \delta_{i+1})(s_i - b_i - \gamma_{i+1})$$ Left hand side expression gets cancelled out so we're left with $$\delta_{i+1}(s_i - b_i) = \gamma_{i+1}(k_i+\delta_{i+1})$$ $$\gamma_{i+1} = \frac{\delta_{i+1}(s_i - b_i)}{k_i+\delta_{i+1}}=\frac{(k_{i+1} - k_i)(t_i - b_i)}{k_{i+1}} = (1 - \frac{k_i}{k_{i+1}})(s_i - b_i)$$ which becomes the target expression after you expand rates and offsets as cumulative sums of adjustments.