Aquila-Consortium/libsharp2
0
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

launch test

Browse source
This commit is contained in:
Andrea Zonca 2017-11-15 07:27:09 -08:00
parent 8c33a5e699
commit cdb11e97d0
1 changed files with 5 additions and 7 deletions
Download patch file Download diff file Expand all files Collapse all files

12
python/libsharp/tests/test_sht.py
View file

@ -1,7 +1,4 @@
import numpy as np import numpy as np
import healpy
from scipy.special import legendre
from scipy.special import p_roots
from numpy.testing import assert_allclose from numpy.testing import assert_allclose
import libsharp import libsharp
@ -28,7 +25,8 @@ def test_basic():
map = libsharp.synthesis(grid, order, np.repeat(alm[None, None, :], 3, 0), comm=MPI.COMM_WORLD) map = libsharp.synthesis(grid, order, np.repeat(alm[None, None, :], 3, 0), comm=MPI.COMM_WORLD)
assert np.all(map[2, :] == map[1, :]) and np.all(map[1, :] == map[0, :]) assert np.all(map[2, :] == map[1, :]) and np.all(map[1, :] == map[0, :])
map = map[0, 0, :] map = map[0, 0, :]
if rank == 0: print(rank, "shape", map.shape)
healpy.mollzoom(map) print(rank, "mean", map.mean())
from matplotlib.pyplot import show
show() if __name__=="__main__":
test_basic()