ELMER SOLVER (v 8.2) STARTED AT: 2017/06/17 16:58:37 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.2 (Rev: 4b428c1, Compiled: 2017-01-20) MAIN: HYPRE library linked in. MAIN: MUMPS library linked in. MAIN: ============================================================= MAIN: MAIN: MAIN: ------------------------------------- MAIN: Reading Model: magnetic_field_with_mortar_only_2_structured_surfaces_V03.sif LoadInputFile: Scanning input file: magnetic_field_with_mortar_only_2_structured_surfaces_V03.sif LoadInputFile: Loading input file: magnetic_field_with_mortar_only_2_structured_surfaces_V03.sif Model Input: Unlisted keyword: [calculate magnetic field density] in section: [solver 3] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 1] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 2] Model Input: Unlisted keyword: [p {e}] in section: [boundary condition 3] Model Input: Unlisted keyword: [magnetic flux density {n}] in section: [boundary condition 4] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_Init0] Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init0] Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init0] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init0] Loading user function library: [SaveData]...[SaveScalars_Init0] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_Init0] LoadMesh: Starting ElmerAsciiMesh: Performing step: 1 LoadMesh: Base mesh name: ./6-4_m_generic_added_one_physical_surface_at_top_and_one_on_bottom_structured_grid LoadMesh: Reading header info from file: ./6-4_m_generic_added_one_physical_surface_at_top_and_one_on_bottom_structured_grid/mesh.header LoadMesh: Number of nodes in mesh: 36747 LoadMesh: Number of bulk elements in mesh: 193890 LoadMesh: Number of boundary elements in mesh: 14048 LoadMesh: Initial number of max element nodes: 5 ElmerAsciiMesh: Performing step: 2 LoadMesh: Reading nodes from file: ./6-4_m_generic_added_one_physical_surface_at_top_and_one_on_bottom_structured_grid/mesh.nodes LoadMesh: Performing coordinate mapping LoadMesh: Dimension of model is: 3 LoadMesh: Dimension of mesh is: 3 ElmerAsciiMesh: Performing step: 3 LoadMesh: Reading bulk elements from file: ./6-4_m_generic_added_one_physical_surface_at_top_and_one_on_bottom_structured_grid/mesh.elements ElmerAsciiMesh: Performing step: 4 LoadMesh: Reading boundary elements from file: ./6-4_m_generic_added_one_physical_surface_at_top_and_one_on_bottom_structured_grid/mesh.boundary LoadMesh: Performing node mapping LoadMesh: Remapping bodies LoadMesh: Minimum initial body index: 1 LoadMesh: Maximum initial body index: 6 LoadMesh: Remapping boundaries LoadMesh: Minimum initial boundary index: 1 LoadMesh: Maximum initial boundary index: 10 ElmerAsciiMesh: Performing step: 5 NonNodalElements: Requested elements require creation of edges FindMeshEdges: Determining faces in 3D mesh FindMeshFaces3D: Number of faces found: 395819 FindMeshEdges: Determining edges in 3D mesh FindMeshEdges3D: Number of edges found: 238675 ElmerAsciiMesh: Performing step: 6 LoadMesh: Loading mesh done LoadMesh: Elapsed time (CPU,REAL): 0.8080 0.9973 (s) MeshStabParams: Computing stabilization parameters MeshStabParams: Elapsed time (CPU,REAL): 0.1080 0.1076 (s) MAIN: ------------------------------------- AddSolvers: Setting up 6 solvers AddSolvers: Setting up solver 1: static current conduction AddEquationBasics: Using procedure: StatCurrentSolve StatCurrentSolver AddEquationBasics: Setting up solver: static current conduction Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_Init] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver_bulk] Loading user function library: [StatCurrentSolve]...[StatCurrentSolver] AddEquationBasics: Creating standard variable: potential OptimizeBandwidth: --------------------------------------------------------- OptimizeBandwidth: Computing matrix structure for: static current conduction...done. OptimizeBandwidth: Half bandwidth without optimization: 667 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 131 OptimizeBandwidth: --------------------------------------------------------- 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: 159474 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 7697 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: 25630 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2057 OptimizeBandwidth: --------------------------------------------------------- AddSolvers: Setting up solver 4: result output AddEquationBasics: Using procedure: ResultOutputSolve ResultOutputSolver AddEquationBasics: Setting up solver: result output Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_Init] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver_bulk] Loading user function library: [ResultOutputSolve]...[ResultOutputSolver] AddSolvers: Setting up solver 5: scalars AddEquationBasics: Using procedure: SaveData SaveScalars AddEquationBasics: Setting up solver: scalars Loading user function library: [SaveData]...[SaveScalars_Init] Loading user function library: [SaveData]...[SaveScalars_bulk] Loading user function library: [SaveData]...[SaveScalars] AddSolvers: Setting up solver 6: never AddEquationBasics: Using procedure: MagnetoDynamics MagnetoDynamics_Dummy AddEquationBasics: Setting up solver: never Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_Init] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy_bulk] Loading user function library: [MagnetoDynamics]...[MagnetoDynamics_Dummy] AddEquationBasics: Creating standard variable: cf_dummy AddMeshCoordinatesAndTime: Setting mesh coordinates and time SetInitialConditions: Setting up initial conditions (if any) MAIN: MAIN: ------------------------------------- MAIN: Steady state iteration: 1 MAIN: ------------------------------------- MAIN: SingleSolver: Attempting to call solver SingleSolver: Solver Equation string is: static current conduction StatCurrentSolve: ------------------------------------- StatCurrentSolve: STAT CURRENT SOLVER: StatCurrentSolve: ------------------------------------- StatElecSolve: Starting Assembly... DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: potential SetDirichletBoundaries: Number of dofs set: 900 DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set StatCurrentSolve: Assembly (s) : 2.4000000000000021E-002 SolveSystem: Solving linear system IterSolver: Using iterative method: bicgstab CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 10765 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 10766 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.2474E-01 2 0.1406E-01 3 0.9820E-02 4 0.6236E-02 5 0.3645E-02 6 0.2470E-02 7 0.1311E-02 8 0.2603E-03 9 0.6497E-04 10 0.1170E-04 11 0.1026E-05 12 0.9705E-07 13 0.9588E-08 13 0.9588E-08 ComputeNorm: Computing norm of solution ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.32110988E-01 2.0000000 ) :: static current conduction StatCurrentSolve: Solve (s) : 4.0000000000000036E-003 StatCurrentSolve: Total Heating Power : 527.97811167748205 StatCurrentSolve: Effective Resistance : 5.7294041800123634E-006 StatCurrentSolve: Total Volume : 3.7948254420513609E-006 StatCurrentSolve: Total Current : 9599.6020304996728 StatCurrentSolve: Control Scaling : 4.1668394036453553E-003 ComputeNorm: Computing norm of solution ComputeChange: SS (ITER=1) (NRM,RELC): ( 0.13380133E-03 2.0000000 ) :: static current conduction GenerateProjectors: Generating mortar projectors PeriodicProjector: ----------------------------------------------------- PeriodicProjector: Creating projector between BCs 5 and 6 CreateInterfaceMeshes: Making a list of elements at interface CreateInterfaceMeshes: Number of interface elements: 1015, 1015 CreateInterfaceMeshes: Number of interface nodes: 1080, 1080 PeriodicProjector: Enforcing > Level Projector = True < with dimensional reduction RadialInterfaceMeshes: Transformed extrema for this boundary (phi,r,z) RadialInterfaceMeshes: Coordinate 1: 5.000E-04 5.090E-02 RadialInterfaceMeshes: Coordinate 2: 9.000E+01 9.000E+01 RadialInterfaceMeshes: Coordinate 3: -7.500E-03 1.833E-02 RadialInterfaceMeshes: Discrepancy from constant angle (degs): 1.421E-14 RadialInterfaceMeshes: Transformed extrema for target boundary (phi,r,z) RadialInterfaceMeshes: Coordinate 1: 5.000E-04 5.090E-02 RadialInterfaceMeshes: Coordinate 2: -9.000E+01 -9.000E+01 RadialInterfaceMeshes: Coordinate 3: -7.500E-03 1.833E-02 RadialInterfaceMeshes: Discrepancy from constant angle (degs): 2.842E-14 RadialInterfaceMeshes: Discrepancy in maximum radius: 0.000E+00 RadialInterfaceMeshes: Discrepancy in minimum radius: 0.000E+00 LevelProjector: Creating projector for a levelized mesh LevelProjector: Checking mesh skew LevelProjector: Maximum skew in this mesh: 0.000E+00 LevelProjector: Checking mesh skew LevelProjector: Maximum skew in target mesh: 3.183E-09 LevelProjector: Initial number of slave nodes 1080 out of 1080 LevelProjector: Number of potential nodes in projector: 1080 LevelProjector: Number of active nodes in projector: 1080 LevelProjector: Number of potential edges in projector: 2094 LevelProjector: Eliminating redundant edges from projector: 0 LevelProjector: Number of active edges in projector: 2094 LevelProjector: Max number of rows in projector: 3174 LevelProjector: Creating strong stride projector for nodal dofs LevelProjector: Creating strong stride projector for edges assuming strides LevelProjector: Minimum absolute sum of edge weights: 1.00000E+00 LevelProjector: Maximum absolute sum of edge weights: 1.00000E+00 LevelProjector: Created strong constraints for edge dofs List_ToCRSMatrix: Number of entries in CRS matrix: 6526 ListToCRSMatrix: Matrix format changed from List to CRS LevelProjector: Projector created PeriodicProjector: Elapsed time (CPU,REAL): 0.0280 0.0420 (s) PeriodicProjector: Projector created, now exiting... 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: 25630 OptimizeBandwidth: OptimizeBandwidth: Bandwidth Optimization ...done. OptimizeBandwidth: Half bandwidth after optimization: 2057 OptimizeBandwidth: --------------------------------------------------------- DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: jfix SetDirichletBoundaries: Number of dofs set: 0 DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set GenerateConstraintMatrix: There are 3174 initial rows in constraint matrices GenerateConstraintMatrix: Adding projector for BC: 5 GenerateConstraintMatrix: Allocating 1080 rows and 2220 nonzeros GenerateConstraintMatrix: Adding projector for BC: 5 SolveSystem: Solving linear system SolveSystem: Solving linear system with constraint matrix SolveWithLinearRestriction: SolveWithLinearRestriction: Adding ConstraintMatrix into CollectionMatrix SolveWithLinearRestriction: Adding Stiffness Matrix into CollectionMatrix SolveWithLinearRestriction: Reverting CollectionMatrix back to CRS matrix List_ToCRSMatrix: Number of entries in CRS matrix: 527737 ListToCRSMatrix: Matrix format changed from List to CRS SolveWithLinearRestriction: CollectionMatrix done SolveWithLinearRestriction: CollectionVector done SolveWithLinearRestriction: Now going for the coupled linear system IterSolver: Using iterative method: bicgstab List_ToCRSMatrix: Number of entries in CRS matrix: 42267 ListToCRSMatrix: Matrix format changed from List to CRS CRS_IncompleteLU: ILU(0) (Real), Starting Factorization: CRS_IncompleteLU: Allocated LU matrix of size: 42267 CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 42268 CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100 CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.00 1 0.2665E+00 2 0.1247E+00 3 0.7439E-01 4 0.4329E-01 5 0.2756E-01 6 0.2064E-01 7 0.1632E-01 8 0.1334E-01 9 0.1110E-01 10 0.9602E-02 11 0.8203E-02 12 0.7163E-02 13 0.6078E-02 14 0.5192E-02 15 0.4384E-02 16 0.3749E-02 17 0.3231E-02 18 0.2822E-02 19 0.2475E-02 20 0.2194E-02 21 0.1927E-02 22 0.1715E-02 23 0.1530E-02 24 0.1374E-02 25 0.1254E-02 26 0.1156E-02 27 0.1063E-02 28 0.9774E-03 29 0.8970E-03 30 0.8220E-03 31 0.7548E-03 32 0.6908E-03 33 0.6338E-03 34 0.5810E-03 35 0.5384E-03 36 0.4997E-03 37 0.4650E-03 38 0.4320E-03 39 0.4037E-03 40 0.3773E-03 41 0.3532E-03 42 0.3278E-03 43 0.3043E-03 44 0.2814E-03 45 0.2622E-03 46 0.2452E-03 47 0.2306E-03 48 0.2167E-03 49 0.2045E-03 50 0.1934E-03 51 0.1842E-03 52 0.1734E-03 53 0.1602E-03 54 0.1609E-03 55 0.1390E-03 56 0.1380E-03 57 0.1356E-03 58 0.1303E-03 59 0.1191E-03 60 0.1177E-03 61 0.1160E-03 62 0.1247E-03 63 0.8974E-04 64 0.7572E-04 65 0.7151E-04 66 0.6573E-04 67 0.6466E-04 68 0.1405E-03 69 0.1453E-03 70 0.5895E-04 71 0.5733E-04 72 0.5516E-04 73 0.5429E-04 74 0.5712E-04 75 0.5262E-04 76 0.5037E-04 77 0.4890E-04 78 0.4540E-04 79 0.4466E-04 80 0.5102E-04 81 0.3551E-04 82 0.3031E-04 83 0.4503E-04 84 0.2753E-04 85 0.2646E-04 86 0.2409E-04 87 0.2376E-04 88 0.2423E-04 89 0.2240E-04 90 0.2073E-04 91 0.2017E-04 92 0.1906E-04 93 0.1829E-04 94 0.1748E-04 95 0.1682E-04 96 0.1529E-04 97 0.1513E-04 98 0.1496E-04 99 0.1338E-04 100 0.1366E-04 101 0.1242E-04 102 0.1130E-04 103 0.1128E-04 104 0.9916E-05 105 0.1001E-04 106 0.9246E-05 107 0.8711E-05 108 0.7881E-05 109 0.6670E-04 110 0.6949E-05 111 0.5884E-05 112 0.6119E-05 113 0.5597E-05 114 0.5371E-05 115 0.5012E-05 116 0.4718E-05 117 0.4847E-05 118 0.4215E-05 119 0.4136E-05 120 0.3829E-05 121 0.3811E-05 122 0.3736E-05 123 0.3565E-05 124 0.3822E-05 125 0.3438E-05 126 0.3390E-05 127 0.3288E-05 128 0.3247E-05 129 0.3158E-05 130 0.3158E-05 131 0.3056E-05 132 0.2993E-05 133 0.2811E-05 134 0.2808E-05 135 0.2677E-05 136 0.2447E-05 137 0.2322E-05 138 0.2351E-05 139 0.2126E-05 140 0.2076E-05 141 0.1938E-05 142 0.1962E-05 143 0.1820E-05 144 0.1764E-05 145 0.1770E-05 146 0.1613E-05 147 0.1455E-05 148 0.1347E-05 149 0.1356E-05 150 0.1267E-05 151 0.1195E-05 152 0.1085E-05 153 0.1015E-05 154 0.9766E-06 155 0.8302E-06 156 0.8191E-06 157 0.8178E-06 158 0.8079E-06 159 0.1635E-03 160 0.3375E-06 161 0.2968E-05 162 0.2079E-05 163 0.1634E-06 164 0.1538E-06 165 0.1289E-06 166 0.1051E-06 167 0.1046E-06 168 0.1016E-06 169 0.1038E-06 170 0.7691E-07 171 0.7184E-07 172 0.6002E-07 173 0.5258E-07 174 0.4165E-07 175 0.4158E-07 176 0.2818E-07 177 0.2744E-07 178 0.2540E-07 179 0.1946E-07 180 0.1927E-07 181 0.1511E-07 182 0.1518E-07 183 0.1163E-07 184 0.1057E-07 185 0.7606E-08 186 0.7776E-08 187 0.6547E-08 188 0.6566E-08 189 0.4831E-08 190 0.3750E-08 191 0.3278E-08 192 0.2560E-08 193 0.1949E-08 194 0.1945E-08 195 0.1250E-08 196 0.1302E-08 197 0.1100E-08 198 0.9174E-09 199 0.7187E-09 200 0.7201E-09 201 0.5742E-09 202 0.5292E-09 203 0.4555E-09 204 0.4183E-09 205 0.3836E-09 206 0.4036E-09 207 0.3608E-09 208 0.2321E-09 209 0.2072E-09 210 0.5873E-09 211 0.1532E-09 212 0.1397E-09 213 0.1672E-09 214 0.1081E-09 215 0.9254E-10 215 0.9254E-10 SolveWithLinearRestriction: Picking solution from collection solution SolveWithLinearRestriction: All done DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions SetNodalLoads: Checking for nodal loads for variable: p SetDirichletBoundaries: Number of dofs set: 0 DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set Program received signal SIGSEGV: Segmentation fault - invalid memory reference. Backtrace for this error: #0 0x7F6482045E08 #1 0x7F6482044F90 #2 0x7F6481C964AF #3 0x7F6472CD0B04 #4 0x7F6472CD217A #5 0x7F6472CE51D7 #6 0x7F648252691E #7 0x7F64825371E1 #8 0x7F6482538331 #9 0x7F648253AABC #10 0x7F64826DBEF2 #11 0x7F64826E1007 #12 0x401075 in MAIN__ at Solver.F90:69 Segmentation fault