Subgridding to Solving Magnetostatics within Discrete Geometric Approach

We propose a recipe to construct a symmetric positive definite and consistent reluctance constitutive matrix to be used within discrete geometric approaches when the primal grid is generated by an enhanced subgridding of a generic hexahedral grid. We focus on a magnetostatic problem as working example.