Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: while the probability density function ( pdf ) of the sum of two independent random variables is easily described as the convolution of their pdfs , the expression for the pdf of the product is significantly more complicated . given two independent random variables @xmath0 and @xmath1 with pdfs @xmath2 and @xmath3 , where @xmath4 and @xmath5 , the pdf @xmath6 of their product , @xmath7 , can be succinctly written as @xmath8 where @xmath9 is the delta function defined as @xmath10 for a class of test functions @xmath11.when written in a more traditional form , @xmath12 the difficulty of computing @xmath6 becomes apparent since the integral in diverges at @xmath13 as long as e.g. @xmath14 and @xmath15 . indeed , in the simplest case when independent random variables are selected from two normal distributions with zero means , the pdf @xmath6 has a logarithmic singularity at @xmath13 ( see below ) . as far as we know , the only universal method currently available for computing the pdf of the product of two independent random variables relies on a monte - carlo approach , where one samples the individual pdfs , computes the products , and collects enough samples to achieve certain accuracy in the computation of @xmath6 in . however , due to the slow convergence of such method ( typically @xmath16 , where @xmath17 is the number of samples ) achieving high accuracy is not feasible .    in this paper we show how to accurately and efficiently compute the pdf @xmath6 in using an approximate multiresolution analysis ( mra ) , where the scaling function is a gaussian . we demonstrate that , in contrast with the sampling method , our new algorithm is both fast and accurate . moreover , while the sampling method only provides a histogram of @xmath6 , we obtain the result in a functional form that can be used in further computations . a representation of a function in gaussian - based mra ( gmra ) can also be viewed as a multiresolution gaussian mixture . approximations by the gaussian mixtures have been used early on , see e.g.  @xcite . we prove that the integral ( [ eq : integral - to - compute-1 ] ) can be efficiently evaluated as a multiresolution gaussian mixture and provide a tight estimate for the accuracy of the result . our approach to constructing a gmra is based on the observation that , for any finite user - selected accuracy @xmath18 , the function @xmath19 satisfies an approximate two scale relation for an appropriately selected parameter @xmath20 . as a consequence , the scaled and shifted functions @xmath21 form a gaussian - based approximate mra . the approximate nature of the gmra does not limit its applicability since any finite accuracy can be selected . for example , in our computations we roughly match the double precision arithmetic of a typical computer by selecting @xmath18 in the range @xmath22 .    the integral ( [ eq : integral - to - compute-1 ] ) can , in principle , be computed using the mellin transform . it has been observed that , for non - negative random variables , the mellin transform of @xmath6 in ( [ eq : integral - to - compute-1 ] ) is equal to the product of the mellin transforms of @xmath2 and @xmath3 ( equivalently , @xmath6 is the so - called mellin convolution of @xmath2 and @xmath3 ) . in e.g.  @xcite , this approach has been extended to random variables taking both positive and negative values . while analytically appealing , a numerical implementation of mellin convolution is problematic and has not resulted in a reliable numerical method .    instead of using the mellin transform to represent the pdfs of non - negative random variables , we can use linear combinations of decaying exponentials with complex exponents constructed via algorithms developed in @xcite . such representation of pdfs of non - negative random variables leads to a universal method for computing their product . we plan to address such approach in a separate paper . however , for random variables that can take both positive and negative values , it is more natural to use multiresolution gaussian mixtures as described in this paper . moreover , a generalization to higher dimensions becomes problematic if one were to split the treatment of positive and negative values of such random variables . computing the pdf of the product of independent random variables using gmra is just the first of several interesting applications of this approximate basis . the idea of using gaussians to approximate functions appeared earlier in @xcite , where the coefficients of the shifted gaussians were associated with function values . it turns out that much more accurate approximations can be achieved if the coefficients are computed differently , as described later in this paper . since many classical operators of mathematical physics can be accurately approximated via a linear combination of gaussians ( see e.g. @xcite ) , using gmra to apply these operators can be reduced to the analytic evaluation of integrals , and hence , suggests a new class of fast multiresolution algorithms . however , constructing such algorithms is not straightforward since integer translates of gaussians are far from being orthogonal ( after all , they are positive functions so that no cancellations are possible ) . therefore , the use of gmra differs from the usual mra with orthogonal or nearly orthogonal scaling functions and we briefly discuss some of these issues in the appendix . we present our basis construction for functions of one variable and note that the generalization to two and three variables is fairly straightforward and we will describe it in appropriate contexts . extending our approach to bases for multivariate functions in higher dimensions is more complicated and we plan to address it later . finally , we note that the quotient of two independent random variables can be evaluated in the same fashion as the product so that gmra allows us to develop a numerical calculus of pdfs . we start by proving a key approximation result and considering its consequences in section  [ sec : approximation - of - shifted ] . in section [ sec : distribution - of - the ] we prove a relative error estimate for the representation of the pdf of the product via gmra when one of the distributions is gaussian . we also present two approaches to approximating pdfs via a linear combination of gaussians . we then turn to numerical examples in section  [ sec : numerical - examples ] . conclusions and further work is discussed in section  [ sec : conclusions - and - further ] . finally , we present a brief discussion of certain features of gmra in the appendix . we start by proving a key estimate ( theorem  [ thm : projection of a gaussian ] below ) that allows us to replace any shifted gaussian by an explicit linear combination of gaussians ( with a larger , fixed exponent ) shifted at integer values . this estimate yields an approximate two - scale relation for gaussians that makes possible to construct an approximate multiresolution basis with gaussians as the scaling functions . since the basis of shifted gaussians is far from being orthogonal , one would expect difficulties in computing projections of functions on such basis . remarkably , theorem  [ thm : projection of a gaussian ] and theorem  [ thm : the - gaussian - on arbitrary scale ] allow us to compute such projections accurately for any gaussian with an arbitrary exponent and shift . since we are interested in applications where all functions and the results of computing integrals can be approximated via linear combinations of gaussians , the following result is the key to using gmra as a numerical tool . [ thm : projection of a gaussian]if @xmath23 then , for @xmath24 , we have    @xmath25    where @xmath26 and @xmath27 is the jacobi theta function . we have @xmath28 and need to estimate @xmath29 where @xmath30 using poissons summation formula , we have @xmath31 where @xmath32 we obtain @xmath33 yielding the estimate . see also ( * ? ? ? * section 2 , eq.5 ) for this type of analysis . [ cor : app_two_scale_eq]defining @xmath19 and assuming that @xmath34 , we obtain from    @xmath35    where @xmath36 also , setting @xmath37 , @xmath38 and replacing @xmath39 by @xmath40 in , we obtain @xmath41 an approximate two - scale relation for gaussians .    consequently , we can use the function @xmath42 as a scaling function for an approximate multiresolution basis . by selecting an appropriate @xmath43 we can reduce the value of @xmath18 in and achieve arbitrary finite precision in . in the examples presented later in the paper , we use @xmath44 so that @xmath45 ( see figures  [ fig : example0101 - 2 ] and [ fig : example010101 - 2 ] ) . we note that selecting @xmath46 yields @xmath47 achieving double precision accuracy .    if @xmath48 in , we recover the approximation described in @xcite . however , if @xmath49 , our approximation is substantially different and significantly more accurate than that obtained in @xcite . the key difference , of course , is in the approach to the selection of coefficients of the shifted gaussians .       while in ( * ? ? ? * section 2 ) we use the step size of the discretization to control the accuracy of approximating integrals by a sum , here we fix a unit step size and select the exponent of the gaussian @xmath42 in order to achieve the desired accuracy . both , in theorem  [ thm : projection of a gaussian ] and corollary  [ cor : app_two_scale_eq ] , the error is described in terms of the quantity @xmath50 for some appropriate @xmath51 . in the next lemma we show that , as a function of @xmath52 , @xmath50 decreases exponentially fast . thus , slightly increasing the value of @xmath52 provides a much improved error . also , using the lemma , for a given @xmath43 and target accuracy @xmath53 we can accurately estimate the maximum value of @xmath54 to use in corollary  [ cor : app_two_scale_eq ] . [ lem : estimate_theta_at_zero ]    for @xmath55 , @xmath56 in particular , for @xmath57 , we have @xmath58 and @xmath59 where @xmath60{\pi}}{\gamma\left(\frac{3}{4}\right)}-1\right)<2.002.\label{constant_c_pi}\ ] ]    consider @xmath61 since @xmath62 , @xmath2 is a decreasing function which implies the first claim . for , use @xcite to write @xmath63{\pi}}{\gamma\left(\frac{3}{4}\right)}-1\right)<2.002 $ ] . for , note that @xmath64 since for @xmath65 we have @xmath66 , we obtain the result using with @xmath52 replaced by @xmath67 . we consider the family of functions @xmath68 normalized so that the @xmath69-norm @xmath70 does not depend on the scale @xmath71 ( _ n.b_. , setting the norm to @xmath72 would have resulted in bulkier formulas ) . in the next theorem , we show that a gaussian with an arbitrary exponent can be projected ( approximated ) at an appropriate scale of gmra . it is remarkable that the projection coefficients are obtained explicitly without solving any equations . [ thm : the - gaussian - on arbitrary scale]the gaussian @xmath73 with exponent @xmath54 and shift @xmath74 @xmath75 where @xmath76    given @xmath54 , we select the scale @xmath71 so that @xmath77 . once @xmath71 is selected , we rescale the exponent @xmath54 and the shift @xmath74 as @xmath78 and @xmath79 and compute coefficients using for the rescaled gaussian @xmath80 , @xmath81 replacing the variable @xmath39 by @xmath82 in the expression above , we arrive at . given two independent gaussian variables @xmath83 , @xmath84 , with @xmath85 , the distribution of their product variable , @xmath86 is given by    @xmath87    where @xmath88 is modified bessel function of the second kind . the pdf @xmath89 has an integrable logarithmic singularity at zero even though the normal distributions do not have any singularities . in general , the contribution of such singularity to the cumulative distribution function ( cdf ) may or may not be significant but it is always present as it follows from the representation of @xmath6 in . specifically , for @xmath90 in , as @xmath91 approaches @xmath92 we have @xmath93 where @xmath52 is the euler constant . for three or more gaussian pdfs with zero means , their product can be expressed via meijer g - functions . however , when the gaussian pdfs have non - zero means , we are not aware of any explicit expression for the pdf of their product in terms of special functions .    in this section , we compute the pdf of the product of two independent random variables , @xmath0 with pdf @xmath94 , @xmath95 , and the gaussian variable @xmath84 with pdf @xmath96 the pdf of their product , @xmath7 , is expressed via the integral @xmath97 we assume that the pdf @xmath2 is a smooth function except at zero and a finite number of points where it may have integrable singularities . specifically , @xmath2 may have an integrable logarithmic singularity at @xmath98 , @xmath99 , for some integer @xmath100 . we also assume a mild decay of @xmath2 as @xmath101 , @xmath102 where @xmath103 and @xmath104 are positive constants . in particular , @xmath2 can be the cauchy distribution .    in the next theorem we obtain a multiresolution approximation of @xmath105 , where the coefficients at all scales are given explicitly by integrals of well - behaved functions over a unit interval . these integrals are evaluated numerically via adaptive integration , see section  [ sub : adaptive - method - for ] below . [ thm : the - pdf - of the product]the pdf of the product of two independent random variables @xmath0 and the gaussian variable @xmath84 in can be approximated as    @xmath106    where @xmath18 is given in , @xmath107    @xmath108 and @xmath109    splitting the integration over @xmath39 in , we have    @xmath110    and changing @xmath111 in the second integral , obtain @xmath112 since @xmath113 for @xmath114 , we have @xmath115 changing variables @xmath116 , we obtain @xmath117 where @xmath118 . splitting the integration in into integrals over unit intervals , we rewrite @xmath119 and changing the variable of integration , @xmath120 , we arrive at @xmath121d\tau.\label{eq : intermediate formula}\end{aligned}\ ] ] for @xmath122 $ ] we have @xmath123 , so that @xmath124 satisfies @xmath125 . consequently , using and , we obtain @xmath126 where @xmath127 and use these expressions to replace the corresponding gaussians in . setting @xmath128 as in , we now estimate @xmath129 . we have @xmath130d\tau.\nonumber \end{aligned}\ ] ] since @xmath2 is non - negative , we obtain @xmath131d\tau.\label{eq : intermediate formula-2}\end{aligned}\ ] ] using and , we have @xmath132 yielding the estimate in . we compute the integrals    @xmath133 using adaptive integration . the integrand may have only integrable singularities in the interior of the interval @xmath134 $ ] due to our assumption on @xmath2 . since @xmath2 may have a logarithmic singularity at @xmath98 , i.e. @xmath135 near @xmath98 , we have @xmath136 so that , if @xmath71 is large , @xmath2 behaves like a polynomial as a function of @xmath137 . such behavior of the integrand explains our choice of a simple and efficient method for adaptive integration . specifically , we subdivide the integral @xmath134 $ ] adaptively in a hierarchical fashion while using a fixed number of legendre nodes to compute the integral on each subinterval . the subdivision terminates when the integral computed on two subintervals matches that computed on the parent interval within an appropriately chosen accuracy . in our examples we used the gauss - legendre quadrature with 10 nodes and set the accuracy parameter of the adaptive integrator to @xmath138 ; the local maximum depth of subdivision remained small in all cases . we consider two approaches to approximating pdfs via a linear combination of gaussians ( a gaussian mixture ) . the first approach is applicable to distributions that can be expressed via integrals involving gaussians . we illustrate this case using the laplace distribution . the second approach relies on a rapidly decaying fourier transform of a pdf so that it can be approximated by a bandlimited function . we illustrate this case using the gumbel distribution . we also refer to @xcite for a way to approximate a smooth pdf via a linear combination of gaussians provided that the accuracy requirements are not too stringent . although the pdf of the laplace distribution has two parameters , it is sufficient to approximate the function in for the particular case @xmath139 and @xmath140 ,    @xmath141 following ( * ? ? 37 ) , we have @xmath142 a linear combination of gaussians is then obtained by discretizing this integral . it is sufficient to use the trapezoidal rule due to the rapid decay of the integrand as @xmath143 . we limit integration to the interval @xmath144 $ ] and use @xmath145 as the step size . setting @xmath146 , and using @xmath147 terms , we arrive at the gaussian mixture approximation , @xmath148 we use this approximation in example  [ sub : the - laplace - and - gumbel ] to compute the pdf of the product of two independent random variables , one with laplace distribution and the other with gumbel distribution . approximation of the latter via a gaussian mixture is described next . an interpolating gaussian - based scaling function spans the same subspace as the function @xmath42 , since it is constructed as @xmath149 using the interpolating property @xmath150 , @xmath151 , we also have @xmath152 multiplying this identity by @xmath153 , and adding over @xmath154 , we obtain @xmath155 by evaluating the fourier transform , @xmath156 and using , we arrive at @xmath157 by poisson summation formula , @xmath158 we finally obtain @xmath159 to illustrate the use of the interpolating scaling functions in the construction of gaussian mixtures , we consider the gumbel distribution in with parameters @xmath160 and @xmath140 , @xmath161 gaussian mixtures for general parameters are then obtained from this case by rescaling . first we select an interval to approximate @xmath3 . for example , outside the interval @xmath162 $ ] the function @xmath3 is less than @xmath163 . we treat @xmath3 on this interval as approximately periodic and , changing variables , @xmath164 , rescale it to the interval @xmath134 $ ] , @xmath165 sampling @xmath166 at @xmath17 points , where @xmath17 is sufficiently large , we write its approximation using the interpolating scaling function @xmath167 as @xmath168 using the scaling function @xmath169 , we seek the coefficients @xmath170 such that @xmath171 computing the fourier transform , we obtain @xmath172 and @xmath173 sampling at @xmath17 points and using , we obtain @xmath174 and compute the coefficients @xmath170 via the fast fourier transform ( fft ) . in this construction , @xmath17 should have small prime factors ( since we use fft ) but @xmath17 does not have to be a power of @xmath175 . we illustrate the result in figure  [ fig : error - of - approximating gumbel ] , where we set @xmath176 . we use this gaussian mixture as an approximation of the gumbel distribution in example  [ sub : the - laplace - and - gumbel ] .    a more efficient representation ( with fewer terms ) can be obtained by a multiresolution version of the approach presented in this section .     of approximating the rescaled gumbel pdf @xmath166 in as a gaussian mixture with @xmath176 terms . ]      for computing the pdf of the product of independent random variables where the original pdfs are in gmra representation , we repeat the derivation in theorem  [ thm : the - pdf - of the product ] using two representative gaussians from our approximate basis , namely @xmath177 we have @xmath178 and employing the same changes of variables as in the proof of theorem  [ thm : the - pdf - of the product ] , arrive at @xmath179 changing variables @xmath180 , @xmath181 we define @xmath182 changing variable @xmath183 , @xmath184 , and splitting the real axis into unit subintervals ( see proof of theorem  [ thm : the - pdf - of the product ] ) , we obtain @xmath185 changing variables one more time , @xmath120 , we arrive at @xmath186 finally , we use theorem  [ thm : the - gaussian - on arbitrary scale ] to obtain    @xmath187 where @xmath188 since the integrand is comprised of well - behaved functions , we can discretize the integral using the gauss - legendre quadrature with @xmath189 nodes directly ( rather than using an adaptive integrator ) . it follows that @xmath190 where @xmath191 , @xmath192 are the weights of the gauss - legendre quadrature , @xmath193 and @xmath194 we note that matrices @xmath195 and @xmath196 can be precomputed and that they are sparse ( banded ) since small entries can be removed . also , these matrices do not depend on the scales @xmath154 and @xmath197 of the original gaussians @xmath2 and @xmath3 in .    [ ordering remark]the integral describing the pdf of the product does not depend on the order of the factors @xmath2 and @xmath3 . however , the computation of relies on the evaluation of the coefficients in and these coefficients do depend on the order of the factors . while the resulting pdf does not depend on the chosen order , one of the two possible representations may require significantly more coefficients . to explain the reason for such behavior , let us assume that both random variables in the product have normal distributions but different means and deviations , @xmath198 , @xmath199 and @xmath200 , @xmath201 . according to theorem  [ thm : the - pdf - of the product ] , we need to evaluate    @xmath202 as @xmath71 becomes large , we have @xmath203 thus , in order to reduce the number of scales , given the ratios , @xmath204 we pick @xmath2 in as the distribution with the largest ratio . the impact of ordering the factors in this way is illustrated in example  [ sub : example6121 ] . this remark implies that the ordering of the gaussians @xmath2 and @xmath3 in - should be such that @xmath205 . we also note that in this particular case the ratios in do not depend on the scales @xmath154 and @xmath197 . [ moments remark]once a pdf @xmath2 is represented via gmra , @xmath206 where @xmath71 and @xmath207 only take a finite number of integer values , it is immediate to compute its moments since the moments of the basis functions @xmath208 , are readily available . to be precise , for any nonnegative integer @xmath154 , we compute @xmath209 where @xmath210 thus we have @xmath211 , @xmath212 , and higher moments of the basis functions @xmath208 can be evaluated using the recurrence , @xmath213 obtained by rewriting as @xmath214 and integrating by parts . for the moments of @xmath2 we arrive at the recurrence @xmath215 where @xmath216 and @xmath217 in order to compute the variance , we have @xmath218=\mathcal{m}^{\left(2\right)}-\left(\mathcal{m}^{\left(1\right)}\right)^{2},\ ] ] with similar formulas for skewness and kurtosis . unless noted otherwise , in all our numerical examples we set the parameter @xmath44 , so that @xmath219 in . for representing the pdf of the product of random variables , we use a gmra with a very large number of scales ( @xmath220 scales ) to demonstrate that we can accurately resolve the singularity at zero . obviously , there is no need to use so many scales in practical applications ; we use them to illustrate that the accuracy of our approach is well beyond the computational feasibility of any monte - carlo method . the test code was written in fortran 90 and run on a laptop with @xmath221 ghz intel processor . while no speed comparison with a monte - carlo method is meaningful , we note the actual time required for our algorithm . computing the coefficients of the product of two normal distributions and using @xmath220 scales takes @xmath222 seconds and evaluation of the resulting function at @xmath223 points takes @xmath224 second . while the computation of the coefficients @xmath128 in is trivially parallelizable , we did not take advantage of it in our code . in all our examples , we display significant coefficients of the gmra representation , @xmath128 in , as a matrix where the horizontal direction corresponds to the shift @xmath207 and the vertical direction to the scale @xmath71 . in interpreting these images ( see figures  [ fig : example0101 - 1 ] , [ fig : example010101 - 2 ] , [ fig : example2111 - 2 ] , [ fig : example6121 - 3 ] , [ fig : cauchy-2 ] , [ fig : laplace2 ] , [ fig : gumbel2 ] and [ fig : the - laplace - and - the - gumbel - figure-2 ] below ) , it is important to remember that the shift @xmath207 of the basis functions @xmath225 indicates the distance to the singularity at @xmath13 , which is measured in scale dependent units , i.e. @xmath226 at the scale @xmath71 .    just as an extra check , in the examples where explicit answers are not available , we verify the computed distributions by evaluating zero , first and second moments . in all cases these moments were within the accuracy of computation . the moments of pdfs in gmra representation were computed as described in remark  [ moments remark ] . the pdfs of the product of two and three normal random variables with zero means are available analytically and we use them to ascertain accuracy of our algorithm . the exact pdf of the product of two random variables @xmath83 , @xmath84 , with @xmath85 , is given in . setting @xmath227 , we compute the pdf of the product @xmath7 via our algorithm and compare it with the exact pdf . the computed pdf @xmath228 and its gmra coefficients are displayed in figure  [ fig : example0101 - 1 ] and the relative errors @xmath229 for different values of the parameter @xmath43 of the gmra , are displayed in figure  [ fig : example0101 - 2 ] . the exact pdf in has a logarithmic singularity at the origin and we use @xmath220 scales to recover the pdf up to the interval @xmath230 $ ] with accuracy @xmath231 . the limit of resolution of this singularity is illustrated in figure  [ fig : example0101 - 2 ] ( again , well beyond what may be needed in applications ) . pdf @xmath228 of the product in example  [ sub : example0101 ] ( left ) and logarithm ( base 10 ) of gmra coefficients of @xmath228 ( right).,title="fig : " ] pdf @xmath228 of the product in example  [ sub : example0101 ] ( left ) and logarithm ( base 10 ) of gmra coefficients of @xmath228 ( right).,title="fig : " ]     the limit of resolution of singularity of the pdf @xmath228 in example [ sub : example0101 ] using gmra with @xmath220 scales and @xmath232 ( left ) and the relative error curves shown over the interval @xmath233 $ ] for different parameters @xmath43 in the gmra . the bottom curve corresponds to @xmath232 , followed by @xmath234 , @xmath235 and @xmath236 , as we go up in the display ( right).,title="fig : " ]   the limit of resolution of singularity of the pdf @xmath228 in example [ sub : example0101 ] using gmra with @xmath220 scales and @xmath232 ( left ) and the relative error curves shown over the interval @xmath233 $ ] for different parameters @xmath43 in the gmra . the bottom curve corresponds to @xmath232 , followed by @xmath234 , @xmath235 and @xmath236 , as we go up in the display ( right).,title="fig : " ]      we compute the pdf of the product of three zero - mean gaussian random variables , starting with the gmra representation of the product of two of them computed in the previous example . given three zero - mean gaussian random variables , @xmath237 , @xmath238 , and @xmath239 , the pdf of their product , @xmath240 , is available analytically as @xmath241 where @xmath242 is a special case of meijer g - function , see ( * ? ? ? * definition 9.301 ) . the computed product pdf @xmath243 is obtained by multiplying @xmath244 by the gmra representation of the pdf of @xmath245 , as computed in example [ sub : example0101 ] . the pdfs of the random variables @xmath246 , @xmath247 , and @xmath248 are illustrated in figure [ fig : example010101 - 1 ] . as in example [ sub : example0101 ] , we use @xmath220 scales and @xmath232 to recover the pdf @xmath243 with accuracy @xmath231 . in figure  [ fig : example010101 - 2 ] we display the gmra coefficients of @xmath243 and its logarithmic singularity resolved to the interval @xmath249 $ ] . we compute the relative error ,    @xmath250    for different values of the parameter @xmath43 in the gmra and display them in figure  [ fig : example010101 - 2 ] . pdfs of the random variables @xmath245 and @xmath247 in example  [ sub : example010101 ] ( left ) and pdf of the product random variable @xmath251 in example  [ sub : example010101 ] ( right).,title="fig : " ] pdfs of the random variables @xmath245 and @xmath247 in example  [ sub : example010101 ] ( left ) and pdf of the product random variable @xmath251 in example  [ sub : example010101 ] ( right).,title="fig : " ]     logarithm ( base 10 ) of gmra coefficients of the pdf @xmath243 in example  [ sub : example010101 ] ( left ) and the relative error curves shown over the interval @xmath233 $ ] for different parameters @xmath43 in the gmra . the bottom curve corresponds to @xmath232 , followed by @xmath234 , @xmath235 and @xmath236 , as we go up in the display , demonstrating resolution up to the interval @xmath249$].,title="fig : " ] logarithm ( base 10 ) of gmra coefficients of the pdf @xmath243 in example  [ sub : example010101 ] ( left ) and the relative error curves shown over the interval @xmath233 $ ] for different parameters @xmath43 in the gmra . the bottom curve corresponds to @xmath232 , followed by @xmath234 , @xmath235 and @xmath236 , as we go up in the display , demonstrating resolution up to the interval @xmath249$].,title="fig : " ]      in the next five examples we use our algorithm to compute the pdfs of the product of two independent normal random variables with non - zero means as well as the product of independent normal random variable and independent random variables with laplace , cauchy and gumbel distributions . we also compute the product of independent random variables with laplace and gumbel distributions where both distributions are approximated via gaussian mixtures as described in section  [ sub : representing - pdfs - in - gmra ] .      for the random variables @xmath252 , @xmath253 with @xmath254 , @xmath255 and @xmath256 , @xmath257 , we compute the pdf @xmath228 of their product , @xmath7 . figure  [ fig : example2111 - 1 ] shows the pdfs of the random variables @xmath0 , @xmath1 and @xmath246 . as far as we know , there is no analytic expression for the result which we compute with accuracy @xmath258 . using @xmath220 scales , we resolve the logarithmic singularity at zero within an interval as small as @xmath259 $ ] ( see figure  [ fig : example2111 - 2 ] ) . we observe that to compute the pdf of @xmath246 to such accuracy and resolution via a monte - carlo method is not feasible . moreover , we obtain the pdf of @xmath246 in a functional form which can be used for further computations . in contrast , a monte - carlo method using @xmath260 samples and @xmath261 bins of approximate size @xmath262 yields just the histogram shown in figure  [ fig : example2111 - 3 ] . pdfs of the random variables @xmath0 and @xmath1 of example  [ sub : example2111 ] ( left ) and the computed pdf @xmath228 of example  [ sub : example2111 ] ( right).,title="fig : " ] pdfs of the random variables @xmath0 and @xmath1 of example  [ sub : example2111 ] ( left ) and the computed pdf @xmath228 of example  [ sub : example2111 ] ( right).,title="fig : " ]     logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : example2111 ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath263$].,title="fig : " ] logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : example2111 ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath263$].,title="fig : " ]    a histogram of the distribution of example  [ sub : example2111 ] computed with @xmath260 samples using @xmath261 bins of approximate size @xmath262 . the logarithmic singularity at the origin is poorly resolved ( cf . the resolution of our method illustrated in figure  [ fig : example2111 - 2 ] ) . ] we demonstrate the effect of the order of factors on the coefficients of the representation of the product as described in remark  [ ordering remark ] . we compute the pdfs of @xmath264 and @xmath265 , where @xmath266 and @xmath267 , and display the results in figure  [ fig : example6121 - 1 ] and [ fig : example6121 - 3 ] . while the values of the pdfs @xmath268 and @xmath269 are the same within the accuracy of computation , the coefficients of @xmath268 and @xmath269 do depend on the order of the factors as illustrated in figure  [ fig : example6121 - 3 ] . pdfs of random variables @xmath0 and @xmath1 in example  [ sub : example6121 ] ( left ) and computed product pdf @xmath268 ( right).,title="fig : " ] pdfs of random variables @xmath0 and @xmath1 in example  [ sub : example6121 ] ( left ) and computed product pdf @xmath268 ( right).,title="fig : " ]    logarithm ( base 10 ) of gmra coefficients of the pdf @xmath268 ( left ) and logarithm ( base 10 ) of gmra coefficients of the pdf @xmath269 ( right).,title="fig:"]logarithm ( base 10 ) of gmra coefficients of the pdf @xmath268 ( left ) and logarithm ( base 10 ) of gmra coefficients of the pdf @xmath269 ( right).,title="fig : " ]      we now turn to computing the pdf of the product of random variables where one of them has a non - gaussian distribution . first , we consider the pdf of the cauchy distribution    @xmath270 where @xmath271 is the location parameter and @xmath52 is the scale parameter . the cauchy distribution is defined on the real line and is an example of a heavy - tailed distribution for which the moments are infinite . in this example , we compute the pdf @xmath228 of @xmath7 , where @xmath272 and @xmath273 . we used @xmath274 `` negative scales '' , i.e. @xmath275 in order to capture the heavy - tail behavior of the product distribution ( which also has no finite moments ) . by using so many coarse scales , the integral of the resulting pdf over the real line differs from @xmath72 by approximately @xmath276 . pdfs of the factors @xmath0 and @xmath1 and the product @xmath246 are illustrated in figure  [ fig : cauchy-1 ] . in figure [ fig : cauchy-2 ] , we display the gmra coefficients of @xmath228 and resolve the singularity at zero of @xmath228 within the interval @xmath249 $ ] . the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : examplecauchy ] ( left ) and computed product pdf @xmath228 ( right).,title="fig : " ] the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : examplecauchy ] ( left ) and computed product pdf @xmath228 ( right).,title="fig : " ]     logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : examplecauchy ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] ( right ) . , title="fig : " ] logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : examplecauchy ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] ( right ) . , title="fig : " ]      we consider the pdf of the laplace distribution    @xmath277    where @xmath278 is the location parameter and @xmath279 is the scale parameter . the density function , @xmath280 , has a cusp at @xmath281 . in this example , we compute the pdf @xmath228 of the product @xmath7 , where @xmath282 and @xmath283 and display these pdfs in figure [ fig : laplace1 ] . in figure [ fig : laplace2 ] we display the gmra coefficients of @xmath228 and resolve , within the interval @xmath249 $ ] , the singularity of @xmath228 at zero . note that the resulting distribution does not have a singularity at @xmath281 but has a singularity at zero as in the other examples . the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : examplelaplace ] ( left ) and computed product pdf @xmath228 ( right ) . , title="fig : " ] the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : examplelaplace ] ( left ) and computed product pdf @xmath228 ( right ) . , title="fig : " ]     logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : examplelaplace ] and ( b ) logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] . , title="fig : " ] logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : examplelaplace ] and ( b ) logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] . , title="fig : " ]      we also consider the gumbel distribution which is a special case of the generalized extreme value distribution . its pdf is not symmetric and is given by    @xmath284    where @xmath278 is the location parameter and @xmath285 is the scale parameter . here , we compute the pdf @xmath228 of the product random variable @xmath7 , where @xmath286 and @xmath287 and display these pdfs in figure  [ fig : gumbel1 ] . in figure [ fig : gumbel2 ] we display the gmra coefficients of @xmath228 and resolve , within the interval @xmath249 $ ] , the singularity of @xmath228 at zero . pdfs of random variables @xmath0 and @xmath1 in example  [ sub : examplegumbel ] ( left ) and computed product pdf @xmath228 ( right).,title="fig : " ] pdfs of random variables @xmath0 and @xmath1 in example  [ sub : examplegumbel ] ( left ) and computed product pdf @xmath228 ( right).,title="fig : " ]     logarithm ( base 10 ) of gmra coefficients of @xmath228 in example  [ sub : examplegumbel ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] ( right ) . , title="fig : " ] logarithm ( base 10 ) of gmra coefficients of @xmath228 in example  [ sub : examplegumbel ] ( left ) and logarithmic singularity of @xmath228 resolved to the interval @xmath249 $ ] ( right ) . , title="fig : " ]      in this example we compute the pdf @xmath228 of the product random variable @xmath7 , where @xmath288 is the laplace distribution and @xmath289 is the gumbel distribution . both distributions are approximated as gaussian mixtures as described in section  [ sub : representing - pdfs - in - gmra ] . in computing @xmath228 we ordered the parameters of gaussians to reduce the total number of coefficients as explained in remark  [ ordering remark ] . the laplace distribution was approximated by a gaussian mixture with @xmath290 terms and the gumbel distribution with 300 terms , see section  [ sub : representing - pdfs - in - gmra ] for details . we display the results in figure  [ fig : the - laplace - and - the - gumbel - figure-1 ] and figure  [ fig : the - laplace - and - the - gumbel - figure-2 ] . the accuracy of the result is @xmath291 as indicated by computing moments of @xmath228 . the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : the - laplace - and - gumbel](left ) and computed product pdf @xmath228 ( right).,title="fig:"]the pdfs of the random variables @xmath0 and @xmath1 in example  [ sub : the - laplace - and - gumbel](left ) and computed product pdf @xmath228 ( right).,title="fig : " ]    logarithm ( base 10 ) of gmra coefficients of the pdf @xmath228 in example  [ sub : the - laplace - and - gumbel ] . ] gmra allows us to develop a numerical calculus of pdfs since we can now compute sums and products of independent random variables ( as well as quotients , _ _ since the integral for the quotient of two independent random variables can be evaluated in a manner similar to that of the product ) . indeed , if the initial pdfs are expressed in the gmra , the resulting pdf is also expressed in the gmra , i.e. as a multiresolution gaussian mixture . similar results hold for cdfs of random variables . we expect that generalizations of the gmra to two and three dimensions can be accomplished in a straightforward manner , perhaps supplemented by employing the singular value decomposition ( svd ) to simplify the representation of coefficients on each scale . we plan to address such generalizations within relevant applications . a practical generalization to dimensions higher than three is more challenging and we plan to work towards it . computing the pdf of the product of independent random variables does not require all of the features of mra since we only compute projections on the scaling functions at appropriate scales . in this appendix we briefly present some of the usual mra constructs . we assume that the reader is familiar with the basic concepts of orthogonal and bi - orthogonal wavelet bases , see e.g.  @xcite . we note that the approximate gmra is in many respects similar to mra based on high order splines , see e.g.  @xcite . however , there are several special properties of gmra that are not shared by spline based mra . for example , as shown in theorem  [ thm : the - gaussian - on arbitrary scale ] , a gaussian with an arbitrary exponent and an arbitrary shift can be represented at an appropriate subspace which does not depend on the shift . another advantage of using gmra is that it becomes possible to apply operators to functions in a semi - analytic fashion which is important in a number of applications we plan to address separately . finally , we would like to emphasize that the fact that the gmra is approximate is not an obstacle in practice since we can adjust the accuracy of the gmra as needed . @xmath293 more precisely , using the poisson summation formula applied to @xmath294 , we have @xmath295 for @xmath296 , we obtain the error estimate @xmath297 where the last equality can be found in ( * ? ? ? 8.180.3 ) . since @xmath298 lemma  [ lem : estimate_theta_at_zero ] shows that the error decays exponentially fast as a function of the parameter @xmath43 . for example , if @xmath46 , the estimate provides an error bound which is approximately @xmath299 . this leads us to define the filter @xmath300 for gmra as @xmath301 and observe that ( for any @xmath302 ) implies @xmath303 where @xmath304 , see . therefore , for the values of @xmath43 of interest , @xmath305 is extremely close to one . for example , for @xmath46 , yields @xmath306 next , we compute    @xmath307 to define the fourier transform of the approximate , orthogonal gaussian - based scaling function as @xmath308 the corresponding approximate quadrature mirror filter ( qmf ) , determined by the identity @xmath309 , is then given by @xmath310 we note that @xmath311 so that , for @xmath46 , @xmath312 as well . the real valued filter @xmath313 satisfies an approximate qmf equation , @xmath314 for @xmath46 , we numerically verify that @xmath315 in order to obtain the filter that is dual to @xmath300 , we rewrite as @xmath316 so that the dual filter @xmath317 is given by @xmath318 we use @xmath317 to obtain a formula for the projection of coefficients onto a coarser scale . it is sufficient to consider a projection from zero scale into the next coarser scale . therefore , starting with @xmath2 of the form @xmath319 so that @xmath320 the coefficients of @xmath2 in the next coarser scale are computed via ( see e.g. ( * ? ? ? * section 8.3.1 ) ) , @xmath321e^{2\pi inp}dp.\nonumber \end{aligned}\ ] ] we would like to point out a numerical problem associated with using the filter @xmath317 . the plot of @xmath317 for @xmath322 in figure  [ fig : filter  for m00 ] shows that one can expect a loss of significant digits due to its large dynamic range .      to avoid the possible loss of significant digits we have several options . first , as is done in this paper , we can avoid using @xmath317 altogether . the second option is to accept a lower ( e.g. single precision ) accuracy in computations with gmra . finally , one can look for an alternative ( nonlinear ) algorithm for computing coefficients @xmath324 in . to illustrate the second option , we choose @xmath236 and observe the relatively small dynamic range of @xmath317 for this value of @xmath43 in figure  [ fig : filter  for m00 ] . using @xmath236 in computing the pdfs results in roughly single precision accuracy as illustrated in figures  [ fig : example0101 - 2 ] and[fig : example010101 - 2 ] . an alternative approach motivated by the construction of gmra is to develop an exact mra such that the resulting scaling function is closely approximated by a single gaussian . although we do not use this approach in this paper , we briefly indicate how to attain such construction . we can demonstrate that the nonnegative filter      is an excellent approximation of the filter @xmath313 in . in fact , for @xmath322 , we verify numerically that @xmath326.\label{diff_m_a_and_m_0}\ ] ] in the next theorem we show that the filter @xmath327 has the necessary properties to generate an exact orthogonal mra . since @xmath332 satisfies the well - known property @xmath333 and @xmath334 , we obtain @xmath335 by construction , we have @xmath331 and , since @xmath336 is an even function , @xmath327 is even . using and observing that @xmath337 is strictly decreasing for @xmath51 , we have @xmath338 hence , denoting by @xmath339 the constant in , we know that @xmath340 . for @xmath341 , we denote by @xmath2 the function @xmath342 and notice that , due to , we have @xmath343 and @xmath344 . since @xmath345 and @xmath346 is strictly decreasing on @xmath347 , to show that @xmath2 is strictly decreasing in this interval , it is enough to show that @xmath3 has the same property . by , we have @xmath348 and , taking derivatives , @xmath349 is equivalent to the inequality @xmath350 from @xcite , we have , for any @xmath51 , @xmath351 and , hence , @xmath352 which implies and , thus , that @xmath2 is strictly decreasing on @xmath347 . therefore , @xmath2 is also positive in that interval . it follows that @xmath353 is strictly decreasing on @xmath347 . since @xmath327 has no zeros in @xmath354 $ ] , using corollary 6.3.2 in @xcite , we know that the qmf @xmath327 generates an orthogonal mra . the fourier transform of the resulting scaling function @xmath355 can now be obtained as @xmath356 due to , if we use a finite number of scales , the exact mra and approximate gmra generated by the function @xmath357 in are exchangeable for an appropriately selected parameter @xmath43 tuned to the user - selected accuracy . g.  beylkin , g.  fann , r.  j. harrison , c.  kurcz , and l.  monzn . multiresolution representation of operators with boundary conditions on simple domains . , 33:109139 , 2012 .
Abstract:

Answer: we introduce a new approximate multiresolution analysis ( mra ) using a single gaussian as the scaling function , which we call gaussian mra ( gmra ) . we employ this new tool to accurately and efficiently compute the probability density function ( pdf ) of the product of two independent random variables . in contrast with monte - carlo ( mc ) type methods ( the only other universal approach known to address this problem ) , our method not only achieves accuracies beyond the reach of mc but also produces a pdf expressed as a gaussian mixture , thus allowing for further efficient computations . we also show that an exact mra corresponding to our gmra can be constructed for a matching user - selected accuracy .