source: trunk/fsource/pyspg.for @ 3396

Last change on this file since 3396 was 3396, checked in by vondreele, 4 years ago

set refList=False for a couple of calls to GetHistogramPhaseData? where reflection list shouldn't be rebuilt
fix OnDrawDist? to report Sym Ops correctly
Move atom labels over 4 spaces to be clear of atom ball

File size: 1.9 KB
Line 
1C Space group access routines for python
2
3      SUBROUTINE SGFORPY(SPG,LAUE,SGINV,SGLATT,SGUNIQ,SGPOL,
4     1  SGNOPS,SGMTRX,SGTRNS,SGGEN,IERR)
5Cf2py intent(in)  SPG
6Cf2py intent(out) LAUE
7Cf2py intent(out) SGINV
8Cf2py intent(out) SGLATT
9Cf2py intent(out) SGUNIQ
10Cf2py intent(out) SGPOL
11Cf2py intent(out) SGNOPS
12Cf2py intent(out) SGMTRX
13Cf2py intent(out) SGTRNS
14Cf2py intent(out) SGGEN
15Cf2py intent(out) IERR
16
17      CHARACTER*(20) SPG
18      INTEGER*4     LAUE,SGINV,SGLATT,SGUNIQ,SGNOPS,IERR,SGNCEN
19      REAL*4        SGMTRX(24,3,3),SGTRNS(24,3),SGGEN(24)
20      REAL*4        RT(5,4,25),CEN(3,4)
21      INTEGER*4     JRT(3,5,24)
22
23
24      CALL SGROUPNP(SPG,LAUE,SGUNIQ,SGINV,SGLATT,SGNOPS,SGPOL,JRT,
25     1  CEN,SGNCEN,RT,IERR)
26      DO K=1,SGNOPS
27        DO I=1,3
28          DO J=1,3
29            SGMTRX(K,I,J) = JRT(I,J,K)
30            SGTRNS(K,I) = JRT(I,4,K)/12.
31          END DO
32        END DO
33        SGGEN(K) = JRT(3,5,K)
34      END DO
35      RETURN
36      END
37
38      SUBROUTINE GENHKLPY(XH,NSYM,SGMTRX,SGTRNS,ICEN,NCV,SGCEN,JHK,
39     1  HKL,IABSNT,MULP)
40Cf2py intent(in)  XH
41Cf2py intent(in)  NSYM
42Cf2py intent(in)  SGMTRX
43Cf2py intent(in)  SGTRNS
44Cf2py depend(NSYM) SGMTRX,SGTRNS
45Cf2py intent(in)  ICEN
46Cf2py intent(in)  NCV
47Cf2py intent(in)  SGCEN
48Cf2py depend(NCV) SGCEN
49Cf2py intent(out) JHK
50Cf2py intent(out) HKL
51Cf2py intent(out) IABSNT
52Cf2py intent(out) MULP
53
54      INTEGER*4     ICEN,NSYM
55      REAL*4        SGMTRX(NSYM,3,3),SGTRNS(NSYM,3),SGCEN(NCV,3)
56      REAL*4        CEN(3,NCV),HKL(4,24),XH(4)
57      INTEGER*4     JRT(3,5,24),JHK,NCV
58
59      DO J=1,NCV
60        DO I=1,3
61          CEN(I,J) = SGCEN(J,I)
62        END DO
63      END DO
64      DO K=1,NSYM
65        DO I=1,3
66          DO J=1,3
67            JRT(I,J,K) = SGMTRX(K,I,J)*1.
68            JRT(I,4,K) = SGTRNS(K,I)*12.
69          END DO
70        END DO
71      END DO
72      CALL GENHKL(XH,NSYM,JRT,ICEN,NCV,CEN,JHK,HKL,IABSNT,MULP)
73      RETURN
74      END
Note: See TracBrowser for help on using the repository browser.