ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ELMER SOLVER (v 8.4) STARTED AT: 2020/04/14 03:39:46 ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= LoadInputFile: Scanning input file: AdvReactDB.sif MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif LoadInputFile: Scanning input file: AdvReactDB.sif ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif LoadInputFile: Scanning input file: AdvReactDB.sif ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= LoadInputFile: Scanning input file: AdvReactDB.sif MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif LoadInputFile: Scanning input file: AdvReactDB.sif ParCommInit: Initialize #PEs: 1 MAIN: MAIN: ============================================================= MAIN: ElmerSolver finite element software, Welcome! MAIN: This program is free software licensed under (L)GPL MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd. MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi MAIN: Version: 8.4 (Rev: bf231bac, Compiled: 2020-04-14) MAIN: Running one task without MPI parallelization. MAIN: Running in parallel with 6 threads per task. MAIN: HYPRE library linked in. MAIN: Trilinos library linked in. MAIN: MUMPS library linked in. MAIN: Intel MKL linked in. MAIN: Lua interpreted linked in. MAIN: Zoltan library linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: AdvReactDB.sif LoadInputFile: Loading input file: AdvReactDB.sif Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 LoadInputFile: Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 LoadInputFile: Number of Bodies: 4 Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] LoadInputFile: Scanning input file: AdvReactDB.sif LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square LoadMesh: Reading header info from file: ./square/mesh.header InitializeMesh: Number of nodes in mesh: 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./square/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 0.0019 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: Completing keywords for mortar BCs MAIN: ------------------------------------- AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 CreateMatrix: Matrix created AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddEquationBasics: Creating null variable with no name AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 ElmerSolver: Number of timesteps to be saved: 1 LoadInputFile: Loading input file: AdvReactDB.sif Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 LoadInputFile: Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 LoadInputFile: Number of Bodies: 4 Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 LoadInputFile: Loading input file: AdvReactDB.sif Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square LoadMesh: Reading header info from file: ./square/mesh.header InitializeMesh: Number of nodes in mesh: 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./square/mesh.boundary LoadInputFile: LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity LoadInputFile: Number of Bodies: 4 FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 0.0008 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: Completing keywords for mortar BCs MAIN: ------------------------------------- AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] LoadInputFile: Loading input file: AdvReactDB.sif AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 CreateMatrix: Matrix created AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 LoadInputFile: Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 LoadInputFile: Number of Bodies: 4 Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] LoadMesh: Reading header info from file: ./square/mesh.header AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] InitializeMesh: Number of nodes in mesh: 240 AddEquationBasics: Creating null variable with no name AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square LoadMesh: Reading header info from file: ./square/mesh.header InitializeMesh: Number of nodes in mesh: 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 LoadMesh: ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes Reading boundary elements from file: ./square/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./square/mesh.boundary Performing step: 5 LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity NonNodalElements: Requested elements require creation of edges FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 0.0008 (s) MeshStabParams: Computing stabilization parameters FindMeshEdges: MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: Completing keywords for mortar BCs MAIN: ------------------------------------- AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...LoadMesh: Preparing mesh done done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 LoadMesh: Elapsed REAL time: 0.0030 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: CreateMatrix: Matrix created Completing keywords for mortar BCs MAIN: ------------------------------------- AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact ElmerSolver: Number of timesteps to be saved: 1 AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 CreateMatrix: Matrix created AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddEquationBasics: Creating null variable with no name AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddEquationBasics: Creating null variable with no name ElmerSolver: Number of timesteps to be saved: 1 AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 ElmerSolver: Number of timesteps to be saved: 1 MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 LoadInputFile: Loading input file: AdvReactDB.sif Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 LoadInputFile: Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 LoadInputFile: Number of Bodies: 4 Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] LoadInputFile: Loading input file: AdvReactDB.sif LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square LoadMesh: Reading header info from file: ./square/mesh.header InitializeMesh: Number of nodes in mesh: 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./square/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 0.0008 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: Completing keywords for mortar BCs MAIN: ------------------------------------- AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix Model Input: Unlisted keyword: [vtu: discontinuous galerkin] in section: [simulation] Model Input: Unlisted keyword: [potential gamma] in section: [material 1] Model Input: Unlisted keyword: [potential source] in section: [body force 1] LoadInputFile: Number of BCs: 1 LoadInputFile: Number of Body Forces: 1 LoadInputFile: Number of Initial Conditions: 0 LoadInputFile: Number of Materials: 1 LoadInputFile: Number of Equations: 1 LoadInputFile: Number of Solvers: 1 LoadInputFile: Number of Bodies: 4 Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init0] OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 CreateMatrix: Matrix created AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./square LoadMesh: Reading header info from file: ./square/mesh.header InitializeMesh: Number of nodes in mesh: 240 InitializeMesh: Number of bulk elements in mesh: 210 InitializeMesh: Number of boundary elements in mesh: 58 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./square/mesh.nodes LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 2 SetMeshDimension: Max dimension of mesh is: 2 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./square/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./square/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Skipping remapping of bodies LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 4 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining edges in 2D mesh FindMeshEdges2D: Allocating edge table of size: 840 FindMeshEdges2D: Creating hash table of size 240 for noto-to-node connectivity FindMeshEdges2D: Number of edges found: 449 FindMeshEdges2D: All done LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 0.0029 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0000 (s) CompleteModelKeywords: Completing keywords for mortar BCs MAIN: ------------------------------------- AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: DBsolution AddVtuOutputSolverHack: Increasing number of solver to: 2 ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous galerkin ListCopyPrefixedKeywords: Copied 1 keywords with prefix: vtu: AddVtuOutputSolverHack: Finished appeding VTU output solver AddSolvers: Setting up 2 solvers AddSolvers: Setting up solver 1: advreact AddEquationBasics: Using procedure: AdvectionReaction AdvectionReactionSolver AddEquationBasics: Setting up solver: advreact Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver] AddEquationBasics: Treating variable string: potential AddEquationBasics: Creating standard variable: potential AddEquationBasics: Computing size of permutation vector AddEquationBasics: Maximum size of permutation vector is: 1072 AddEquationBasics: Creating solver matrix topology CreateMatrix: creating initial permutation InitialPermutation: Group 1 starts from index 0 InitialPermutation: Group 1 has 64 db dofs InitialPermutation: Group 2 starts from index 64 InitialPermutation: Group 2 has 64 db dofs InitialPermutation: Group 3 starts from index 128 InitialPermutation: Group 3 has 72 db dofs InitialPermutation: Group 4 starts from index 200 InitialPermutation: Group 4 has 72 db dofs InitialPermutation: Numbered 272 db nodes from bulk hits CreateMatrix: Creating list matrix for equation: advreact MakeListMatrix: Creating list matrix AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddEquationBasics: Creating null variable with no name AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: advreact...done. OptimizeBandwidth: Half bandwidth without optimization: 202 OptimizeBandwidth: --------------------------------------------------------- CreateMatrix: Initializing list matrix for equation CRS_CreateMatrix: Creating CRS Matrix of size: 272 CRS_CreateMatrix: Creating CRS Matrix with nofs: 2796 CreateMatrix: Matrix created AddEquationBasics: Number of rows in CRS matrix: 272 AddEquationBasics: Creating solver variable VariableAdd: Adding variable > potential < of size 272 AddEquationBasics: Creating exported variable: nodal result VariableAdd: Adding variable > nodal result < of size 272 AddEquationBasics: Succesfully created variable: nodal result AddSolvers: Setting up solver 2: internalvtuoutputsolver AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: internalvtuoutputsolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] AddEquationBasics: Checking for _init solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] AddEquationBasics: Checking for _bulk solver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddEquationBasics: Creating null variable with no name AddSolvers: Setting up solvers done AddMeshCoordinatesAndTime: Setting mesh coordinates and time VariableAdd: Adding variable > Coordinate 1 < of size 240 VariableAdd: Adding variable > Coordinate 2 < of size 240 VariableAdd: Adding variable > Coordinate 3 < of size 240 VariableAdd: Adding variable > Time < of size 1 VariableAdd: Adding variable > Periodic Time < of size 1 VariableAdd: Adding variable > Timestep < of size 1 VariableAdd: Adding variable > Timestep size < of size 1 VariableAdd: Adding variable > Timestep interval < of size 1 VariableAdd: Adding variable > nonlin iter < of size 1 VariableAdd: Adding variable > coupled iter < of size 1 VariableAdd: Adding variable > Partition < of size 1 SetInitialConditions: Setting up initial conditions (if any) InitCond mesh:square 2 240 ElmerSolver: Number of timesteps to be saved: 1 ElmerSolver: Number of timesteps to be saved: 1 MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SolveEquations: Solvers before timestep SolveEquations: Solvers in main iteration loop SolveEquations: Performing set of solvers in sequence SetActiveElementsTable: Creating active element table for: advreact SetActiveElementsTable: Number of active elements found : 210 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: advreact AdvectionReaction (potential): AdvectionReactionSolver for variable potential AdvectionReaction (potential): Memory allocation done WARNING:: AdvectionReaction (potential): No >Nonlinear System Min Iterations< found. Setting 1 AdvectionReaction (potential): No keyword > Limit Solution < found. Solution will not be limited AdvectionReaction (potential): AdvectionReaction (potential): AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): Nonlinear iteration no. 1 of max 1 AdvectionReaction (potential): ------------------------------------- AdvectionReaction (potential): AdvectionReaction (potential): Starting Assembly... InitializeToZero: Initializing the linear system to zero GetNOFActive: Number of active elements: 210 AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.1045E+01 2 0.7378E+00 3 0.2108E+01 4 0.4189E+00 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 34 0.8240E-05 35 0.6736E-05 36 0.6808E-05 37 0.1171E-04 38 0.5715E-05 39 0.2418E-05 40 0.1746E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.15 0.15 AdvectionReaction (potential): iter: 1 Solve: (s) 0.02 0.02 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.1045E+01 2 0.7378E+00 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges 3 0.2108E+01 ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists 4 0.4189E+00 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 34 0.8240E-05 35 0.6736E-05 36 0.6808E-05 37 0.1171E-04 38 0.5715E-05 39 0.2418E-05 40 0.1746E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.16 0.16 AdvectionReaction (potential): iter: 1 Solve: (s) 0.02 0.02 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.1045E+01 2 0.7378E+00 3 0.2108E+01 4 0.4189E+00 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 34 0.8240E-05 35 0.6736E-05 36 0.6808E-05 37 0.1171E-04 38 0.5715E-05 39 0.2418E-05 40 0.1746E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.15 0.15 AdvectionReaction (potential): iter: 1 Solve: (s) 0.02 0.02 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ComputeNorm: Computing norm of solution ComputeChange: ElmerSolver: The end SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact SOLVER TOTAL TIME(CPU,REAL): 0.42 0.87 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47 ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.1045E+01 2 0.7378E+00 3 0.2108E+01 ElmerSolver: The end SOLVER TOTAL TIME(CPU,REAL): 0.41 0.89 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47 4 0.4189E+00 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 34 0.8240E-05 35 0.6736E-05 36 0.6808E-05 37 0.1171E-04 38 0.5715E-05 39 0.2418E-05 40 0.1746E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.16 0.16 AdvectionReaction (potential): iter: 1 Solve: (s) 0.01 0.01 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists ElmerSolver: The end SOLVER TOTAL TIME(CPU,REAL): 0.35 0.82 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47 AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.1045E+01 2 0.7378E+00 3 0.2108E+01 4 0.4189E+00 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 AdvectionReaction (potential): Assembly done DefaultSolve: Solving linear system with default routines GenerateConstraintMatrix: Nothing to do for now DefaultSolve: Calling SolveSystem for linear solution range b0: 1.3234893494237605E-003 0.20638599743213587 9.4808268267976636 9.4808268267976636 range A0: -0.13605442176849725 6.9786563799379059E-002 6.0000000000002158 55.713794986593030 SolveSystem: Solving linear system SolveSystem: Solving linear system without constraint matrix SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity ScaleLinearSystem: Computing > DiagScaling < vector ScaleLinearSystem: Scaling matrix values ScaleLinearSystem: Scaling Rhs vector SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab IterSolver: Matrix is real valued CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 2796 CRS_IncompleteLU: Performing incomplete LU ILU range: -12.765847660953263 34.326590513362945 1.6498214042595511 2.1322003800851173 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 2797 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 34 0.8240E-05 1 0.1045E+01 35 0.6736E-05 2 0.7378E+00 36 0.6808E-05 3 0.2108E+01 37 0.1171E-04 4 0.4189E+00 38 0.5715E-05 5 0.4342E+00 6 0.1939E+01 7 0.2904E+00 8 0.2781E+00 9 0.2399E+00 10 0.2543E+00 11 0.2707E+00 12 0.1954E+00 13 0.3222E+00 14 0.2679E+00 39 0.2418E-05 ElmerSolver: The end SOLVER TOTAL TIME(CPU,REAL): 0.33 0.79 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47 15 0.6317E-01 16 0.6245E-01 17 0.2500E-01 18 0.9198E-02 19 0.8049E-02 20 0.9932E-02 21 0.4210E-02 22 0.5439E-02 40 0.1746E-05 23 0.3920E-02 24 0.4408E-02 25 0.4145E-02 26 0.3735E-03 27 0.6058E-03 28 0.2315E-03 29 0.2875E-03 30 0.9641E-03 31 0.9501E-03 32 0.2139E-04 33 0.1254E-04 34 0.8240E-05 35 0.6736E-05 36 0.6808E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.16 0.16 AdvectionReaction (potential): iter: 1 Solve: (s) 0.02 0.02 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! 37 0.1171E-04 ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields 38 0.5715E-05 AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists 39 0.2418E-05 40 0.1746E-05 41 0.1687E-06 41 0.1687E-06 BackScaleLinearSystem: Scaling back to original scale ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact UpdateDependentObjects: Updating objects depending on primary field in nonlinear system SolveSystem: System solved SolveSystem: Finished solving the system range x0: 0.99632790241631930 2.0076405490417035 367.45093174344083 367.45093174344083 AdvectionReaction (potential): iter: 1 Assembly: (s) 0.16 0.16 AdvectionReaction (potential): iter: 1 Solve: (s) 0.02 0.02 AdvectionReaction (potential): Using >Exported Variable 1< for the nodal output field Loading user function library: [AdvectionReaction]...[AdvectionReactionSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 1.3927519 2.0000000 ) :: advreact ListPopNameSpace: No namespace entry to delete SolveEquations: Solvers after timestep SetActiveElementsTable: Creating active element table for: internalvtuoutputsolver SetActiveElementsTable: No active elements found SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: internalvtuoutputsolver ResultOutputSolver: ------------------------------------- ResultOutputSolver: Saving with prefix: DBsolution ResultOutputSolver: Number of output meshes: 1 ResultOutputSolver: Working on mesh: square ResultOutputSolver: Dimension of mesh is: 2 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: nodal result ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format CheckAndCreateDGIndexes: Creating DG indexes! CheckAndCreateDGIndexes: Creating DG 840 indexes VtuOutputSolver: Saving results in VTK XML format with prefix: DBsolution VtuOutputSolver: Saving number of partitions: 1 SolverOutputDirectory: Creating directory: square VtuOutputSolver: Full filename base is: square/DBsolution VtuOutputSolver: Number of active elements 268 out of 268 VtuOutputSolver: Number of geometry nodes 240 out of 240 VtuOutputSolver: Saving results as discontinuous DG fields VtuOutputSolver: Maximum number of dofs in DG: 840 VtuOutputSolver: Independent dofs in discontinuous mesh: 840 VtuOutputSolver: Total number of geometry nodes to save: 240 VtuOutputSolver: Total number of dof nodes to save: 840 VtuOutputSolver: Total number of elements to save: 268 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: square/DBsolution_t0001.vtu AscBinWriteInit: Initializing buffered ascii/binary writing AscBinWriteInit: Writing in binary AscBinWriteInit: Writing in double precision AscBinWriteInit: Writing to unit number: 58 AscBinWriteInit: Size of buffer is: 1000 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Number of nodal fields written: 2 VtuOutputSolver: Writing elemental fields VtuOutputSolver: Number of elemental fields written: 0 VtuOutputSolver: Writing entity IDs for bodies and boundaries VtuOutputSolver: Writing coordinates for each used node VtuOutputSolver: Writing the elemental connectivity data VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential WriteVTUFile: Setting field type to discontinuous VtuOutputSolver: Saving variable: nodal result WriteVTUFile: Setting field type to nodal VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing WriteVtuFile: Finished writing file VtuOutputSolver: Finished writing vtu file VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] UpdateDependentObjects: Updating objects depending on primary field in steady state DerivateExportedVariables: Derivating variables, if any! ListPopNameSpace: No namespace entry to delete ListPopNameSpace: No namespace entry to delete CompareToReferenceSolution: Solver 1 PASSED: Norm = 1.39275190E+00 RefNorm = 1.39275188E+00 CompareToReferenceSolution: Relative Error to reference norm: 1.106709E-08 ReloadInputFile: Realoading input file LoadInputFile: Loading input file: CompareToReferenceSolution: PASSED all 1 tests! ElmerSolver: *** Elmer Solver: ALL DONE *** FreeModel: Freeing meshes ReleaseMesh: Releasing mesh variables ReleaseMesh: Releasing mesh nodes ReleaseMesh: Releasing mesh edges ReleaseMesh: Releasing mesh elements ReleaseMesh: Releasing mesh finished FreeModel: Freeing constants list FreeModel: Freeing simulation list FreeModel: Freeing boundary lists FreeModel: Freeing solvers FreeModel: Solver: 1 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Solver: 2 FreeSolver: Free solver matrix FreeSolver: Free solver miscallenous FreeModel: Freeing initial conditions lists FreeModel: Freeing body lists FreeModel: Freeing equations lists FreeModel: Freeing body forces lists ElmerSolver: The end ElmerSolver: The end SOLVER TOTAL TIME(CPU,REAL): 0.31 0.87 SOLVER TOTAL TIME(CPU,REAL): 0.35 0.90 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47 ELMER SOLVER FINISHED AT: 2020/04/14 03:39:47