Révision 9dca524b

b/modules/screening.py
878 878 
    """
879 879 
    from copy import deepcopy
880 880 
    slab_ads_list = []
881 
    # Rotation over bond angle # TODO Check sampling
882 
    for alpha in np.arange(90, 180, 90 / min(num_pts, 2)):
881 
    # Rotation over bond angle
882 
    for alpha in np.arange(90, 180+1, 90 / max(1, num_pts-1))[:num_pts]:
883 883 
        # Rotation over surf-adsorbate dihedral angle
884 884 
        for beta in np.arange(0, max_hel, max_hel / num_pts):
885 885 
            # Rotation over adsorbate bond dihedral angle
886 
            for gamma in np.arange(0, 360, 360 / num_pts):
886 
            for gamma in np.arange(90, 270+1, 180/max(1, num_pts-1))[:num_pts]:
887 
                # Avoid duplicates as gamma rotation has no effect on plane
888 
                # angles.
889 
                if alpha == 180 and gamma > 90:
890 
                    continue
887 891 
                new_molec = deepcopy(orig_molec)
888 892 
                chemcat_rotate(new_molec, slab, ctr1_mol, ctr2_mol, ctr3_mol,
889 893 
                               ctr1_surf, ctr2_surf, norm_vect, alpha,

Formats disponibles : Unified diff