eic-opticks/ana/mm.py

0001 #!/usr/bin/env python
0002 """
0003 This demonstrates that volume_meshes.npy can be eliminated, as it matches volume_identity[:,1] for all mm
0004 """
0005 import os, numpy as np
0006
0007 m = np.load("volume_meshes.npy")
0008 i = np.load("volume_identity.npy")
0009 n = np.load("volume_nodeinfo.npy")
0010
0011 print("volume_meshes mesh_index\n",m)
0012 print("volume_identity node_index/mesh_index/boundary/identity_index(copyNumber) GVolume::getIdentity\n",i)
0013 print("volume_nodeinfo num_face/num_vert/node_index/parent_index\n",n)
0014
0015 assert np.all( i[:,1] ==  m[:,0] )
0016