1
+ from __future__ import unicode_literals
2
+
1
3
from pyslalib import slalib as S
2
4
import numpy as N
3
5
import numpy .testing as T
@@ -216,13 +218,13 @@ def testcc62s(self):
216
218
217
219
def testcd2tf (self ):
218
220
s , ihmsf = S .sla_cd2tf (4 , - 0.987654321 )
219
- T .assert_equal (s , '-' , 'sla_cd2tf, s' )
221
+ T .assert_equal (s , b '-' , 'sla_cd2tf, s' )
220
222
T .assert_equal (ihmsf [0 ], 23 , 'sla_cd2tf, (1)' )
221
223
T .assert_equal (ihmsf [1 ], 42 , 'sla_cd2tf, (2)' )
222
224
T .assert_equal (ihmsf [2 ], 13 , 'sla_cd2tf, (3)' )
223
225
T .assert_almost_equal (ihmsf [3 ], 3333 , - 3 , 'sla_cd2tf, (4)' )
224
226
s , ihmsf = S .sla_dd2tf (4 , - 0.987654321 )
225
- T .assert_equal (s , '-' , 'sla_dd2tf, s' )
227
+ T .assert_equal (s , b '-' , 'sla_dd2tf, s' )
226
228
T .assert_equal (ihmsf [0 ], 23 , 'sla_dd2tf, (1)' )
227
229
T .assert_equal (ihmsf [1 ], 42 , 'sla_dd2tf, (2)' )
228
230
T .assert_equal (ihmsf [2 ], 13 , 'sla_dd2tf, (3)' )
@@ -235,27 +237,27 @@ def testcldj(self):
235
237
236
238
def testcr2af (self ):
237
239
s , idmsf = S .sla_cr2af (4 , 2.345 )
238
- T .assert_equal (s , '+' , 'sla_cr2af, s' )
240
+ T .assert_equal (s , b '+' , 'sla_cr2af, s' )
239
241
T .assert_equal (idmsf [0 ], 134 , 'sla_cr2af, (1)' )
240
242
T .assert_equal (idmsf [1 ], 21 , 'sla_cr2af, (2)' )
241
243
T .assert_equal (idmsf [2 ], 30 , 'sla_cr2af, (3)' )
242
244
T .assert_almost_equal (idmsf [3 ], 9706 , - 3 , 'sla_cr2af, (4)' )
243
245
s , idmsf = S .sla_dr2af (4 , 2.345 )
244
- T .assert_equal (s , '+' , 'sla_dr2af, s' )
246
+ T .assert_equal (s , b '+' , 'sla_dr2af, s' )
245
247
T .assert_equal (idmsf [0 ], 134 , 'sla_dr2af, (1)' )
246
248
T .assert_equal (idmsf [1 ], 21 , 'sla_dr2af, (2)' )
247
249
T .assert_equal (idmsf [2 ], 30 , 'sla_dr2af, (3)' )
248
250
T .assert_equal (idmsf [3 ], 9706 , 'sla_dr2af, (4)' )
249
251
250
252
def testcr2tf (self ):
251
253
s , ihmsf = S .sla_cr2tf (4 , - 3.01234 )
252
- T .assert_equal (s , '-' , 'sla_cr2tf, s' )
254
+ T .assert_equal (s , b '-' , 'sla_cr2tf, s' )
253
255
T .assert_equal (ihmsf [0 ], 11 , 'sla_cr2tf, (1)' )
254
256
T .assert_equal (ihmsf [1 ], 30 , 'sla_cr2tf, (2)' )
255
257
T .assert_equal (ihmsf [2 ], 22 , 'sla_cr2tf, (3)' )
256
258
T .assert_almost_equal (ihmsf [3 ], 6484 , - 3 , 'sla_cr2tf, (4)' )
257
259
s , ihmsf = S .sla_dr2tf (4 , - 3.01234 )
258
- T .assert_equal (s , '-' , 'sla_dr2tf, s' )
260
+ T .assert_equal (s , b '-' , 'sla_dr2tf, s' )
259
261
T .assert_equal (ihmsf [0 ], 11 , 'sla_dr2tf, (1)' )
260
262
T .assert_equal (ihmsf [1 ], 30 , 'sla_dr2tf, (2)' )
261
263
T .assert_equal (ihmsf [2 ], 22 , 'sla_dr2tf, (3)' )
@@ -493,13 +495,13 @@ def testfitxy(self):
493
495
45 , 32.5 , 8.6 , - 17 ,
494
496
15.3 , 10 , 121.7 , - 3.8 ])
495
497
# The reshaping is necessary due to C/Fortran order
496
- xye = N .reshape (xye , (2 ,8 ), order = 'Fortran ' )
498
+ xye = N .reshape (xye , (2 ,8 ), order = 'F ' )
497
499
xym = N .asarray ([- 23.41 , 12.12 , 32.03 , 15.34 ,
498
500
10.93 , - 23.72 , - 3.01 , - 16.10 ,
499
501
44.90 , - 32.46 , 8.55 , 17.02 ,
500
502
15.31 , - 10.07 , 120.92 , 3.81 ])
501
503
# The reshaping is necessary due to C/Fortran order
502
- xym = N .reshape (xym , (2 ,8 ), order = 'Fortran ' )
504
+ xym = N .reshape (xym , (2 ,8 ), order = 'F ' )
503
505
coeffs , j = S .sla_fitxy (4 , xye , xym )
504
506
ans1 = N .asarray ([- 7.938263381515947e-3 , 1.004640925187200 ,
505
507
3.976948048238268e-4 , - 2.501031681585021e-2 ,
@@ -522,7 +524,7 @@ def testfitxy(self):
522
524
15.348618307280820 , 10.07063070741086835 ,
523
525
121.5833272936291482 , - 3.788442308260240 ])
524
526
# And another reshaping...
525
- ans3 = N .reshape (ans3 , (2 ,8 ), order = 'Fortran ' )
527
+ ans3 = N .reshape (ans3 , (2 ,8 ), order = 'F ' )
526
528
T .assert_array_almost_equal (xyp , ans3 , 12 , 'sla_fitxy, xyp' )
527
529
T .assert_almost_equal (xrms ,0.1087247110488075 , 13 , 'sla_pxy, xrms' )
528
530
T .assert_almost_equal (yrms , 0.03224481175794666 , 13 , 'sla_pxy, yrms' )
@@ -721,29 +723,29 @@ def testkbj(self):
721
723
k = '?'
722
724
e = 1950
723
725
k , j = S .sla_kbj (- 1 , e )
724
- T .assert_equal (k , ' ' , 'sla_kbj, jb1' )
726
+ T .assert_equal (k , b ' ' , 'sla_kbj, jb1' )
725
727
T .assert_equal (j , 1 , 'sla_kbj, j1' )
726
728
k , j = S .sla_kbj (0 , e )
727
- T .assert_equal (k , 'B' , 'sla_kbj, jb2' )
729
+ T .assert_equal (k , b 'B' , 'sla_kbj, jb2' )
728
730
T .assert_equal (j , 0 , 'sla_kbj, j2' )
729
731
k , j = S .sla_kbj (1 , e )
730
- T .assert_equal (k , 'B' , 'sla_kbj, jb3' )
732
+ T .assert_equal (k , b 'B' , 'sla_kbj, jb3' )
731
733
T .assert_equal (j , 0 , 'sla_kbj, j3' )
732
734
k , j = S .sla_kbj (2 , e )
733
- T .assert_equal (k , 'J' , 'sla_kbj, jb4' )
735
+ T .assert_equal (k , b 'J' , 'sla_kbj, jb4' )
734
736
T .assert_equal (j , 0 , 'sla_kbj, j4' )
735
737
k , j = S .sla_kbj (3 , e )
736
- T .assert_equal (k , ' ' , 'sla_kbj, jb5' )
738
+ T .assert_equal (k , b ' ' , 'sla_kbj, jb5' )
737
739
T .assert_equal (j , 1 , 'sla_kbj, j5' )
738
740
e = 2000
739
741
k , j = S .sla_kbj (0 , e )
740
- T .assert_equal (k , 'J' , 'sla_kbj, jb6' )
742
+ T .assert_equal (k , b 'J' , 'sla_kbj, jb6' )
741
743
T .assert_equal (j , 0 , 'sla_kbj, j6' )
742
744
k , j = S .sla_kbj (1 , e )
743
- T .assert_equal (k , 'B' , 'sla_kbj, jb7' )
745
+ T .assert_equal (k , b 'B' , 'sla_kbj, jb7' )
744
746
T .assert_equal (j , 0 , 'sla_kbj, j7' )
745
747
k , j = S .sla_kbj (2 , e )
746
- T .assert_equal (k , 'J' , 'sla_kbj, jb8' )
748
+ T .assert_equal (k , b 'J' , 'sla_kbj, jb8' )
747
749
T .assert_equal (j , 0 , 'sla_kbj, j8' )
748
750
749
751
def testmap (self ):
@@ -808,31 +810,31 @@ def testobs(self):
808
810
n = 0
809
811
c = 'MMT'
810
812
c , name , w , p , h = S .sla_obs (n , c )
811
- T .assert_equal (c .strip (), 'MMT' , 'sla_obs, 1/c' )
812
- T .assert_equal (name .strip (), 'MMT 6.5m, Mt Hopkins' , 'sla_obs, 1/name' )
813
+ T .assert_equal (c .strip (), b 'MMT' , 'sla_obs, 1/c' )
814
+ T .assert_equal (name .strip (), b 'MMT 6.5m, Mt Hopkins' , 'sla_obs, 1/name' )
813
815
T .assert_almost_equal (w , 1.935300584055477 , 8 , 'sla_obs, 1/w' )
814
816
T .assert_almost_equal (p , 0.5530735081550342238 , 10 , 'sla_obs, 1/p' )
815
817
T .assert_almost_equal (h , 2608 , 10 , 'sla_obs, 1/h' )
816
818
n = 61
817
819
c = 20 * ' '
818
820
c , name , w , p , h = S .sla_obs (n , c )
819
- T .assert_equal (c .strip (), 'KECK1' , 'sla_obs, 2/c' )
820
- T .assert_equal (name .strip (), 'Keck 10m Telescope #1' , 'sla_obs, 2/name' )
821
+ T .assert_equal (c .strip (), b 'KECK1' , 'sla_obs, 2/c' )
822
+ T .assert_equal (name .strip (), b 'Keck 10m Telescope #1' , 'sla_obs, 2/name' )
821
823
T .assert_almost_equal (w , 2.713545757918895 , 8 , 'sla_obs, 2/w' )
822
824
T .assert_almost_equal (p , 0.3460280563536619 , 8 , 'sla_obs, 2/p' )
823
825
T .assert_almost_equal (h , 4160 , 10 , 'sla_obs, 2/h' )
824
826
n = 83
825
827
c = 20 * ' '
826
828
c , name , w , p , h = S .sla_obs (n , c )
827
- T .assert_equal (c .strip (), 'MAGELLAN2' , 'sla_obs, 3/c' )
828
- T .assert_equal (name .strip (), 'Magellan 2, 6.5m, Las Campanas' , 'sla_obs, 3/name' )
829
+ T .assert_equal (c .strip (), b 'MAGELLAN2' , 'sla_obs, 3/c' )
830
+ T .assert_equal (name .strip (), b 'Magellan 2, 6.5m, Las Campanas' , 'sla_obs, 3/name' )
829
831
T .assert_almost_equal (w , 1.233819305534497 , 8 , 'sla_obs, 3/w' )
830
832
T .assert_almost_equal (p , - 0.506389344359954 , 8 , 'sla_obs, 3/p' )
831
833
T .assert_almost_equal (h , 2408 , 10 , 'sla_obs, 3/h' )
832
834
n = 84
833
835
c = 20 * ' '
834
836
c , name , w , p , h = S .sla_obs (n , c )
835
- T .assert_equal (name .strip (), '?' , 'sla_obs, 4/name' )
837
+ T .assert_equal (name .strip (), b '?' , 'sla_obs, 4/name' )
836
838
837
839
def testpa (self ):
838
840
T .assert_almost_equal (S .sla_pa (- 1.567 , 1.5123 , 0.987 ),
@@ -1192,13 +1194,13 @@ def testsmat(self):
1192
1194
a = N .asarray ([[2.22 , 1.6578 , 1.380522 ],
1193
1195
[1.6578 , 1.380522 , 1.22548578 ],
1194
1196
[1.380522 , 1.22548578 , 1.1356276122 ]],
1195
- order = 'Fortran ' )
1197
+ order = 'F ' )
1196
1198
v = N .asarray ([2.28625 , 1.7128825 , 1.429432225 ])
1197
1199
a , v , d , j , iw = S .sla_smat (a , v , N .empty (3 , dtype = N .float32 ))
1198
1200
ans = N .asarray ([[18.02550629769198 , - 52.16386644917280607 , 34.37875949717850495 ],
1199
1201
[- 52.16386644917280607 , 168.1778099099805627 , - 118.0722869694232670 ],
1200
1202
[34.37875949717850495 , - 118.0722869694232670 , 86.50307003740151262 ]],
1201
- order = 'Fortran ' )
1203
+ order = 'F ' )
1202
1204
T .assert_array_almost_equal (a , ans , 2 , 'sla_smat, a(3,3)' )
1203
1205
ans = N .asarray ([1.002346480763383 , 0.03285594016974583489 , 0.004760688414885247309 ])
1204
1206
T .assert_array_almost_equal (v , ans , 4 , 'sla_smat, v(3)' )
@@ -1216,7 +1218,7 @@ def testsvd(self):
1216
1218
nc = 7
1217
1219
m = 5
1218
1220
n = 4
1219
- a = N .empty ((m , n ), order = 'Fortran ' )
1221
+ a = N .empty ((m , n ), order = 'F ' )
1220
1222
b = N .empty (m )
1221
1223
for ii in range (m ):
1222
1224
val = 0.5 * (ii + 1.0 )
0 commit comments