From 3481fc9ebee180f7583aba9f283dfc15c9c47950 Mon Sep 17 00:00:00 2001 From: "David L. F. Gaden" Date: Mon, 4 Jul 2011 16:47:59 -0500 Subject: [PATCH] Added equationReader tutorial files to tutorials/ --- .../equationReaderDemo/expectedOutput | 55 + .../equationReaderDemo/system/controlDict | 47 + .../equationReaderDemo/testDict | 47 + .../equationReaderDemo/testDict2 | 25 + .../equationReaderDemo/testDict3 | 34 + .../equationReader/equationReaderTester/0/U | 41 + .../equationReaderTester/0/eqns | 201 ++ .../equationReader/equationReaderTester/0/p | 39 + .../constant/polyMesh/blockMeshDict | 63 + .../constant/polyMesh/boundary | 40 + .../constant/transportProperties | 19 + .../equationReaderTester/expectedOutput | 2968 +++++++++++++++++ .../equationReaderTester/system/controlDict | 33 + .../equationReaderTester/system/fvSchemes | 58 + .../equationReaderTester/system/fvSolution | 44 + .../equationReaderTester/testDict1 | 49 + .../equationReaderTester/testDict10 | 27 + .../equationReaderTester/testDict2 | 27 + .../equationReaderTester/testDict3 | 30 + .../equationReaderTester/testDict4 | 27 + .../equationReaderTester/testDict5 | 27 + .../equationReaderTester/testDict6 | 27 + .../equationReaderTester/testDict7 | 27 + .../equationReaderTester/testDict8 | 27 + .../equationReaderTester/testDict9 | 27 + 25 files changed, 4009 insertions(+) create mode 100644 tutorials/equationReader/equationReaderDemo/expectedOutput create mode 100644 tutorials/equationReader/equationReaderDemo/system/controlDict create mode 100644 tutorials/equationReader/equationReaderDemo/testDict create mode 100644 tutorials/equationReader/equationReaderDemo/testDict2 create mode 100644 tutorials/equationReader/equationReaderDemo/testDict3 create mode 100644 tutorials/equationReader/equationReaderTester/0/U create mode 100644 tutorials/equationReader/equationReaderTester/0/eqns create mode 100644 tutorials/equationReader/equationReaderTester/0/p create mode 100644 tutorials/equationReader/equationReaderTester/constant/polyMesh/blockMeshDict create mode 100644 tutorials/equationReader/equationReaderTester/constant/polyMesh/boundary create mode 100644 tutorials/equationReader/equationReaderTester/constant/transportProperties create mode 100644 tutorials/equationReader/equationReaderTester/expectedOutput create mode 100644 tutorials/equationReader/equationReaderTester/system/controlDict create mode 100644 tutorials/equationReader/equationReaderTester/system/fvSchemes create mode 100644 tutorials/equationReader/equationReaderTester/system/fvSolution create mode 100644 tutorials/equationReader/equationReaderTester/testDict1 create mode 100644 tutorials/equationReader/equationReaderTester/testDict10 create mode 100644 tutorials/equationReader/equationReaderTester/testDict2 create mode 100644 tutorials/equationReader/equationReaderTester/testDict3 create mode 100644 tutorials/equationReader/equationReaderTester/testDict4 create mode 100644 tutorials/equationReader/equationReaderTester/testDict5 create mode 100644 tutorials/equationReader/equationReaderTester/testDict6 create mode 100644 tutorials/equationReader/equationReaderTester/testDict7 create mode 100644 tutorials/equationReader/equationReaderTester/testDict8 create mode 100644 tutorials/equationReader/equationReaderTester/testDict9 diff --git a/tutorials/equationReader/equationReaderDemo/expectedOutput b/tutorials/equationReader/equationReaderDemo/expectedOutput new file mode 100644 index 000000000..cc33b1625 --- /dev/null +++ b/tutorials/equationReader/equationReaderDemo/expectedOutput @@ -0,0 +1,55 @@ +/*---------------------------------------------------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM Extend Project: Open source CFD | +| \\ / O peration | Version: 1.6-ext | +| \\ / A nd | Web: www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +Build : 1.6-ext-1f09367282cf +Exec : equationReaderDemo +Date : Apr 07 2011 +Time : 09:03:55 +Host : Bruce +PID : 4063 +Case : /home/dave/OpenFOAM/dave-1.6-ext/run/tutorials/equationReaderDemo +nProcs : 1 +SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +Reading testDict dictionary + +Begin stand-alone operation... +Reading a scalar... done. Result = 1.94281 +Reading a dimensionedScalar ... done. Result = standAloneDScalar [0 0 0 0 0 0 0] 0.080125 + +Creating the equationReader object + +Creating data sources: dictionary ptrs... scalars... dimensionedScalar ptrs... output dimensionedScalar ptrs... done. +Linking in the data sources: dictionary ptrs... dimensionedScalar ptrs... scalar ptrs... done. + +Reading equation a from testDict with no output variable +Evaluating equation a ... done. Result = a [0 0 0 0 0 0 0] 0.625 + +Reading equation b from testDict, linking an output variable +Output variable before update() = activeOutB [0 0 0 0 0 0 0] 0 +Begining .update() - this evaluates all equations with active output... +Done. Output variable after update() = activeOutB [0 0 0 0 0 0 0] 0.3235 + +Equation c depends on equation a. Reading it will link them. +Reading equation c from testDict... done. +Evaluating c will force an evaluate of a. +Evaluating ... done. Result = c [0 0 0 0 0 0 0] 0.805821 + +Equation d depends on equation e, but equation e is never explicitly +read by equationReaderDemo. Reading equation d will automatically +create equation e on-the-fly. +Reading equation d from testDict ... done. +Again, evaluating d will force an evaluate of e. +Evaluating d ... done. The result is = d [0 0 0 0 0 0 0] 2.04124 + +Equations can draw from any sources added to equationReader. +Equation f is very complex, drawing from numerous sources. +Reading equation f ... done. Evaluating equation f ... done. +The result is: f [0 0 0 0 0 0 0] -86.1698 + +Creating output... diff --git a/tutorials/equationReader/equationReaderDemo/system/controlDict b/tutorials/equationReader/equationReaderDemo/system/controlDict new file mode 100644 index 000000000..89d60a63c --- /dev/null +++ b/tutorials/equationReader/equationReaderDemo/system/controlDict @@ -0,0 +1,47 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application equationReaderDemo; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 1; + +deltaT 1; + +writeControl timeStep; + +writeInterval 1; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderDemo/testDict b/tutorials/equationReader/equationReaderDemo/testDict new file mode 100644 index 000000000..b211daaec --- /dev/null +++ b/tutorials/equationReader/equationReaderDemo/testDict @@ -0,0 +1,47 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +// For demonstrating stand-alone operation +// readSa and readDSa equations can only use constants +standAloneScalar "1 + 1/3 * max(2 * sqrt(2), e_)"; + +standAloneDScalar standAloneDScalar [0 0 0 0 0 0 0] "sin(pi_ / 6) * 1 / 3^(1+2/3)"; + +// This equation is handled 'passively' +a "5 / 2^3"; + +// This one is handled 'actively' +b "cos(5 * theta + 8)"; + +// This one demonstrates variable dependence - a is another equation; DSa is a +// dimensionedScalar explicitly added to the equationReader; and Sb is a scalar +// added to the equationReader +c "a^2 + DSa - cbrt(Sb)"; + +// This one demonstrates on-the-fly equation parsing - it depends on equation +// e, but equation e is never explicitly read in equationTester.C. +d "2 * exp(inv(e)^2)"; + +// equationReader will automatically find and parse equation e when evaluating +// equation d. +//e "1 / (a + SMALL_)"; +e 7; + +// There shouldn't be any limits to the parenthesis depth, variable dependence, +// or number of sources from which the equation draws. +//f "beta"; + +f "DSa*(Sa+(DSb^Sc-(min(dictTwoA * dictThreeB, dictTwoB * dictThreeA)^DSb*cos(pi_/6)^arbitrary+7)))"; diff --git a/tutorials/equationReader/equationReaderDemo/testDict2 b/tutorials/equationReader/equationReaderDemo/testDict2 new file mode 100644 index 000000000..b997b7601 --- /dev/null +++ b/tutorials/equationReader/equationReaderDemo/testDict2 @@ -0,0 +1,25 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dictTwoA dictTwoA [0 0 0 0 0 0 0] 0.9; + +dictTwoB 0.58; + +arbitrary -1; + +beta "gamma * kappa"; + +kappa 0.87; diff --git a/tutorials/equationReader/equationReaderDemo/testDict3 b/tutorials/equationReader/equationReaderDemo/testDict3 new file mode 100644 index 000000000..a9908448c --- /dev/null +++ b/tutorials/equationReader/equationReaderDemo/testDict3 @@ -0,0 +1,34 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dictThreeA "alpha + beta"; + +dictThreeB "dictThreeA + alpha"; + +// Uh oh - a duplicate variable. equationReader doesn't check for these. +// It's anyone's guess which one will be used... +arbitrary 1; + +// This one will still be parsed, even though it is a single constant +alpha "7.0"; + +gamma "epsilon + whatever"; + +epsilon 0.001; + +theta 0.665; + +whatever 8.5;//elseYouWant [0 0 0 0 0 0 0] 8.5; diff --git a/tutorials/equationReader/equationReaderTester/0/U b/tutorials/equationReader/equationReaderTester/0/U new file mode 100644 index 000000000..fab2c48c0 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/0/U @@ -0,0 +1,41 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + object U; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (0 0 0); + +boundaryField +{ + movingWall + { + type fixedValue; + value uniform (1 0 0); + } + + fixedWalls + { + type fixedValue; + value uniform (0 0 0); + } + + frontAndBack + { + type empty; + } +} + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/0/eqns b/tutorials/equationReader/equationReaderTester/0/eqns new file mode 100644 index 000000000..a4cc74d58 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/0/eqns @@ -0,0 +1,201 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5-dev | +| \\ / A nd | Revision: exported | +| \\/ M anipulation | Web: http://www.OpenFOAM.org | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class equationReader; + location "0.1"; + object eqns; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + + +{ + equations + { + Pa + { + Pa "twoa"; + lastResult Pa [ 0 0 0 0 0 0 0 ] 1; + } + twoa + { + twoa "1"; + lastResult twoa [ 0 0 0 0 0 0 0 ] 1; + } + Pb + { + Pb "sign(Sa)*cosh(twoa)"; + lastResult Pb [ 0 0 0 0 0 0 0 ] 1.54308; + } + Pc + { + Pc "DStime * 1e2 / sinh(foura + fivea / sixa^(sevena-eighta/ninea^inv(tena)))"; + lastResult Pc [ 0 0 0 0 0 0 0 ] 2.75721; + } + foura + { + foura "1"; + lastResult foura [ 0 0 0 0 0 0 0 ] 1; + } + ninea + { + ninea "nineb"; + lastResult ninea [ 0 0 0 0 0 0 0 ] 1; + } + Pd + { + Pd "Pa+Pb+Pc+Pe+Pf+Aa+Ab+Ac+Ad+Ae+Af"; + lastResult Pd [ 0 0 0 0 0 0 0 ] 586.516; + } + Pe + { + Pe "Sa+(Sb*Sc^(Sd-Se/Sf^(Sg+Sh^Si)))"; + lastResult Pe [ 0 0 0 0 0 0 0 ] 0.411732; + } + Pf + { + Pf "DSa+(DSb*DSc^(DSd-DSe/DSf^inv(DSg+DSh^DSi)))"; + lastResult Pf [ 0 0 0 0 0 0 0 ] 1.66667; + } + Aa + { + Aa "twob/log10(twoc*10)+cosh(threef) * DStime/Sa+fiveg^max(sixf,sevenj)"; + lastResult Aa [ 0 0 0 0 0 0 0 ] 3.54308; + outputVar Aa [ 0 0 0 0 0 0 0 ] 3.54308; + } + Ab + { + Ab "fourf*log(pos(sevend))*inv(stabilise(Af, SMALL_))"; + lastResult Ab [ 0 0 0 0 0 0 0 ] 0; + outputVar Ab [ 0 0 0 0 0 0 0 ] 0; + } + Ac + { + Ac "DSo + DSm * erf(e_^0) + Ad*Ae+Ad"; + lastResult Ac [ 0 0 0 0 0 0 0 ] 523.335; + outputVar Ac [ 0 0 0 0 0 0 0 ] 523.335; + } + Ad + { + Ad "Ae+Ae*Af+Pb*inv(stabilise(Af, SMALL_))^Pe+DStime^DSa*twog+sevenh"; + lastResult Ad [ 0 0 0 0 0 0 0 ] 38.26; + outputVar Ad [ 0 0 0 0 0 0 0 ] 38.26; + } + Ae + { + Ae "threej + eightg / nineb*DSk - lgamma(ninej)"; + lastResult Ae [ 0 0 0 0 0 0 0 ] 12; + outputVar Ae [ 0 0 0 0 0 0 0 ] 12; + } + Af + { + Af "fiveg^fivec+eighti"; + lastResult Af [ 0 0 0 0 0 0 0 ] 2; + outputVar Af [ 0 0 0 0 0 0 0 ] 2; + } + twob + { + twob "1"; + lastResult twob [ 0 0 0 0 0 0 0 ] 1; + } + twoc + { + twoc "1"; + lastResult twoc [ 0 0 0 0 0 0 0 ] 1; + } + threef + { + threef "1"; + lastResult threef [ 0 0 0 0 0 0 0 ] 1; + } + fourf + { + fourf "1"; + lastResult fourf [ 0 0 0 0 0 0 0 ] 1; + } + twog + { + twog "1"; + lastResult twog [ 0 0 0 0 0 0 0 ] 1; + } + threej + { + threej "1"; + lastResult threej [ 0 0 0 0 0 0 0 ] 1; + } + nineb + { + nineb "ninec"; + lastResult nineb [ 0 0 0 0 0 0 0 ] 1; + } + nu + { + nu [ 0 2 -1 0 0 0 0 ] "0.01 + nuAdd"; + lastResult nu [ 0 2 -1 0 0 0 0 ] 0.01001; + outputVar nu [ 0 2 -1 0 0 0 0 ] 0.01001; + } + nuAdd + { + nuAdd [ 0 2 -1 0 0 0 0 ] "DStime / 10000"; + lastResult nuAdd [ 0 2 -1 0 0 0 0 ] 1e-05; + } + ninec + { + ninec "nined"; + lastResult ninec [ 0 0 0 0 0 0 0 ] 1; + } + nined + { + nined "ninee"; + lastResult nined [ 0 0 0 0 0 0 0 ] 1; + } + ninee + { + ninee "ninef"; + lastResult ninee [ 0 0 0 0 0 0 0 ] 1; + } + ninef + { + ninef "nineg"; + lastResult ninef [ 0 0 0 0 0 0 0 ] 1; + } + nineg + { + nineg "nineh"; + lastResult nineg [ 0 0 0 0 0 0 0 ] 1; + } + nineh + { + nineh "ninei"; + lastResult nineh [ 0 0 0 0 0 0 0 ] 1; + } + ninei + { + ninei "ninej"; + lastResult ninei [ 0 0 0 0 0 0 0 ] 1; + } + } +} + +{ + dataSources + { + dictSources 10 ( "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict1" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict2" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict3" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict4" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict5" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict6" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict7" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict8" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict9" "/home/dave/OpenFOAM/dave-1.5-dev/run/tutorials/equationReaderTester/testDict10" ); + dictLookups 14 ( fivea sixa sevena eighta tena fiveg sixf sevenj sevend sevenh eightg ninej fivec eighti ); + externalDScalars 16 ( DSa [ 0 0 0 0 0 0 0 ] 1 DSb [ 0 0 0 0 0 0 0 ] 2 DSc [ 0 0 0 0 0 0 0 ] 3 DSd [ 0 0 0 0 0 0 0 ] 4 DSe [ 0 0 0 0 0 0 0 ] 5 DSf [ 0 0 0 0 0 0 0 ] 6 DSg [ 0 0 0 0 0 0 0 ] 7 DSh [ 0 0 0 0 0 0 0 ] 8 DSi [ 0 0 0 0 0 0 0 ] 9 DSj [ 0 0 0 0 0 0 0 ] 10 DSk [ 0 0 0 0 0 0 0 ] 11 DSl [ 0 0 0 0 0 0 0 ] 12 DSm [ 0 0 0 0 0 0 0 ] 13 DSn [ 0 0 0 0 0 0 0 ] 14 DSo [ 0 0 0 0 0 0 0 ] 15 DStime [ 0 0 0 0 0 0 0 ] 0.1 ); + externalScalars 12 ( 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.11 0.12 ); + externalScalarNames 12 ( Sa Sb Sc Se Sd Sf Sg Sh Si Sj Sk Sl ); + outputScalars 7 ( Aa [ 0 0 0 0 0 0 0 ] 3.54308 Ab [ 0 0 0 0 0 0 0 ] 0 Ac [ 0 0 0 0 0 0 0 ] 523.335 Ad [ 0 0 0 0 0 0 0 ] 38.26 Ae [ 0 0 0 0 0 0 0 ] 12 Af [ 0 0 0 0 0 0 0 ] 2 nu [ 0 2 -1 0 0 0 0 ] 0.01001 ); + internalScalars 8 ( 100 10 1e-15 2.71828 0 0.01 1 10000 ); + } +} + + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/0/p b/tutorials/equationReader/equationReaderTester/0/p new file mode 100644 index 000000000..8c5250182 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/0/p @@ -0,0 +1,39 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + movingWall + { + type zeroGradient; + } + + fixedWalls + { + type zeroGradient; + } + + frontAndBack + { + type empty; + } +} + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/constant/polyMesh/blockMeshDict b/tutorials/equationReader/equationReaderTester/constant/polyMesh/blockMeshDict new file mode 100644 index 000000000..d89701b42 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/constant/polyMesh/blockMeshDict @@ -0,0 +1,63 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 0.1; + +vertices +( + (0 0 0) + (1 0 0) + (1 1 0) + (0 1 0) + (0 0 0.1) + (1 0 0.1) + (1 1 0.1) + (0 1 0.1) +); + +blocks +( + hex (0 1 2 3 4 5 6 7) (20 20 1) simpleGrading (1 1 1) +); + +edges +( +); + +patches +( + wall movingWall + ( + (3 7 6 2) + ) + wall fixedWalls + ( + (0 4 7 3) + (2 6 5 1) + (1 5 4 0) + ) + empty frontAndBack + ( + (0 3 2 1) + (4 5 6 7) + ) +); + +mergePatchPairs +( +); + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/constant/polyMesh/boundary b/tutorials/equationReader/equationReaderTester/constant/polyMesh/boundary new file mode 100644 index 000000000..ba0763bac --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/constant/polyMesh/boundary @@ -0,0 +1,40 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5-dev | +| \\ / A nd | Revision: exported | +| \\/ M anipulation | Web: http://www.OpenFOAM.org | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class polyBoundaryMesh; + location "constant/polyMesh"; + object boundary; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +3 +( + movingWall + { + type wall; + nFaces 20; + startFace 760; + } + fixedWalls + { + type wall; + nFaces 60; + startFace 780; + } + frontAndBack + { + type empty; + nFaces 800; + startFace 840; + } +) + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/constant/transportProperties b/tutorials/equationReader/equationReaderTester/constant/transportProperties new file mode 100644 index 000000000..cdba3d784 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/constant/transportProperties @@ -0,0 +1,19 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object transportProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +nu nu [0 2 -1 0 0 0 0] 0.01; + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/expectedOutput b/tutorials/equationReader/equationReaderTester/expectedOutput new file mode 100644 index 000000000..e4fb25117 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/expectedOutput @@ -0,0 +1,2968 @@ +/*---------------------------------------------------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM Extend Project: Open source CFD | +| \\ / O peration | Version: 1.6-ext | +| \\ / A nd | Web: www.extend-project.de | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +Build : 1.6-ext-1f09367282cf +Exec : equationReaderTester +Date : Apr 07 2011 +Time : 09:03:35 +Host : Bruce +PID : 4052 +Case : /home/dave/OpenFOAM/dave-1.6-ext/run/tutorials/equationReaderTester +nProcs : 1 +SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +Create time + +Create mesh for time = 0 + +Reading transportProperties + +Reading field p + +Creating dummy field + +Dummy name dummy +Reading field U + +Reading/calculating face flux field phi + +Creating the equationReader object + +Reading testDict dictionaries + +Stand-alone test: +saA = 1 +saB = 2 +saC = 3 +saD = 4 +saE = 5 +saF = 6 +dsaA = dsaA [0.1 0 0 0 0 0 0] 1e-05 +dsaB = dsaB [0 0.1 0 0 0 0 0] 20 +dsaC = dsaC [0 0 0.1 0 0 0 0] 30 +dsaD = dsaD [0 0 0 0.1 0 0 0] 40 +dsaE = dsaE [0 0 0 0 0.1 0 0] 50 +dsaF = dsaF [0 0 0 0 0 0.1 0] 60 +Pa is at index 0 +Pb is at index 2 +Pc is at index 3 +Pd is at index 6 +Pe is at index 7 +Pf is at index 8 +Aa is at index 9 +Ab is at index 10 +Ac is at index 11 +Ad is at index 12 +Ae is at index 13 +Af is at index 14 +nu is at index 22 + +Starting time loop + +Time = 0.01 + +Moving p index to 1... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.15431 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.165 +Ad = Ad [0 0 0 0 0 0 0] 38.17 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010001 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 0.275721 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 581.386 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0 max: 0 velocity magnitude: 0 +DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 4.03296e-06, No Iterations 10 +DILUPBiCG: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0 +DICPCG: Solving for p, Initial residual = 1, Final residual = 6.19585e-07, No Iterations 35 +time step continuity errors : sum local = 1.19691e-08, global = -1.21761e-19, cumulative = -1.21761e-19 +DICPCG: Solving for p, Initial residual = 0.591076, Final residual = 8.95649e-07, No Iterations 34 +time step continuity errors : sum local = 2.68811e-08, global = -7.49095e-19, cumulative = -8.70856e-19 +ExecutionTime = 0.08 s ClockTime = 1 s + +Time = 0.02 + +Moving p index to 2... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.30862 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.295 +Ad = Ad [0 0 0 0 0 0 0] 38.18 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010002 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 0.551441 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 581.956 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.128578 max: 1.29396 velocity magnitude: 0.64698 +DILUPBiCG: Solving for Ux, Initial residual = 0.129651, Final residual = 4.71277e-06, No Iterations 8 +DILUPBiCG: Solving for Uy, Initial residual = 0.262473, Final residual = 4.09295e-06, No Iterations 9 +DICPCG: Solving for p, Initial residual = 0.498818, Final residual = 3.08891e-07, No Iterations 34 +time step continuity errors : sum local = 7.65506e-09, global = -3.00101e-19, cumulative = -1.17096e-18 +DICPCG: Solving for p, Initial residual = 0.392883, Final residual = 6.42421e-07, No Iterations 33 +time step continuity errors : sum local = 1.77384e-08, global = 1.77362e-18, cumulative = 6.02663e-19 +ExecutionTime = 0.1 s ClockTime = 1 s + +Time = 0.03 + +Moving p index to 3... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.46292 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.425 +Ad = Ad [0 0 0 0 0 0 0] 38.19 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010003 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 0.827162 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 582.526 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.185741 max: 1.6461 velocity magnitude: 0.823048 +DILUPBiCG: Solving for Ux, Initial residual = 0.0294556, Final residual = 4.04557e-06, No Iterations 8 +DILUPBiCG: Solving for Uy, Initial residual = 0.0544908, Final residual = 2.35694e-06, No Iterations 8 +DICPCG: Solving for p, Initial residual = 0.229358, Final residual = 4.74348e-07, No Iterations 33 +time step continuity errors : sum local = 8.32516e-09, global = 1.15144e-19, cumulative = 7.17807e-19 +DICPCG: Solving for p, Initial residual = 0.187652, Final residual = 3.33756e-07, No Iterations 33 +time step continuity errors : sum local = 5.79388e-09, global = -9.85668e-19, cumulative = -2.67862e-19 +ExecutionTime = 0.12 s ClockTime = 1 s + +Time = 0.04 + +Moving p index to 4... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.61723 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.555 +Ad = Ad [0 0 0 0 0 0 0] 38.2 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010004 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 1.10288 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 583.096 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.204428 max: 1.65755 velocity magnitude: 0.828775 +DILUPBiCG: Solving for Ux, Initial residual = 0.0120538, Final residual = 1.74808e-06, No Iterations 8 +DILUPBiCG: Solving for Uy, Initial residual = 0.0339442, Final residual = 6.40872e-06, No Iterations 7 +DICPCG: Solving for p, Initial residual = 0.076503, Final residual = 5.37689e-07, No Iterations 32 +time step continuity errors : sum local = 9.36536e-09, global = 2.3793e-18, cumulative = 2.11144e-18 +DICPCG: Solving for p, Initial residual = 0.0620413, Final residual = 3.54663e-07, No Iterations 32 +time step continuity errors : sum local = 6.24356e-09, global = -1.38172e-18, cumulative = 7.29718e-19 +ExecutionTime = 0.14 s ClockTime = 1 s + +Time = 0.05 + +Moving p index to 5... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.77154 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.685 +Ad = Ad [0 0 0 0 0 0 0] 38.21 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010005 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 1.3786 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 583.666 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.216232 max: 1.68509 velocity magnitude: 0.842547 +DILUPBiCG: Solving for Ux, Initial residual = 0.0068727, Final residual = 2.45261e-06, No Iterations 7 +DILUPBiCG: Solving for Uy, Initial residual = 0.0159984, Final residual = 3.46637e-06, No Iterations 7 +DICPCG: Solving for p, Initial residual = 0.0371471, Final residual = 5.5855e-07, No Iterations 30 +time step continuity errors : sum local = 9.54531e-09, global = 7.30897e-19, cumulative = 1.46061e-18 +DICPCG: Solving for p, Initial residual = 0.0305971, Final residual = 3.92368e-07, No Iterations 30 +time step continuity errors : sum local = 6.69483e-09, global = -1.73377e-18, cumulative = -2.73156e-19 +ExecutionTime = 0.17 s ClockTime = 1 s + +Time = 0.06 + +Moving p index to 6... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 2.92585 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.815 +Ad = Ad [0 0 0 0 0 0 0] 38.22 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010006 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 1.65432 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 584.236 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.222634 max: 1.69185 velocity magnitude: 0.845924 +DILUPBiCG: Solving for Ux, Initial residual = 0.00382478, Final residual = 9.07759e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.0100601, Final residual = 4.60975e-06, No Iterations 6 +DICPCG: Solving for p, Initial residual = 0.0152903, Final residual = 7.61663e-07, No Iterations 29 +time step continuity errors : sum local = 1.30035e-08, global = 1.71475e-18, cumulative = 1.44159e-18 +DICPCG: Solving for p, Initial residual = 0.0127069, Final residual = 5.48597e-07, No Iterations 29 +time step continuity errors : sum local = 9.36373e-09, global = -1.8595e-19, cumulative = 1.25564e-18 +ExecutionTime = 0.19 s ClockTime = 1 s + +Time = 0.07 + +Moving p index to 7... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.08016 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 522.945 +Ad = Ad [0 0 0 0 0 0 0] 38.23 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010007 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 1.93004 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 584.806 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.226715 max: 1.69731 velocity magnitude: 0.848657 +DILUPBiCG: Solving for Ux, Initial residual = 0.00239759, Final residual = 5.46641e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.00612175, Final residual = 2.74341e-06, No Iterations 6 +DICPCG: Solving for p, Initial residual = 0.0104525, Final residual = 9.9741e-07, No Iterations 28 +time step continuity errors : sum local = 1.70227e-08, global = -1.31853e-19, cumulative = 1.12379e-18 +DICPCG: Solving for p, Initial residual = 0.00866341, Final residual = 6.43813e-07, No Iterations 28 +time step continuity errors : sum local = 1.09749e-08, global = 3.8497e-19, cumulative = 1.50876e-18 +ExecutionTime = 0.21 s ClockTime = 1 s + +Time = 0.08 + +Moving p index to 8... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.23446 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.075 +Ad = Ad [0 0 0 0 0 0 0] 38.24 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010008 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 2.20576 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 585.376 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.229256 max: 1.69995 velocity magnitude: 0.849973 +DILUPBiCG: Solving for Ux, Initial residual = 0.00150719, Final residual = 3.3473e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.00392926, Final residual = 1.85486e-06, No Iterations 6 +DICPCG: Solving for p, Initial residual = 0.00595047, Final residual = 5.19658e-07, No Iterations 28 +time step continuity errors : sum local = 8.93036e-09, global = -2.94145e-19, cumulative = 1.21461e-18 +DICPCG: Solving for p, Initial residual = 0.00489196, Final residual = 4.48192e-07, No Iterations 28 +time step continuity errors : sum local = 7.71646e-09, global = -2.96627e-19, cumulative = 9.17984e-19 +ExecutionTime = 0.24 s ClockTime = 1 s + +Time = 0.09 + +Moving p index to 9... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.38877 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.205 +Ad = Ad [0 0 0 0 0 0 0] 38.25 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010009 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 2.48149 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 585.946 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.230882 max: 1.70162 velocity magnitude: 0.850809 +DILUPBiCG: Solving for Ux, Initial residual = 0.00097323, Final residual = 2.43569e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.00252728, Final residual = 8.35714e-06, No Iterations 7 +DICPCG: Solving for p, Initial residual = 0.00399438, Final residual = 8.08472e-07, No Iterations 27 +time step continuity errors : sum local = 1.38743e-08, global = 1.1344e-18, cumulative = 2.05238e-18 +DICPCG: Solving for p, Initial residual = 0.00327854, Final residual = 5.10516e-07, No Iterations 25 +time step continuity errors : sum local = 8.75289e-09, global = -7.85822e-20, cumulative = 1.9738e-18 +ExecutionTime = 0.26 s ClockTime = 1 s + +Time = 0.1 + +Moving p index to 10... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.54308 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.335 +Ad = Ad [0 0 0 0 0 0 0] 38.26 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01001 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 2.75721 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 586.516 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.231926 max: 1.70261 velocity magnitude: 0.851303 +DILUPBiCG: Solving for Ux, Initial residual = 0.000632732, Final residual = 1.5485e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.0016409, Final residual = 7.81803e-06, No Iterations 5 +DICPCG: Solving for p, Initial residual = 0.00250763, Final residual = 5.76385e-07, No Iterations 27 +time step continuity errors : sum local = 9.84437e-09, global = 1.20735e-18, cumulative = 3.18115e-18 +DICPCG: Solving for p, Initial residual = 0.00205286, Final residual = 5.73723e-07, No Iterations 27 +time step continuity errors : sum local = 9.81059e-09, global = 8.62088e-19, cumulative = 4.04324e-18 +ExecutionTime = 0.29 s ClockTime = 1 s + +Time = 0.11 + +Moving p index to 11... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.69739 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.465 +Ad = Ad [0 0 0 0 0 0 0] 38.27 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010011 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 3.03293 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 587.086 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.232602 max: 1.70322 velocity magnitude: 0.851609 +DILUPBiCG: Solving for Ux, Initial residual = 0.000414761, Final residual = 1.02601e-06, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.00106888, Final residual = 4.44886e-06, No Iterations 5 +DICPCG: Solving for p, Initial residual = 0.00166151, Final residual = 9.28221e-07, No Iterations 26 +time step continuity errors : sum local = 1.59185e-08, global = 3.98701e-20, cumulative = 4.08311e-18 +DICPCG: Solving for p, Initial residual = 0.00135473, Final residual = 7.66074e-07, No Iterations 26 +time step continuity errors : sum local = 1.31088e-08, global = -6.42719e-19, cumulative = 3.44039e-18 +ExecutionTime = 0.31 s ClockTime = 1 s + +Time = 0.12 + +Moving p index to 12... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 3.8517 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.595 +Ad = Ad [0 0 0 0 0 0 0] 38.28 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010012 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 3.30865 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 587.656 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23304 max: 1.7036 velocity magnitude: 0.851799 +DILUPBiCG: Solving for Ux, Initial residual = 0.000273734, Final residual = 6.65266e-07, No Iterations 5 +DILUPBiCG: Solving for Uy, Initial residual = 0.000698747, Final residual = 6.73378e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.0010952, Final residual = 6.86661e-07, No Iterations 26 +time step continuity errors : sum local = 1.17782e-08, global = -5.29396e-19, cumulative = 2.91099e-18 +DICPCG: Solving for p, Initial residual = 0.000890157, Final residual = 5.29018e-07, No Iterations 26 +time step continuity errors : sum local = 9.05946e-09, global = -1.99896e-18, cumulative = 9.12029e-19 +ExecutionTime = 0.34 s ClockTime = 1 s + +Time = 0.13 + +Moving p index to 13... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.006 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.725 +Ad = Ad [0 0 0 0 0 0 0] 38.29 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010013 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 3.58437 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 588.226 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233324 max: 1.70384 velocity magnitude: 0.85192 +DILUPBiCG: Solving for Ux, Initial residual = 0.000182181, Final residual = 6.35526e-06, No Iterations 4 +DILUPBiCG: Solving for Uy, Initial residual = 0.000457311, Final residual = 4.77931e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.000720863, Final residual = 5.11064e-07, No Iterations 26 +time step continuity errors : sum local = 8.75327e-09, global = -2.93649e-19, cumulative = 6.18379e-19 +DICPCG: Solving for p, Initial residual = 0.000583231, Final residual = 9.69214e-07, No Iterations 25 +time step continuity errors : sum local = 1.6556e-08, global = -1.1157e-18, cumulative = -4.97322e-19 +ExecutionTime = 0.36 s ClockTime = 1 s + +Time = 0.14 + +Moving p index to 14... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.16031 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.855 +Ad = Ad [0 0 0 0 0 0 0] 38.3 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010014 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 3.86009 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 588.796 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233508 max: 1.70399 velocity magnitude: 0.851996 +DILUPBiCG: Solving for Ux, Initial residual = 0.000123966, Final residual = 9.28349e-06, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 0.000300651, Final residual = 2.97326e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.000482526, Final residual = 4.53165e-07, No Iterations 26 +time step continuity errors : sum local = 7.81137e-09, global = -4.00025e-19, cumulative = -8.97346e-19 +DICPCG: Solving for p, Initial residual = 0.000386241, Final residual = 5.82688e-07, No Iterations 16 +time step continuity errors : sum local = 1.35169e-08, global = 1.40207e-18, cumulative = 5.04725e-19 +ExecutionTime = 0.38 s ClockTime = 1 s + +Time = 0.15 + +Moving p index to 15... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.31462 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 523.985 +Ad = Ad [0 0 0 0 0 0 0] 38.31 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010015 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 4.13581 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 589.367 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233626 max: 1.70409 velocity magnitude: 0.852045 +DILUPBiCG: Solving for Ux, Initial residual = 8.60485e-05, Final residual = 5.39774e-06, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 0.000200679, Final residual = 2.0296e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.000350351, Final residual = 8.00625e-07, No Iterations 18 +time step continuity errors : sum local = 1.72981e-08, global = 1.42771e-18, cumulative = 1.93244e-18 +DICPCG: Solving for p, Initial residual = 0.000277604, Final residual = 8.62603e-07, No Iterations 15 +time step continuity errors : sum local = 1.74024e-08, global = 7.84167e-20, cumulative = 2.01086e-18 +ExecutionTime = 0.4 s ClockTime = 1 s + +Time = 0.16 + +Moving p index to 16... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.46893 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.115 +Ad = Ad [0 0 0 0 0 0 0] 38.32 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010016 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 4.41153 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 589.937 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233704 max: 1.70415 velocity magnitude: 0.852076 +DILUPBiCG: Solving for Ux, Initial residual = 6.01462e-05, Final residual = 3.03685e-06, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 0.000135901, Final residual = 1.46879e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.000257985, Final residual = 6.23369e-07, No Iterations 25 +time step continuity errors : sum local = 1.06746e-08, global = 2.12089e-19, cumulative = 2.22294e-18 +DICPCG: Solving for p, Initial residual = 0.000203616, Final residual = 9.12131e-07, No Iterations 17 +time step continuity errors : sum local = 1.89831e-08, global = 2.39188e-18, cumulative = 4.61482e-18 +ExecutionTime = 0.42 s ClockTime = 1 s + +Time = 0.17 + +Moving p index to 17... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.62324 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.245 +Ad = Ad [0 0 0 0 0 0 0] 38.33 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010017 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 4.68725 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 590.507 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233756 max: 1.70419 velocity magnitude: 0.852096 +DILUPBiCG: Solving for Ux, Initial residual = 4.28138e-05, Final residual = 1.72046e-06, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 9.3866e-05, Final residual = 1.12685e-06, No Iterations 4 +DICPCG: Solving for p, Initial residual = 0.000200357, Final residual = 7.45204e-07, No Iterations 24 +time step continuity errors : sum local = 1.27699e-08, global = -6.45201e-20, cumulative = 4.5503e-18 +DICPCG: Solving for p, Initial residual = 0.000157745, Final residual = 5.67908e-07, No Iterations 24 +time step continuity errors : sum local = 9.73194e-09, global = -1.09932e-18, cumulative = 3.45098e-18 +ExecutionTime = 0.44 s ClockTime = 1 s + +Time = 0.18 + +Moving p index to 18... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.77755 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.375 +Ad = Ad [0 0 0 0 0 0 0] 38.34 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010018 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 4.96297 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 591.077 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233791 max: 1.70422 velocity magnitude: 0.852109 +DILUPBiCG: Solving for Ux, Initial residual = 3.15539e-05, Final residual = 1.1144e-06, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 6.58072e-05, Final residual = 6.48167e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000164491, Final residual = 5.87135e-07, No Iterations 24 +time step continuity errors : sum local = 1.00558e-08, global = -2.29625e-19, cumulative = 3.22135e-18 +DICPCG: Solving for p, Initial residual = 0.000130599, Final residual = 9.81702e-07, No Iterations 23 +time step continuity errors : sum local = 1.67585e-08, global = -1.06888e-18, cumulative = 2.15247e-18 +ExecutionTime = 0.46 s ClockTime = 1 s + +Time = 0.19 + +Moving p index to 19... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 4.93185 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.505 +Ad = Ad [0 0 0 0 0 0 0] 38.35 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010019 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 5.23869 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 591.647 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233814 max: 1.70424 velocity magnitude: 0.852118 +DILUPBiCG: Solving for Ux, Initial residual = 2.40382e-05, Final residual = 7.26458e-07, No Iterations 3 +DILUPBiCG: Solving for Uy, Initial residual = 4.64734e-05, Final residual = 8.34047e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000137099, Final residual = 9.18422e-07, No Iterations 18 +time step continuity errors : sum local = 2.14501e-08, global = -1.07848e-18, cumulative = 1.07399e-18 +DICPCG: Solving for p, Initial residual = 0.000108256, Final residual = 9.6813e-07, No Iterations 16 +time step continuity errors : sum local = 1.7575e-08, global = -9.67801e-20, cumulative = 9.7721e-19 +ExecutionTime = 0.48 s ClockTime = 1 s + +Time = 0.2 + +Moving p index to 20... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.08616 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.635 +Ad = Ad [0 0 0 0 0 0 0] 38.36 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01002 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 5.51441 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 592.217 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233829 max: 1.70425 velocity magnitude: 0.852124 +DILUPBiCG: Solving for Ux, Initial residual = 1.91765e-05, Final residual = 3.14833e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 3.44636e-05, Final residual = 6.16516e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000142509, Final residual = 6.7399e-07, No Iterations 24 +time step continuity errors : sum local = 1.15301e-08, global = -2.70819e-19, cumulative = 7.06391e-19 +DICPCG: Solving for p, Initial residual = 0.000114909, Final residual = 5.05187e-07, No Iterations 24 +time step continuity errors : sum local = 8.65079e-09, global = -2.88355e-19, cumulative = 4.18036e-19 +ExecutionTime = 0.51 s ClockTime = 1 s + +Time = 0.21 + +Moving p index to 21... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.24047 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.765 +Ad = Ad [0 0 0 0 0 0 0] 38.37 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010021 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 5.79013 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 592.787 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23384 max: 1.70426 velocity magnitude: 0.852128 +DILUPBiCG: Solving for Ux, Initial residual = 1.49147e-05, Final residual = 1.40638e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.91358e-05, Final residual = 4.51047e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000112109, Final residual = 5.52778e-07, No Iterations 24 +time step continuity errors : sum local = 9.46135e-09, global = 8.36941e-19, cumulative = 1.25498e-18 +DICPCG: Solving for p, Initial residual = 8.7367e-05, Final residual = 4.14559e-07, No Iterations 24 +time step continuity errors : sum local = 7.10321e-09, global = -1.70896e-18, cumulative = -4.53977e-19 +ExecutionTime = 0.54 s ClockTime = 1 s + +Time = 0.22 + +Moving p index to 22... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.39478 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 524.895 +Ad = Ad [0 0 0 0 0 0 0] 38.38 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010022 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 6.06585 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 593.357 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233846 max: 1.70426 velocity magnitude: 0.852131 +DILUPBiCG: Solving for Ux, Initial residual = 1.35049e-05, Final residual = 1.32916e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.43295e-05, Final residual = 3.16036e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000110479, Final residual = 5.8499e-07, No Iterations 24 +time step continuity errors : sum local = 1.0035e-08, global = -1.89987e-18, cumulative = -2.35385e-18 +DICPCG: Solving for p, Initial residual = 8.59878e-05, Final residual = 4.3554e-07, No Iterations 24 +time step continuity errors : sum local = 7.47543e-09, global = -1.3665e-18, cumulative = -3.72035e-18 +ExecutionTime = 0.56 s ClockTime = 1 s + +Time = 0.23 + +Moving p index to 23... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.54909 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.025 +Ad = Ad [0 0 0 0 0 0 0] 38.39 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010023 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 6.34157 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 593.927 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23385 max: 1.70426 velocity magnitude: 0.852132 +DILUPBiCG: Solving for Ux, Initial residual = 1.21645e-05, Final residual = 1.35615e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.27789e-05, Final residual = 2.00117e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000105374, Final residual = 5.56894e-07, No Iterations 24 +time step continuity errors : sum local = 9.54349e-09, global = 1.0439e-18, cumulative = -2.67645e-18 +DICPCG: Solving for p, Initial residual = 8.10236e-05, Final residual = 4.12499e-07, No Iterations 24 +time step continuity errors : sum local = 7.07586e-09, global = 2.106e-19, cumulative = -2.46585e-18 +ExecutionTime = 0.58 s ClockTime = 1 s + +Time = 0.24 + +Moving p index to 24... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.70339 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.155 +Ad = Ad [0 0 0 0 0 0 0] 38.4 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010024 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 6.61729 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 594.497 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233853 max: 1.70427 velocity magnitude: 0.852133 +DILUPBiCG: Solving for Ux, Initial residual = 1.14376e-05, Final residual = 1.38145e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.23563e-05, Final residual = 2.18639e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000102807, Final residual = 5.43654e-07, No Iterations 24 +time step continuity errors : sum local = 9.31635e-09, global = 1.05681e-18, cumulative = -1.40904e-18 +DICPCG: Solving for p, Initial residual = 7.84452e-05, Final residual = 9.84023e-07, No Iterations 23 +time step continuity errors : sum local = 1.68273e-08, global = 6.68858e-19, cumulative = -7.40182e-19 +ExecutionTime = 0.6 s ClockTime = 1 s + +Time = 0.25 + +Moving p index to 25... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 5.8577 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.285 +Ad = Ad [0 0 0 0 0 0 0] 38.41 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010025 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 6.89301 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 595.067 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233855 max: 1.70427 velocity magnitude: 0.852134 +DILUPBiCG: Solving for Ux, Initial residual = 1.09681e-05, Final residual = 1.41882e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.26966e-05, Final residual = 4.42931e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000101097, Final residual = 5.13595e-07, No Iterations 24 +time step continuity errors : sum local = 8.79767e-09, global = -2.51231e-18, cumulative = -3.25249e-18 +DICPCG: Solving for p, Initial residual = 7.63743e-05, Final residual = 9.38463e-07, No Iterations 23 +time step continuity errors : sum local = 1.60407e-08, global = -7.91943e-19, cumulative = -4.04444e-18 +ExecutionTime = 0.62 s ClockTime = 1 s + +Time = 0.26 + +Moving p index to 26... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.01201 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.415 +Ad = Ad [0 0 0 0 0 0 0] 38.42 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010026 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 7.16873 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 595.637 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233856 max: 1.70427 velocity magnitude: 0.852135 +DILUPBiCG: Solving for Ux, Initial residual = 1.06896e-05, Final residual = 1.45615e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.34574e-05, Final residual = 7.22144e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 0.000100651, Final residual = 4.93519e-07, No Iterations 24 +time step continuity errors : sum local = 8.45355e-09, global = -3.77178e-18, cumulative = -7.81622e-18 +DICPCG: Solving for p, Initial residual = 7.56163e-05, Final residual = 9.0138e-07, No Iterations 23 +time step continuity errors : sum local = 1.54055e-08, global = -1.35492e-19, cumulative = -7.95171e-18 +ExecutionTime = 0.64 s ClockTime = 1 s + +Time = 0.27 + +Moving p index to 27... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.16632 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.545 +Ad = Ad [0 0 0 0 0 0 0] 38.43 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010027 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 7.44446 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 596.207 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233857 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.06061e-05, Final residual = 9.96445e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.45016e-05, Final residual = 1.34963e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000139181, Final residual = 7.44327e-07, No Iterations 24 +time step continuity errors : sum local = 1.27595e-08, global = -5.1401e-19, cumulative = -8.46572e-18 +DICPCG: Solving for p, Initial residual = 0.000111618, Final residual = 5.8629e-07, No Iterations 24 +time step continuity errors : sum local = 1.00614e-08, global = 3.91918e-19, cumulative = -8.0738e-18 +ExecutionTime = 0.66 s ClockTime = 1 s + +Time = 0.28 + +Moving p index to 28... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.32063 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.675 +Ad = Ad [0 0 0 0 0 0 0] 38.44 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010028 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 7.72018 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 596.777 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233858 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 9.9815e-06, Final residual = 9.9815e-06, No Iterations 0 +DILUPBiCG: Solving for Uy, Initial residual = 2.56066e-05, Final residual = 7.16223e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 7.69103e-05, Final residual = 8.31322e-07, No Iterations 12 +time step continuity errors : sum local = 2.10947e-08, global = -7.6233e-19, cumulative = -8.83613e-18 +DICPCG: Solving for p, Initial residual = 4.63569e-05, Final residual = 8.92534e-07, No Iterations 11 +time step continuity errors : sum local = 1.64036e-08, global = -2.84385e-19, cumulative = -9.12051e-18 +ExecutionTime = 0.68 s ClockTime = 1 s + +Time = 0.29 + +Moving p index to 29... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.47493 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.805 +Ad = Ad [0 0 0 0 0 0 0] 38.45 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010029 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 7.9959 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 597.347 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233857 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.23244e-05, Final residual = 7.98139e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.43003e-05, Final residual = 6.9121e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000154198, Final residual = 6.68257e-07, No Iterations 24 +time step continuity errors : sum local = 1.14736e-08, global = -1.99963e-18, cumulative = -1.11201e-17 +DICPCG: Solving for p, Initial residual = 0.000118736, Final residual = 5.0479e-07, No Iterations 24 +time step continuity errors : sum local = 8.68328e-09, global = -9.52085e-19, cumulative = -1.20722e-17 +ExecutionTime = 0.71 s ClockTime = 1 s + +Time = 0.3 + +Moving p index to 30... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.62924 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 525.935 +Ad = Ad [0 0 0 0 0 0 0] 38.46 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01003 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 8.27162 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 597.917 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.00212e-05, Final residual = 1.31625e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.51639e-05, Final residual = 7.92632e-06, No Iterations 2 +DICPCG: Solving for p, Initial residual = 9.51516e-05, Final residual = 4.33567e-07, No Iterations 24 +time step continuity errors : sum local = 7.42502e-09, global = 2.4104e-19, cumulative = -1.18312e-17 +DICPCG: Solving for p, Initial residual = 7.2261e-05, Final residual = 7.80594e-07, No Iterations 23 +time step continuity errors : sum local = 1.33379e-08, global = -6.5612e-19, cumulative = -1.24873e-17 +ExecutionTime = 0.74 s ClockTime = 1 s + +Time = 0.31 + +Moving p index to 31... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.78355 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.065 +Ad = Ad [0 0 0 0 0 0 0] 38.47 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010031 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 8.54734 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 598.487 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.06935e-05, Final residual = 8.41088e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.55959e-05, Final residual = 9.53806e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000140761, Final residual = 6.02978e-07, No Iterations 24 +time step continuity errors : sum local = 1.03339e-08, global = 1.8514e-18, cumulative = -1.06359e-17 +DICPCG: Solving for p, Initial residual = 0.000109401, Final residual = 4.55244e-07, No Iterations 24 +time step continuity errors : sum local = 7.81418e-09, global = -1.46461e-18, cumulative = -1.21005e-17 +ExecutionTime = 0.76 s ClockTime = 1 s + +Time = 0.32 + +Moving p index to 32... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 6.93786 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.195 +Ad = Ad [0 0 0 0 0 0 0] 38.48 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010032 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 8.82306 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 599.057 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.01805e-05, Final residual = 1.39733e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.76532e-05, Final residual = 1.5962e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.00010114, Final residual = 5.63167e-07, No Iterations 24 +time step continuity errors : sum local = 9.65348e-09, global = -2.30188e-18, cumulative = -1.44024e-17 +DICPCG: Solving for p, Initial residual = 7.93505e-05, Final residual = 9.5282e-07, No Iterations 23 +time step continuity errors : sum local = 1.63039e-08, global = -6.64557e-19, cumulative = -1.5067e-17 +ExecutionTime = 0.78 s ClockTime = 1 s + +Time = 0.33 + +Moving p index to 33... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.09217 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.325 +Ad = Ad [0 0 0 0 0 0 0] 38.49 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010033 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 9.09878 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 599.627 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.055e-05, Final residual = 7.45796e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.48608e-05, Final residual = 8.75212e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129175, Final residual = 5.08864e-07, No Iterations 24 +time step continuity errors : sum local = 8.72486e-09, global = -1.46295e-18, cumulative = -1.65299e-17 +DICPCG: Solving for p, Initial residual = 9.73802e-05, Final residual = 3.88992e-07, No Iterations 24 +time step continuity errors : sum local = 6.67783e-09, global = 2.05108e-18, cumulative = -1.44788e-17 +ExecutionTime = 0.8 s ClockTime = 1 s + +Time = 0.34 + +Moving p index to 34... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.24647 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.455 +Ad = Ad [0 0 0 0 0 0 0] 38.5 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010034 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 9.3745 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 600.197 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03495e-05, Final residual = 1.39857e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73181e-05, Final residual = 1.53211e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.00010263, Final residual = 5.93843e-07, No Iterations 24 +time step continuity errors : sum local = 1.01783e-08, global = 2.14289e-18, cumulative = -1.23359e-17 +DICPCG: Solving for p, Initial residual = 8.05563e-05, Final residual = 9.95561e-07, No Iterations 23 +time step continuity errors : sum local = 1.70372e-08, global = 2.77767e-19, cumulative = -1.20582e-17 +ExecutionTime = 0.82 s ClockTime = 1 s + +Time = 0.35 + +Moving p index to 35... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.40078 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.585 +Ad = Ad [0 0 0 0 0 0 0] 38.51 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010035 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 9.65022 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 600.767 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.05243e-05, Final residual = 7.63003e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.4688e-05, Final residual = 9.0102e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129707, Final residual = 5.03164e-07, No Iterations 24 +time step continuity errors : sum local = 8.62691e-09, global = 5.60001e-19, cumulative = -1.14982e-17 +DICPCG: Solving for p, Initial residual = 9.79878e-05, Final residual = 3.87893e-07, No Iterations 24 +time step continuity errors : sum local = 6.65755e-09, global = -5.16822e-19, cumulative = -1.2015e-17 +ExecutionTime = 0.85 s ClockTime = 1 s + +Time = 0.36 + +Moving p index to 36... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.55509 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.715 +Ad = Ad [0 0 0 0 0 0 0] 38.52 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010036 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 9.92594 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 601.337 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03409e-05, Final residual = 1.3919e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73717e-05, Final residual = 1.54077e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102473, Final residual = 5.97345e-07, No Iterations 24 +time step continuity errors : sum local = 1.02378e-08, global = -2.37053e-18, cumulative = -1.43855e-17 +DICPCG: Solving for p, Initial residual = 8.04774e-05, Final residual = 9.9965e-07, No Iterations 23 +time step continuity errors : sum local = 1.71068e-08, global = 6.41727e-19, cumulative = -1.37438e-17 +ExecutionTime = 0.87 s ClockTime = 1 s + +Time = 0.37 + +Moving p index to 37... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.7094 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.845 +Ad = Ad [0 0 0 0 0 0 0] 38.53 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010037 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 10.2017 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 601.907 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.05225e-05, Final residual = 7.63311e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46815e-05, Final residual = 8.9904e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129707, Final residual = 5.04275e-07, No Iterations 24 +time step continuity errors : sum local = 8.64563e-09, global = 2.43853e-18, cumulative = -1.13053e-17 +DICPCG: Solving for p, Initial residual = 9.79724e-05, Final residual = 3.89046e-07, No Iterations 24 +time step continuity errors : sum local = 6.67693e-09, global = -7.7755e-21, cumulative = -1.1313e-17 +ExecutionTime = 0.89 s ClockTime = 1 s + +Time = 0.38 + +Moving p index to 38... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 7.86371 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 526.975 +Ad = Ad [0 0 0 0 0 0 0] 38.54 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010038 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 10.4774 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 602.477 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03385e-05, Final residual = 1.3904e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73616e-05, Final residual = 1.54145e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102427, Final residual = 5.97549e-07, No Iterations 24 +time step continuity errors : sum local = 1.02411e-08, global = 1.33838e-19, cumulative = -1.11792e-17 +DICPCG: Solving for p, Initial residual = 8.04475e-05, Final residual = 9.99861e-07, No Iterations 23 +time step continuity errors : sum local = 1.71102e-08, global = -3.50344e-18, cumulative = -1.46826e-17 +ExecutionTime = 0.92 s ClockTime = 1 s + +Time = 0.39 + +Moving p index to 39... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.01801 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.105 +Ad = Ad [0 0 0 0 0 0 0] 38.55 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010039 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 10.7531 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 603.047 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.05211e-05, Final residual = 7.6389e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46745e-05, Final residual = 8.97564e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129743, Final residual = 5.0487e-07, No Iterations 24 +time step continuity errors : sum local = 8.65579e-09, global = -6.73656e-19, cumulative = -1.53563e-17 +DICPCG: Solving for p, Initial residual = 9.80089e-05, Final residual = 3.89566e-07, No Iterations 24 +time step continuity errors : sum local = 6.68582e-09, global = 1.03083e-18, cumulative = -1.43255e-17 +ExecutionTime = 0.94 s ClockTime = 1 s + +Time = 0.4 + +Moving p index to 40... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.17232 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.235 +Ad = Ad [0 0 0 0 0 0 0] 38.56 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01004 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 11.0288 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 603.617 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03359e-05, Final residual = 1.38976e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73554e-05, Final residual = 1.54186e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.0001024, Final residual = 5.9741e-07, No Iterations 24 +time step continuity errors : sum local = 1.02386e-08, global = -9.501e-19, cumulative = -1.52756e-17 +DICPCG: Solving for p, Initial residual = 8.04296e-05, Final residual = 9.99607e-07, No Iterations 23 +time step continuity errors : sum local = 1.71058e-08, global = -1.65304e-18, cumulative = -1.69286e-17 +ExecutionTime = 0.97 s ClockTime = 1 s + +Time = 0.41 + +Moving p index to 41... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.32663 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.365 +Ad = Ad [0 0 0 0 0 0 0] 38.57 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010041 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 11.3045 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 604.187 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852136 +DILUPBiCG: Solving for Ux, Initial residual = 1.05195e-05, Final residual = 7.64105e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46701e-05, Final residual = 8.96839e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.00012975, Final residual = 5.05058e-07, No Iterations 24 +time step continuity errors : sum local = 8.65902e-09, global = 1.61399e-18, cumulative = -1.53146e-17 +DICPCG: Solving for p, Initial residual = 9.80187e-05, Final residual = 3.89728e-07, No Iterations 24 +time step continuity errors : sum local = 6.6886e-09, global = -5.82335e-19, cumulative = -1.58969e-17 +ExecutionTime = 0.99 s ClockTime = 1 s + +Time = 0.42 + +Moving p index to 42... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.48094 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.495 +Ad = Ad [0 0 0 0 0 0 0] 38.58 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010042 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 11.5803 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 604.757 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03339e-05, Final residual = 1.38949e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.7351e-05, Final residual = 1.54221e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102382, Final residual = 5.97283e-07, No Iterations 24 +time step continuity errors : sum local = 1.02365e-08, global = -2.92326e-19, cumulative = -1.61893e-17 +DICPCG: Solving for p, Initial residual = 8.04176e-05, Final residual = 9.99385e-07, No Iterations 23 +time step continuity errors : sum local = 1.7102e-08, global = -2.4832e-19, cumulative = -1.64376e-17 +ExecutionTime = 1 s ClockTime = 1 s + +Time = 0.43 + +Moving p index to 43... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.63525 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.625 +Ad = Ad [0 0 0 0 0 0 0] 38.59 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010043 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 11.856 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 605.327 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05178e-05, Final residual = 7.64237e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46665e-05, Final residual = 8.96536e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129747, Final residual = 5.05093e-07, No Iterations 24 +time step continuity errors : sum local = 8.65963e-09, global = -1.48032e-18, cumulative = -1.79179e-17 +DICPCG: Solving for p, Initial residual = 9.80202e-05, Final residual = 3.89769e-07, No Iterations 24 +time step continuity errors : sum local = 6.68931e-09, global = -3.19953e-19, cumulative = -1.82379e-17 +ExecutionTime = 1.02 s ClockTime = 2 s + +Time = 0.44 + +Moving p index to 44... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.78955 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.755 +Ad = Ad [0 0 0 0 0 0 0] 38.6 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010044 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 12.1317 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 605.897 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.03319e-05, Final residual = 1.38934e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73477e-05, Final residual = 1.5426e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102366, Final residual = 5.97179e-07, No Iterations 24 +time step continuity errors : sum local = 1.02347e-08, global = 4.58258e-20, cumulative = -1.8192e-17 +DICPCG: Solving for p, Initial residual = 8.04074e-05, Final residual = 9.992e-07, No Iterations 23 +time step continuity errors : sum local = 1.70988e-08, global = 9.82691e-19, cumulative = -1.72093e-17 +ExecutionTime = 1.05 s ClockTime = 2 s + +Time = 0.45 + +Moving p index to 45... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 8.94386 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 527.885 +Ad = Ad [0 0 0 0 0 0 0] 38.61 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010045 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 12.4074 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 606.467 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.0516e-05, Final residual = 7.6434e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46633e-05, Final residual = 8.96398e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129742, Final residual = 5.05075e-07, No Iterations 24 +time step continuity errors : sum local = 8.65933e-09, global = 8.73337e-19, cumulative = -1.6336e-17 +DICPCG: Solving for p, Initial residual = 9.80188e-05, Final residual = 3.8977e-07, No Iterations 24 +time step continuity errors : sum local = 6.68934e-09, global = 7.73745e-19, cumulative = -1.55623e-17 +ExecutionTime = 1.07 s ClockTime = 2 s + +Time = 0.46 + +Moving p index to 46... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.09817 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.015 +Ad = Ad [0 0 0 0 0 0 0] 38.62 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010046 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 12.6831 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 607.037 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.033e-05, Final residual = 1.38923e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73449e-05, Final residual = 1.54301e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102351, Final residual = 5.97085e-07, No Iterations 24 +time step continuity errors : sum local = 1.02331e-08, global = -1.99599e-18, cumulative = -1.75583e-17 +DICPCG: Solving for p, Initial residual = 8.03978e-05, Final residual = 9.99031e-07, No Iterations 23 +time step continuity errors : sum local = 1.70959e-08, global = 1.13638e-18, cumulative = -1.64219e-17 +ExecutionTime = 1.09 s ClockTime = 2 s + +Time = 0.47 + +Moving p index to 47... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.25248 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.145 +Ad = Ad [0 0 0 0 0 0 0] 38.63 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010047 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 12.9589 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 607.607 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05143e-05, Final residual = 7.64432e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46603e-05, Final residual = 8.96323e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129735, Final residual = 5.05039e-07, No Iterations 24 +time step continuity errors : sum local = 8.65872e-09, global = -1.1023e-18, cumulative = -1.75242e-17 +DICPCG: Solving for p, Initial residual = 9.80165e-05, Final residual = 3.89758e-07, No Iterations 24 +time step continuity errors : sum local = 6.68915e-09, global = 1.26691e-18, cumulative = -1.62573e-17 +ExecutionTime = 1.11 s ClockTime = 2 s + +Time = 0.48 + +Moving p index to 48... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.40679 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.275 +Ad = Ad [0 0 0 0 0 0 0] 38.64 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010048 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 13.2346 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 608.177 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03282e-05, Final residual = 1.38913e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73424e-05, Final residual = 1.54344e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102336, Final residual = 5.96995e-07, No Iterations 24 +time step continuity errors : sum local = 1.02316e-08, global = 1.34417e-18, cumulative = -1.49131e-17 +DICPCG: Solving for p, Initial residual = 8.03884e-05, Final residual = 9.98869e-07, No Iterations 23 +time step continuity errors : sum local = 1.70932e-08, global = 3.07364e-18, cumulative = -1.18395e-17 +ExecutionTime = 1.14 s ClockTime = 2 s + +Time = 0.49 + +Moving p index to 49... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.5611 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.405 +Ad = Ad [0 0 0 0 0 0 0] 38.65 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010049 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 13.5103 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 608.747 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05125e-05, Final residual = 7.64521e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46573e-05, Final residual = 8.96269e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129728, Final residual = 5.04996e-07, No Iterations 24 +time step continuity errors : sum local = 8.65801e-09, global = -1.9131e-18, cumulative = -1.37526e-17 +DICPCG: Solving for p, Initial residual = 9.80139e-05, Final residual = 3.89741e-07, No Iterations 24 +time step continuity errors : sum local = 6.68887e-09, global = 9.35541e-19, cumulative = -1.2817e-17 +ExecutionTime = 1.16 s ClockTime = 2 s + +Time = 0.5 + +Moving p index to 50... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.7154 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.535 +Ad = Ad [0 0 0 0 0 0 0] 38.66 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01005 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 13.786 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 609.318 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03263e-05, Final residual = 1.38904e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73399e-05, Final residual = 1.54386e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102322, Final residual = 5.96907e-07, No Iterations 24 +time step continuity errors : sum local = 1.02301e-08, global = -1.81649e-19, cumulative = -1.29987e-17 +DICPCG: Solving for p, Initial residual = 8.03791e-05, Final residual = 9.98709e-07, No Iterations 23 +time step continuity errors : sum local = 1.70904e-08, global = 7.84167e-20, cumulative = -1.29203e-17 +ExecutionTime = 1.19 s ClockTime = 2 s + +Time = 0.51 + +Moving p index to 51... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 9.86971 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.665 +Ad = Ad [0 0 0 0 0 0 0] 38.67 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010051 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 14.0617 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 609.888 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233859 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05107e-05, Final residual = 7.64608e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46545e-05, Final residual = 8.96224e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.00012972, Final residual = 5.04952e-07, No Iterations 24 +time step continuity errors : sum local = 8.65726e-09, global = -2.01286e-18, cumulative = -1.49331e-17 +DICPCG: Solving for p, Initial residual = 9.80111e-05, Final residual = 3.89723e-07, No Iterations 24 +time step continuity errors : sum local = 6.68858e-09, global = -1.96075e-18, cumulative = -1.68939e-17 +ExecutionTime = 1.21 s ClockTime = 2 s + +Time = 0.52 + +Moving p index to 52... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.024 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.795 +Ad = Ad [0 0 0 0 0 0 0] 38.68 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010052 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 14.3375 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 610.458 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03244e-05, Final residual = 1.38895e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73375e-05, Final residual = 1.54429e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102307, Final residual = 5.9682e-07, No Iterations 24 +time step continuity errors : sum local = 1.02286e-08, global = 6.51984e-19, cumulative = -1.62419e-17 +DICPCG: Solving for p, Initial residual = 8.03699e-05, Final residual = 9.9855e-07, No Iterations 23 +time step continuity errors : sum local = 1.70877e-08, global = -8.48026e-19, cumulative = -1.70899e-17 +ExecutionTime = 1.23 s ClockTime = 2 s + +Time = 0.53 + +Moving p index to 53... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.1783 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 528.925 +Ad = Ad [0 0 0 0 0 0 0] 38.69 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010053 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 14.6132 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 611.028 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05089e-05, Final residual = 7.64694e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46516e-05, Final residual = 8.9618e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129713, Final residual = 5.04906e-07, No Iterations 24 +time step continuity errors : sum local = 8.65649e-09, global = -3.83481e-19, cumulative = -1.74734e-17 +DICPCG: Solving for p, Initial residual = 9.80083e-05, Final residual = 3.89705e-07, No Iterations 24 +time step continuity errors : sum local = 6.68827e-09, global = 1.63285e-18, cumulative = -1.58405e-17 +ExecutionTime = 1.26 s ClockTime = 2 s + +Time = 0.54 + +Moving p index to 54... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.3326 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.055 +Ad = Ad [0 0 0 0 0 0 0] 38.7 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010054 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 14.8889 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 611.598 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03226e-05, Final residual = 1.38886e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73351e-05, Final residual = 1.54472e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102292, Final residual = 5.96733e-07, No Iterations 24 +time step continuity errors : sum local = 1.02271e-08, global = 4.62559e-19, cumulative = -1.5378e-17 +DICPCG: Solving for p, Initial residual = 8.03606e-05, Final residual = 9.98391e-07, No Iterations 23 +time step continuity errors : sum local = 1.7085e-08, global = -9.5953e-19, cumulative = -1.63375e-17 +ExecutionTime = 1.28 s ClockTime = 2 s + +Time = 0.55 + +Moving p index to 55... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.4869 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.185 +Ad = Ad [0 0 0 0 0 0 0] 38.71 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010055 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 15.1646 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 612.168 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05071e-05, Final residual = 7.64781e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46488e-05, Final residual = 8.96138e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129706, Final residual = 5.0486e-07, No Iterations 24 +time step continuity errors : sum local = 8.65572e-09, global = 5.47594e-20, cumulative = -1.62827e-17 +DICPCG: Solving for p, Initial residual = 9.80055e-05, Final residual = 3.89686e-07, No Iterations 24 +time step continuity errors : sum local = 6.68797e-09, global = -2.7694e-19, cumulative = -1.65597e-17 +ExecutionTime = 1.3 s ClockTime = 2 s + +Time = 0.56 + +Moving p index to 56... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.6413 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.315 +Ad = Ad [0 0 0 0 0 0 0] 38.72 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010056 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 15.4404 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 612.738 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03207e-05, Final residual = 1.38877e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73327e-05, Final residual = 1.54514e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102278, Final residual = 5.96646e-07, No Iterations 24 +time step continuity errors : sum local = 1.02256e-08, global = 9.53905e-19, cumulative = -1.56058e-17 +DICPCG: Solving for p, Initial residual = 8.03513e-05, Final residual = 9.98232e-07, No Iterations 23 +time step continuity errors : sum local = 1.70823e-08, global = 4.40391e-19, cumulative = -1.51654e-17 +ExecutionTime = 1.32 s ClockTime = 2 s + +Time = 0.57 + +Moving p index to 57... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.7956 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.445 +Ad = Ad [0 0 0 0 0 0 0] 38.73 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010057 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 15.7161 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 613.308 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05053e-05, Final residual = 7.64867e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46459e-05, Final residual = 8.96096e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129699, Final residual = 5.04815e-07, No Iterations 24 +time step continuity errors : sum local = 8.65495e-09, global = 1.44442e-18, cumulative = -1.3721e-17 +DICPCG: Solving for p, Initial residual = 9.80028e-05, Final residual = 3.89668e-07, No Iterations 24 +time step continuity errors : sum local = 6.68766e-09, global = -4.67357e-19, cumulative = -1.41883e-17 +ExecutionTime = 1.34 s ClockTime = 2 s + +Time = 0.58 + +Moving p index to 58... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 10.9499 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.575 +Ad = Ad [0 0 0 0 0 0 0] 38.74 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010058 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 15.9918 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 613.878 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03188e-05, Final residual = 1.38868e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73303e-05, Final residual = 1.54557e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102263, Final residual = 5.96559e-07, No Iterations 24 +time step continuity errors : sum local = 1.02241e-08, global = -1.31803e-18, cumulative = -1.55063e-17 +DICPCG: Solving for p, Initial residual = 8.03421e-05, Final residual = 9.98074e-07, No Iterations 23 +time step continuity errors : sum local = 1.70796e-08, global = 1.56321e-18, cumulative = -1.39431e-17 +ExecutionTime = 1.37 s ClockTime = 2 s + +Time = 0.59 + +Moving p index to 59... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.1042 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.705 +Ad = Ad [0 0 0 0 0 0 0] 38.75 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010059 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 16.2675 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 614.448 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05035e-05, Final residual = 7.64953e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46431e-05, Final residual = 8.96054e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129691, Final residual = 5.04769e-07, No Iterations 24 +time step continuity errors : sum local = 8.65418e-09, global = 3.92414e-19, cumulative = -1.35507e-17 +DICPCG: Solving for p, Initial residual = 9.8e-05, Final residual = 3.89649e-07, No Iterations 24 +time step continuity errors : sum local = 6.68735e-09, global = 4.30134e-20, cumulative = -1.35077e-17 +ExecutionTime = 1.39 s ClockTime = 2 s + +Time = 0.6 + +Moving p index to 60... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.2585 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.835 +Ad = Ad [0 0 0 0 0 0 0] 38.76 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01006 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 16.5432 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 615.018 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.0317e-05, Final residual = 1.38859e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73279e-05, Final residual = 1.546e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102249, Final residual = 5.96472e-07, No Iterations 24 +time step continuity errors : sum local = 1.02226e-08, global = -4.64876e-19, cumulative = -1.39726e-17 +DICPCG: Solving for p, Initial residual = 8.03329e-05, Final residual = 9.97916e-07, No Iterations 23 +time step continuity errors : sum local = 1.70769e-08, global = 3.52544e-19, cumulative = -1.362e-17 +ExecutionTime = 1.42 s ClockTime = 2 s + +Time = 0.61 + +Moving p index to 61... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.4128 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 529.965 +Ad = Ad [0 0 0 0 0 0 0] 38.77 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010061 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 16.819 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 615.588 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70427 velocity magnitude: 0.852137 +DILUPBiCG: Solving for Ux, Initial residual = 1.05017e-05, Final residual = 7.6504e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46403e-05, Final residual = 8.96012e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129684, Final residual = 5.04723e-07, No Iterations 24 +time step continuity errors : sum local = 8.65341e-09, global = -2.28881e-18, cumulative = -1.59089e-17 +DICPCG: Solving for p, Initial residual = 9.79972e-05, Final residual = 3.8963e-07, No Iterations 24 +time step continuity errors : sum local = 6.68704e-09, global = 1.26906e-18, cumulative = -1.46398e-17 +ExecutionTime = 1.44 s ClockTime = 2 s + +Time = 0.62 + +Moving p index to 62... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.5671 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.095 +Ad = Ad [0 0 0 0 0 0 0] 38.78 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010062 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 17.0947 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 616.158 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03151e-05, Final residual = 1.3885e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73255e-05, Final residual = 1.54643e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102234, Final residual = 5.96385e-07, No Iterations 24 +time step continuity errors : sum local = 1.02211e-08, global = 7.44628e-19, cumulative = -1.38952e-17 +DICPCG: Solving for p, Initial residual = 8.03237e-05, Final residual = 9.97757e-07, No Iterations 23 +time step continuity errors : sum local = 1.70742e-08, global = -6.29981e-19, cumulative = -1.45251e-17 +ExecutionTime = 1.46 s ClockTime = 2 s + +Time = 0.63 + +Moving p index to 63... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.7214 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.225 +Ad = Ad [0 0 0 0 0 0 0] 38.79 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010063 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 17.3704 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 616.728 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04999e-05, Final residual = 7.65126e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46374e-05, Final residual = 8.9597e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129677, Final residual = 5.04677e-07, No Iterations 24 +time step continuity errors : sum local = 8.65264e-09, global = 8.96498e-19, cumulative = -1.36286e-17 +DICPCG: Solving for p, Initial residual = 9.79944e-05, Final residual = 3.89612e-07, No Iterations 24 +time step continuity errors : sum local = 6.68673e-09, global = 6.96652e-19, cumulative = -1.2932e-17 +ExecutionTime = 1.49 s ClockTime = 2 s + +Time = 0.64 + +Moving p index to 64... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 11.8757 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.355 +Ad = Ad [0 0 0 0 0 0 0] 38.8 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010064 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 17.6461 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 617.298 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.03133e-05, Final residual = 1.38841e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73231e-05, Final residual = 1.54686e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102219, Final residual = 5.96299e-07, No Iterations 24 +time step continuity errors : sum local = 1.02197e-08, global = -4.62725e-19, cumulative = -1.33947e-17 +DICPCG: Solving for p, Initial residual = 8.03145e-05, Final residual = 9.97599e-07, No Iterations 23 +time step continuity errors : sum local = 1.70715e-08, global = -5.66122e-19, cumulative = -1.39608e-17 +ExecutionTime = 1.51 s ClockTime = 2 s + +Time = 0.65 + +Moving p index to 65... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.03 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.485 +Ad = Ad [0 0 0 0 0 0 0] 38.81 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010065 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 17.9218 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 617.868 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04982e-05, Final residual = 7.65212e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46346e-05, Final residual = 8.95928e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.00012967, Final residual = 5.04631e-07, No Iterations 24 +time step continuity errors : sum local = 8.65187e-09, global = -3.24867e-18, cumulative = -1.72095e-17 +DICPCG: Solving for p, Initial residual = 9.79916e-05, Final residual = 3.89593e-07, No Iterations 24 +time step continuity errors : sum local = 6.68642e-09, global = -6.32462e-19, cumulative = -1.7842e-17 +ExecutionTime = 1.53 s ClockTime = 2 s + +Time = 0.66 + +Moving p index to 66... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.1843 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.615 +Ad = Ad [0 0 0 0 0 0 0] 38.82 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010066 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 18.1976 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 618.438 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03114e-05, Final residual = 1.38832e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73208e-05, Final residual = 1.54729e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102205, Final residual = 5.96212e-07, No Iterations 24 +time step continuity errors : sum local = 1.02182e-08, global = -7.28581e-19, cumulative = -1.85706e-17 +DICPCG: Solving for p, Initial residual = 8.03052e-05, Final residual = 9.9744e-07, No Iterations 23 +time step continuity errors : sum local = 1.70688e-08, global = -3.99363e-19, cumulative = -1.89699e-17 +ExecutionTime = 1.55 s ClockTime = 2 s + +Time = 0.67 + +Moving p index to 67... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.3386 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.745 +Ad = Ad [0 0 0 0 0 0 0] 38.83 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010067 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 18.4733 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 619.008 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04964e-05, Final residual = 7.65298e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46318e-05, Final residual = 8.95886e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129662, Final residual = 5.04586e-07, No Iterations 24 +time step continuity errors : sum local = 8.6511e-09, global = -1.57776e-18, cumulative = -2.05477e-17 +DICPCG: Solving for p, Initial residual = 9.79888e-05, Final residual = 3.89574e-07, No Iterations 24 +time step continuity errors : sum local = 6.68611e-09, global = -7.40988e-19, cumulative = -2.12887e-17 +ExecutionTime = 1.58 s ClockTime = 2 s + +Time = 0.68 + +Moving p index to 68... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.4929 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 530.875 +Ad = Ad [0 0 0 0 0 0 0] 38.84 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010068 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 18.749 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 619.578 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03095e-05, Final residual = 1.38823e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73184e-05, Final residual = 1.54771e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.00010219, Final residual = 5.96125e-07, No Iterations 24 +time step continuity errors : sum local = 1.02167e-08, global = -1.8954e-18, cumulative = -2.31841e-17 +DICPCG: Solving for p, Initial residual = 8.0296e-05, Final residual = 9.97282e-07, No Iterations 23 +time step continuity errors : sum local = 1.70661e-08, global = -4.19215e-19, cumulative = -2.36033e-17 +ExecutionTime = 1.6 s ClockTime = 2 s + +Time = 0.69 + +Moving p index to 69... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.6473 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.005 +Ad = Ad [0 0 0 0 0 0 0] 38.85 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010069 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 19.0247 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 620.148 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04946e-05, Final residual = 7.65384e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46289e-05, Final residual = 8.95844e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129655, Final residual = 5.0454e-07, No Iterations 24 +time step continuity errors : sum local = 8.65033e-09, global = -1.62061e-18, cumulative = -2.52239e-17 +DICPCG: Solving for p, Initial residual = 9.79861e-05, Final residual = 3.89556e-07, No Iterations 24 +time step continuity errors : sum local = 6.6858e-09, global = -2.39469e-18, cumulative = -2.76186e-17 +ExecutionTime = 1.62 s ClockTime = 2 s + +Time = 0.7 + +Moving p index to 70... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.8016 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.135 +Ad = Ad [0 0 0 0 0 0 0] 38.86 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01007 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 19.3004 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 620.718 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03077e-05, Final residual = 1.38814e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.7316e-05, Final residual = 1.54814e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102175, Final residual = 5.96038e-07, No Iterations 24 +time step continuity errors : sum local = 1.02152e-08, global = 1.58819e-18, cumulative = -2.60304e-17 +DICPCG: Solving for p, Initial residual = 8.02868e-05, Final residual = 9.97124e-07, No Iterations 23 +time step continuity errors : sum local = 1.70634e-08, global = -1.92799e-18, cumulative = -2.79584e-17 +ExecutionTime = 1.65 s ClockTime = 2 s + +Time = 0.71 + +Moving p index to 71... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 12.9559 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.265 +Ad = Ad [0 0 0 0 0 0 0] 38.87 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010071 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 19.5762 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 621.288 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04928e-05, Final residual = 7.65471e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46261e-05, Final residual = 8.95802e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129648, Final residual = 5.04494e-07, No Iterations 24 +time step continuity errors : sum local = 8.64956e-09, global = -6.56947e-19, cumulative = -2.86153e-17 +DICPCG: Solving for p, Initial residual = 9.79833e-05, Final residual = 3.89537e-07, No Iterations 24 +time step continuity errors : sum local = 6.6855e-09, global = -1.01015e-18, cumulative = -2.96255e-17 +ExecutionTime = 1.67 s ClockTime = 2 s + +Time = 0.72 + +Moving p index to 72... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.1102 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.395 +Ad = Ad [0 0 0 0 0 0 0] 38.88 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010072 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 19.8519 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 621.858 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03058e-05, Final residual = 1.38805e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73136e-05, Final residual = 1.54857e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102161, Final residual = 5.95951e-07, No Iterations 24 +time step continuity errors : sum local = 1.02137e-08, global = 6.35109e-19, cumulative = -2.89904e-17 +DICPCG: Solving for p, Initial residual = 8.02776e-05, Final residual = 9.96966e-07, No Iterations 23 +time step continuity errors : sum local = 1.70607e-08, global = 6.85567e-19, cumulative = -2.83048e-17 +ExecutionTime = 1.7 s ClockTime = 2 s + +Time = 0.73 + +Moving p index to 73... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.2645 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.525 +Ad = Ad [0 0 0 0 0 0 0] 38.89 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010073 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 20.1276 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 622.428 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.0491e-05, Final residual = 7.65557e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46233e-05, Final residual = 8.9576e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129641, Final residual = 5.04448e-07, No Iterations 24 +time step continuity errors : sum local = 8.64879e-09, global = 1.18518e-18, cumulative = -2.71196e-17 +DICPCG: Solving for p, Initial residual = 9.79805e-05, Final residual = 3.89518e-07, No Iterations 24 +time step continuity errors : sum local = 6.68519e-09, global = -1.1607e-18, cumulative = -2.82803e-17 +ExecutionTime = 1.72 s ClockTime = 2 s + +Time = 0.74 + +Moving p index to 74... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.4188 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.655 +Ad = Ad [0 0 0 0 0 0 0] 38.9 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010074 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 20.4033 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 622.998 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.0304e-05, Final residual = 1.38796e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73112e-05, Final residual = 1.549e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102146, Final residual = 5.95865e-07, No Iterations 24 +time step continuity errors : sum local = 1.02123e-08, global = -1.4249e-18, cumulative = -2.97052e-17 +DICPCG: Solving for p, Initial residual = 8.02684e-05, Final residual = 9.96808e-07, No Iterations 23 +time step continuity errors : sum local = 1.7058e-08, global = 1.45815e-18, cumulative = -2.82471e-17 +ExecutionTime = 1.74 s ClockTime = 2 s + +Time = 0.75 + +Moving p index to 75... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.5731 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.785 +Ad = Ad [0 0 0 0 0 0 0] 38.91 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010075 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 20.679 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 623.568 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04892e-05, Final residual = 7.65643e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46205e-05, Final residual = 8.95718e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129633, Final residual = 5.04403e-07, No Iterations 24 +time step continuity errors : sum local = 8.64802e-09, global = 1.32564e-18, cumulative = -2.69214e-17 +DICPCG: Solving for p, Initial residual = 9.79777e-05, Final residual = 3.895e-07, No Iterations 24 +time step continuity errors : sum local = 6.68488e-09, global = 1.66594e-19, cumulative = -2.67548e-17 +ExecutionTime = 1.76 s ClockTime = 2 s + +Time = 0.76 + +Moving p index to 76... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.7274 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 531.915 +Ad = Ad [0 0 0 0 0 0 0] 38.92 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010076 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 20.9548 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 624.138 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03021e-05, Final residual = 1.38787e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73089e-05, Final residual = 1.54943e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102132, Final residual = 5.95778e-07, No Iterations 24 +time step continuity errors : sum local = 1.02108e-08, global = 3.19176e-18, cumulative = -2.35631e-17 +DICPCG: Solving for p, Initial residual = 8.02592e-05, Final residual = 9.96649e-07, No Iterations 23 +time step continuity errors : sum local = 1.70553e-08, global = -2.05108e-18, cumulative = -2.56142e-17 +ExecutionTime = 1.78 s ClockTime = 2 s + +Time = 0.77 + +Moving p index to 77... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 13.8817 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.045 +Ad = Ad [0 0 0 0 0 0 0] 38.93 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010077 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 21.2305 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 624.708 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04874e-05, Final residual = 7.65729e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46176e-05, Final residual = 8.95675e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129626, Final residual = 5.04357e-07, No Iterations 24 +time step continuity errors : sum local = 8.64725e-09, global = 8.25195e-19, cumulative = -2.4789e-17 +DICPCG: Solving for p, Initial residual = 9.7975e-05, Final residual = 3.89481e-07, No Iterations 24 +time step continuity errors : sum local = 6.68457e-09, global = -1.06011e-18, cumulative = -2.58491e-17 +ExecutionTime = 1.81 s ClockTime = 2 s + +Time = 0.78 + +Moving p index to 78... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.036 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.175 +Ad = Ad [0 0 0 0 0 0 0] 38.94 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010078 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 21.5062 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 625.278 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.03003e-05, Final residual = 1.38778e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73065e-05, Final residual = 1.54985e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102117, Final residual = 5.95691e-07, No Iterations 24 +time step continuity errors : sum local = 1.02093e-08, global = 1.19081e-18, cumulative = -2.46583e-17 +DICPCG: Solving for p, Initial residual = 8.025e-05, Final residual = 9.96491e-07, No Iterations 23 +time step continuity errors : sum local = 1.70527e-08, global = -1.43284e-18, cumulative = -2.60911e-17 +ExecutionTime = 1.83 s ClockTime = 2 s + +Time = 0.79 + +Moving p index to 79... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.1903 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.305 +Ad = Ad [0 0 0 0 0 0 0] 38.95 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010079 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 21.7819 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 625.848 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04857e-05, Final residual = 7.65816e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46148e-05, Final residual = 8.95633e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129619, Final residual = 5.04311e-07, No Iterations 24 +time step continuity errors : sum local = 8.64648e-09, global = -1.25765e-18, cumulative = -2.73488e-17 +DICPCG: Solving for p, Initial residual = 9.79722e-05, Final residual = 3.89463e-07, No Iterations 24 +time step continuity errors : sum local = 6.68426e-09, global = 4.97963e-19, cumulative = -2.68508e-17 +ExecutionTime = 1.85 s ClockTime = 2 s + +Time = 0.8 + +Moving p index to 80... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.3446 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.435 +Ad = Ad [0 0 0 0 0 0 0] 38.96 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01008 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 22.0576 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 626.418 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.02984e-05, Final residual = 1.38769e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73041e-05, Final residual = 1.55028e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102103, Final residual = 5.95605e-07, No Iterations 24 +time step continuity errors : sum local = 1.02078e-08, global = 1.99185e-19, cumulative = -2.66516e-17 +DICPCG: Solving for p, Initial residual = 8.02408e-05, Final residual = 9.96333e-07, No Iterations 23 +time step continuity errors : sum local = 1.705e-08, global = -9.87654e-19, cumulative = -2.76393e-17 +ExecutionTime = 1.88 s ClockTime = 2 s + +Time = 0.81 + +Moving p index to 81... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.499 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.565 +Ad = Ad [0 0 0 0 0 0 0] 38.97 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010081 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 22.3334 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 626.988 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852138 +DILUPBiCG: Solving for Ux, Initial residual = 1.04839e-05, Final residual = 7.65902e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.4612e-05, Final residual = 8.95591e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129612, Final residual = 5.04266e-07, No Iterations 24 +time step continuity errors : sum local = 8.64571e-09, global = -8.09314e-19, cumulative = -2.84486e-17 +DICPCG: Solving for p, Initial residual = 9.79694e-05, Final residual = 3.89444e-07, No Iterations 24 +time step continuity errors : sum local = 6.68396e-09, global = 5.35682e-19, cumulative = -2.79129e-17 +ExecutionTime = 1.91 s ClockTime = 2 s + +Time = 0.82 + +Moving p index to 82... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.6533 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.695 +Ad = Ad [0 0 0 0 0 0 0] 38.98 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010082 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 22.6091 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 627.558 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.85214 +DILUPBiCG: Solving for Ux, Initial residual = 1.02966e-05, Final residual = 1.38759e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.73017e-05, Final residual = 1.55071e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102088, Final residual = 5.95518e-07, No Iterations 24 +time step continuity errors : sum local = 1.02063e-08, global = 1.72434e-18, cumulative = -2.61886e-17 +DICPCG: Solving for p, Initial residual = 8.02316e-05, Final residual = 9.96175e-07, No Iterations 23 +time step continuity errors : sum local = 1.70473e-08, global = -7.68616e-19, cumulative = -2.69572e-17 +ExecutionTime = 1.93 s ClockTime = 2 s + +Time = 0.83 + +Moving p index to 83... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.8076 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.825 +Ad = Ad [0 0 0 0 0 0 0] 38.99 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010083 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 22.8848 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 628.128 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.23386 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04821e-05, Final residual = 7.65988e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46091e-05, Final residual = 8.95548e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129604, Final residual = 5.0422e-07, No Iterations 24 +time step continuity errors : sum local = 8.64495e-09, global = 7.84167e-20, cumulative = -2.68788e-17 +DICPCG: Solving for p, Initial residual = 9.79667e-05, Final residual = 3.89425e-07, No Iterations 24 +time step continuity errors : sum local = 6.68365e-09, global = -1.4747e-18, cumulative = -2.83535e-17 +ExecutionTime = 1.95 s ClockTime = 2 s + +Time = 0.84 + +Moving p index to 84... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 14.9619 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 532.955 +Ad = Ad [0 0 0 0 0 0 0] 39 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010084 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 23.1605 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 628.699 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02947e-05, Final residual = 1.3875e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72993e-05, Final residual = 1.55114e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102073, Final residual = 5.95431e-07, No Iterations 24 +time step continuity errors : sum local = 1.02049e-08, global = -2.20328e-18, cumulative = -3.05567e-17 +DICPCG: Solving for p, Initial residual = 8.02224e-05, Final residual = 9.96017e-07, No Iterations 23 +time step continuity errors : sum local = 1.70446e-08, global = 1.40687e-18, cumulative = -2.91499e-17 +ExecutionTime = 1.97 s ClockTime = 2 s + +Time = 0.85 + +Moving p index to 85... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.1162 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.085 +Ad = Ad [0 0 0 0 0 0 0] 39.01 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010085 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 23.4362 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 629.269 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04803e-05, Final residual = 7.66074e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46063e-05, Final residual = 8.95506e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129597, Final residual = 5.04174e-07, No Iterations 24 +time step continuity errors : sum local = 8.64418e-09, global = -4.12929e-19, cumulative = -2.95628e-17 +DICPCG: Solving for p, Initial residual = 9.79639e-05, Final residual = 3.89407e-07, No Iterations 24 +time step continuity errors : sum local = 6.68334e-09, global = -8.28835e-20, cumulative = -2.96457e-17 +ExecutionTime = 1.99 s ClockTime = 2 s + +Time = 0.86 + +Moving p index to 86... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.2705 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.215 +Ad = Ad [0 0 0 0 0 0 0] 39.02 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010086 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 23.712 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 629.839 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02928e-05, Final residual = 1.38741e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.7297e-05, Final residual = 1.55157e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102059, Final residual = 5.95345e-07, No Iterations 24 +time step continuity errors : sum local = 1.02034e-08, global = -5.69762e-19, cumulative = -3.02154e-17 +DICPCG: Solving for p, Initial residual = 8.02132e-05, Final residual = 9.95859e-07, No Iterations 23 +time step continuity errors : sum local = 1.70419e-08, global = 5.99044e-19, cumulative = -2.96164e-17 +ExecutionTime = 2.02 s ClockTime = 2 s + +Time = 0.87 + +Moving p index to 87... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.4248 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.345 +Ad = Ad [0 0 0 0 0 0 0] 39.03 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010087 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 23.9877 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 630.409 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04785e-05, Final residual = 7.6616e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46035e-05, Final residual = 8.95464e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.00012959, Final residual = 5.04129e-07, No Iterations 24 +time step continuity errors : sum local = 8.64341e-09, global = -1.59596e-18, cumulative = -3.12124e-17 +DICPCG: Solving for p, Initial residual = 9.79612e-05, Final residual = 3.89388e-07, No Iterations 24 +time step continuity errors : sum local = 6.68303e-09, global = 5.12025e-19, cumulative = -3.07003e-17 +ExecutionTime = 2.04 s ClockTime = 3 s + +Time = 0.88 + +Moving p index to 88... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.5791 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.475 +Ad = Ad [0 0 0 0 0 0 0] 39.04 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010088 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 24.2634 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 630.979 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.0291e-05, Final residual = 1.38732e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72946e-05, Final residual = 1.55199e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102044, Final residual = 5.95258e-07, No Iterations 24 +time step continuity errors : sum local = 1.02019e-08, global = -1.00751e-19, cumulative = -3.08011e-17 +DICPCG: Solving for p, Initial residual = 8.0204e-05, Final residual = 9.95701e-07, No Iterations 23 +time step continuity errors : sum local = 1.70392e-08, global = 1.83005e-18, cumulative = -2.8971e-17 +ExecutionTime = 2.06 s ClockTime = 3 s + +Time = 0.89 + +Moving p index to 89... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.7334 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.605 +Ad = Ad [0 0 0 0 0 0 0] 39.05 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010089 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 24.5391 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 631.549 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04767e-05, Final residual = 7.66247e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.46007e-05, Final residual = 8.95421e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129583, Final residual = 5.04083e-07, No Iterations 24 +time step continuity errors : sum local = 8.64264e-09, global = -1.43003e-18, cumulative = -3.04011e-17 +DICPCG: Solving for p, Initial residual = 9.79584e-05, Final residual = 3.8937e-07, No Iterations 24 +time step continuity errors : sum local = 6.68272e-09, global = -2.12023e-18, cumulative = -3.25213e-17 +ExecutionTime = 2.08 s ClockTime = 3 s + +Time = 0.9 + +Moving p index to 90... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 15.8877 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.735 +Ad = Ad [0 0 0 0 0 0 0] 39.06 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.01009 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 24.8149 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 632.119 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02891e-05, Final residual = 1.38723e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72922e-05, Final residual = 1.55242e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.00010203, Final residual = 5.95171e-07, No Iterations 24 +time step continuity errors : sum local = 1.02004e-08, global = -8.84587e-19, cumulative = -3.34059e-17 +DICPCG: Solving for p, Initial residual = 8.01949e-05, Final residual = 9.95543e-07, No Iterations 23 +time step continuity errors : sum local = 1.70365e-08, global = -5.92096e-19, cumulative = -3.3998e-17 +ExecutionTime = 2.11 s ClockTime = 3 s + +Time = 0.91 + +Moving p index to 91... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.042 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.865 +Ad = Ad [0 0 0 0 0 0 0] 39.07 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010091 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 25.0906 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 632.689 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.0475e-05, Final residual = 7.66333e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.45978e-05, Final residual = 8.95379e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129575, Final residual = 5.04038e-07, No Iterations 24 +time step continuity errors : sum local = 8.64188e-09, global = -1.53359e-19, cumulative = -3.41513e-17 +DICPCG: Solving for p, Initial residual = 9.79557e-05, Final residual = 3.89351e-07, No Iterations 24 +time step continuity errors : sum local = 6.68242e-09, global = 1.02984e-18, cumulative = -3.31215e-17 +ExecutionTime = 2.14 s ClockTime = 3 s + +Time = 0.92 + +Moving p index to 92... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.1963 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 533.995 +Ad = Ad [0 0 0 0 0 0 0] 39.08 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010092 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 25.3663 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 633.259 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02873e-05, Final residual = 1.38714e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72898e-05, Final residual = 1.55285e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102015, Final residual = 5.95085e-07, No Iterations 24 +time step continuity errors : sum local = 1.0199e-08, global = 1.6001e-18, cumulative = -3.15214e-17 +DICPCG: Solving for p, Initial residual = 8.01857e-05, Final residual = 9.95385e-07, No Iterations 23 +time step continuity errors : sum local = 1.70338e-08, global = 7.32717e-19, cumulative = -3.07887e-17 +ExecutionTime = 2.16 s ClockTime = 3 s + +Time = 0.93 + +Moving p index to 93... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.3506 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.125 +Ad = Ad [0 0 0 0 0 0 0] 39.09 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010093 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 25.642 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 633.829 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04732e-05, Final residual = 7.66419e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.4595e-05, Final residual = 8.95336e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129568, Final residual = 5.03992e-07, No Iterations 24 +time step continuity errors : sum local = 8.64111e-09, global = 1.73361e-18, cumulative = -2.90551e-17 +DICPCG: Solving for p, Initial residual = 9.79529e-05, Final residual = 3.89333e-07, No Iterations 24 +time step continuity errors : sum local = 6.68211e-09, global = -8.3082e-19, cumulative = -2.98859e-17 +ExecutionTime = 2.18 s ClockTime = 3 s + +Time = 0.94 + +Moving p index to 94... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.505 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.255 +Ad = Ad [0 0 0 0 0 0 0] 39.1 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010094 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 25.9177 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 634.399 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02854e-05, Final residual = 1.38705e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72875e-05, Final residual = 1.55328e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000102001, Final residual = 5.94998e-07, No Iterations 24 +time step continuity errors : sum local = 1.01975e-08, global = 1.02934e-18, cumulative = -2.88565e-17 +DICPCG: Solving for p, Initial residual = 8.01765e-05, Final residual = 9.95228e-07, No Iterations 23 +time step continuity errors : sum local = 1.70311e-08, global = -2.87743e-18, cumulative = -3.1734e-17 +ExecutionTime = 2.2 s ClockTime = 3 s + +Time = 0.95 + +Moving p index to 95... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.6593 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.385 +Ad = Ad [0 0 0 0 0 0 0] 39.11 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010095 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 26.1935 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 634.969 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04714e-05, Final residual = 7.66505e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.45922e-05, Final residual = 8.95294e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129561, Final residual = 5.03946e-07, No Iterations 24 +time step continuity errors : sum local = 8.64034e-09, global = 2.29956e-18, cumulative = -2.94344e-17 +DICPCG: Solving for p, Initial residual = 9.79501e-05, Final residual = 3.89314e-07, No Iterations 24 +time step continuity errors : sum local = 6.6818e-09, global = 1.00221e-18, cumulative = -2.84322e-17 +ExecutionTime = 2.22 s ClockTime = 3 s + +Time = 0.96 + +Moving p index to 96... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.8136 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.515 +Ad = Ad [0 0 0 0 0 0 0] 39.12 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010096 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 26.4692 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 635.539 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02836e-05, Final residual = 1.38696e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72851e-05, Final residual = 1.55371e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000101986, Final residual = 5.94912e-07, No Iterations 24 +time step continuity errors : sum local = 1.0196e-08, global = 2.3707e-19, cumulative = -2.81951e-17 +DICPCG: Solving for p, Initial residual = 8.01673e-05, Final residual = 9.9507e-07, No Iterations 23 +time step continuity errors : sum local = 1.70284e-08, global = -2.51033e-18, cumulative = -3.07055e-17 +ExecutionTime = 2.25 s ClockTime = 3 s + +Time = 0.97 + +Moving p index to 97... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 16.9679 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.645 +Ad = Ad [0 0 0 0 0 0 0] 39.13 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010097 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 26.7449 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 636.109 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04696e-05, Final residual = 7.66591e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.45894e-05, Final residual = 8.95251e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129554, Final residual = 5.03901e-07, No Iterations 24 +time step continuity errors : sum local = 8.63958e-09, global = 1.47569e-19, cumulative = -3.05579e-17 +DICPCG: Solving for p, Initial residual = 9.79474e-05, Final residual = 3.89296e-07, No Iterations 24 +time step continuity errors : sum local = 6.6815e-09, global = 8.44551e-19, cumulative = -2.97133e-17 +ExecutionTime = 2.27 s ClockTime = 3 s + +Time = 0.98 + +Moving p index to 98... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 17.1222 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.775 +Ad = Ad [0 0 0 0 0 0 0] 39.14 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010098 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 27.0206 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 636.679 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02817e-05, Final residual = 1.38687e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72827e-05, Final residual = 1.55414e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000101972, Final residual = 5.94825e-07, No Iterations 24 +time step continuity errors : sum local = 1.01945e-08, global = -6.46359e-19, cumulative = -3.03597e-17 +DICPCG: Solving for p, Initial residual = 8.01581e-05, Final residual = 9.94912e-07, No Iterations 23 +time step continuity errors : sum local = 1.70257e-08, global = -3.45761e-20, cumulative = -3.03943e-17 +ExecutionTime = 2.29 s ClockTime = 3 s + +Time = 0.99 + +Moving p index to 99... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 17.2765 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 534.905 +Ad = Ad [0 0 0 0 0 0 0] 39.15 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.010099 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 27.2963 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 637.249 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233861 max: 1.70428 velocity magnitude: 0.852139 +DILUPBiCG: Solving for Ux, Initial residual = 1.04678e-05, Final residual = 7.66678e-06, No Iterations 1 +DILUPBiCG: Solving for Uy, Initial residual = 2.45866e-05, Final residual = 8.95208e-06, No Iterations 1 +DICPCG: Solving for p, Initial residual = 0.000129547, Final residual = 5.03855e-07, No Iterations 24 +time step continuity errors : sum local = 8.63881e-09, global = -2.39816e-18, cumulative = -3.27924e-17 +DICPCG: Solving for p, Initial residual = 9.79447e-05, Final residual = 3.89277e-07, No Iterations 24 +time step continuity errors : sum local = 6.68119e-09, global = 8.78135e-19, cumulative = -3.19143e-17 +ExecutionTime = 2.31 s ClockTime = 3 s + +Time = 1 + +Moving p index to 100... +Passive reread... +Active reread... +Updating active equations... +Aa = Aa [0 0 0 0 0 0 0] 17.4308 +Ab = Ab [0 0 0 0 0 0 0] 0 +Ac = Ac [0 0 0 0 0 0 0] 535.035 +Ad = Ad [0 0 0 0 0 0 0] 39.16 +Ae = Ae [0 0 0 0 0 0 0] 12 +Af = Af [0 0 0 0 0 0 0] 2 +nu = nu [0 2 -1 0 0 0 0] 0.0101 +Evaluating passive equations: Pa, Pb, Pc, Pd, Pe, Pf. +Pa = Pa [0 0 0 0 0 0 0] 1 +Pb = Pb [0 0 0 0 0 0 0] 1.54308 +Pc = Pc [0 0 0 0 0 0 0] 27.5721 +Pd = Pe [0 0 0 0 0 0 0] 0.411732 +Pe = Pd [0 0 0 0 0 0 0] 637.819 +Pf = Pf [0 0 0 0 0 0 0] 1.66667 +Courant Number mean: 0.233862 max: 1.70428 velocity magnitude: 0.852141 +DILUPBiCG: Solving for Ux, Initial residual = 1.02799e-05, Final residual = 1.38678e-06, No Iterations 2 +DILUPBiCG: Solving for Uy, Initial residual = 2.72803e-05, Final residual = 1.55456e-06, No Iterations 3 +DICPCG: Solving for p, Initial residual = 0.000101957, Final residual = 5.94739e-07, No Iterations 24 +time step continuity errors : sum local = 1.0193e-08, global = -9.94602e-19, cumulative = -3.29089e-17 +DICPCG: Solving for p, Initial residual = 8.0149e-05, Final residual = 9.94754e-07, No Iterations 23 +time step continuity errors : sum local = 1.7023e-08, global = 1.6585e-18, cumulative = -3.12504e-17 +ExecutionTime = 2.34 s ClockTime = 3 s + +End + diff --git a/tutorials/equationReader/equationReaderTester/system/controlDict b/tutorials/equationReader/equationReaderTester/system/controlDict new file mode 100644 index 000000000..a845aa339 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/system/controlDict @@ -0,0 +1,33 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application equationReaderTester; +startFrom startTime; +startTime 0; +stopAt endTime; +endTime "2/2"; +deltaT 0.01; +writeControl timeStep; +writeInterval 10; +purgeWrite 0; +writeFormat ascii; +writePrecision 6; +writeCompression uncompressed; +timeFormat general; +timePrecision 6; +runTimeModifiable yes; + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/system/fvSchemes b/tutorials/equationReader/equationReaderTester/system/fvSchemes new file mode 100644 index 000000000..712cc1a60 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/system/fvSchemes @@ -0,0 +1,58 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; + grad(p) Gauss linear; +} + +divSchemes +{ + default none; + div(phi,U) Gauss linear; +} + +laplacianSchemes +{ + default none; + laplacian(nu,U) Gauss linear corrected; + laplacian((1|A(U)),p) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; + interpolate(HbyA) linear; +} + +snGradSchemes +{ + default corrected; +} + +fluxRequired +{ + default no; + p; +} + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/system/fvSolution b/tutorials/equationReader/equationReaderTester/system/fvSolution new file mode 100644 index 000000000..f29239916 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/system/fvSolution @@ -0,0 +1,44 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + p + { + solver PCG; + preconditioner DIC; + tolerance 1e-06; + relTol 0; + }; + + U + { + solver PBiCG; + preconditioner DILU; + tolerance 1e-05; + relTol 0; + }; +} + +PISO +{ + nCorrectors 2; + nNonOrthogonalCorrectors 0; + pRefCell 0; + pRefValue 0; +} + +// ************************************************************************* // diff --git a/tutorials/equationReader/equationReaderTester/testDict1 b/tutorials/equationReader/equationReaderTester/testDict1 new file mode 100644 index 000000000..1d576a857 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict1 @@ -0,0 +1,49 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +// Scalar standalone test +saA "5^0"; +saB "1+1"; +saC "1+2"; +saD "2^2"; +saE "sqrt(25)"; +saF "sin(pi_) + cbrt(6^2 * (3 + 2 + 1))"; + +// Dimensioned scalar standalone test +dsaA dsaA [0.1 0 0 0 0 0 0] "1e-5"; +dsaB dsaB [0 0.1 0 0 0 0 0] "sqrt(400)"; +dsaC dsaC [0 0 0.1 0 0 0 0] "(5+5*2)*2"; +dsaD dsaD [0 0 0 0.1 0 0 0] "5*2^3"; +dsaE dsaE [0 0 0 0 0.1 0 0] "log10(1e5)*2*(2^2+1)"; +dsaF dsaF [0 0 0 0 0 0.1 0] "60"; + +// Passive equations +Pa "twoa"; +//Pb "2* DSa + sign(Sa) * cosh(twoa^threea) / (Pc + Aa)"; +Pb "sign(Sa)*cosh(twoa)";//^threea)"; +Pc "DStime * 1e2 / sinh(foura + fivea / sixa^(sevena-eighta/ninea^inv(tena)))"; +Pd "Pa+Pb+Pc+Pe+Pf+Aa+Ab+Ac+Ad+Ae+Af"; +Pe "Sa+(Sb*Sc^(Sd-Se/Sf^(Sg+Sh^Si)))"; +Pf "DSa+(DSb*DSc^(DSd-DSe/DSf^inv(DSg+DSh^DSi)))"; + +// Active equations +Aa "twob/log10(twoc*10)+cosh(threef) * DStime/Sa+fiveg^max(sixf,sevenj) * dummy"; +Ab "fourf*log(pos(sevend))*inv(stabilise(Af, SMALL_))"; +Ac "DSo + DSm * erf(e_^0) + Ad*Ae+Ad + pdimless"; +Ad "Ae+Ae*Af+Pb*inv(stabilise(Af, SMALL_))^Pe+DStime^DSa*twog+sevenh"; +Ae "threej + eightg / nineb*DSk - lgamma(ninej)"; +Af "fiveg^fivec+eighti"; +nu [0 2 -1 0 0 0 0] "0.01 + nuAdd"; diff --git a/tutorials/equationReader/equationReaderTester/testDict10 b/tutorials/equationReader/equationReaderTester/testDict10 new file mode 100644 index 000000000..7a055725b --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict10 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict10; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +tena 1; +tenb 1; +tenc 1; +tend 1; +tene 1; +tenf 1; +teng 1; +tenh 1; +teni 1; +tenj 1; diff --git a/tutorials/equationReader/equationReaderTester/testDict2 b/tutorials/equationReader/equationReaderTester/testDict2 new file mode 100644 index 000000000..c5c15ee06 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict2 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict2; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +twoa "1"; +twob "1"; +twoc "1"; +twod "1"; +twoe "1"; +twof "1"; +twog "1"; +twoh "1"; +twoi "1"; +twoj "1"; diff --git a/tutorials/equationReader/equationReaderTester/testDict3 b/tutorials/equationReader/equationReaderTester/testDict3 new file mode 100644 index 000000000..2a95edde0 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict3 @@ -0,0 +1,30 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict3; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +threea "1"; +threeb "1"; +threec "1"; +threed "1"; +threee "1"; +threef "1"; +threeg "1"; +threeh "1"; +threei "1"; +threej "1"; +nuAdd [0 2 -1 0 0 0 0] "DStime / 10000"; + +pdimless [0 0 0 0 0 0 0] "p"; diff --git a/tutorials/equationReader/equationReaderTester/testDict4 b/tutorials/equationReader/equationReaderTester/testDict4 new file mode 100644 index 000000000..5f0afb8cc --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict4 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict4; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +foura "1"; +fourb "1"; +fourc "1"; +fourd "1"; +foure "1"; +fourf "1"; +fourg "1"; +fourh "1"; +fouri "1"; +fourj "1"; diff --git a/tutorials/equationReader/equationReaderTester/testDict5 b/tutorials/equationReader/equationReaderTester/testDict5 new file mode 100644 index 000000000..3286bf8d8 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict5 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict5; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +fivea 1.0; +fiveb 1.0; +fivec 1.0; +fived 1.0; +fivee 1.0; +fivef 1.0; +fiveg 1.0; +fiveh 1.0; +fivei 1.0; +fivej 1.0; diff --git a/tutorials/equationReader/equationReaderTester/testDict6 b/tutorials/equationReader/equationReaderTester/testDict6 new file mode 100644 index 000000000..428b19759 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict6 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict6; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +sixa 1.0; +sixb 1.0; +sixc 1.0; +sixd 1.0; +sixe 1.0; +sixf 1.0; +sixg 1.0; +sixh 1.0; +sixi 1.0; +sixj 1.0; diff --git a/tutorials/equationReader/equationReaderTester/testDict7 b/tutorials/equationReader/equationReaderTester/testDict7 new file mode 100644 index 000000000..624a19e24 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict7 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict7; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +sevena sevena [0 0 0 0 0 0 0] 1.0; +sevenb sevenb [0 0 0 0 0 0 0] 1.0; +sevenc sevenc [0 0 0 0 0 0 0] 1.0; +sevend sevend [0 0 0 0 0 0 0] 1.0; +sevene sevene [0 0 0 0 0 0 0] 1.0; +sevenf sevenf [0 0 0 0 0 0 0] 1.0; +seveng seveng [0 0 0 0 0 0 0] 1.0; +sevenh sevenh [0 0 0 0 0 0 0] 1.0; +seveni seveni [0 0 0 0 0 0 0] 1.0; +sevenj sevenj [0 0 0 0 0 0 0] 1.0; diff --git a/tutorials/equationReader/equationReaderTester/testDict8 b/tutorials/equationReader/equationReaderTester/testDict8 new file mode 100644 index 000000000..fc3e698d3 --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict8 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict8; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +eighta eighta [0 0 0 0 0 0 0] 1.0; +eightb eightb [0 0 0 0 0 0 0] 1.0; +eightc eightc [0 0 0 0 0 0 0] 1.0; +eightd eightd [0 0 0 0 0 0 0] 1.0; +eighte eighte [0 0 0 0 0 0 0] 1.0; +eightf eightf [0 0 0 0 0 0 0] 1.0; +eightg eightg [0 0 0 0 0 0 0] 1.0; +eighth eighth [0 0 0 0 0 0 0] 1.0; +eighti eighti [0 0 0 0 0 0 0] 1.0; +eightj eightj [0 0 0 0 0 0 0] 1.0; diff --git a/tutorials/equationReader/equationReaderTester/testDict9 b/tutorials/equationReader/equationReaderTester/testDict9 new file mode 100644 index 000000000..3a53afcad --- /dev/null +++ b/tutorials/equationReader/equationReaderTester/testDict9 @@ -0,0 +1,27 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object testDict9; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +// Sa..Sl; DSa..DSo; Aa..Af; Pa..Pf + +ninea "nineb"; +nineb "ninec"; +ninec "nined"; +nined "ninee"; +ninee "ninef"; +ninef "nineg"; +nineg "nineh"; +nineh "ninei"; +ninei "ninej"; +ninej 1.0;