Checking the stability regions of omegas for the parametrized Cumulant Just omega(2) is given in input. omega(2)=-1 means omega2=omegaBulk. Limiters read from input. lim(N)=10^10 means unlimited. lim(N) is for omega(N+2). Just omega(3:5) are limited as in the paper. omega(6:10) = 1._rk

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

real(kind=rk), | intent(in) | :: | omegaVisc(:) | vector of omegas in the level |
||

real(kind=rk), | intent(in) | :: | omegaBulk | omega bulk value in the level |
||

real(kind=rk), | intent(in) | :: | omegaIn(:) | vector of omegas as given in musubi.lua |
||

integer, | intent(in) | :: | nSolve | number of elements solved in kernel |
||

integer, | intent(in) | :: | level | current level |

Type | Visibility | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|---|

real(kind=rk), | private | :: | omega(10) | ||||

real(kind=rk), | private | :: | omegaPrev | ||||

integer, | private | :: | iElem | ||||

logical, | private | :: | test(5) |