ELMER SOLVER (v 8.4) STARTED AT: 2019/01/15 17:17:47 ParCommInit: Initialize #PEs: 1 MAIN: OMP_NUM_THREADS not set. Using only 1 thread per task. 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: dbe89bf, Compiled: 2019-01-15) MAIN: Running one task without MPI parallelization. MAIN: Running with just one thread per task. MAIN: HYPRE library linked in. MAIN: MUMPS library linked in. MAIN: Lua interpreted linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: case.sif LoadInputFile: Scanning input file: case.sif LoadInputFile: Loading input file: case.sif Model Input: Unlisted keyword: [calculate magnetic field density] in section: [solver 3] Model Input: Unlisted keyword: [p {e}] in section: [boundary condition 1] Model Input: Unlisted keyword: [current density normalize by area] in section: [boundary condition 2] Model Input: Unlisted keyword: [current density normalize by area] in section: [boundary condition 4] Loading user function library: [StatCurrentSolveVec]...[StatCurrentSolver_Init0] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init0] WhitneyAVSolver_Init0: Ignoring scalar potential in AV equation! Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init0] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0] Loading user function library: [SaveData]...[SaveScalars_Init0] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./adjacent_open_coil_V01 LoadMesh: Reading header info from file: ./adjacent_open_coil_V01/mesh.header InitializeMesh: Number of nodes in mesh: 15854 InitializeMesh: Number of bulk elements in mesh: 95282 InitializeMesh: Number of boundary elements in mesh: 1810 InitializeMesh: Initial number of max element nodes: 4 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./adjacent_open_coil_V01/mesh.nodes LoadMesh: Performing coordinate mapping SetMeshDimension: Dimension of mesh is: 3 SetMeshDimension: Max dimension of mesh is: 3 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./adjacent_open_coil_V01/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./adjacent_open_coil_V01/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Remapping bodies LoadMesh: Minimum initial body index: 1 LoadMesh: Maximum initial body index: 3 LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 5 ElmerAsciiMesh: Performing step: 5 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining faces in 3D mesh FindMeshFaces3D: Number of faces found: 191455 FindMeshEdges: Determining edges in 3D mesh FindMeshEdges3D: Number of edges found: 112026 LoadMesh: Preparing mesh done LoadMesh: Elapsed REAL time: 1.1362 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed REAL time: 0.0303 (s) MAIN: ------------------------------------- AddSolvers: Setting up 5 solvers AddSolvers: Setting up solver 1: solving currents AddEquationBasics: Using procedure: StatCurrentSolveVec StatCurrentSolver AddEquationBasics: Setting up solver: solving currents Loading user function library: [StatCurrentSolveVec]...[StatCurrentSolver_Init] Loading user function library: [StatCurrentSolveVec]...[StatCurrentSolver_bulk] Loading user function library: [StatCurrentSolveVec]...[StatCurrentSolver] AddEquationBasics: Creating standard variable: potential OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: solving currents...done. OptimizeBandwidth: Half bandwidth without optimization: 1523 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 36 OptimizeBandwidth: --------------------------------------------------------- CreateDGPerm: Created permutation for DG nodes: 30712 AddSolvers: Setting up solver 2: mgdynamics AddEquationBasics: Using procedure: MagnetoDynamics WhitneyAVSolver AddEquationBasics: Setting up solver: mgdynamics Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_bulk] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver] AddEquationBasics: Creating standard variable: p OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamics...done. OptimizeBandwidth: Half bandwidth without optimization: 111793 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 9546 OptimizeBandwidth: --------------------------------------------------------- AddSolvers: Setting up solver 3: mgdynamicscalc AddEquationBasics: Using procedure: MagnetoDynamics MagnetoDynamicsCalcFields AddEquationBasics: Setting up solver: mgdynamicscalc Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_bulk] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields] AddEquationBasics: Creating standard variable: hr_dummy OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamicscalc...done. OptimizeBandwidth: Half bandwidth without optimization: 15584 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2643 OptimizeBandwidth: --------------------------------------------------------- CreateDGPerm: Created permutation for DG nodes: 381128 AddSolvers: Setting up solver 4: AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddSolvers: Setting up solver 5: AddEquationBasics: Using procedure: SaveData SaveScalars Loading user function library: [SaveData]...[SaveScalars_Init] Loading user function library: [SaveData]...[SaveScalars] AddMeshCoordinatesAndTime: Setting mesh coordinates and time SetInitialConditions: Setting up initial conditions (if any) MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: ComputeNodalWeights: Computing weights for the mesh entities BC weight: 1 0.58415932318820074 BC weight: 2 1.0945640753648100E-005 BC weight: 3 1.0945640755712418E-005 BC weight: 4 1.0945640754496450E-005 BC weight: 5 1.0945640754494295E-005 BF weight: 1 7.5355051262507254E-006 Body weight: 1 4.1921374419380965E-002 Body weight: 2 4.0837387672323649E-006 Body weight: 3 3.4517663590184210E-006 Mat weight: 1 4.1921374419380965E-002 Mat weight: 2 7.5355051262507254E-006 CalculateEntityWeights: All done ListSetCoefficients: Normalizing > current density < by 9.13606E+04 ListSetCoefficients: Normalizing > current density < by 9.13606E+04 SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: solving currents StatCurrentSolver: ------------------------------------------------ StatCurrentSolver: Solving static current conduction solver DefaultStart: Starting solver: solving currents StatCurrentSolver: Performing non-vectorized bulk element assembly ListInitElementKeyword: Treating keyword: Current Source ListInitElementKeyword: Initiated handle for: > Current Source < of type: -1 ListInitElementKeyword: Treating keyword: Electric Conductivity ListInitElementKeyword: Initiated handle for: > Electric Conductivity < of type: 4 ListInitElementKeyword: Treating keyword: Relative Permittivity ListInitElementKeyword: Initiated handle for: > Relative Permittivity < of type: 4 StatCurrentSolverBulkAssembly: Elapsed REAL time: 0.0924 (s) ModelPDEthreaded: Assembly of boundary colour: 1 ListInitElementKeyword: Treating keyword: Current Density ListInitElementKeyword: Initiated handle for: > Current Density < of type: 4 ListInitElementKeyword: Treating keyword: Electric Resistivity ListInitElementKeyword: Initiated handle for: > Electric Resistivity < of type: -1 ListInitElementKeyword: Treating keyword: External Potential ListInitElementKeyword: Initiated handle for: > External Potential < of type: -1 ListInitElementKeyword: Treating keyword: Farfield Potential ListInitElementKeyword: Initiated handle for: > Farfield Potential < of type: -1 StatCurrentSolverBCAssembly: Elapsed REAL time: 0.0008 (s) DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions EnforceDirichletConditions: Applying Dirichlet conditions using scaled diagonal ScaleLinearSystem: Scaling diagonal entries to unity DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu1 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab CRS_IncompleteLU: ILU(1) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 42596 CRS_IncompleteLU: ILU(1) (Real), NOF nonzeros: 42597 CRS_IncompleteLU: ILU(1) (Real), filling (%) : 141 CRS_IncompleteLU: ILU(1) (Real), Factorization ready at (s): 0.00 10 0.2012E+00 20 0.2625E+00 30 0.1340E+00 40 0.1436E+00 50 0.1017E+00 60 0.4874E-01 70 0.1530E-01 80 0.1079E-01 90 0.6120E-02 100 0.1295E-02 110 0.1032E-02 120 0.7929E-03 130 0.4063E-03 140 0.1771E-06 150 0.7158E-09 160 0.1331E-09 161 0.7802E-10 ComputeNorm: Computing norm of solution ComputedNorm: 9.7283914152921230E-003 2 ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.97283914E-02 2.0000000 ) :: solving currents DefaultFinish: Finished solver: solving currents Loading user function library: [StatCurrentSolveVec]...[StatCurrentSolver_post] StatCurrentSolver_post: ------------------------------------------------ StatCurrentSolver_post: Calculating postprocessing fields StatCurrentSolver_post: Number of 1 postprocessing fields ListInitElementKeyword: Treating keyword: Current Source ListInitElementKeyword: Initiated handle for: > Current Source < of type: -1 ListInitElementKeyword: Treating keyword: Electric Conductivity ListInitElementKeyword: Initiated handle for: > Electric Conductivity < of type: 4 ListInitElementKeyword: Treating keyword: Relative Permittivity ListInitElementKeyword: Initiated handle for: > Relative Permittivity < of type: 4 ComputeNorm: Computing norm of solution ComputedNorm: 9.7283914152921230E-003 2 ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.97283914E-02 2.0000000 ) :: solving currents SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: mgdynamics WhitneyAVSolver: ------------------------------------------- WhitneyAVSolver: Solving the AV equations with edge elements OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: mgdynamics...done. OptimizeBandwidth: Half bandwidth without optimization: 15584 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2643 OptimizeBandwidth: --------------------------------------------------------- DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions EnforceDirichletConditions: No Dirichlet conditions to enforce, exiting! DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstab CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 239906 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 239907 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01 10 0.1813E-03 20 0.1262E-05 30 0.6590E-08 39 0.4052E-11 DefaultStart: Starting solver: mgdynamics MGDynAssembly: Elapsed REAL time: 5.0452 (s) DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions EnforceDirichletConditions: Applying Dirichlet conditions using scaled diagonal ScaleLinearSystem: Scaling diagonal entries to unity DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: none SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: bicgstabl 10 0.3089E-01 0.3089E-01 20 0.2555E-02 0.2555E-02 30 0.2217E-03 0.2217E-03 40 0.1679E-04 0.1679E-04 50 0.3109E-05 0.3109E-05 60 0.4036E-06 0.4036E-06 68 0.9452E-08 0.9452E-08 ComputeNorm: Computing norm of solution ComputedNorm: 3.1223674370985976E-008 2 ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.31223674E-07 2.0000000 ) :: mgdynamics DefaultFinish: Finished solver: mgdynamics WhitneyAVSolver: All done WhitneyAVSolver: ------------------------------------------- Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_post] ComputeNorm: Computing norm of solution ComputedNorm: 3.1223674370985976E-008 2 ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.31223674E-07 2.0000000 ) :: mgdynamics SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: mgdynamicscalc MagnetoDynamicsCalcFields: ------------------------------ MagnetoDynamicsCalcFields: Computing postprocessed fields MagnetoDynamicsCalcFields: Solving for field: magnetic flux density[magnetic flux density:3] DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: cg CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 239906 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 239907 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.01 1 0.3243E-01 2 0.2242E-02 3 0.1279E-03 4 0.6803E-05 5 0.4278E-06 6 0.2675E-07 7 0.1620E-08 7 0.1620E-08 ComputeNorm: Computing norm of solution ComputedNorm: 5.3432418127216370E-004 2 ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.53432418E-03 2.0000000 ) :: mgdynamicscalc DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: cg 1 0.3221E-01 2 0.2217E-02 3 0.1245E-03 4 0.6489E-05 5 0.3975E-06 6 0.2453E-07 7 0.1489E-08 7 0.1489E-08 ComputeNorm: Computing norm of solution ComputedNorm: 5.2937937488501275E-004 2 ComputeChange: NS (ITER=2) (NRM,RELC): ( 0.52937937E-03 0.92973392E-02 ) :: mgdynamicscalc DefaultSolve: Solving linear system with default routines SolveSystem: Solving linear system SolveLinearSystem: Assuming real valued linear system ScaleLinearSystem: Scaling diagonal entries to unity SolveSystem: Linear System Solver: iterative SolveSystem: Linear System Preconditioning: ilu0 SolveSystem: Serial linear System Solver: iterative IterSolver: Using iterative method: cg 1 0.3283E-01 2 0.2153E-02 3 0.1223E-03 4 0.6218E-05 5 0.3767E-06 6 0.2341E-07 7 0.1448E-08 7 0.1448E-08 ComputeNorm: Computing norm of solution ComputedNorm: 7.8873071483073385E-004 2 ComputeChange: NS (ITER=3) (NRM,RELC): ( 0.78873071E-03 0.39352000 ) :: mgdynamicscalc MagnetoDynamicsCalcFields: Eddy current power: 0.0000000000000000 MagnetoDynamicsCalcFields: (Electro)Magnetic Field Energy: 2.1380630991181497E-004 Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_post] ComputeNorm: Computing norm of solution ComputedNorm: 7.8873071483073385E-004 2 ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.78873071E-03 2.0000000 ) :: mgdynamicscalc SingleSolver: Attempting to call solver ResultOutputSolver: ------------------------------------- ResultOutputSolve: Saving with prefix: results ResultOutputSolver: Working on mesh: adjacent_open_coil_V01 ResultOutputSolver: Dimension of mesh is: 3 ResultOutputSolver: Creating list for saving - if not present CreateListForSaving: Field Variables for Saving CreateListForSaving: Scalar Field 1: potential CreateListForSaving: Scalar Field 2: jfix CreateListForSaving: Vector Field 1: volume current e CreateListForSaving: Vector Field 2: magnetic flux density CreateListForSaving: Vector Field 3: magnetic flux density e ResultOutputSolver: Saving in unstructured VTK XML (.vtu) format VtuOutputSolver: Saving results in VTK XML format with prefix: results VtuOutputSolver: Saving number of partitions: 1 VtuOutputSolver: Full filename base is: ./adjacent_open_coil_V01/results VtuOutputSolver: Number of active elements 97092 out of 97092 VtuOutputSolver: Number of geometry nodes 15854 out of 15854 VtuOutputSolver: Total number of geometry nodes to save: 15854 VtuOutputSolver: Total number of dof nodes to save: 15854 VtuOutputSolver: Total number of elements to save: 97092 VtuOutputSolver: Setting offset for boundary entities: 100 VtuOutputSolver: Writing the vtu file: ./adjacent_open_coil_V01/results0001.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: 97092 VtuOutputSolver: Writing nodal fields VtuOutputSolver: Saving variable: potential VtuOutputSolver: Saving variable: jfix VtuOutputSolver: Saving variable: volume current e VtuOutputSolver: Saving variable: magnetic flux density VtuOutputSolver: Saving variable: magnetic flux density e VtuOutputSolver: Number of nodal fields written: 3 VtuOutputSolver: Writing elemental fields WriteVtuFile: Writing variable: volume current e WriteVtuFile: Writing variable: magnetic flux density e VtuOutputSolver: Number of elemental fields written: 2 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 VtuOutputSolver: Saving variable: jfix VtuOutputSolver: Saving variable: volume current e VtuOutputSolver: Saving variable: magnetic flux density VtuOutputSolver: Saving variable: magnetic flux density e VtuOutputSolver: Writing elemental fields AscBinWriteFree: Terminating buffered ascii/binary writing VtuOutputSolver: All done for now ResultOutputSolver: ------------------------------------- Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_post] SingleSolver: Attempting to call solver SaveScalars: ----------------------------------------- SaveScalars: Saving scalar values of various kinds SaveScalars: Performed 0 reduction operations SaveScalars: Found 2 result scalars in simulation section SaveScalars: Found 2 values to save in total SaveScalars: Saving names of values to file: ./data_scalars.dat.names SaveScalars: Saving values to file: ./data_scalars.dat SaveScalars: ----------------------------------------- Loading user function library: [SaveData]...[SaveScalars_post] ReloadInputFile: Realoading input file LoadInputFile: Loading input file: ElmerSolver: *** Elmer Solver: ALL DONE *** ElmerSolver: The end SOLVER TOTAL TIME(CPU,REAL): 18.46 18.27 ELMER SOLVER FINISHED AT: 2019/01/15 17:18:05