@@ 1067-1084 (lines=18) @@ | ||
1064 | assert_allclose(e3, e3_) |
|
1065 | ||
1066 | ||
1067 | def test_basevectors_apex_scalar_qd(): |
|
1068 | apex_out = Apex(date=2000, refh=300) |
|
1069 | ||
1070 | (f1, f2, f3, g1, g2, g3, d1, d2, d3, e1, e2, |
|
1071 | e3) = apex_out.basevectors_apex(60, 15, 100, coords='qd', precision=1e-2) |
|
1072 | ||
1073 | glat, glon, _ = apex_out.qd2geo(60, 15, 100, precision=1e-2) |
|
1074 | (_, _, _, _, f1_, f2_, _, d1_, d2_, d3_, _, e1_, e2_, |
|
1075 | e3_) = apex_out._geo2apexall(glat, glon, 100) |
|
1076 | ||
1077 | assert_allclose(f1, f1_) |
|
1078 | assert_allclose(f2, f2_) |
|
1079 | assert_allclose(d1, d1_) |
|
1080 | assert_allclose(d2, d2_) |
|
1081 | assert_allclose(d3, d3_) |
|
1082 | assert_allclose(e1, e1_) |
|
1083 | assert_allclose(e2, e2_) |
|
1084 | assert_allclose(e3, e3_) |
|
1085 | ||
1086 | ||
1087 | # test shapes and vectorization of arguments |
|
@@ 1047-1064 (lines=18) @@ | ||
1044 | assert_allclose(e3, e3_) |
|
1045 | ||
1046 | ||
1047 | def test_basevectors_apex_scalar_apex(): |
|
1048 | apex_out = Apex(date=2000, refh=300) |
|
1049 | ||
1050 | (f1, f2, f3, g1, g2, g3, d1, d2, d3, e1, e2, |
|
1051 | e3) = apex_out.basevectors_apex(60, 15, 100, coords='apex', precision=1e-2) |
|
1052 | ||
1053 | glat, glon, _ = apex_out.apex2geo(60, 15, 100, precision=1e-2) |
|
1054 | (_, _, _, _, f1_, f2_, _, d1_, d2_, d3_, _, e1_, e2_, |
|
1055 | e3_) = apex_out._geo2apexall(glat, glon, 100) |
|
1056 | ||
1057 | assert_allclose(f1, f1_) |
|
1058 | assert_allclose(f2, f2_) |
|
1059 | assert_allclose(d1, d1_) |
|
1060 | assert_allclose(d2, d2_) |
|
1061 | assert_allclose(d3, d3_) |
|
1062 | assert_allclose(e1, e1_) |
|
1063 | assert_allclose(e2, e2_) |
|
1064 | assert_allclose(e3, e3_) |
|
1065 | ||
1066 | ||
1067 | def test_basevectors_apex_scalar_qd(): |