This file is indexed.

/usr/bin/freefoam-graphResUVWP is in freefoam 0.1.0+dfsg+1-3.

This file is owned by root:root, with mode 0o755.

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
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
#! /usr/bin/python
#-------------------------------------------------------------------------------
#               ______                _     ____          __  __
#              |  ____|             _| |_  / __ \   /\   |  \/  |
#              | |__ _ __ ___  ___ /     \| |  | | /  \  | \  / |
#              |  __| '__/ _ \/ _ ( (| |) ) |  | |/ /\ \ | |\/| |
#              | |  | | |  __/  __/\_   _/| |__| / ____ \| |  | |
#              |_|  |_|  \___|\___|  |_|   \____/_/    \_\_|  |_|
#
#                   FreeFOAM: The Cross-Platform CFD Toolkit
#
# Copyright (C) 2008-2012 Michael Wild <themiwi@users.sf.net>
#                         Gerber van der Graaf <gerber_graaf@users.sf.net>
#-------------------------------------------------------------------------------
# License
#   This file is part of FreeFOAM.
#
#   FreeFOAM is free software: you can redistribute it and/or modify it
#   under the terms of the GNU General Public License as published by the
#   Free Software Foundation, either version 3 of the License, or (at your
#   option) any later version.
#
#   FreeFOAM 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 General Public License
#   for more details.
#
#   You should have received a copy of the GNU General Public License
#   along with FreeFOAM.  If not, see <http://www.gnu.org/licenses/>.
#
# Script
#     freefoam-graphResUVWP
#
# Description
#     The utility extracts the residuals of the velocity components 'U', 'V'
#     and 'W' and the pressure 'p' at each time step and writes them to
#     'residualUVWP.dat' in TecPlot format.
#------------------------------------------------------------------------------

"""Usage: freefoam graphResUVWP [-h, -help] <logFile>

The utility extracts the residuals of the velocity components 'U', 'V' and 'W'
and the pressure 'p' at each time step and writes them to 'residualUVWP.dat' in
TecPlot format.

Options
-------
<logFile> The log file from which to extract the data.
-h, -help Display this help message.

"""

import sys
import os.path
import re
# want to be future proof
sys.path.insert(0, '/usr/lib/python2.7/site-packages')
from FreeFOAM.compat import *

if len(sys.argv[1:]) != 1:
   sys.stderr.write('Error: Require a log file as argument\n')
   sys.stderr.write(__doc__+'\n')
   sys.exit(1)

if sys.argv[1] == '-h' or sys.argv[1] == '-help':
   print(__doc__)
   sys.exit(0)

logFile = sys.argv[1]
if not os.path.isfile(logFile):
   sys.stderr.write('Error: no such file')
   sys.exit(1)

# parse the residuals from the log file
regex = re.compile(r'Solving for (?P<var>U[xyz]|p),\s+.*Initial residual\s+=\s+(?P<res>\S+),')
res = {}
for v in 'Ux Uy Uz p'.split():
   res[v] = []
for l in open(logFile, 'rt'):
   m = regex.search(l)
   if m:
      res[m.group('var')].append(float(m.group('res')))

# write file header
ngraphs = 0
for v in res.values():
   if len(v) > 0:
      ngraphs += 1
outFile = open('residualUVWP.dat', 'wt')
outFile.write("""Solver Performance
number of iterations
residual
0 0
0 0
%d
"""%ngraphs)

# write data
for n, r in res.items():
   if len(r) > 0:
      outFile.write("""%s
0
%d
"""%(n, len(r)))
      for i, v in enumerate(r):
         outFile.write('%d  %g\n'%(i+1, v))

outFile.close()

# ------------------- vim: set sw=3 sts=3 ft=python et: ------------ end-of-file