Normalized Solutions and Numerical Approximation for a Critical p-Laplacian Schrödinger–Bopp–Podolsky System with Variable Mass and Logarithmic Nonlinearity