def read_cyc(this, fn, conv=<NUM_LIT:1.0>):
f = paropen(fn, "<STR_LIT:r>")<EOL>f.readline()<EOL>f.readline()<EOL>f.readline()<EOL>f.readline()<EOL>cell = np.array( [ [ <NUM_LIT:0.0>, <NUM_LIT:0.0>, <NUM_LIT:0.0> ], [ <NUM_LIT:0.0>, <NUM_LIT:0.0>, <NUM_LIT:0.0> ], [ <NUM_LIT:0.0>, <NUM_LIT:0.0>, <NUM_LIT:0.0> ] ] )<EOL>l = f.readline()<EOL>s = map(float, l.split())<EOL>cell[<NUM_LIT:0>, <NUM_LIT:0>] = s[<NUM_LIT:0>]*conv<EOL>cell[<NUM_LIT:1>, <NUM_LIT:0>] = s[<NUM_LIT:1>]*conv<EOL>cell[<NUM_LIT:2>, <NUM_LIT:0>] = s[<NUM_LIT:2>]*conv<EOL>l = f.readline()<EOL>s = map(float, l.split())<EOL>cell[<NUM_LIT:0>, <NUM_LIT:1>] = s[<NUM_LIT:0>]*conv<EOL>cell[<NUM_LIT:1>, <NUM_LIT:1>] = s[<NUM_LIT:1>]*conv<EOL>cell[<NUM_LIT:2>, <NUM_LIT:1>] = s[<NUM_LIT:2>]*conv<EOL>l = f.readline()<EOL>s = map(float, l.split())<EOL>cell[<NUM_LIT:0>, <NUM_LIT:2>] = s[<NUM_LIT:0>]*conv<EOL>cell[<NUM_LIT:1>, <NUM_LIT:2>] = s[<NUM_LIT:1>]*conv<EOL>cell[<NUM_LIT:2>, <NUM_LIT:2>] = s[<NUM_LIT:2>]*conv<EOL>this.set_cell(cell)<EOL>this.set_pbc(True)<EOL>f.close()<EOL>
Read the lattice information from a cyc.dat file (i.e., tblmd input file)