def METHOD_NAME(permutation):
    n_atoms_max = 5
    system = get_simple_finite()
    desc = CoulombMatrix(n_atoms_max=n_atoms_max, permutation=permutation)
    cm = desc.create(system)
    cm_assumed = cm_python(system, n_atoms_max, permutation, False)
    assert np.allclose(cm, cm_assumed)