First of all, for dimensionality, you should be able to convince yourself that $\dim \mathfrak{s} = \frac{n(n+1)}{2}$. This can be done easily through counting the number of independent parameters in a symmetric matrix. To prove that $S$ is a manifold, one can either show that $S$ is an open subset of $\mathfrak{s}$ or that $S_n \approx GL_n(\mathbb{R})/O(n)$ using the homogeneous space construction theorem. This also proves that $\dim S = \dim GL_n(\mathbb{R}) - \dim O(n) = n^2 - \frac{n(n-1)}{2} = \frac{n(n+1)}{2}$. Therefore they are of the same dimension.
To further show how $S$ and $\mathfrak{s}$ are related, you should recognize that $\mathfrak{s} = T_IS$ since if $\gamma:(-\epsilon, \epsilon)\to GL_n(\mathbb{R})$ with $\gamma(0) = I$ and $\gamma'(0) = X$, then $\eta(t) = \gamma(t)\gamma^T(t)$ defines a smooth curve in $S$. Then, differentiating we find that
\begin{eqnarray*}
\dot{\eta}(0) & = & \dot{\gamma}(0)\gamma^T(0) + \gamma(0)\dot{\gamma}^T(0) \\
& = & X + X^T
\end{eqnarray*}
which is symmetric, proving that $T_IS \subseteq \mathfrak{s}$. Then by dimensionality equality holds.
What's remarkable is that the exponential map $\exp:\mathfrak{s} \to S$ is actually a global diffeomorphism. We can see this actually just through the eigendecomposition of symmetric matrices. If $A \in \mathfrak{s}$ then $A = PDP^T$ where $D$ is diagonal with arbitrary entries, and $P$ is orthogonal. Then basic properties of the Taylor series of $\exp$ show that
$$
\exp(A) \;\; =\;\; P\exp(D)P^T
$$
where $\exp(D)$ will just have the diagonal entries $e^{d_i}$ which will necessarily be positive and therefore force $\exp(A) \in S$. Similarly, the matrix logarithm, defined by the Taylor series
$$
\log(B) \;\; =\;\; \sum_{k=1}^\infty \frac{(-1)^{k+1}}{k} (I - B)^k
$$
serves as an inverse to $\exp$ since the logarithm carries the same computational property: $\log(B) = \log(QEQ^T) = Q\log(E)Q^T$ with the diagonal entries of $\log(E)$ just being $\log(\lambda_i)$ where $\lambda_i$ are the diagonal entries of $E$.
This answer should hopefully give you some insight into the relationship between symmetric matrices and the symmetric positive definite matrices, however in your prompt you seem to be further curious about the pushforward of the exponential. In other words, while $\exp:\mathfrak{s} \to S$ is a diffeomorphism, we still have the pushforward map
$$
d\exp_X: T_X\mathfrak{s} \to T_{\exp(X)}S.
$$
Computing the purshforward is a much more involved computation, however the fact that $\exp$ is a global diffeomorphism should convince you that $d\exp_X$ is an isomorphism between the tangent spaces of the two manifolds.