def METHOD_NAME(hexa_x, hexa_y, hexa_z, xsi, eta, zet, mesh):
    dxdxsi, dxdeta, dxdzet, dydxsi, dydeta, dydzet, dzdxsi, dzdeta, dzdzet = dxdxsi3D_lin(hexa_x, hexa_y, hexa_z, xsi, eta, zet, mesh)
    jac = dxdxsi * (dydeta*dzdzet - dzdeta*dydzet)\
        - dxdeta * (dydxsi*dzdzet - dzdxsi*dydzet)\
        + dxdzet * (dydxsi*dzdeta - dzdxsi*dydeta)
    return jac