This function performs the integration near resonances as described in Section 3.1 of the code paper for a relativistic calculation. It is only called if resonances are present in or near the integration domain.

Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|

integer, | intent(in) | :: | sproc_rel |
is_rel of the current process. |
||

double complex, | intent(in) | :: | om |
Complex wave frequency . |
||

integer, | intent(in) | :: | nn |
Order of the Bessel function. |
||

integer, | intent(in) | :: | mode |
Index of the entries in the T-tensor of Eq. (2.10). |
||

integer, | intent(in) | :: | igamma |
Index to loop over . |
||

integer, | intent(in) | :: | ipparbar_res |
Index of the nearest relativistic parallel momentum to the resonance. |
||

integer, | intent(in) | :: | upperlimit |
Index of upper limit for integration according to Eq. (3.5). |