case.sif

! Vectorized version of StatElecSolve where computation of capacitance
! matrix is done using the old keywords and constraint modes library
! functionality.
!
! This one sets the modes using Body Force section instead of
! BC section. Results should be exactly the same. 
!
! P.R. 30.5.2024

Check Keywords Warn

Header
  Mesh DB "." "plates"
End

Simulation
  Max Output Level = 7
  Coordinate System = Cartesian
  Simulation Type = Steady State
  Steady State Max Iterations = 1

!  Mesh Levels = 2
  Serendipity P Elements = False
End

Constants
!  Permittivity Of Vacuum = 8.8542e-12
!  This one is changed to one to allow easier verification
  Permittivity Of Vacuum = 1.0
End

Body 1
  Equation = 1
  Material = 1
  Body Force = 1
End
Body 2
  Body Force = 2
End
Body 3
  Body Force = 3
End
Body 4
  Body Force = 4
End
Body 5
  Body Force = 5
End
Body 6
  Body Force = 6
End

Body Force 1
  Calc This = Logical True
End 

Body Force 2
  Capacitance Body = Integer 1
End 

Body Force 3
  Capacitance Body = Integer 2
End 
Body Force 4
  Capacitance Body = Integer 3
End 
Body Force 5
  Capacitance Body = Integer 4
End 
Body Force 6
  Capacitance Body = Integer 5
End 

Equation 1
  Active Solvers(2) = 1 2
  Calculate Electric Energy = True
End

Solver 1
  Equation = Stat Elec Solver
  Variable = Potential
  Procedure = "StatElecSolveVec" "StatElecSolver"

  Linear System Solver = iterative
  Linear System Iterative Method = BiCGStabl
  Linear System Max Iterations = 1000
  Linear System Convergence Tolerance = 1.0e-08
  Linear System Preconditioning = ILU1
  Linear System Residual Output = 10

  Nonlinear System Max Iterations = 1
  Nonlinear System Consistent Norm = True

  Calculate Capacitance Matrix = True
  Element = p:2
End

Solver 2
  Exec Solver = never
  Equation = "result output"
  Procedure = "ResultOutputSolve" "ResultOutputSolver"
  Output File Name = case
  Vtu Format = Logical True

! This enforces the constraints to be written to individual files 
  Constraint Modes Analysis = Logical True
End

Solver 3
  Equation = SaveScalars
  Procedure = "SaveData" "SaveScalars"
  Filename = f.dat
  Show Norm = True
  Show Norm Index = 1

  Variable 1 = Potential
  Coefficient 1 = "Relative Permittivity"
  Operator 1 = "diffusive energy"
  Mask Name 1 = "calc this"

  Parallel Reduce = True
End


Material 1
  Relative Permittivity = 1
End

Boundary Condition 1
  Target Boundaries = 1
  Potential = 0.0
End

Solver 1 :: Reference Norm = 1.37719121E-01
Solver 3 :: Reference Norm = 5.08932089E+00