refactor(engine): update workspace to use earth property matrices

Author: amaurigmartins

src/engine/workspace.jl

@@ -43,10 +43,12 @@ $(TYPEDFIELDS)
 	phase_map::Vector{Int}
 	"Vector of cable mapping indices."
 	cable_map::Vector{Int}
-	"Effective earth parameters as a vector of NamedTuples."
-	earth::Vector{
-		NamedTuple{(:rho_g, :eps_g, :mu_g), Tuple{Vector{T}, Vector{T}, Vector{T}}},
-	}
+	"Effective earth resistivity (layers × freq)."
+	rho_g::Matrix{T}
+	"Effective earth permittivity (layers × freq)."
+	eps_g::Matrix{T}
+	"Effective earth permeability (layers × freq)."
+	mu_g::Matrix{T}
 	"Operating temperature [°C]."
 	temp::T
 	"Number of frequency samples."
@@ -135,10 +137,10 @@ function init_workspace(
 		end
 	end
 
-	earth = _get_earth_data(
+	(rho_g, eps_g, mu_g) = _get_earth_data(
 		formulation.equivalent_earth,
 		problem.earth_props,
-		problem.frequencies,
+		freq,
 		T,
 	)
 
@@ -152,7 +154,8 @@ function init_workspace(
 		r_ins_in = r_ins_in, r_ins_ext = r_ins_ext,
 		rho_cond = rho_cond, alpha_cond = alpha_cond, mu_cond = mu_cond,
 		eps_cond = eps_cond, rho_ins = rho_ins, mu_ins = mu_ins, eps_ins = eps_ins,
-		tan_ins = tan_ins, phase_map = phase_map, cable_map = cable_map, earth = earth,
+		tan_ins = tan_ins, phase_map = phase_map, cable_map = cable_map, rho_g = rho_g,
+		eps_g = eps_g, mu_g = mu_g,
 		temp = temp, n_frequencies = n_frequencies, n_phases = n_phases,
 		n_cables = system.num_cables, Zprim = Zprim, Yprim = Yprim,
 	)