/usr/lib/python2.7/dist-packages/FIAT/asci2vtk3d.py is in python-fiat 1.6.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 | # Copyright (C) 2008-2012 Robert C. Kirby (Texas Tech University)
#
# This file is part of FIAT.
#
# FIAT is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# FIAT is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with FIAT. If not, see <http://www.gnu.org/licenses/>.
#!/usr/bin/env python
# 3d mode: x y z f, f = f(x,y,z)
import sys
if len(sys.argv) > 1:
filename = sys.argv[1]
print(filename)
base = filename.split(".")[0]
output = "%s.vtk" % (base,)
print("output to %s" % (output,))
else:
print("python asci2vtk3d.py foo")
sys.exit(0)
fin = open( filename, "r" )
coords = [ ]
for line in fin:
coords.append( line.split() )
fin.close()
n = len( coords )
print("%s points" % (str(n),))
fout = open( output, "w" )
fout.write("""# vtk DataFile Version 2.0
points
ASCII
DATASET UNSTRUCTURED_GRID
POINTS %s float\n""" % (str(n),))
for c in coords:
fout.write("%s %s %s\n" % (c[0], c[1], c[2]))
fout.write("CELLS %s %s\n" % (n, 2*n))
for i in range( n ):
fout.write("1 %s\n" % (i,))
fout.write("CELL_TYPES %s\n" % (n,))
for i in range( n ):
fout.write("1\n")
fout.write("POINT_DATA %s\n" % (n,))
fout.write("""SCALARS Z float 1
LOOKUP_TABLE default\n""")
for i in range( n ):
fout.write("%s\n", ncoords[i][3])
fout.close()
|