Fast Computation of Singular Oscillatory Fourier Transforms

<table>Illustration of the integration paths for the integrand <svg height="21.1875" id="M47" style="vertical-align:-2.34499pt" version="1.1" viewbox="0 0 166.53751 21.1875" width="166.53751" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.7)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,5.944,16.7)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,19.22,16.7)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,32.972,16.7)"><path d="M483 97q-42 -50 -88.5 -79.5t-68.5 -29.5q-37 0 -17 93l22 102h-2q-54 -79 -144 -149q-59 -46 -100 -46q-24 0 -43 29t-19 86q0 78 34.5 153t94.5 120q41 31 94 51.5t98 20.5q29 0 72 -9q26 -6 39 -6l2 -4q-30 -117 -67 -323q-8 -41 2 -41q16 0 79 58zM374 387
q-32 15 -73 15q-52 0 -83 -23q-48 -36 -78 -108.5t-30 -152.5q0 -33 8.5 -50.5t20.5 -17.5q31 0 107 79t99 132q15 40 29 126z" id="x1D44E"></path></g><g transform="matrix(.017,-0,0,-.017,41.573,16.7)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
<g transform="matrix(.012,-0,0,-.012,47.462,8.537)"><path d="M545 106q-67 -118 -134 -118q-24 0 -40 37.5t-30 129.5h-2q-47 -72 -103 -119.5t-108 -47.5q-47 0 -76 45.5t-29 119.5q0 113 85 204t174 91q47 0 70 -33.5t43 -119.5h3q32 47 80 140l55 13l10 -9q-47 -80 -138 -201q17 -99 27.5 -136t22.5 -37q23 0 69 61zM333 204
q-14 98 -31 149.5t-50 51.5q-49 0 -94 -70t-45 -164q0 -55 15.5 -86t40.5 -31q70 0 164 150z" id="x1D6FC"></path></g>
<g transform="matrix(.017,-0,0,-.017,54.838,16.7)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,60.719,16.7)"><path d="M452 333q0 -82 -46 -164t-116 -128q-81 -53 -158 -53q-26 0 -51 13t-37 32q-26 41 -6 134l91 431q7 40 2.5 47t-34.5 7h-33l2 25q36 4 72.5 13t58.5 15.5t28 6.5q10 0 4 -29l-91 -391h2q65 77 130 116.5t105 39.5q36 0 56.5 -32t20.5 -83zM365 316q0 76 -35 76
q-31 0 -93.5 -49t-107.5 -110q-12 -37 -17 -70q-9 -65 12 -96.5t59 -31.5q31 0 56 14q53 29 89.5 105t36.5 162z" id="x1D44F"></path></g><g transform="matrix(.017,-0,0,-.017,72.465,16.7)"><path d="M535 230h-483v50h483v-50z" id="x2212"></path></g><g transform="matrix(.017,-0,0,-.017,86.217,16.7)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,95.719,16.7)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g>
<g transform="matrix(.012,-0,0,-.012,101.612,8.537)"><path d="M558 587q0 -32 -14 -61t-40 -53.5t-48.5 -41t-54.5 -36.5q144 -51 144 -174q0 -55 -43.5 -108t-104.5 -87q-77 -42 -131 -42q-31 0 -54 20t-31 47l11 18q48 -29 108 -29q79 0 119.5 43t40.5 109t-44.5 107.5t-119.5 50.5l22 47q34 1 65 21q96 61 96 157q0 42 -24 67.5
t-62 25.5q-24 0 -43.5 -9t-35 -29.5t-27 -44t-22.5 -63t-19.5 -75.5t-18.5 -91q-57 -294 -68 -380q-26 -190 -35 -200q-26 -31 -97 -37l-4 26q19 9 48 170l77 413q23 121 52.5 187.5t83.5 114.5q70 62 148 62q51 0 88.5 -34t37.5 -91z" id="x1D6FD"></path></g>
<g transform="matrix(.017,-0,0,-.017,109.15,16.7)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g><g transform="matrix(.017,-0,0,-.017,120.063,16.7)"><path d="M300 -147l-18 -23q-106 71 -159 185.5t-53 254.5v1q0 139 53 252.5t159 186.5l18 -24q-74 -62 -115.5 -173.5t-41.5 -242.5q0 -130 41.5 -242.5t115.5 -174.5z" id="x28"></path></g><g transform="matrix(.017,-0,0,-.017,125.945,16.7)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g><g transform="matrix(.017,-0,0,-.017,135.447,16.7)"><path d="M275 270q0 -296 -211 -440l-19 23q75 62 116.5 174t41.5 243t-42 243t-116 173l19 24q211 -144 211 -440z" id="x29"></path></g><g transform="matrix(.017,-0,0,-.017,141.329,16.7)"><path d="M391 364q0 -28 -20.5 -53.5t-50 -43.5t-70 -34t-73.5 -25t-68 -17v-29q0 -48 20 -81t65 -33q73 0 157 76l16 -23q-107 -113 -221 -113q-21 0 -40.5 6.5t-39 22t-31.5 47t-12 75.5q0 69 33.5 139.5t95.5 114.5q77 55 143 55q42 0 69 -24.5t27 -59.5zM313 350
q0 26 -14.5 40.5t-37.5 14.5q-32 0 -44 -7q-82 -46 -104 -175q200 49 200 127z" id="x1D452"></path></g>
<g transform="matrix(.012,-0,0,-.012,148.362,8.537)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g><g transform="matrix(.012,-0,0,-.012,151.611,8.537)"><path d="M615 290q0 -148 -114 -245q-65 -57 -125 -57q-39 0 -62 27t-29 70q-31 -41 -72 -69t-74 -28q-49 0 -82.5 43t-33.5 124q0 97 55 174.5t148 124.5l15 -24q-59 -45 -94 -113t-35 -146q0 -58 19.5 -88.5t48.5 -30.5q33 0 62 34.5t45 101.5q16 73 20 135l80 17l5 -6
q-11 -37 -23.5 -100t-12.5 -92q0 -47 16.5 -71.5t44.5 -24.5q62 0 95.5 67t33.5 160q0 107 -63 152l1 9q14 14 45 14q29 0 57.5 -44t28.5 -114z" id="x1D714"></path></g><g transform="matrix(.012,-0,0,-.012,159.202,8.537)"><path d="M536 404q0 -17 -13.5 -31.5t-26.5 -14.5q-8 0 -15 10q-11 14 -25 14q-22 0 -67 -50q-47 -52 -68 -82l37 -102q31 -88 55 -88t78 59l16 -23q-32 -48 -68.5 -78t-65.5 -30q-19 0 -37.5 20t-29.5 53l-41 116q-72 -106 -114.5 -147.5t-79.5 -41.5q-21 0 -34.5 14t-13.5 37
q0 16 13.5 31.5t28.5 15.5q12 0 17 -11q5 -10 25 -10q22 0 57.5 36t89.5 111l-40 108q-22 58 -36 58q-21 0 -67 -57l-19 20q81 107 125 107q17 0 30 -22t39 -88l22 -55q68 92 108.5 128.5t74.5 36.5q20 0 32.5 -14t12.5 -30z" id="x1D465"></path></g>
</svg>.</table>

Abstract and Applied Analysis

Fast Computation of Singular Oscillatory Fourier Transforms

Figure 1