In the SVI-JW parametrization, we have
$$ w(k; a, b, \rho, m, \sigma) = a + b \left [ \rho(k-m) + \sqrt{(k-m)^{2} + \sigma^{2}} \right ] $$
Which gives us
$$ \begin{align*} \sigma_{BS}(k) &= \frac{1}{\sqrt{t}}\sqrt{a + b \left [ \rho(k-m) + \sqrt{(k-m)^{2} + \sigma^{2}} \right ]} \\ \\ \\ \frac{\partial \sigma_{BS}}{\partial k} &= \frac{b\left [\rho + \frac{(k - m)}{\sqrt{(k-m)^{2} + v^2}}\right ]}{2\sqrt{t}\sqrt{a + b \left [ \rho(k-m) + \sqrt{(k-m)^{2} + \sigma^{2}} \right ]}} \end{align*} $$
We can evaluate ATM variance $v_{t}$ by setting $k=0$ in $w(k; a, b, \rho, m, \sigma)$, we can evaluate ATM skew $\psi_{t}$ by evaluating $\frac{\partial \sigma_{BS}}{\partial k}|_{k=0}$ and we can evaluate minimum implied variance $\tilde{v}_{t}$ by setting $\frac{\partial \sigma_{BS}}{\partial k} = 0$ and plugging $k$ into $w(k)$.
How can we find the put and call slopes $p_{t}$ and $c_{t}$? I assumed it would be the limit of the variance $\frac{w(k)}{t}$ when $k \rightarrow \pm \infty$ but this gives me $\frac{b}{\sqrt{t}}(\rho \pm 1)$ which does not match Gatheral's results in his original paper.
Link to original paper Arbitrage-free SVI volatility surfaces by Jim Gatheral, Antoine Jacquier here.