Let $X$ and $Y$ be independent non-negative random variables, with density functions $f_X(x)$ and $f_Y(y)$. Let $Z=XY$.
$$P(Z\le z)=\iint_D f_X(x)f_Y(y)\,dx\,dy.$$
Here $D$ is the region in the first quadrant which is "below" the hyperbola $xy=z$.
Evaluate the integral, and differentiate the result with respect to $z$ to get the density function of $Z$. We can usually arrange to do the differentiation under the integral sign, but that still leaves one integral that may, like most integrals, not be expressible in terms of standard functions.
An alternative approach is to find the density functions of the random variables $\ln X$ and $\ln Y$, by using standard methods. Then use the usual "convolution" formula to find the density function of $U$, where $U=\ln X +\ln Y$. Finally, find the density function of $\exp(U)$.
Added: To my surprise, an integral not far from what is necessary for the first approach can be expressed in terms of modified Bessel functions and modified Struve functions (whatever those are). So says Wolfram Alpha. Thus there may be a sort of closed form for your density function. There would be, anyway, if what is called $v$ in the problem was a plain exponential. I would suggest also trying the second approach.