The non-negativity condition of a discretized mean-reverting Heston model with stochastic violatilities

Question

I happened to encounter the following discretized mean-reverting Heston model with stochastic volatilities in a paper $$ P(t) = P(t-1) + v_1(u_1-P(t-1))+\sqrt{\sigma(t)}\cdot \epsilon_1(t) \\ \sigma(t) = \sigma(t-1) + v_2(u_2-\sigma(t-1))+\sqrt{\sigma(t-1)}\cdot \epsilon_2(t) $$ where $v_1=v_2=0.1,u_1=100,u_2=0.01$ are pre-set parameters, and $\epsilon_1,\epsilon_2 \sim N(0,1)$ follow the normal distribution IID. Recall that in the original Heston Model formulation, there is a condition (known as the Feller condition) to make sure that the values under the square root is positive. See wiki for more info. But in this case, how can I ensure that the value of $\sigma$ to be positive?