Calculate Omega matrices

calc_omega(eval, D_l)

## Arguments

eval vector of eigenvalues from decomposition of relatedness matrix vector of length d_size

## Value

list of length 2. First entry in the list is the symmetric matrix OmegaU. Second entry in the list is the symmetric matrix OmegaE.

## Examples

calc_omega(eval = 50:1, D_l = runif(2))
#> [[1]]
#>            [,1]       [,2]       [,3]       [,4]       [,5]       [,6]
#> [1,] 0.01814778 0.01848321 0.01883127 0.01919269 0.01956826 0.01995882
#> [2,] 0.01960022 0.01999206 0.02039990 0.02082472 0.02126762 0.02172976
#>            [,7]       [,8]       [,9]      [,10]      [,11]      [,12]
#> [1,] 0.02036529 0.02078865 0.02122999 0.02169048 0.02217139 0.02267411
#> [2,] 0.02221243 0.02271703 0.02324509 0.02379828 0.02437844 0.02498760
#>           [,13]      [,14]      [,15]      [,16]      [,17]      [,18]
#> [1,] 0.02320015 0.02375118 0.02432903 0.02493569 0.02557338 0.02624454
#> [2,] 0.02562798 0.02630205 0.02701254 0.02776247 0.02855524 0.02939461
#>           [,19]      [,20]      [,21]      [,22]      [,23]      [,24]
#> [1,] 0.02695188 0.02769840 0.02848746 0.02932279 0.03020859 0.03114957
#> [2,] 0.03028482 0.03123063 0.03223743 0.03331130 0.03445918 0.03568899
#>           [,25]     [,26]      [,27]      [,28]      [,29]      [,30]
#> [1,] 0.03215107 0.0332191 0.03436052 0.03558318 0.03689606 0.03830953
#> [2,] 0.03700983 0.0384322 0.03996827 0.04163224 0.04344078 0.04541358
#>           [,31]      [,32]      [,33]      [,34]      [,35]      [,36]
#> [1,] 0.03983561 0.04148832 0.04328411 0.04524238 0.04738625 0.04974341
#> [2,] 0.04757409 0.04995043 0.05257666 0.05549437 0.05875494 0.06242257
#>           [,37]      [,38]      [,39]      [,40]      [,41]      [,42]
#> [1,] 0.05234734 0.05523895 0.05846870 0.06209958 0.06621128 0.07090606
#> [2,] 0.06657858 0.07132746 0.07680582 0.08319574 0.09074537 0.09980193
#>           [,43]     [,44]      [,45]      [,46]     [,47]     [,48]     [,49]
#> [1,] 0.07631743 0.0826230 0.09006439 0.09897886 0.1098519 0.1234085 0.1407823
#> [2,] 0.11086663 0.1246907 0.14245326 0.16611719 0.1992093 0.2487657 0.3311427
#>          [,50]
#> [1,] 0.1638493
#> [2,] 0.4950871
#>
#> [[2]]
#>           [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
#> [1,] 0.9073889 0.9056771 0.9039009 0.9020565 0.9001399 0.8981468 0.8960725
#> [2,] 0.9800108 0.9796112 0.9791952 0.9787620 0.9783103 0.9778390 0.9773467
#>           [,8]      [,9]     [,10]     [,11]     [,12]     [,13]     [,14]
#> [1,] 0.8939120 0.8916598 0.8893098 0.8868557 0.8842902 0.8816057 0.8787937
#> [2,] 0.9768321 0.9762936 0.9757294 0.9751377 0.9745165 0.9738634 0.9731759
#>          [,15]     [,16]     [,17]     [,18]     [,19]     [,20]     [,21]
#> [1,] 0.8758449 0.8727490 0.8694948 0.8660697 0.8624601 0.8586504 0.8546237
#> [2,] 0.9724514 0.9716865 0.9708780 0.9700220 0.9691141 0.9681495 0.9671228
#>          [,22]     [,23]     [,24]     [,25]     [,26]     [,27]     [,28]
#> [1,] 0.8503609 0.8458405 0.8410385 0.8359277 0.8304774 0.8246525 0.8184131
#> [2,] 0.9660276 0.9648569 0.9636027 0.9622556 0.9608051 0.9592385 0.9575415
#>          [,29]     [,30]     [,31]     [,32]     [,33]     [,34]    [,35]
#> [1,] 0.8117132 0.8045001 0.7967122 0.7882781 0.7791140 0.7691205 0.758180
#> [2,] 0.9556971 0.9536851 0.9514817 0.9490582 0.9463799 0.9434043 0.940079
#>          [,36]     [,37]     [,38]     [,39]     [,40]     [,41]     [,42]
#> [1,] 0.7461511 0.7328628 0.7181064 0.7016244 0.6830954 0.6621128 0.6381545
#> [2,] 0.9363386 0.9321001 0.9272569 0.9216698 0.9151531 0.9074537 0.8982173
#>          [,43]     [,44]     [,45]     [,46]     [,47]     [,48]     [,49]
#> [1,] 0.6105394 0.5783610 0.5403863 0.4948943 0.4394075 0.3702256 0.2815645
#> [2,] 0.8869330 0.8728347 0.8547195 0.8305860 0.7968371 0.7462971 0.6622853
#>          [,50]
#> [1,] 0.1638493
#> [2,] 0.4950871
#>