/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
|