I want to convert SI to Gaussian units given vectors:

Subscript[Ε, si] = {Ex, Ey, Ez}; Subscript[B, si] = {Bx, By, Bz}; Subscript[S, si] = 1/Subscript[μ, 0] (Subscript[Ε, si]/[Cross]Subscript[B, si]) Subscript[S, si] /. {Subscript[B, si] -> Subscript[B, si]/c} why doesn't the change apply?
Also, I have this horrible vector I'm trying to simplify in a similar manner (and apply certain conditions to it such as $B_z E_x = B_x E_z$):

Subscript[Ey, g] = γ (Ey + v/c Bz); Subscript[Ez, g] = γ (Ez - v/c By); Subscript[By, g] = γ (By - v/c Ez); Subscript[Bz, g] = γ (Bz + v/c Ey); Subscript[S, p] = c/(4 π) (Subscript[Ε, si]\[Cross]Subscript[B, si]) /. {Ey -> Subscript[Ey, g], Ez -> Subscript[Ez, g], By -> Subscript[By, g], Bz -> Subscript[Bz, g]}; Subscript[S, p] /. {v -> c Tanh[ζ], Bz Ex -> Bx Ez, Tanh[ζ]/Sqrt[Sech[ζ]^2] -> Sinh[ζ]} // FullSimplify I see that Rule doesn't work on combinations of variables, at least not in this expression for some reason. I tried inputting Bz Ex - Bx Ez /.Bz Ex -> Bx Ez on a separate line and it returned 0, like it should. Is there some trick to it? Why is it being finicky? As far as trig functions, I tried TrigReduce, but it didn't do anything. The square root of a square thing always bugs me. I also tried using assumptions with $\zeta>0$ and $\text{sech}(\zeta)>0$, but again to no avail. Is there a way to reduce that as well? And is there a general way to change certain chunks of expressions like that?
ReplaceAllis that the expression is evaluated before replacing, so after evaluation MMA can't find any instances of what you're trying to replace $\endgroup$