I tried running the simulation with the given .sif file but my simulation diverged. Are there any setting that needs to be set outside of the .sif file?. Can you please have a look at my solver log.
Solver log:
Code: Select all
ELMER SOLVER (v 9.0) STARTED AT: 2021/03/07 20:50:38
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: 9.0 (Rev: Release, Compiled: 2020-11-10)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: Lua interpreted linked in.
MAIN: =============================================================
LoadInputFile: Reading only "Run Control" section
MAIN:
MAIN:
MAIN: -------------------------------------
MAIN: Reading Model: case.sif
LoadInputFile: Scanning input file: case.sif
LoadInputFile: Scanning only size info
LoadInputFile: First time visiting
LoadInputFile: Reading base load of sif file
LoadInputFile: Loading input file: case.sif
LoadInputFile: Reading base load of sif file
CheckKeyword: Unlisted keyword: [vtu: save bulk only] in section: [simulation]
CheckKeyword: Unlisted keyword: [vtu: discontinuous bodies] in section: [simulation]
LoadInputFile: Number of BCs: 1
LoadInputFile: Number of Body Forces: 1
LoadInputFile: Number of Initial Conditions: 0
LoadInputFile: Number of Materials: 3
LoadInputFile: Number of Equations: 2
LoadInputFile: Number of Solvers: 3
LoadInputFile: Number of Bodies: 3
WARNING:: LoadInputFile: > Material 2 < not used in any Body!
Loading user function library: [CoilSolver]...[CoilSolver_Init0]
Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init0]
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init0]
MagnetoDynamicsCalcFields_Init0: The target solver seems to be real valued
MagnetoDynamicsCalcFields_Init0: Target Variable Solver Index: 2
LoadMesh: Starting
ElmerAsciiMesh: Performing step: 1
ElmerAsciiMesh: Base mesh name: ./.
ReadHeaderFile: Reading header info from file: ././mesh.header
InitializeMesh: Number of nodes in mesh: 83561
InitializeMesh: Number of bulk elements in mesh: 466992
InitializeMesh: Number of boundary elements in mesh: 50256
InitializeMesh: Initial number of max element nodes: 4
ElmerAsciiMesh: Performing step: 2
ReadNodesFile: Reading nodes from file: ././mesh.nodes
MapCoordinates: Performing coordinate mapping
SetMeshDimension: Dimension of mesh is: 3
SetMeshDimension: Max dimension of mesh is: 3
MapCoordinates: Scaling coordinates: 1.000E-03 1.000E-03 1.000E-03
ElmerAsciiMesh: Performing step: 3
ReadElementsFile: Reading bulk elements from file: ././mesh.elements
ElmerAsciiMesh: Performing step: 4
ReadBoundaryFile: Reading boundary elements from file: ././mesh.boundary
PermuteNodeNumbering: Performing node mapping
MapBodiesAndBCs: Remapping bodies
MapBodiesAndBCs: Minimum initial body index: 1
MapBodiesAndBCs: Maximum initial body index: 3
MapBodiesAndBCs: Remapping boundaries
MapBodiesAndBCs: Minimum initial boundary index: 1
MapBodiesAndBCs: Maximum initial boundary index: 3
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: Finding mesh faces in 3D mesh
FindMeshFaces3D: Number of faces found: 951640
FindMeshFaces3D: All done
FindMeshEdges: Determining edges in 3D mesh
FindMeshEdges3D: Finding mesh edges in 3D mesh
FindMeshEdges3D: Number of edges found: 568208
FindMeshEdges3D: All done
LoadMesh: Preparing mesh done
LoadMesh: Elapsed REAL time: 5.6750 (s)
MeshStabParams: Computing stabilization parameters
MeshStabParams: Elapsed REAL time: 0.2350 (s)
CompleteModelKeywords: Completing keywords for mortar BCs
MAIN: -------------------------------------
AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: case
AddVtuOutputSolverHack: Increasing number of solver to: 4
ListCopyPrefixedKeywords: Prefix: vtu: Keyword: save bulk only
ListCopyPrefixedKeywords: Prefix: vtu: Keyword: discontinuous bodies
ListCopyPrefixedKeywords: Copied 2 keywords with prefix: vtu:
AddVtuOutputSolverHack: Finished appending VTU output solver
AddSolvers: Setting up 4 solvers
AddSolvers: Setting up solver 1: coilsolver
AddEquationBasics: Using procedure: CoilSolver CoilSolver
AddEquationBasics: Setting up solver: coilsolver
Loading user function library: [CoilSolver]...[CoilSolver_Init]
AddEquationBasics: Checking for _init solver
Loading user function library: [CoilSolver]...[CoilSolver_bulk]
AddEquationBasics: Checking for _bulk solver
Loading user function library: [CoilSolver]...[CoilSolver]
AddEquationBasics: Creating standard variable: coiltmp
AddEquationBasics: Computing size of permutation vector
AddEquationBasics: Maximum size of permutation vector is: 651769
AddEquationBasics: Creating solver matrix topology
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: coilsolver...done.
OptimizeBandwidth: Half bandwidth without optimization: 7901
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 579
OptimizeBandwidth: ---------------------------------------------------------
CRS_CreateMatrix: Creating CRS Matrix of size: 13261
AddEquationBasics: Number of rows in CRS matrix: 13261
AddEquationBasics: Creating solver variable
AddEquationBasics: Creating exported variable: coilpot
AddEquationBasics: Succesfully created variable: coilpot
AddEquationBasics: Creating exported variable: coilpotb
AddEquationBasics: Succesfully created variable: coilpotb
AddEquationBasics: Creating exported variable: potselect
AddEquationBasics: Succesfully created variable: potselect
AddEquationBasics: Creating exported variable: -dg coilcurrent e[coilcurrent e:3]
CreateDGPerm: Creating permutation for DG variable
CreateDGPerm: Size of DgPerm table: 1867968
CreateDGPerm: Created permutation for DG nodes: 234212
AddEquationBasics: Succesfully created variable: coilcurrent e
AddEquationBasics: Succesfully created variable: coilcurrent e 1
AddEquationBasics: Succesfully created variable: coilcurrent e 2
AddEquationBasics: Succesfully created variable: coilcurrent e 3
AddEquationBasics: Creating exported variable: coilcurrent[coilcurrent:3]
AddEquationBasics: Succesfully created variable: coilcurrent
AddEquationBasics: Succesfully created variable: coilcurrent 1
AddEquationBasics: Succesfully created variable: coilcurrent 2
AddEquationBasics: Succesfully created variable: coilcurrent 3
AddSolvers: Setting up solver 2: mgdynamics
AddEquationBasics: Using procedure: MagnetoDynamics WhitneyAVSolver
AddEquationBasics: Setting up solver: mgdynamics
Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_Init]
AddEquationBasics: Checking for _init solver
Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver_bulk]
AddEquationBasics: Checking for _bulk solver
Loading user function library: [MagnetoDynamics]...[WhitneyAVSolver]
AddEquationBasics: Creating standard variable: av
AddEquationBasics: Computing size of permutation vector
AddEquationBasics: Maximum size of permutation vector is: 651769
AddEquationBasics: Creating solver matrix topology
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: mgdynamics...done.
OptimizeBandwidth: Half bandwidth without optimization: 558599
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 25081
OptimizeBandwidth: ---------------------------------------------------------
CRS_CreateMatrix: Creating CRS Matrix of size: 568208
AddEquationBasics: Number of rows in CRS matrix: 568208
AddEquationBasics: Creating solver variable
AddSolvers: Setting up solver 3: mgdynamicscalc
AddEquationBasics: Using procedure: MagnetoDynamics MagnetoDynamicsCalcFields
AddEquationBasics: Setting up solver: mgdynamicscalc
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_Init]
AddEquationBasics: Checking for _init solver
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields_bulk]
AddEquationBasics: Checking for _bulk solver
Loading user function library: [MagnetoDynamics]...[MagnetoDynamicsCalcFields]
AddEquationBasics: Creating standard variable: hr_dummy
AddEquationBasics: Computing size of permutation vector
AddEquationBasics: Maximum size of permutation vector is: 651769
AddEquationBasics: Creating solver matrix topology
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: mgdynamicscalc...done.
OptimizeBandwidth: Half bandwidth without optimization: 77605
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 5361
OptimizeBandwidth: ---------------------------------------------------------
CRS_CreateMatrix: Creating CRS Matrix of size: 77774
AddEquationBasics: Number of rows in CRS matrix: 77774
AddEquationBasics: Creating solver variable
AddEquationBasics: Creating exported variable: magnetic flux density[magnetic flux density:3]
AddEquationBasics: Succesfully created variable: magnetic flux density
AddEquationBasics: Succesfully created variable: magnetic flux density 1
AddEquationBasics: Succesfully created variable: magnetic flux density 2
AddEquationBasics: Succesfully created variable: magnetic flux density 3
AddEquationBasics: Creating exported variable: magnetic field strength[magnetic field strength:3]
AddEquationBasics: Succesfully created variable: magnetic field strength
AddEquationBasics: Succesfully created variable: magnetic field strength 1
AddEquationBasics: Succesfully created variable: magnetic field strength 2
AddEquationBasics: Succesfully created variable: magnetic field strength 3
AddEquationBasics: Creating exported variable: maxwell stress[maxwell stress:6]
AddEquationBasics: Succesfully created variable: maxwell stress
AddEquationBasics: Succesfully created variable: maxwell stress 1
AddEquationBasics: Succesfully created variable: maxwell stress 2
AddEquationBasics: Succesfully created variable: maxwell stress 3
AddEquationBasics: Succesfully created variable: maxwell stress 4
AddEquationBasics: Succesfully created variable: maxwell stress 5
AddEquationBasics: Succesfully created variable: maxwell stress 6
AddEquationBasics: Creating exported variable: current density[current density:3]
AddEquationBasics: Succesfully created variable: current density
AddEquationBasics: Succesfully created variable: current density 1
AddEquationBasics: Succesfully created variable: current density 2
AddEquationBasics: Succesfully created variable: current density 3
AddEquationBasics: Creating exported variable: -dg magnetic flux density e[magnetic flux density e:3]
CreateDGPerm: Creating permutation for DG variable
CreateDGPerm: Size of DgPerm table: 1867968
CreateDGPerm: Created permutation for DG nodes: 1633756
AddEquationBasics: Succesfully created variable: magnetic flux density e
AddEquationBasics: Succesfully created variable: magnetic flux density e 1
AddEquationBasics: Succesfully created variable: magnetic flux density e 2
AddEquationBasics: Succesfully created variable: magnetic flux density e 3
AddEquationBasics: Creating exported variable: -dg magnetic field strength e[magnetic field strength e:3]
CreateDGPerm: Creating permutation for DG variable
CreateDGPerm: Size of DgPerm table: 1867968
CreateDGPerm: Created permutation for DG nodes: 1633756
AddEquationBasics: Succesfully created variable: magnetic field strength e
AddEquationBasics: Succesfully created variable: magnetic field strength e 1
AddEquationBasics: Succesfully created variable: magnetic field strength e 2
AddEquationBasics: Succesfully created variable: magnetic field strength e 3
AddEquationBasics: Creating exported variable: -dg maxwell stress e[maxwell stress e:6]
CreateDGPerm: Creating permutation for DG variable
CreateDGPerm: Size of DgPerm table: 1867968
CreateDGPerm: Created permutation for DG nodes: 1633756
AddEquationBasics: Succesfully created variable: maxwell stress e
AddEquationBasics: Succesfully created variable: maxwell stress e 1
AddEquationBasics: Succesfully created variable: maxwell stress e 2
AddEquationBasics: Succesfully created variable: maxwell stress e 3
AddEquationBasics: Succesfully created variable: maxwell stress e 4
AddEquationBasics: Succesfully created variable: maxwell stress e 5
AddEquationBasics: Succesfully created variable: maxwell stress e 6
AddEquationBasics: Creating exported variable: -dg current density e[current density e:3]
CreateDGPerm: Creating permutation for DG variable
CreateDGPerm: Size of DgPerm table: 1867968
CreateDGPerm: Created permutation for DG nodes: 1633756
AddEquationBasics: Succesfully created variable: current density e
AddEquationBasics: Succesfully created variable: current density e 1
AddEquationBasics: Succesfully created variable: current density e 2
AddEquationBasics: Succesfully created variable: current density e 3
AddSolvers: Setting up solver 4: 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]
AddSolvers: Setting up solvers done
AddMeshCoordinatesAndTime: Setting mesh coordinates and time
SetInitialConditions: Setting up initial conditions (if any)
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: coilsolver
SetActiveElementsTable: Number of active elements found : 58553
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: coilsolver
CoilSolver: --------------------------------------
CoilSolver: Solving current distribution in a coil
CoilSolver: --------------------------------------
CoilSolver: Treating coil in Component: 1
CoilSolver: Coil volume: 9.5686E-11
CoilSolver: Coil center: 2.4788E-11 -1.2958E-10 5.0000E-04
DefineCoilParametes: Fitting the coil by maximizing inertia
CoilSolver: Coil inertia eigenvalues: -1.0000E+00 -1.0000E+00 -1.0000E+00
CoilSolver: Coil axis normal: -5.9976E-02 -9.9820E-01 0.0000E+00
CoilSolver: Coil tangent1: 9.9820E-01 -5.9976E-02 0.0000E+00
CoilSolver: Coil tangent2: 0.0000E+00 0.0000E+00 1.0000E+00
CoilSolver: Choosing fixing nodes for set: 1
CoilSolver: Choosing fixing nodes for set: 2
CoilSolver: Coil system consists of 1 coils
CoilSolver: Set1 : 662 +nodes and 658 -nodes
CoilSolver: Set2 : 658 +nodes and 662 -nodes
CoilSolver: Computing the dummy potential field
InitializeToZero: Initializing the linear system to zero
DefaultFinishBulkAssembly: Saving bulk values for: coilsolver
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
SetNodalLoads: Checking for nodal loads for variable: coiltmp
SetNodalLoads: Finished checking for nodal loads
SetDirichletBoundaries: Number of dofs set for coiltmp: 0
EnforceDirichletConditions: Applying Dirichlet conditions using scaled diagonal
ScaleLinearSystem: Scaling diagonal entries to unity
ScaleLinearSystem: Computing > DiagScaling < vector
EnforceDirichletConditions: Dirichlet boundary conditions enforced
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Number of nonzeros: 348891
CRS_IncompleLU: Recursive round: 1
CRS_IncompleteLU: Number of nonzeros: 850521
CRS_IncompleteLU: Allocated LU matrix of size: 592635
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.10
1 0.5158E-02 0.5158E-02
2 0.3192E-03 0.3192E-03
3 0.1364E-03 0.1364E-03
4 0.3459E-04 0.3459E-04
SolveLinearSystem: Calculating nodal loads
ComputeNorm: Computing norm of solution
ComputeChange: NS (ITER=1) (NRM,RELC): ( 0.50146751 2.0000000 ) :: coilsolver
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Positive coil currents: 365 4.7527E-05
CoilSolver: Negative coil currents: 360 -4.7801E-05
CoilSolver: Discrepancy of start and end coil currents: 5.7490E-03
CoilSolver: Initial coil current for coil 1: 4.7664E-05
CoilSolver: Coil potential multiplier: 2.0980E+04
InitializeToZero: Initializing the linear system to zero
DefaultFinishBulkAssembly: Saving bulk values for: coilsolver
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
SetNodalLoads: Checking for nodal loads for variable: coiltmp
SetNodalLoads: Finished checking for nodal loads
SetDirichletBoundaries: Number of dofs set for coiltmp: 0
EnforceDirichletConditions: Applying Dirichlet conditions using scaled diagonal
ScaleLinearSystem: Scaling diagonal entries to unity
ScaleLinearSystem: Computing > DiagScaling < vector
EnforceDirichletConditions: Dirichlet boundary conditions enforced
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.1403E-02 0.1403E-02
2 0.1153E-03 0.1153E-03
3 0.6762E-04 0.6762E-04
SolveLinearSystem: Calculating nodal loads
ComputeNorm: Computing norm of solution
ComputeChange: NS (ITER=2) (NRM,RELC): ( 0.49766008 0.76215044E-02 ) :: coilsolver
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Positive coil currents: 360 4.9401E-05
CoilSolver: Negative coil currents: 365 -4.9292E-05
CoilSolver: Discrepancy of start and end coil currents: 2.2100E-03
CoilSolver: Initial coil current for coil 1: 4.9346E-05
CoilSolver: Coil potential multiplier: 2.0265E+04
CoilSolver: Computing current component: 1
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.7517E-06 0.7517E-06
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Computing current component: 2
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.2485E-05 0.2485E-05
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Computing current component: 3
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.8996E-06 0.8996E-06
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Minimum current components: -2.3752E+08 -2.9558E+08 -3.4483E+08
CoilSolver: Maximum current components: 2.5393E+08 1.4103E+08 6.3508E+08
CoilSolver: Computing current component: 1
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.1261E-05 0.1261E-05
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Computing current component: 2
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.2485E-05 0.2485E-05
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Computing current component: 3
InitializeToZero: Initializing the linear system to zero
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
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
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu2
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(2) (Real), Performing Factorization:
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(2) (Real), NOF nonzeros: 592636
CRS_IncompleteLU: ILU(2) (Real), filling (%) : 344
CRS_IncompleteLU: ILU(2) (Real), Factorization ready at (s): 0.03
1 0.8996E-06 0.8996E-06
SolveSystem: System solved
SolveSystem: Finished solving the system
CoilSolver: Minimum current components: -2.3752E+08 -2.9558E+08 -3.4483E+08
CoilSolver: Maximum current components: 2.5393E+08 1.4103E+08 6.3508E+08
CoilSolver: Normalizing current density to a constant value
CoilSolver: All done
CoilSolver: --------------------------------------
Loading user function library: [CoilSolver]...[CoilSolver_post]
ComputeNorm: Computing norm of solution
ComputeChange: SS (ITER=1) (NRM,RELC): ( 10520.825 2.0000000 ) :: coilsolver
SetActiveElementsTable: Creating active element table for: mgdynamics
SetActiveElementsTable: Number of active elements found : 466992
SingleSolver: Attempting to call solver
SingleSolver: Solver Equation string is: mgdynamics
WhitneyAVSolver: -------------------------------------------
WhitneyAVSolver: Solving the AV equations with edge elements
JfixPotentialSolver: Computing fixing potential for given current density
ListPushNameSpace: Adding name space: jfix:
OptimizeBandwidth: ---------------------------------------------------------
OptimizeBandwidth: Computing matrix structure for: mgdynamics...done.
OptimizeBandwidth: Half bandwidth without optimization: 82683
OptimizeBandwidth:
OptimizeBandwidth: Bandwidth Optimization ...done.
OptimizeBandwidth: Half bandwidth after optimization: 6645
OptimizeBandwidth: ---------------------------------------------------------
CRS_CreateMatrix: Creating CRS Matrix of size: 83561
ListPopNameSpace: Deleting entry from name space:
JfixPotentialSolver: Finished creating matrix equation
InitializeToZero: Initializing the linear system to zero
MarkOuterNodes: Marking outer nodes on outer boundary
MarkOuterNodes: Total number of surface nodes: 17658
DefaultStart: Starting solver: mgdynamics
InitializeToZero: Initializing the linear system to zero
WhitneyAVSolver: Solving the fixing potential
JfixPotentialSolver: Computing fixing potential for given current density
ListPushNameSpace: Adding name space: jfix:
JfixPotentialSolver: Setting all out-of-plane source nodes to zero
JfixBCs: Maximum source term on boundaries: 0.000E+00
JfixBCs: Using jfix epsilon for flux: 2.220E-16
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
SetNodalLoads: Checking for nodal loads for variable: jfix
SetNodalLoads: Finished checking for nodal loads
SetDirichletBoundaries: Number of dofs set for jfix: 0
EnforceDirichletConditions: No Dirichlet conditions to enforce, exiting!
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
WARNING:: JfixPotentialSolver: No Dirichlet conditions used to define Jfix level!
JfixPotentialSolver: Solving for Jfix
SolveSystem: Solving linear system
SolveSystem: Solving linear system without constraint matrix
SolveLinearSystem: Assuming real valued linear system
RowEquilibration: Unscaled matrix norm: 9.2567977426944725E-003
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: ilu
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstab
IterSolver: Matrix is real valued
CRS_IncompleteLU: ILU(0) (Real), Performing Factorization:
CRS_IncompleteLU: Allocated LU matrix of size: 1219977
CRS_IncompleteLU: Performing incomplete LU
CRS_IncompleteLU: ILU(0) (Real), NOF nonzeros: 1219978
CRS_IncompleteLU: ILU(0) (Real), filling (%) : 100
CRS_IncompleteLU: ILU(0) (Real), Factorization ready at (s): 0.05
20 0.6278E-04
40 0.5604E-05
60 0.1341E-07
72 0.5870E-11
SolveSystem: System solved
SolveSystem: Finished solving the system
JfixPotentialSolver: Norm for Jfix computation: 4.287E-02
ListPopNameSpace: Deleting entry from name space:
WhitneyAVSolver: Adding the fixing potential to the r.h.s. of AV equation
WhitneyAVSolver: Finished adding the fixing potential
MGDynAssembly: Elapsed REAL time: 17.7730 (s)
DefUtils::DefaultDirichletBCs: Setting Dirichlet boundary conditions
SetNodalLoads: Checking for nodal loads for variable: av
SetNodalLoads: Finished checking for nodal loads
SetDirichletBoundaries: Number of dofs set for av: 0
EnforceDirichletConditions: Applying Dirichlet conditions using scaled diagonal
ScaleLinearSystem: Scaling diagonal entries to unity
ScaleLinearSystem: Computing > DiagScaling < vector
EnforceDirichletConditions: Dirichlet boundary conditions enforced
DefUtils::DefaultDirichletBCs: Dirichlet boundary conditions set
DefaultSolve: Solving linear system with default routines
GenerateConstraintMatrix: Nothing to do for now
SolveSystem: Solving linear system
SolveSystem: Solving linear system without constraint matrix
SolveLinearSystem: Assuming real valued linear system
RowEquilibration: Unscaled matrix norm: 7020936616764.4375
SolveLinearSystem: Linear System Solver: iterative
SolveLinearSystem: Linear System Preconditioning: none
SolveLinearSystem: Serial linear System Solver: iterative
IterSolver: Using iterative method: bicgstabl
IterSolver: Matrix is real valued
1 0.5943E-08 0.1069E+12
2 0.1154E-08 0.2075E+11
3 0.3969E-09 0.7138E+10
4 0.8903E-09 0.1601E+11
5 0.1024E-09 0.1842E+10
6 0.6924E-10 0.1245E+10
7 0.6456E-10 0.1161E+10
8 0.2737E-10 0.4922E+09
9 0.2131E-10 0.3832E+09
10 0.3867E-10 0.6954E+09
11 0.1393E-09 0.2504E+10
12 0.9767E-10 0.1756E+10
13 0.7941E-10 0.1428E+10
14 0.6952E-10 0.1250E+10
15 0.4786E-10 0.8607E+09
16 0.5440E-11 0.9783E+08
17 0.3613E-11 0.6498E+08
18 0.5010E-11 0.9009E+08
19 0.8447E-12 0.1519E+08
20 0.7908E-12 0.1422E+08
21 0.6967E-12 0.1253E+08
22 0.5773E-12 0.1038E+08
23 0.4413E-12 0.7937E+07
24 0.3729E-12 0.6707E+07
25 0.3418E-12 0.6147E+07
26 0.2935E-12 0.5278E+07
27 0.2942E-12 0.5290E+07
28 0.2309E-12 0.4153E+07
29 0.1766E-12 0.3176E+07
30 0.1270E-12 0.2284E+07
31 0.9861E-13 0.1773E+07
32 0.8380E-13 0.1507E+07
33 0.5549E-13 0.9979E+06
34 0.6571E-14 0.1182E+06
35 0.1203E-13 0.2163E+06
36 0.1102E-13 0.1981E+06
37 0.2513E-14 0.4519E+05
38 0.3196E-14 0.5747E+05
39 0.9869E-15 0.1775E+05
40 0.1223E-14 0.2199E+05
41 0.4923E-15 0.8853E+04
42 0.5348E-15 0.9618E+04
43 0.4994E-15 0.8981E+04
44 0.4678E-15 0.8413E+04
45 0.8691E-15 0.1563E+05
46 0.1820E-15 0.3274E+04
47 0.1530E-15 0.2752E+04
48 0.7550E-16 0.1358E+04
49 0.3089E-16 0.5555E+03
50 0.2667E-16 0.4796E+03
51 0.2113E-16 0.3800E+03
52 0.3191E-16 0.5739E+03
53 0.4357E-16 0.7835E+03
54 0.1461E-16 0.2627E+03
55 0.1152E-16 0.2072E+03
56 0.1013E-16 0.1821E+03
57 0.9438E-17 0.1697E+03
58 0.9073E-17 0.1632E+03
59 0.4052E-17 0.7286E+02
60 0.4823E-17 0.8674E+02
61 0.1273E-16 0.2289E+03
62 0.7897E-18 0.1420E+02
63 0.8400E-18 0.1511E+02
64 0.8278E-18 0.1489E+02
65 0.5508E-18 0.9906E+01
66 0.4819E-18 0.8667E+01
67 0.4779E-18 0.8595E+01
68 0.4614E-18 0.8297E+01
69 0.4317E-18 0.7763E+01
70 0.3766E-17 0.6773E+02
71 0.1417E-18 0.2547E+01
72 0.1709E-18 0.3073E+01
73 0.1889E-17 0.3397E+02
74 0.2099E-18 0.3774E+01
75 0.4627E-18 0.8320E+01
76 0.3930E-18 0.7067E+01
77 0.5143E-18 0.9248E+01
78 0.5663E-18 0.1018E+02
79 0.1955E-18 0.3515E+01
80 0.5660E-19 0.1018E+01
81 0.1535E-16 0.2761E+03
82 0.6443E-17 0.1159E+03
83 0.1448E-16 0.2605E+03
84 0.6538E-17 0.1176E+03
85 0.2222E-17 0.3996E+02
86 0.1500E-17 0.2698E+02
87 0.1356E-17 0.2439E+02
88 0.1259E-17 0.2265E+02
89 0.1257E-17 0.2261E+02
90 0.1112E-17 0.1999E+02
91 0.1471E-17 0.2645E+02
92 0.1257E-17 0.2260E+02
93 0.1272E-17 0.2287E+02
94 0.1572E-17 0.2828E+02
95 0.2084E-17 0.3747E+02
96 0.1394E-17 0.2508E+02
97 0.1391E-17 0.2502E+02
98 0.3411E-17 0.6135E+02
99 0.1750E-16 0.3148E+03
100 0.7976E-17 0.1434E+03
101 0.3049E-17 0.5483E+02
102 0.3788E-17 0.6813E+02
103 0.3167E-17 0.5695E+02
104 0.3729E-17 0.6706E+02
105 0.4355E-18 0.7832E+01
106 0.5452E-17 0.9804E+02
107 0.1922E-15 0.3457E+04
108 0.5235E-17 0.9415E+02
109 0.1685E-16 0.3031E+03
110 0.2571E-16 0.4623E+03
111 0.2680E-16 0.4820E+03
112 0.2659E-16 0.4783E+03
113 0.2229E-16 0.4008E+03
114 0.4729E-17 0.8505E+02
115 0.6809E-17 0.1225E+03
116 0.6222E-17 0.1119E+03
117 0.1261E-16 0.2267E+03
118 0.6036E-17 0.1085E+03
119 0.5372E-17 0.9660E+02
120 0.6167E-17 0.1109E+03
121 0.5623E-17 0.1011E+03
122 0.1056E-16 0.1898E+03
123 0.5944E-17 0.1069E+03
124 0.6734E-17 0.1211E+03
125 0.7265E-17 0.1306E+03
126 0.7000E-17 0.1259E+03
127 0.7110E-17 0.1279E+03
128 0.6905E-17 0.1242E+03
129 0.7405E-17 0.1332E+03
130 0.1152E-16 0.2072E+03
131 0.1650E-16 0.2967E+03
132 0.8850E-17 0.1592E+03
133 0.1301E-16 0.2340E+03
134 0.1513E-16 0.2721E+03
135 0.5559E-17 0.9998E+02
136 0.8295E-17 0.1492E+03
137 0.5513E-17 0.9914E+02
138 0.5587E-17 0.1005E+03
139 0.5981E-17 0.1076E+03
140 0.5796E-17 0.1042E+03
141 0.3898E-17 0.7009E+02
142 0.2832E-17 0.5094E+02
143 0.2285E-17 0.4110E+02
144 0.4499E-17 0.8090E+02
145 0.5449E-17 0.9800E+02
146 0.7832E-17 0.1409E+03
147 0.5714E-17 0.1028E+03
148 0.4106E-17 0.7384E+02
149 0.3341E-17 0.6008E+02
150 0.3698E-17 0.6650E+02
151 0.3089E-17 0.5554E+02
152 0.4997E-17 0.8987E+02
153 0.4397E-17 0.7907E+02
154 0.4820E-17 0.8668E+02
155 0.4259E-17 0.7660E+02
156 0.6510E-17 0.1171E+03
157 0.5368E-17 0.9653E+02
158 0.1116E-17 0.2007E+02
159 0.5829E-17 0.1048E+03
160 0.3816E-17 0.6863E+02
161 0.3563E-17 0.6407E+02
162 0.3688E-17 0.6633E+02
163 0.1536E-17 0.2762E+02
164 0.1720E-18 0.3093E+01
165 0.9107E-18 0.1638E+02
166 0.2938E-17 0.5283E+02
167 0.1528E-17 0.2748E+02
168 0.2955E-17 0.5315E+02
169 0.4167E-16 0.7493E+03
170 0.5865E-17 0.1055E+03
171 0.1738E-16 0.3125E+03
172 0.8458E-16 0.1521E+04
173 0.1567E-15 0.2817E+04
174 0.1409E-15 0.2533E+04
175 0.1754E-15 0.3155E+04
176 0.2615E-16 0.4703E+03
177 0.3375E-15 0.6070E+04
178 0.9327E-14 0.1677E+06
179 0.2069E-15 0.3721E+04
180 0.4704E-14 0.8460E+05
181 0.4731E-14 0.8508E+05
182 0.2704E-14 0.4863E+05
183 0.2135E-14 0.3840E+05
184 0.1512E-13 0.2719E+06
185 0.6151E-14 0.1106E+06
186 0.5656E-14 0.1017E+06
187 0.7829E-14 0.1408E+06
188 0.7106E-14 0.1278E+06
189 0.1578E-14 0.2837E+05
190 0.7895E-14 0.1420E+06
191 0.6881E-14 0.1237E+06
192 0.7002E-14 0.1259E+06
193 0.6058E-14 0.1089E+06
194 0.5527E-14 0.9939E+05
195 0.5483E-14 0.9860E+05
196 0.5796E-14 0.1042E+06
197 0.6085E-14 0.1094E+06
198 0.1644E-13 0.2956E+06
199 0.1287E-14 0.2314E+05
200 0.6186E-14 0.1113E+06
201 0.1124E-13 0.2021E+06
202 0.1517E-13 0.2727E+06
203 0.3557E-13 0.6396E+06
204 0.1935E-13 0.3480E+06
205 0.2226E-13 0.4003E+06
206 0.2986E-13 0.5371E+06
207 0.3759E-13 0.6760E+06
208 0.3718E-13 0.6686E+06
209 0.3572E-13 0.6424E+06
210 0.2920E-13 0.5251E+06
211 0.2270E-13 0.4082E+06
212 0.2261E-13 0.4066E+06
213 0.1634E-13 0.2939E+06
214 0.1143E-13 0.2056E+06
215 0.1041E-13 0.1872E+06
216 0.1055E-13 0.1897E+06
217 0.1057E-13 0.1901E+06
218 0.1400E-13 0.2518E+06
219 0.8041E-14 0.1446E+06
220 0.7022E-14 0.1263E+06
221 0.4270E-14 0.7679E+05
222 0.4904E-14 0.8818E+05
223 0.3813E-14 0.6857E+05
224 0.2374E-14 0.4270E+05
225 0.1439E-14 0.2587E+05
226 0.1327E-14 0.2386E+05
227 0.1045E-14 0.1879E+05
228 0.8873E-15 0.1596E+05
229 0.8156E-15 0.1467E+05
230 0.1346E-14 0.2421E+05
231 0.8083E-15 0.1454E+05
232 0.8651E-15 0.1556E+05
233 0.9217E-15 0.1657E+05
234 0.1791E-14 0.3221E+05
235 0.8182E-16 0.1471E+04
236 0.5725E-15 0.1029E+05
237 0.6209E-15 0.1117E+05
238 0.2548E-15 0.4582E+04
239 0.1561E-15 0.2808E+04
240 0.9827E-16 0.1767E+04
241 0.1071E-15 0.1925E+04
242 0.1168E-15 0.2100E+04
243 0.1108E-15 0.1993E+04
244 0.1048E-15 0.1885E+04
245 0.7296E-16 0.1312E+04
246 0.7693E-16 0.1384E+04
247 0.7918E-16 0.1424E+04
248 0.7686E-16 0.1382E+04
249 0.7702E-16 0.1385E+04
250 0.9015E-16 0.1621E+04
251 0.7404E-16 0.1331E+04
252 0.7484E-16 0.1346E+04
253 0.7556E-16 0.1359E+04
254 0.7337E-16 0.1320E+04
255 0.7334E-16 0.1319E+04
256 0.1064E-15 0.1914E+04
257 0.2556E-15 0.4596E+04
258 0.1634E-15 0.2939E+04
259 0.8726E-16 0.1569E+04
260 0.8065E-16 0.1450E+04
261 0.7777E-16 0.1399E+04
262 0.7816E-16 0.1406E+04
263 0.7790E-16 0.1401E+04
264 0.7684E-16 0.1382E+04
265 0.7687E-16 0.1382E+04
266 0.7396E-16 0.1330E+04
267 0.7282E-16 0.1309E+04
268 0.8910E-16 0.1602E+04
269 0.1883E-15 0.3386E+04
270 0.3160E-15 0.5682E+04
271 0.2371E-14 0.4264E+05
272 0.1224E-15 0.2200E+04
273 0.1074E-14 0.1931E+05
274 0.2232E-14 0.4014E+05
275 0.1804E-14 0.3244E+05
276 0.9121E-15 0.1640E+05
277 0.3651E-13 0.6565E+06
278 0.1641E-14 0.2950E+05
279 0.3000E-13 0.5396E+06
280 0.1678E-12 0.3018E+07
281 0.2609E-13 0.4691E+06
282 0.1505E-12 0.2706E+07
283 0.4026E-12 0.7239E+07
284 0.3332E-12 0.5992E+07
285 0.3364E-12 0.6049E+07
286 0.1816E-12 0.3266E+07
287 0.9244E-13 0.1662E+07
288 0.8248E-13 0.1483E+07
289 0.1175E-12 0.2114E+07
290 0.1159E-12 0.2085E+07
291 0.2964E-13 0.5330E+06
292 0.2241E-13 0.4031E+06
293 0.3891E-13 0.6996E+06
294 0.3824E-13 0.6877E+06
295 0.3394E-13 0.6103E+06
296 0.4859E-13 0.8739E+06
297 0.2356E-13 0.4237E+06
298 0.2011E-13 0.3617E+06
299 0.8437E-14 0.1517E+06
300 0.5248E-14 0.9438E+05
301 0.5298E-14 0.9528E+05
302 0.4886E-14 0.8787E+05
303 0.5327E-14 0.9580E+05
304 0.5329E-14 0.9583E+05
305 0.5216E-14 0.9379E+05
306 0.5238E-14 0.9420E+05
307 0.5248E-14 0.9438E+05
308 0.1333E-13 0.2397E+06
309 0.5292E-14 0.9516E+05
310 0.5046E-14 0.9075E+05
311 0.4917E-14 0.8842E+05
312 0.5275E-14 0.9486E+05
313 0.5319E-14 0.9565E+05
314 0.5035E-14 0.9055E+05
315 0.5279E-14 0.9494E+05
316 0.5253E-14 0.9446E+05
317 0.5224E-14 0.9394E+05
318 0.5614E-14 0.1010E+06
319 0.5252E-14 0.9445E+05
320 0.5225E-14 0.9396E+05
321 0.5220E-14 0.9386E+05
322 0.5262E-14 0.9462E+05
323 0.5215E-14 0.9378E+05
324 0.4977E-14 0.8950E+05
325 0.5314E-14 0.9556E+05
326 0.5143E-14 0.9248E+05
327 0.5305E-14 0.9539E+05
328 0.5094E-14 0.9160E+05
329 0.5276E-14 0.9488E+05
330 0.5264E-14 0.9466E+05
331 0.5031E-14 0.9047E+05
332 0.5319E-14 0.9566E+05
333 0.4923E-14 0.8853E+05
334 0.5158E-14 0.9275E+05
335 0.5148E-14 0.9258E+05
336 0.5754E-14 0.1035E+06
337 0.5383E-14 0.9680E+05
338 0.5315E-14 0.9558E+05
339 0.5247E-14 0.9436E+05
340 0.5263E-14 0.9465E+05
341 0.5327E-14 0.9580E+05
342 0.5325E-14 0.9577E+05
343 0.4797E-14 0.8626E+05
344 0.5330E-14 0.9585E+05
345 0.1703E-12 0.3063E+07
346 0.5145E-14 0.9252E+05
347 0.5052E-14 0.9086E+05
348 0.5396E-14 0.9704E+05
349 0.5218E-14 0.9384E+05
350 0.5316E-14 0.9560E+05
351 0.5231E-14 0.9407E+05
352 0.5316E-14 0.9561E+05
353 0.5008E-14 0.9006E+05
354 0.5025E-14 0.9037E+05
355 0.5291E-14 0.9516E+05
356 0.5283E-14 0.9500E+05
357 0.6476E-14 0.1165E+06
358 0.5321E-14 0.9570E+05
359 0.5183E-14 0.9321E+05
360 0.8582E-14 0.1543E+06
361 0.1159E-13 0.2084E+06
362 0.6341E-14 0.1140E+06
363 0.6798E-14 0.1222E+06
364 0.8672E-14 0.1560E+06
365 0.8222E-14 0.1479E+06
366 0.7325E-14 0.1317E+06
367 0.1177E-13 0.2117E+06
368 0.9906E-14 0.1781E+06
369 0.6238E-14 0.1122E+06
370 0.5824E-14 0.1047E+06
371 0.8233E-14 0.1481E+06
372 0.5214E-14 0.9377E+05
373 0.5253E-14 0.9447E+05
374 0.5174E-14 0.9304E+05
375 0.5307E-14 0.9544E+05
376 0.5099E-14 0.9169E+05
377 0.5119E-14 0.9206E+05
378 0.5147E-14 0.9256E+05
379 0.5304E-14 0.9539E+05
380 0.5291E-14 0.9516E+05
381 0.5124E-14 0.9215E+05
382 0.4840E-14 0.8704E+05
383 0.5261E-14 0.9461E+05
384 0.1997E-13 0.3591E+06
385 0.5256E-14 0.9453E+05
386 0.5106E-14 0.9182E+05
387 0.5273E-14 0.9482E+05
388 0.5137E-14 0.9237E+05
389 0.5295E-14 0.9521E+05
390 0.5273E-14 0.9483E+05
391 0.5227E-14 0.9399E+05
392 0.5297E-14 0.9526E+05
393 0.5611E-14 0.1009E+06
394 0.5092E-14 0.9158E+05
395 0.5210E-14 0.9370E+05
396 0.5258E-14 0.9456E+05
397 0.8544E-14 0.1536E+06
398 0.5157E-14 0.9274E+05
399 0.4868E-14 0.8755E+05
400 0.5220E-14 0.9388E+05
401 0.5515E-14 0.9918E+05
402 0.5550E-14 0.9981E+05
403 0.5815E-14 0.1046E+06
404 0.5827E-14 0.1048E+06
405 0.5835E-14 0.1049E+06
406 0.5677E-14 0.1021E+06
407 0.5757E-14 0.1035E+06
408 0.5850E-14 0.1052E+06
409 0.5850E-14 0.1052E+06
410 0.5798E-14 0.1043E+06
411 0.5303E-14 0.9536E+05
412 0.5760E-14 0.1036E+06
413 0.5825E-14 0.1047E+06
414 0.5297E-14 0.9525E+05
415 0.8321E-14 0.1496E+06
416 0.5712E-14 0.1027E+06
417 0.6494E-14 0.1168E+06
418 0.6129E-14 0.1102E+06
419 0.3236E-13 0.5819E+06
420 0.1940E-13 0.3489E+06
421 0.2057E-13 0.3699E+06
422 0.2723E-13 0.4896E+06
423 0.2584E-13 0.4647E+06
424 0.2478E-13 0.4456E+06
425 0.2753E-13 0.4952E+06
426 0.2205E-13 0.3966E+06
427 0.2086E-13 0.3752E+06
428 0.1791E-13 0.3221E+06
429 0.1656E-13 0.2977E+06
430 0.4232E-13 0.7610E+06
431 0.4080E-13 0.7337E+06
432 0.9808E-14 0.1764E+06
433 0.3807E-13 0.6847E+06
434 0.2360E-13 0.4244E+06
435 0.2224E-13 0.4000E+06
436 0.1742E-13 0.3132E+06
437 0.2295E-13 0.4128E+06
438 0.1296E-13 0.2330E+06
439 0.7086E-14 0.1274E+06
440 0.7592E-14 0.1365E+06
441 0.7435E-14 0.1337E+06
442 0.6985E-14 0.1256E+06
443 0.7260E-14 0.1305E+06
444 0.1187E-13 0.2135E+06
445 0.8442E-14 0.1518E+06
446 0.2708E-13 0.4870E+06
447 0.1177E-13 0.2117E+06
448 0.2136E-13 0.3841E+06
449 0.2319E-13 0.4170E+06
450 0.8212E-13 0.1477E+07
451 0.1505E-12 0.2707E+07
452 0.2014E-12 0.3622E+07
453 0.1495E-12 0.2688E+07
454 0.1474E-12 0.2650E+07
455 0.1337E-12 0.2404E+07
456 0.1294E-12 0.2326E+07
457 0.1085E-12 0.1950E+07
458 0.1074E-12 0.1931E+07
459 0.1221E-12 0.2195E+07
460 0.1112E-12 0.2000E+07
461 0.8983E-13 0.1615E+07
462 0.8721E-13 0.1568E+07
463 0.6224E-13 0.1119E+07
464 0.4142E-13 0.7449E+06
465 0.4045E-13 0.7274E+06
466 0.4175E-13 0.7508E+06
467 0.5651E-13 0.1016E+07
468 0.3660E-13 0.6581E+06
469 0.4779E-13 0.8594E+06
470 0.4324E-13 0.7776E+06
471 0.2351E-13 0.4228E+06
472 0.1607E-13 0.2889E+06
473 0.1467E-13 0.2638E+06
474 0.1291E-13 0.2322E+06
475 0.1160E-13 0.2087E+06
476 0.1183E-13 0.2128E+06
477 0.6318E-12 0.1136E+08
478 0.5782E-12 0.1040E+08
479 0.6298E-12 0.1133E+08
480 0.4826E-12 0.8678E+07
481 0.4576E-12 0.8229E+07
482 0.4412E-12 0.7935E+07
483 0.4393E-12 0.7900E+07
484 0.3897E-12 0.7008E+07
485 0.1771E-12 0.3185E+07
486 0.1062E-11 0.1910E+08
487 0.2690E-12 0.4838E+07
488 0.1820E-11 0.3274E+08
489 0.5314E-11 0.9556E+08
490 0.1052E-11 0.1892E+08
491 0.3430E-10 0.6168E+09
492 0.2005E-10 0.3605E+09
493 0.8450E-09 0.1519E+11
494 0.9580E-10 0.1723E+10
495 0.8443E-10 0.1518E+10
496 0.8473E-10 0.1524E+10
497 0.8854E-10 0.1592E+10
498 0.1023E-09 0.1840E+10
499 0.9680E-10 0.1741E+10
500 0.1031E-09 0.1854E+10
501 0.1023E-09 0.1839E+10
502 0.1198E-09 0.2154E+10
503 0.1166E-09 0.2097E+10
504 0.1153E-09 0.2074E+10
505 0.1117E-09 0.2009E+10
506 0.1082E-09 0.1946E+10
507 0.1000E-09 0.1799E+10
508 0.9808E-10 0.1764E+10
509 0.1034E-09 0.1859E+10
510 0.8452E-10 0.1520E+10
511 0.8790E-10 0.1581E+10
512 0.8766E-10 0.1576E+10
513 0.1114E-09 0.2004E+10
514 0.1082E-09 0.1945E+10
515 0.1032E-09 0.1856E+10
516 0.1039E-09 0.1868E+10
517 0.1040E-09 0.1870E+10
518 0.9771E-10 0.1757E+10
519 0.9896E-10 0.1780E+10
520 0.9652E-10 0.1736E+10
521 0.1006E-09 0.1810E+10
522 0.8259E-09 0.1485E+11
523 0.1069E-09 0.1922E+10
524 0.1574E-08 0.2831E+11
525 0.8474E-08 0.1524E+12
526 0.4373E-08 0.7864E+11
527 0.9836E-08 0.1769E+12
528 0.7772E-08 0.1398E+12
529 0.9603E-08 0.1727E+12
530 0.9394E-08 0.1689E+12
531 0.7786E-08 0.1400E+12
532 0.8859E-08 0.1593E+12
533 0.8539E-08 0.1536E+12
534 0.8559E-08 0.1539E+12
535 0.1024E-07 0.1842E+12
536 0.9762E-08 0.1756E+12
537 0.1253E-07 0.2254E+12
538 0.1900E-07 0.3418E+12
539 0.1455E-07 0.2616E+12
540 0.1127E-07 0.2026E+12
541 0.9876E-08 0.1776E+12
542 0.7413E-08 0.1333E+12
543 0.2753E-07 0.4951E+12
544 0.1421E-06 0.2556E+13
545 0.2048E-06 0.3683E+13
546 0.6678E-06 0.1201E+14
547 0.4116E-06 0.7401E+13
548 0.3658E-06 0.6579E+13
549 0.3308E-06 0.5949E+13
550 0.2646E-06 0.4758E+13
551 0.3465E-06 0.6231E+13
552 0.7739E-06 0.1392E+14
553 0.6279E-06 0.1129E+14
554 0.4422E-06 0.7953E+13
555 0.4377E-06 0.7872E+13
556 0.3969E-06 0.7138E+13
557 0.4358E-06 0.7837E+13
558 0.3883E-06 0.6983E+13
559 0.4282E-06 0.7701E+13
560 0.3450E-06 0.6204E+13
561 0.1441E-06 0.2592E+13
562 0.4867E-05 0.8753E+14
563 0.6552E-05 0.1178E+15
564 0.4674E-05 0.8405E+14
565 0.4739E-05 0.8521E+14
566 0.2265E-03 0.4073E+16
567 0.2192E-04 0.3942E+15
568 0.7906E-05 0.1422E+15
569 0.6352E-05 0.1142E+15
570 0.8829E-05 0.1588E+15
571 0.3461E-04 0.6224E+15
572 0.2160E-03 0.3885E+16
573 0.2746E-03 0.4937E+16
574 0.2159E-03 0.3882E+16
575 0.4249E-03 0.7641E+16
576 0.6762E-02 0.1216E+18
577 0.6263E-02 0.1126E+18
578 0.1897E-01 0.3411E+18
579 0.2712E-01 0.4877E+18
580 0.1473E-01 0.2649E+18
581 0.2903E-02 0.5220E+17
582 0.6895E-01 0.1240E+19
583 0.9198E-02 0.1654E+18
584 0.8402E-01 0.1511E+19
585 0.5654E-01 0.1017E+19
586 0.2743E-01 0.4933E+18
587 0.3707E-01 0.6667E+18
588 0.1361E-01 0.2448E+18
589 0.2204E+00 0.3964E+19
590 0.1542E+00 0.2772E+19
591 0.2616E+00 0.4705E+19
592 0.2160E+00 0.3884E+19
593 0.1613E+00 0.2901E+19
594 0.1439E+00 0.2589E+19
595 0.1391E+00 0.2502E+19
596 0.2518E+00 0.4528E+19
597 0.1758E+00 0.3161E+19
598 0.1886E+00 0.3392E+19
599 0.1896E+00 0.3409E+19
600 0.1808E+00 0.3252E+19
601 0.1626E+00 0.2925E+19
602 0.1587E+00 0.2854E+19
603 0.1517E+00 0.2728E+19
604 0.4189E+00 0.7534E+19
605 0.1374E+00 0.2471E+19
606 0.2357E+00 0.4238E+19
607 0.5520E-01 0.9927E+18
608 0.2557E+00 0.4599E+19
609 0.2936E+00 0.5280E+19
610 0.1846E+00 0.3320E+19
611 0.1893E+00 0.3404E+19
612 0.2024E+00 0.3639E+19
613 0.7428E-01 0.1336E+19
614 0.9510E+00 0.1710E+20
615 0.5656E+00 0.1017E+20
616 0.1814E+01 0.3262E+20
617 0.1427E+02 0.2566E+21
617 0.1427E+02 0.2566E+21
Best norm better than final one: 80 1.0177769197390787 2.5656546662490971E+020 617
ERROR:: IterSolve: Numerical Error: System diverged over maximum tolerance.
STOP 1
Best Regards,
Mandeep