|
| enum class | MatrixFlag : index_t {
M = 0
, D = 1
, K = 2
, F = 3
,
dMdX_times_x = 4
, dMdX_times_h = 5
, dKdX_times_x = 6
, dFdX = 7
} |
| | Flags indicating which matrices/vectors have been set by the IWG. More...
|
| enum class | SideSetDofLinkMode {
FacetOnly = 0
, FacetAndMaster = 1
, FacetAndSlave = 2
, MasterAndSlave = 3
,
Cut = 4
, Inactive = 5
, UNDEFINED = 6
} |
| | how the dofs of a sideset element are linked to its facet, master and slave elements ( six modes plus UNDEFINED ) More...
|
| enum class | BoundaryConditionImposing {
Free = 0
, Dirichlet = 1
, Neumann = 2
, Alpha = 3
,
Lambda = 4
, Weak = 5
, UNDEFINED = 6
} |
| enum class | BoundaryConditionType {
Neumann
, Dirichlet
, Bearing
, Gauge
,
Current
, Voltage
, CircuitCurrent
, CircuitVoltage
,
Background
, BackgroundDirichlet
, UNDEFINED
} |
| enum class | GroupActivationMode { GeometryAndDofs = 0
, GeometryOnly = 1
, Inactive = 2
} |
| | Defines how a fem::Block or fem::SideSet is activated based on its DomainType. More...
|
| enum class | SolverAlgorithm { Direct = 0
, NewtonRaphson = 1
, Picard = 2
, UNDEFINED = 3
} |
| enum class | MaxwellPostprocessorType {
Air
, Ferro
, Conductor
, SuperConductor
,
ThinShellConductor
, ThinShellSuperConductor
, SideConnector
, UNDEFINED
} |
|
| void | initialize_integration_points (const ElementType &aElementType, Vector< real > &aWeights, Matrix< real > &aPoints, const uint aIntegrationOrder, const IntegrationScheme aIntegrationScheme) |
| void | initialize_integration_points (const GeometryType &aGeometryType, Vector< real > &aWeights, Matrix< real > &aPoints, const uint aIntegrationOrder, const IntegrationScheme aIntegrationScheme) |
| void | initialize_integration_points_on_facet (const ElementType aElementType, const uint aSideIndex, Vector< real > &aWeights, Matrix< real > &aPoints, const uint aIntegrationOrder, const IntegrationScheme aIntegrationScheme) |
| void | initialize_integration_points_on_facet (const ElementType aElementType, const uint aSideIndex, const uint aOrientation, Vector< real > &aWeights, Matrix< real > &aPoints, const uint aIntegrationOrder, const IntegrationScheme aIntegrationScheme) |
| InterpolationFunction * | initialize_shape_function (const ElementType &aElementType, const Matrix< real > &aXi, Cell< Matrix< real > > &aN, Cell< Matrix< real > > &adNdXi, Cell< Matrix< real > > &ad2NdXi2) |
| uint | num_nedelec_dofs (const ElementType aElementType) |
| void | compute_h_trace (Calculator *aNormalCalc, const bool aMaster, const bool aIsConductor, const Vector< real > &aPhi, Vector< real > &aScratch, Vector< real > &aH) |
| | Full field trace of one volume side of a thin-shell facet, for linear elements ( the k = 0 contract of compute_hn ):
|
| const Vector< real > & | compute_hn (Calculator *aCalc, const uint k) |
| | Computes the purely-normal magnetic field hn at the thin-shell facet: the average of the master and slave volume traces ( compute_h_trace, per-side dispatch phi-region / h-conductor ), projected onto the facet normal.
|
| bool | pin_dirichlet_dof (const DofManagerBase *aParent, Dof *aDof, const real aValue, const id_t aGroupID, const id_t aNodeID, index_t &aNumFirstFlips) |
| | shared Dirichlet pin for a single node dof, used by SideSet::impose_dirichlet and the Maxwell factory pre-fix loop so both sites apply ONE shape test.
|
| string | to_string (const BoundaryConditionType aBoundaryConditiontype) |
| BoundaryConditionType | boundary_condition_type (const string &aString) |
| bool | has_block_global (const BoundaryConditionType aBoundaryConditionType) |
| | true if a condition of this type publishes its imposed value as a mesh global named after its deck block.
|
| uint | anderson_mixing_step (const Vector< real > &aX, const Vector< real > &aR, const ShiftRegister< Vector< real > > &aXHistory, const ShiftRegister< Vector< real > > &aRHistory, const real aBeta, Matrix< real > &aDeltaR, Vector< real > &aRhs, Vector< real > &aWork, Vector< real > &aGamma, Vector< real > &aColNorm, Vector< real > &aXNew) |
| | one type-II Anderson mixing step (Walker & Ni 2011) for the fixed-point iteration x <- G( x ) with residual r = G( x ) - x :
|
| const input::Section * | winning_nonlinear_section (const input::Section *aSolver) |
| | The nonlinear section the penalty blocks are read from: the same alias rule Controller::set_params applies to every key in it – nonlinear magnetic wins over nonlinear, and a block under the losing section is silently ignored.
|
| real | read_ghost_eta (const input::Section *aSolver) |
| | The Nitsche ghost switch as the deck states it, read from the solver section:
|
| bool | ghost_facets_requested (const input::Section *aSolver) |
| | true if the deck asks for duplicate interface dofs and ghost facets
|
| bool | ghost_facets_requested (const InputFile &aInputFile) |
| bool | recovery_block_inverse (const Matrix< real > &aGram, const real aMaxCondition, Matrix< real > &aInverse, real &aCondition, Matrix< real > &aSchur, Vector< real > &aCoupling, Vector< real > &aSchurCoupling, Matrix< real > &aSvdU, Matrix< real > &aSvdVT, Vector< real > &aSvdS, Vector< real > &aSvdWork) |
| | Inverse of the Gram matrix of a patch-recovery fit whose basis is one in slot 0 at the node ( postprocessor_recovery_theory.md §3.2 ).
|
| void | mises_planestress (DofManagerBase *aField) |
| void | T_h_picard (Calculator *aCalc, TimestepMatrices *aMatrices) |
| | Collapsed thermal Picard kernel: all per-point material math is delegated to aCalc->maxwell() ( thermal-side instance; the Maxwell peer element is linked automatically by Calculator::link ), whose constructor-time dispatch covers the bulk/thin-shell, metal/alloy/ HTS, defect, and piecewise axes.
|
| void | T_h_newton (Calculator *aCalc, TimestepMatrices *aMatrices) |
| void | T_phi (Calculator *aCalc, TimestepMatrices *aMatrices) |
| uint belfem::fem::anderson_mixing_step |
( |
const Vector< real > & | aX, |
|
|
const Vector< real > & | aR, |
|
|
const ShiftRegister< Vector< real > > & | aXHistory, |
|
|
const ShiftRegister< Vector< real > > & | aRHistory, |
|
|
const real | aBeta, |
|
|
Matrix< real > & | aDeltaR, |
|
|
Vector< real > & | aRhs, |
|
|
Vector< real > & | aWork, |
|
|
Vector< real > & | aGamma, |
|
|
Vector< real > & | aColNorm, |
|
|
Vector< real > & | aXNew ) |
|
inline |
one type-II Anderson mixing step (Walker & Ni 2011) for the fixed-point iteration x <- G( x ) with residual r = G( x ) - x :
gamma = argmin || r_k - DR * gamma ||_2 x_{k+1} = x_k + beta r_k - ( DX + beta DR ) gamma
where the difference columns DR / DX are formed from the committed history (newest aXHistory.size() snapshots) chained with the current pair ( aX, aR ). The least-squares solve runs on column-normalized DR via QR ( gels ); gamma is un-scaled afterwards.
Failure policy: a rank-deficient solve ( info > 0 ), a non-finite or oversized gamma ( gAndersonGammaMax ), or a vanishing difference column drops the OLDEST column and retries once; if the retry fails too, the plain relaxed Picard step x + beta r is written instead. An illegal-argument error ( info < 0 ) is a programming error and always aborts.
With an empty history ( or on fallback ) the result is exactly the relaxed Picard step: x_{k+1} = ( 1 - beta ) x + beta G( x ).
All scratch is caller-owned and must be presized: aDeltaR to n x depth, aRhs to n, aGamma and aColNorm to depth. aWork is grown to the optimal gels size on first use and reused thereafter.
- Returns
- number of history columns used ( 0 = plain step )
shared Dirichlet pin for a single node dof, used by SideSet::impose_dirichlet and the Maxwell factory pre-fix loop so both sites apply ONE shape test.
If the dof is condensed onto a single source with unit weight, the SOURCE is fixed to aValue / weight and the function returns true ( incrementing aNumFirstFlips when the source was still free ); a dof that is not hanging returns false and the CALLER fixes it directly. BELFEM_ERRORs on a hanging dof with more than one source ( the reconstruction would overwrite a fix on the alias after every solve ), a non-unit single-source weight, a broken source, or a first free-to-fixed flip after aParent->is_initialized(). Call only for non-duplicate nodes: a cut duplicate hangs on its original plus the cut currents at unit weight, and its caller fixes the alias deliberately, where the fix is inert: enforcing the value would constrain the sum of those currents, to zero when the original is held at the same value.
The Nitsche ghost switch as the deck states it, read from the solver section:
-1 no nitsche ghost penalty block -> ghost OFF ( opt-in ) 0 eta : 0 -> ghost OFF >0 eta : <value> -> ghost ON with that eta
A present block without eta is a setup error: k_reg alone cannot say whether the ghost is wanted, and an empty block that silently meant "off" would be indistinguishable from a typo.
ONE reader for three consumers – the thin-shell factory ( which decides whether duplicate interface dofs and ghost facets are created at all ), the controller ( which writes eta into the IWG ) and the mesh cache tag ( which must miss the cache when the switch flips ) – so the three cannot disagree about what "off" means. Setup-path code, called a handful of times per run
Collapsed thermal Picard kernel: all per-point material math is delegated to aCalc->maxwell() ( thermal-side instance; the Maxwell peer element is linked automatically by Calculator::link ), whose constructor-time dispatch covers the bulk/thin-shell, metal/alloy/ HTS, defect, and piecewise axes.
Dispatch in IWG_MaxwellThermal::link_to_group() is by domain type only.
The Newton companion T_h_newton adds the consistent tangent blocks ( dcp/dT with qhist contraction, the dlambda/dT mixed-operator outer product, and the drho/dT quench feedback — B and beta are frozen per thermal solve, not thermal dofs ); derivation in thermal_matrices_cleanup_and_newton_plan.md §3.
1.18.0