def METHOD_NAME(Grid):
    grid = Grid((4, 5))
    value_at_node = np.full(grid.number_of_nodes, 1.0)
    coef_at_link = np.full(grid.number_of_links, 1.0)
    mat, rhs = get_core_node_matrix(grid, value_at_node, coef_at_link=coef_at_link)
    actual, actual_rhs = get_core_node_matrix(grid, value_at_node * 2, coef_at_link)
    assert_array_equal(mat.toarray(), actual.toarray())
    assert_array_equal(rhs, actual_rhs / 2.0)