The Riemann distance function $d(p,q)$ is usually defined as the infimum of the lengths of all piecewise smooth paths between $p$ and $q$.
Does it change if we take the infimum only over smooth paths? (Note that if a smooth manifold is connected, then it is smoothly path connected).
I am quite certain the distance does not change. I think that every piecewise smooth path can be approximated by a smooth path.
Around any singular point of the original path, we can take a coordinate ball, and create smomehow a smoothing of a relevant segment of the path which is not much longer than the original.
An explicit construction such as this can be found here. However, the point there is only to show smooth path connectivity, and we also need some bound on the "added length".
Partial Result (Reduction to the case of Euclidean metric):
I show that the specific Riemannian metric does not matter. That is, if we can create a smoothing with small elongation measured by one metric $g_1$ then we can do the same for any other metric $g_2$.
Hence it is enough to prove the claim for $\mathbb{R}^n$ with the standard metric.
Proof:
Since the question is local (we focus around some point $p$ of non-smoothness of the original piecewise-smooth path) we can take an orthonormal frame for $g_1$, denoted by $E_i$. write $g_{ij}=g_2(E_i,E_j)$, I want to find $\text{max} \{g_2(v,v)|v\in \mathbb{S}^{n-1}_{g_1}\} = \text{max} \{g_2(v,v)|v=x^iE_i , x=(x^1,...,x^n) \in \mathbb{S}^{n-1}_{Euclidean}\} = \text{max} \{g_{ij}x^ix^j| \sum(x^i)^2=1 \} = \text{max} \{x^T \cdot G \cdot x | \|x\|=1 \} = \text{max}{\lambda(G)}$.
Since the roots of a polynomial are continuous in in terms of its coefficients, and the coefficients of the charactersitic polynomial of a matrix depends continuously on the matrix entries, it follows that the eigenvalues of a matrix depends continuously on the matrix entries. Hence, since the matrix $g_{ij}(q)$ is a continuous function of $q$, it follows that if we restrict to a compact small enough neigbourhood of $p$ we the function $f(q)= \text{max}{\lambda(g_{ij}(q))}$ is continuous and in particularly bounded by some constant $C$. Hence for any path $\gamma$ which is contained in a small enough neighbourhood of $p$ $L_{g_2}(\gamma) \le \sqrt C L_{g_1}(\gamma)$.
In particular we can take $g_1$ to be the pullback metric of the standrad Euclidean metric via some coordinate ball around $p$. Now solving the problem for the Euclidean case (which implies solving it for $g_1$), we obtain a solution for an arbitrary $g_2$ as required.