Physicists sometimes get the lagrangian $$\mathcal{L}=-\frac{1}{2}\mbox dA \wedge \star \mbox dA - A \wedge \star J$$ define a functional given by $$S(A)=\int_{N_4} \mathcal{L}= \int_{N_4}-\frac{1}{2}\mbox dA \wedge \star \mbox dA - A \wedge \star J$$ and calculate variation of such functional $$\delta S(A)(\delta A)=\delta \int_{N_4} \mathcal L = \int_{N_4} \delta \mathcal{L}=\int_{N_4}-\frac{1}{2}\delta(\mbox dA\wedge \star \mbox dA)-\delta A\wedge \star J$$ etc. and they derive Maxwell's equation $d\star F=-\star J$. The functional is well defined - it's simply an integral of a 4-differential form over an oriented manifold (submanifold). But how is the variation defined? If $M$ is Riemannian manifold we've got Riemannian metric and induced inner product $(\cdot, \cdot)$ given by $$(\omega, \eta):=\int_M \left<\omega, \eta \right>\mathrm{dvol}$$ Hence we've got induced norm of a differential form and can define variation of a functional and of a differential form as its Frechet derivative. Alas, in physics, manifold is psudoriemannian manifold hence $$(\omega, \eta)=\int_M \left<\omega, \eta \right>\mathrm{dvol}$$ is no longer an inner product. We have no norm of a differential form, and without norm we cannot define variations as Frechet derivatives. Moreover! We cannot even define a local extremum of a functional.
Can you tell me how is defined variation of such functional, and a differential form? Does it even make sence what physicists are doing?