Header
    CHECK KEYWORDS Warn
    Mesh DB "elmer_mesh" "."
    Include Path ""
    Results Directory "."
End

Simulation
    Use Mesh Names = Logical True
    Max Output Level = 5
    Coordinate System = Cartesian
    Coordinate Mapping(3) = 1 2 3
    Simulation Type = Steady state
    Steady State Max Iterations = 1
    Output Intervals = 1
    Coordinate Scaling = 1e-3
End

Constants
    Stefan Boltzmann = 5.670374419e-08
    Permittivity of Vacuum = 8.85418781e-12
    Permeability of Vacuum = 1.25663706e-6
    Boltzmann Constant = 1.380649e-23
End

Body 1
    Name = air
    Equation = 1
    Material = 1
End

Body 2
    Name = magnet-N
    Equation = 1
    Material = 2
End

Body 3
    Name = magnet-S
    Equation = 1
    Material = 3
End

Body 4
    Name = iron
    Equation = 1
    Material = 4
End

Solver 1
    Equation = MgDyn
    Procedure = "MagnetoDynamics" "WhitneyAVSolver"
    Optimize Bandwidth = True

    Nonlinear System Convergence Tolerance = 1.0e-10
    Nonlinear System Max Iterations = 50
    Nonlinear System Newton After Iterations = 3
    Nonlinear System Newton After Tolerance = 1.0e-3
    Nonlinear System Relaxation Factor = 1

    Linear System Solver = Iterative
    ! Bicgstabl is usually performing well for AV solver
    Linear System Iterative Method = BiCGStabl
    Linear System Max Iterations = 500
    Linear System Convergence Tolerance = 1.0e-8
    BiCGstabl polynomial degree = 4
    ! We cannot use ILU since vector potential is not uniquely defined
    Linear System Preconditioning = none
    Linear System ILUT Tolerance = 1.0e-3
    Linear System Abort Not Converged = False
    Linear System Residual Output = 10
    Linear System Precondition Recompute = 1
End

Solver 2
    Equation = MgDynPost
    Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields"
    Optimize Bandwidth = True

    Calculate Nodal Fields = False
    Calculate Elemental Fields = True

    Calculate Magnetic Field Strength = True

    ! Perform the averaging within this solver
    Average Within Materials = Logical True
End

Solver 3
  Exec Solver = After Timestep
  Equation = "VtuOutput"
  Procedure = "ResultOutputSolve" "ResultOutputSolver"
  Output File Name = static_mover
  Vtu Format = True
  Binary Output = True
  Single Precision = True
  Save Geometry Ids = True
  Show Variables = True
! This looks much better!
  Discontinuous Bodies = True
End

Equation 1
    Name = "Equation 1"
    Active Solvers(3) = 1 2 3
End

Material 1
    Name = "Air"
    Relative Permeability = Real 1
    Relative Permittivity = Real 1

End

Material 2
    Name = "Neodymium"
    Relative Permeability = Real 1.05
    Relative Permittivity = Real 1
    Magnetization 3 = Real 1178e3
End

Material 3
    Name = "Neodymium"
    Relative Permeability = Real 1.05
    Relative Permittivity = Real 1
    Magnetization 3 = Real -1178e3
End

Material 4
    Name = "Iron"
    Relative Permeability = Real 4000
    Relative Permittivity = Real 1
End


Boundary Condition 1
    Name = Farfield
    AV {e} = 0.0
End