The Hodge star operator is customarily denoted by * and maps p-forms on an n-dimensional manifold to (n - p)-forms.  To define the Hodge star operator, let omega_1, omega_2, ..., omega_n be an orthonormal basis of 1-forms with respect to the metric g.  Pick integers i_1, i_2, ..., i_r with 1 <= i_1 < i_2 < ... <i_r <= n and let 1 <= j_1 < j_2 < ... < j_s <= n be the complementary set of integers (r = n - s).  Let

The program HodgeStar(g, omega) returns the Hodge star of the form omega with respect to the metric g.

The program HodgeStar(g, omega) calculates *(omega) as follows. [i] Convert omega to skew-symmetric, covariant rank r tensor T. [ii] Fully contract T with the contravariant permutation symbol and multiply by the weight 1 scalar density defined by the metric g.  The result is a skew-symmetric contravariant rank n-r tensor S. [iii] Lower all the indices of S with the metric g and convert the resulting skew-symmetric, contravariant rank n-r tensor to a differential form.

This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form HodgeStar(...) only after executing the command with(DifferentialGeometry) and with(Tensor) in that order.  It can always be used in the long form DifferentialGeometry:-Tensor:-HodgeStar.