“…In data preparation stage, the following are determined (a) geometry and material physical properties of the subregions, (b) a nodal distribution for discretizing the full domain, (c) the domain of influence for each node, (d) the number and size of the integration cells, and the number of integration points inside each integration cell, (e) support domain for each integration point, (f) weights, shape functions, and shape function derivatives, for each node in the entire region, (g) nodes and integration points along electrodes, and (h) computation of the second term on the right hand side of Eq. (22) and the terms on the right hand side of Eqs. (23) and (24) for each electrode.…”