def get_cube_points(origin, axes, nrep):
points = np.zeros((nrep[<NUM_LIT:0>], nrep[<NUM_LIT:1>], nrep[<NUM_LIT:2>], <NUM_LIT:3>), float)<EOL>points[:] = origin<EOL>points[:] += np.outer(np.arange(nrep[<NUM_LIT:0>], dtype=float), axes[<NUM_LIT:0>]).reshape((-<NUM_LIT:1>,<NUM_LIT:1>,<NUM_LIT:1>,<NUM_LIT:3>))<EOL>points[:] += np.outer(np.arange(nrep[<NUM_LIT:1>], dtype=float), axes[<NUM_LIT:1>]).reshape((<NUM_LIT:1>,-<NUM_LIT:1>,<NUM_LIT:1>,<NUM_LIT:3>))<EOL>points[:] += np.outer(np.arange(nrep[<NUM_LIT:2>], dtype=float), axes[<NUM_LIT:2>]).reshape((<NUM_LIT:1>,<NUM_LIT:1>,-<NUM_LIT:1>,<NUM_LIT:3>))<EOL>return points<EOL>
Generate the Cartesian coordinates of the points in a cube file

       *Arguemnts:*

       origin
            The cartesian coordinate for the origin of the grid.

       axes
            The 3 by 3 array with the grid spacings as rows.

       nrep
            The number of grid points along each axis.