def METHOD_NAME(self):
    # Here we will make up some j vectors that vary in space
    # j = [j_r, j_z] - to test face inner products
    fcts = FaceInnerProductFctsIsotropic()
    sig, jv = fcts.vectors(self.mesh)
    MfSig = self.mesh.get_face_inner_product(sig)
    numeric_ans = jv.T.dot(MfSig.dot(jv))
    ans = fcts.sol()
    print("------ Testing Face Inner Product-----------")
    print(
        " Analytic: {analytic}, Numeric: {numeric}, "
        "ratio (num/ana): {ratio}".format(
            analytic=ans, numeric=numeric_ans, ratio=float(numeric_ans) / ans
        )
    )
    assert np.abs(ans - numeric_ans) < TOL