GPE |ψ|² and v(ψ) reconstruct the Unruh acoustic metric; matches Painlevé–Gullstrand PG within declared tolerances.

Close the gap between 'visual analogy' and an auditable numerical pipeline: from the same GPE background $\psi$ as sim_03, rebuild the Unruh hydrodynamic acoustic metric and compare — component-wise — to the Painlevé–Gullstrand coefficients and radial null slopes.

On the $y=0$ slice, using $|\psi|^2$ for $\rho$ and the analytic PG radial velocity $v_r$, the Unruh form $ds^2 = (\rho/c_s)[-(c_s^2-v^2)dt^2 - 2\mathbf{v}\!\cdot d\mathbf{x}\,dt + d\mathbf{x}^2]$ matches the metric reconstructed from Madelung velocities from $\psi$ to sub-$10^{-2}$ absolute tolerance in $g_{tt}$ and sub-$4\cdot 10^{-3}$ in $g_{tx}$ outside the core and away from the ergosurface band; outgoing radial null speeds agree with sim_03 to small RMS relative error; the $|v|/c_s=1$ locus matches $r_h$.

sim_03 demonstrates trapping and PG geometry; sim_58 shows the **same** wavefunction defines a Lorentzian $g_{\mu\nu}^{\mathrm{eff}}$ in the standard analogue-gravity normalization, strengthening the 'gravity = acoustic metric' claim for reviewers who want a line-by-line numerical check.