def test_abp2latlon_2():
'''
abp2latlon(): check (ne=120, ei=84, ej=79, panel=2)
'''
from cs_ll import ij2ab, abp2latlon, latlon2abp
ne, ngq = 120, 4
panel = 2
gi1, gj1, ei1, ej1 = 4, 4, 84, 79
gi2, gj2, ei2, ej2 = 1, 4, 85, 79
gi3, gj3, ei3, ej3 = 1, 1, 85, 80
gi4, gj4, ei4, ej4 = 4, 1, 84, 80
a1, b1 = ij2ab(ne, ngq, ei1, ej1, gi1, gj1)
a2, b2 = ij2ab(ne, ngq, ei2, ej2, gi2, gj2)
a3, b3 = ij2ab(ne, ngq, ei3, ej3, gi3, gj3)
a4, b4 = ij2ab(ne, ngq, ei4, ej4, gi4, gj4)
lat1, lon1 = abp2latlon(a1,b1,panel)
lat2, lon2 = abp2latlon(a2,b2,panel)
lat3, lon3 = abp2latlon(a3,b3,panel)
lat4, lon4 = abp2latlon(a4,b4,panel)
aa_equal([lat1,lon1], [lat2,lon2], 15)
aa_equal([lat3,lon3], [lat2,lon2], 15)
aa_equal([lat3,lon3], [lat4,lon4], 15)
aa_equal([lat1,lon1], [lat4,lon4], 15)
def test_abp2latlon():
'''
abp2latlon(): MVP low accuracy
'''
from cs_ll import ij2ab, abp2latlon
ne, ngq = 120, 4
panel = 2
gi1, gj1, ei1, ej1 = 4, 4, 84, 79
gi2, gj2, ei2, ej2 = 1, 4, 85, 79
gi3, gj3, ei3, ej3 = 1, 1, 85, 80
gi4, gj4, ei4, ej4 = 4, 1, 84, 80
a1, b1 = ij2ab(ne, ngq, ei1, ej1, gi1, gj1)
a2, b2 = ij2ab(ne, ngq, ei2, ej2, gi2, gj2)
a3, b3 = ij2ab(ne, ngq, ei3, ej3, gi3, gj3)
a4, b4 = ij2ab(ne, ngq, ei4, ej4, gi4, gj4)
lat1, lon1 = abp2latlon(a1,b1,panel)
lat2, lon2 = abp2latlon(a2,b2,panel)
lat3, lon3 = abp2latlon(a3,b3,panel)
lat4, lon4 = abp2latlon(a4,b4,panel)
'''
print('')
print('{:.18f}, {:.18f}'.format(a1, b1))
print('{:.18f}, {:.18f}'.format(a2, b2))
#print('{:.18f}, {:.18f}'.format(a3, b3))
#print('{:.18f}, {:.18f}'.format(a4, b4))
print('')
print('{:.18f}, {:.18f}'.format(lat1, lon1))
print('{:.18f}, {:.18f}'.format(lat2, lon2))
#print('{:.18f}, {:.18f}'.format(lat3, lon3))
#print('{:.18f}, {:.18f}'.format(lat4, lon4))
'''
aa_equal([lat1,lon1], [lat2,lon2], 15)
aa_equal([lat2,lon2], [lat3,lon3], 15)
aa_equal([lat3,lon3], [lat4,lon4], 15)
aa_equal([lat1,lon1], [lat4,lon4], 15)
def test_ij2ab():
'''
ij2ab(): center of panel, at panel border
'''
from cs_ll import ij2ab
alpha, beta = ij2ab(ne=16, ngq=4, ei=1, ej=1, gi=1, gj=1)
a_equal([alpha,beta], [-pi/4,-pi/4])
alpha, beta = ij2ab(ne=16, ngq=4, ei=8, ej=8, gi=4, gj=4)
a_equal([alpha,beta], [0,0])
#------------------------------------------------
# MVP accuracy test
#------------------------------------------------
ne, ngq = 120, 4
panel = 2
gi1, gj1, ei1, ej1 = 4, 4, 84, 79
gi2, gj2, ei2, ej2 = 1, 4, 85, 79
gi3, gj3, ei3, ej3 = 1, 1, 85, 80
gi4, gj4, ei4, ej4 = 4, 1, 84, 80
a1, b1 = ij2ab(ne, ngq, ei1, ej1, gi1, gj1)
a2, b2 = ij2ab(ne, ngq, ei2, ej2, gi2, gj2)
a3, b3 = ij2ab(ne, ngq, ei3, ej3, gi3, gj3)
a4, b4 = ij2ab(ne, ngq, ei4, ej4, gi4, gj4)
#print('')
#print('{:.15f}, {:.15f}'.format(a1, b1))
#print('{:.15f}, {:.15f}'.format(a2, b2))
aa_equal([a1,b1], [a2,b2], 15)
aa_equal([a2,b2], [a3,b3], 15)
aa_equal([a3,b3], [a4,b4], 15)
aa_equal([a1,b1], [a4,b4], 15)
