bin/foamGraphResKE

   1 #!/bin/csh -f
   2 #------------------------------------------------------------------------------
   3 # =========                 |
   4 # \\      /  F ield         | foam-extend: Open Source CFD
   5 #  \\    /   O peration     |
   6 #   \\  /    A nd           | For copyright notice see file Copyright
   7 #    \\/     M anipulation  |
   8 #------------------------------------------------------------------------------
   9 # License
  10 #     This file is part of foam-extend.
  11 #
  12 #     foam-extend is free software: you can redistribute it and/or modify it
  13 #     under the terms of the GNU General Public License as published by the
  14 #     Free Software Foundation, either version 3 of the License, or (at your
  15 #     option) any later version.
  16 #
  17 #     foam-extend is distributed in the hope that it will be useful, but
  18 #     WITHOUT ANY WARRANTY; without even the implied warranty of
  19 #     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  20 #     General Public License for more details.
  21 #
  22 #     You should have received a copy of the GNU General Public License
  23 #     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  24 #
  25 # Script
  26 #     foamGraphResKE
  27 #
  28 # Description
  29 #
  30 #------------------------------------------------------------------------------
  31
  32 set tmpK = newK$$.dat
  33 cat $1 | grep "Solving for k" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpK
  34
  35 set tmpEpsilon = newEpsilon$$.dat
  36 cat $1 | grep "Solving for epsilon" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpEpsilon
  37
  38 set linesK = `cat $tmpK | wc -l`
  39 set linesEpsilon = `cat $tmpEpsilon | wc -l`
  40
  41 set ngraphs = 0
  42
  43 if ($linesK > 1) then
  44 @ ngraphs = $ngraphs + 1
  45 endif
  46
  47 if ($linesEpsilon > 1) then
  48 @ ngraphs = $ngraphs + 1
  49 endif
  50
  51 cat << EOF > residualKE.dat
  52 Solver Performance
  53 number of iterations
  54 residual
  55 0 0
  56 0 0
  57 $ngraphs
  58 EOF
  59
  60 if ($linesK > 1) then
  61 cat << EOF >> residualKE.dat
  62 k
  63 0
  64 $linesK
  65 EOF
  66 cat $tmpK >> residualKE.dat
  67 endif
  68
  69 if ($linesEpsilon > 1) then
  70 cat << EOF >> residualKE.dat
  71 epsilon
  72 0
  73 $linesEpsilon
  74 EOF
  75 cat $tmpEpsilon >> residualKE.dat
  76 endif
  77
  78 rm $tmpK
  79 rm $tmpEpsilon
  80
  81 #------------------------------------------------------------------------------