Created
April 3, 2020 04:40
-
-
Save YingboMa/f8c52532eff56b1506a9d0e4aa0796ab to your computer and use it in GitHub Desktop.
Radau IIA transformed coefficients
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
A = [11//45-7sqrt(6)/360 37//225-169sqrt(6)/1800 -2//225+sqrt(6)/75 | |
37//225+169sqrt(6)/1800 11//45+7sqrt(6)/360 -2//225-sqrt(6)/75 | |
4//9-sqrt(6)/36 4//9+sqrt(6)/36 1//9] | |
T = Float64 | |
T11 = convert(T, 9.1232394870892942792e-02) | |
T12 = convert(T, -0.14125529502095420843e0) | |
T13 = convert(T, -3.0029194105147424492e-02) | |
T21 = convert(T, 0.24171793270710701896e0) | |
T22 = convert(T, 0.20412935229379993199e0) | |
T23 = convert(T, 0.38294211275726193779e0) | |
T31 = convert(T, 0.96604818261509293619e0) | |
T32 = 1 | |
T33 = 0 | |
TI11 = convert(T, 4.3255798900631553510e0) | |
TI12 = convert(T, 0.33919925181580986954e0) | |
TI13 = convert(T, 0.54177053993587487119e0) | |
TI21 = convert(T, -4.1787185915519047273e0) | |
TI22 = convert(T, -0.32768282076106238708e0) | |
TI23 = convert(T, 0.47662355450055045196e0) | |
TI31 = convert(T, -0.50287263494578687595e0) | |
TI32 = convert(T, 2.5719269498556054292e0) | |
TI33 = convert(T, -0.59603920482822492497e0) | |
T = [T11 T12 T13 | |
T21 T22 T23 | |
T31 T32 T33] | |
foreach(x->normalize!(x), eachcol(T)) | |
V = eigvecs(inv(A)) | |
(V[:, 1] - conj(V[:, 1])) |> normalize | |
T | |
TI = [TI11 TI12 TI13 | |
TI21 TI22 TI23 | |
TI31 TI32 TI33] |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment