def METHOD_NAME(base_box, expected_matrix):
    assert np.allclose(base_box.to_matrix(), expected_matrix)
    base_box.xy *= 2
    assert isclose(base_box.to_matrix()[0, 1], 2 * expected_matrix[0, 1])
    base_box.yz *= 0.5
    assert isclose(base_box.to_matrix()[1, 2], 0.5 * expected_matrix[1, 2])
    base_box.Lx *= 3
    assert isclose(base_box.to_matrix()[0, 0], 3 * expected_matrix[0, 0])