def Hex8_shape_function_loc(xi_vec,node_local_coords): #Test function written
    """Compute the value of the shape function at xi_vec
    for a node with local coordinates node_local_coords"""
    xi1,xi2,xi3 = xi_vec
    nc1,nc2,nc3 = node_local_coords
    return 0.125*(1+xi1*nc1)*(1+xi2*nc2)*(1+xi3*nc3)