eic-opticks/ana/NRngDiffuseTest.py

0001 #!/usr/bin/env python
0002 #
0003 # Copyright (c) 2019 Opticks Team. All Rights Reserved.
0004 #
0005 # This file is part of Opticks
0006 # (see https://bitbucket.org/simoncblyth/opticks).
0007 #
0008 # Licensed under the Apache License, Version 2.0 (the "License");
0009 # you may not use this file except in compliance with the License.
0010 # You may obtain a copy of the License at
0011 #
0012 #   http://www.apache.org/licenses/LICENSE-2.0
0013 #
0014 # Unless required by applicable law or agreed to in writing, software
0015 # distributed under the License is distributed on an "AS IS" BASIS,
0016 # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
0017 # See the License for the specific language governing permissions and
0018 # limitations under the License.
0019 #
0020
0021 """
0022 NRngDiffuseTest.py
0023 =======================
0024
0025 See npy-/tests/NRngDiffuseTest.cc
0026
0027 """
0028 import os, sys, logging
0029 import numpy as np
0030 log = logging.getLogger(__name__)
0031 import matplotlib.pyplot as plt
0032 from mpl_toolkits.mplot3d import Axes3D
0033
0034 np.set_printoptions(suppress=True, precision=3)
0035
0036
0037 if __name__ == '__main__':
0038
0039
0040     plt.ion()
0041     fig = plt.figure()
0042
0043
0044
0045     d = np.load(os.path.expandvars("$TMP/NRngDiffuseTest_diffuse.npy"))
0046     s = np.load(os.path.expandvars("$TMP/NRngDiffuseTest_sphere.npy"))
0047
0048
0049     ax = fig.add_subplot(111, projection='3d')
0050     ax.scatter(d[:,0], d[:,1], d[:,2])
0051
0052
0053