eic-opticks/ana/cfgeocache.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 See also ana/geocache.bash
0023 """
0024 import os, logging, numpy as np
0025 log = logging.getLogger(__name__)
0026
0027 idp_ = lambda _:os.path.expandvars("$IDPATH/%s" % _ )
0028 idp2_ = lambda _:os.path.expandvars("$IDPATH2/%s" % _ )
0029
0030
0031
0032 def cflib(aa, bb):
0033     """
0034     Compare buffers between two geocache
0035     """
0036     assert aa.shape == bb.shape
0037     print aa.shape
0038
0039     for i in range(len(aa)):
0040         a = aa[i]
0041         b = bb[i]
0042         assert len(a) == 2
0043         assert len(b) == 2
0044
0045         g0 = a[0] - b[0]
0046         g1 = a[1] - b[1]
0047
0048         assert g0.shape == g1.shape
0049
0050         print i, g0.shape, "g0max: ", np.max(g0), "g1max: ", np.max(g1)
0051     pass
0052
0053
0054 def cflib_test():
0055     #rel = "GMaterialLib/GMaterialLib.npy"
0056     rel = "GSurfaceLib/GSurfaceLib.npy"
0057
0058     aa = np.load(idp_(rel))
0059     bb = np.load(idp2_(rel))
0060     cflib(aa,bb)
0061
0062
0063
0064 if __name__ == '__main__':
0065     pass
0066     logging.basicConfig(level=logging.INFO)
0067
0068     assert 0, "ancient code depending in IDPATH"
0069