The definition of orientation on a vector space I'm using is as an equivalence class of ordered bases, where two bases are related if and only if the determinant of the change of basis is positive. Now, one has the following theorem
Theorem.
Let $(V, g, \mathcal{Or})$ be an oriented $n$-dimensional inner product space over $\Bbb{R}$. Then, there is a unique volume form $\omega$ on $V$, such that for every positively-oriented, ordered, orthonormal basis $\{e_1, \dots, e_n\}$ of $V$, we have \begin{align} \omega(e_1, \dots, e_n) &= 1. \end{align} In fact for any positively-oriented basis $\{f_1, \dots, f_n\}$, if we let $\{\varphi^1, \dots, \varphi^n\}$ be the dual basis for $V^*$, then \begin{align} \omega &= \sqrt{\det \left[ g(f_i, f_j)\right]}\,\,\varphi^1 \wedge \dots \wedge \varphi^n \end{align}
The sketch of the proof is: pick some positively-oriented, ordered, orthonormal basis $\{e_1, \dots, e_n\}$, and denote $\{\epsilon^1, \dots, \epsilon^n\}$ to be the dual basis. Next, let $\{f_1, \dots, f_n\}$ be any positvely-oriented basis, with $\{\varphi^1, \dots, \varphi^n\}$ its dual basis, and let $A$ be the change of basis matrix $(f_j = \sum_{i}A_{ij}e_i)$. Then, \begin{align} \epsilon^1 \wedge \dots \wedge \epsilon^n &= \det(A) \cdot \varphi^1 \wedge \dots \wedge \varphi^n. \end{align} In this case, it's relatively straight-forward to show that $A^tA = [g(f_i,f_j)]$. From this, it is easily deduced that we can define $\omega = \epsilon^1 \wedge \dots \wedge \epsilon^n$, this definition is basis-independent, and that $\omega$ is unique.
My question is whether there is an analogous statement in the pseudo-Riemannian case; i.e when we require $g:V \times V \to \Bbb{R}$ to only be bilinear, symmetric, and non-degenerate (as opposed to positive-definite). I tried to adapt the same proof to this situation. However, here's the trouble I faced: let $\{e_1, \dots, e_n\}$ be positively-oriented and orthonormal with respect to $g$, so that $g(e_i,e_j) = \pm \delta_{ij}$. Next, let $\{f_1, \dots, f_n\}$ be any positively-oriented basis. After some calculations, I found that \begin{align} (A^tA)_{ij} &= \sum_{k, \beta, \alpha}g(e_k, e_{\alpha}) g(e_{\alpha}, f_i) g(e_k, e_{\beta}) g(e_{\beta}, f_j). \end{align} In the Riemannian case, the $g(e_{(\cdot)}, e_{(\cdot)})$ would just become a Kronecker delta with a $+$ sign, and then (making use of the formula for expressing vectors in terms of an orthonormal basis) it immediately reduces to $g(f_i,f_j)$. However, in general, the minus signs cause trouble, and I'm not sure how to further simplify this.
So, in the general case, is there even a unique volume element determined only from $g$ and the orientation? If so, how does one go about proving it? What other modifications (if any) to the theorem are necessary?
Having had some more time to think about this, I realized that I've unnecessarily confused myself, and that I'd been trying to prove $A^tA = [g(f_i,f_j)]$ in the pseudo-Riemannian case, which is actually incorrect. The fix however, is very simple. For fun, I'll include the modified statement, along with two proofs.
The first proof makes use of the fact that a (pseudo) inner product, $g$, on $V$ induces for each $q \in \Bbb{N}$, a pseudo-inner-product $g_q$ on the subspace $\mathcal{A}^q(V)$ of alternating $q$-tensors over $V$.
The second proof establishes more in the process:
I found proof $1$ more conceptually clear: the pseudo-inner product $g$ induces one on every space of alternating tensors, and by using the fact $\dim \mathcal{A}^n(V) = 1$, it follows that there are exactly two volume forms which are "normalized", to $(-1)^{\#}$, and these volume forms only differ by a sign. Then, the orientation $\mathcal{Or}$ helps us to pick one of these two volume forms.
The second proof is more "down-to-earth", but necessarily more computationally tedious, because it doesn't require us to define the inner product on each space of alternating tensors (that by itself requires a little bit of work to show everything is well-defined etc). It only requires basic linear algebra, and knowing the relationship between the determinant and wedge products to prove that $\omega$ is well-defined and unique.
One advantage though is that for computational purposes, it gives us an explicit formula for the volume element, in terms of any positively-oriented basis. Such a formula could potentially be useful in certain concrete situations, for example, when performing integration on manifolds, and choosing a chart which is particularly suited to the specific problem at hand. In such a case, orthogonalizing, or normalizing the tangent vectors in the might be more tedious than simply calculating a determinant (for example, the tangent vectors in the spherical-coordinate chart are not normalized).