qudarap/tests/QPMT_vs_SPMT_Test.py

0001 #!/usr/bin/env python
0002
0003 import numpy as np, textwrap
0004 from opticks.ana.fold import Fold
0005
0006 if __name__ == '__main__':
0007
0008     s = Fold.Load("$SFOLD", symbol="s")
0009     q = Fold.Load("$QFOLD", symbol="q")
0010
0011     print(repr(s))
0012     print(repr(q))
0013
0014     exprs = textwrap.dedent(r"""
0015     s.test.get_stackspec.shape
0016     q.stackspec_interp.shape
0017     np.abs(s.test.get_stackspec-q.stackspec_interp).max()
0018     """).split("\n")
0019
0020     for expr in exprs:
0021         print("[%s]"%expr)
0022         if len(expr) == 0: continue
0023         print(eval(expr))
0024     pass
0025
0026
0027
0028
0029