TABLE OF CONTENTS


ABINIT/ptg_Td [ Functions ]

[ Top ] [ Functions ]

NAME

 ptg_Td

FUNCTION

COPYRIGHT

 Copyright (C) 2010-2018 ABINIT group (MG)
 This file is distributed under the terms of the
 GNU General Public License, see ~abinit/COPYING
 or http://www.gnu.org/copyleft/gpl.txt .
 For the initials of contributors, see ~abinit/doc/developers/contributors.txt .

INPUTS

OUTPUT

PARENTS

      m_ptgroups

CHILDREN

SOURCE

 27 !********************************************************************************
 28 ! This include file has been automatically generated by the script ptg.py
 29 ! Do not edit! Change the script source instead.
 30 !********************************************************************************
 31 
 32 ! Point group name  Td (-43m)
 33 
 34 #if defined HAVE_CONFIG_H
 35 #include "config.h"
 36 #endif
 37 
 38 #include "abi_common.h"
 39 
 40 
 41 module m_ptg_Td
 42 contains
 43 !!**
 44 
 45 
 46 
 47  subroutine ptg_Td (nsym,nclass,sym,class_ids,class_names,Irr)
 48 
 49 !DEC$ NOOPTIMIZE
 50 
 51  use defs_basis
 52  use m_abicore
 53  use m_defs_ptgroups,  only : irrep_t
 54 
 55 !This section has been created automatically by the script Abilint (TD).
 56 !Do not modify the following lines by hand.
 57 #undef ABI_FUNC
 58 #define ABI_FUNC 'ptg_Td'
 59 !End of the abilint section
 60 
 61  implicit none
 62 !Arguments ------------------------------------
 63  integer,intent(out) :: nclass,nsym 
 64  !arrays
 65  integer,allocatable,intent(out) :: sym(:,:,:), class_ids(:,:)
 66  character(len=5),allocatable,intent(out) :: class_names(:)
 67  type(irrep_t),allocatable,intent(out) :: Irr(:)
 68  !Local variables-------------------------------
 69  complex(dpc) :: j=(0.0_dp,1.0_dp) 
 70  ! ********************************************************************************
 71 ! List of symmetries packed in classes
 72  nsym = 24
 73  ABI_MALLOC(sym, (3,3,nsym))
 74  sym(:,:,1) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
 75  sym(:,:,2) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
 76  sym(:,:,3) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
 77  sym(:,:,4) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
 78  sym(:,:,5) = RESHAPE( (/0, 1, 0, 0, 0, 1, 1, 0, 0/) ,(/3,3/) )
 79  sym(:,:,6) = RESHAPE( (/0, -1, 0, 0, 0, -1, 1, 0, 0/) ,(/3,3/) )
 80  sym(:,:,7) = RESHAPE( (/0, 1, 0, 0, 0, -1, -1, 0, 0/) ,(/3,3/) )
 81  sym(:,:,8) = RESHAPE( (/0, -1, 0, 0, 0, 1, -1, 0, 0/) ,(/3,3/) )
 82  sym(:,:,9) = RESHAPE( (/0, 0, 1, 1, 0, 0, 0, 1, 0/) ,(/3,3/) )
 83  sym(:,:,10) = RESHAPE( (/0, 0, -1, -1, 0, 0, 0, 1, 0/) ,(/3,3/) )
 84  sym(:,:,11) = RESHAPE( (/0, 0, 1, -1, 0, 0, 0, -1, 0/) ,(/3,3/) )
 85  sym(:,:,12) = RESHAPE( (/0, 0, -1, 1, 0, 0, 0, -1, 0/) ,(/3,3/) )
 86  sym(:,:,13) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 87  sym(:,:,14) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 88  sym(:,:,15) = RESHAPE( (/1, 0, 0, 0, 0, 1, 0, 1, 0/) ,(/3,3/) )
 89  sym(:,:,16) = RESHAPE( (/1, 0, 0, 0, 0, -1, 0, -1, 0/) ,(/3,3/) )
 90  sym(:,:,17) = RESHAPE( (/0, 0, 1, 0, 1, 0, 1, 0, 0/) ,(/3,3/) )
 91  sym(:,:,18) = RESHAPE( (/0, 0, -1, 0, 1, 0, -1, 0, 0/) ,(/3,3/) )
 92  sym(:,:,19) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 93  sym(:,:,20) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 94  sym(:,:,21) = RESHAPE( (/-1, 0, 0, 0, 0, -1, 0, 1, 0/) ,(/3,3/) )
 95  sym(:,:,22) = RESHAPE( (/-1, 0, 0, 0, 0, 1, 0, -1, 0/) ,(/3,3/) )
 96  sym(:,:,23) = RESHAPE( (/0, 0, -1, 0, -1, 0, 1, 0, 0/) ,(/3,3/) )
 97  sym(:,:,24) = RESHAPE( (/0, 0, 1, 0, -1, 0, -1, 0, 0/) ,(/3,3/) )
 98 
 99 ! Number of classes and corresponding indeces
100  nclass = 5
101  ABI_MALLOC(class_ids, (2,nclass))
102  class_ids(1,1) = 1
103  class_ids(2,1) = 1
104  class_ids(1,2) = 2
105  class_ids(2,2) = 4
106  class_ids(1,3) = 5
107  class_ids(2,3) = 12
108  class_ids(1,4) = 13
109  class_ids(2,4) = 18
110  class_ids(1,5) = 19
111  class_ids(2,5) = 24
112 
113 ABI_MALLOC(class_names,(5))
114  class_names(1) = "1+" 
115  class_names(2) = "2+" 
116  class_names(3) = "3+" 
117  class_names(4) = "-2+" 
118  class_names(5) = "-4+" 
119 
120 ! List of irreducible representations.
121  ABI_DT_MALLOC(Irr, (5))
122  Irr(1)%name = "A1"
123  Irr(1)%dim = 1
124  Irr(1)%nsym = 24
125  ABI_MALLOC(Irr(1)%mat, (1,1,24))
126  Irr(1)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
127  Irr(1)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
128  Irr(1)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
129  Irr(1)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
130  Irr(1)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
131  Irr(1)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
132  Irr(1)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
133  Irr(1)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
134  Irr(1)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
135  Irr(1)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
136  Irr(1)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
137  Irr(1)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
138  Irr(1)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
139  Irr(1)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
140  Irr(1)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
141  Irr(1)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
142  Irr(1)%mat(:,:,17) =  RESHAPE( (/1.0/), (/1, 1/) )
143  Irr(1)%mat(:,:,18) =  RESHAPE( (/1.0/), (/1, 1/) )
144  Irr(1)%mat(:,:,19) =  RESHAPE( (/1.0/), (/1, 1/) )
145  Irr(1)%mat(:,:,20) =  RESHAPE( (/1.0/), (/1, 1/) )
146  Irr(1)%mat(:,:,21) =  RESHAPE( (/1.0/), (/1, 1/) )
147  Irr(1)%mat(:,:,22) =  RESHAPE( (/1.0/), (/1, 1/) )
148  Irr(1)%mat(:,:,23) =  RESHAPE( (/1.0/), (/1, 1/) )
149  Irr(1)%mat(:,:,24) =  RESHAPE( (/1.0/), (/1, 1/) )
150 
151  Irr(2)%name = "A2"
152  Irr(2)%dim = 1
153  Irr(2)%nsym = 24
154  ABI_MALLOC(Irr(2)%mat, (1,1,24))
155  Irr(2)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
156  Irr(2)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
157  Irr(2)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
158  Irr(2)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
159  Irr(2)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
160  Irr(2)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
161  Irr(2)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
162  Irr(2)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
163  Irr(2)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
164  Irr(2)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
165  Irr(2)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
166  Irr(2)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
167  Irr(2)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
168  Irr(2)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
169  Irr(2)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
170  Irr(2)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
171  Irr(2)%mat(:,:,17) =  RESHAPE( (/-1.0/), (/1, 1/) )
172  Irr(2)%mat(:,:,18) =  RESHAPE( (/-1.0/), (/1, 1/) )
173  Irr(2)%mat(:,:,19) =  RESHAPE( (/-1.0/), (/1, 1/) )
174  Irr(2)%mat(:,:,20) =  RESHAPE( (/-1.0/), (/1, 1/) )
175  Irr(2)%mat(:,:,21) =  RESHAPE( (/-1.0/), (/1, 1/) )
176  Irr(2)%mat(:,:,22) =  RESHAPE( (/-1.0/), (/1, 1/) )
177  Irr(2)%mat(:,:,23) =  RESHAPE( (/-1.0/), (/1, 1/) )
178  Irr(2)%mat(:,:,24) =  RESHAPE( (/-1.0/), (/1, 1/) )
179 
180  Irr(3)%name = "E"
181  Irr(3)%dim = 2
182  Irr(3)%nsym = 24
183  ABI_MALLOC(Irr(3)%mat, (2,2,24))
184  Irr(3)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
185  Irr(3)%mat(:,:,2) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
186  Irr(3)%mat(:,:,3) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
187  Irr(3)%mat(:,:,4) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
188  Irr(3)%mat(:,:,5) =  RESHAPE( (/-0.5+0.86603*j, 0*j, 0*j, -0.5-0.86603*j/), (/2, 2/) )
189  Irr(3)%mat(:,:,6) =  RESHAPE( (/-0.5+0.86603*j, 0*j, 0*j, -0.5-0.86603*j/), (/2, 2/) )
190  Irr(3)%mat(:,:,7) =  RESHAPE( (/-0.5+0.86603*j, 0*j, 0*j, -0.5-0.86603*j/), (/2, 2/) )
191  Irr(3)%mat(:,:,8) =  RESHAPE( (/-0.5+0.86603*j, 0*j, 0*j, -0.5-0.86603*j/), (/2, 2/) )
192  Irr(3)%mat(:,:,9) =  RESHAPE( (/-0.5-0.86603*j, 0*j, 0*j, -0.5+0.86603*j/), (/2, 2/) )
193  Irr(3)%mat(:,:,10) =  RESHAPE( (/-0.5-0.86603*j, 0*j, 0*j, -0.5+0.86603*j/), (/2, 2/) )
194  Irr(3)%mat(:,:,11) =  RESHAPE( (/-0.5-0.86603*j, 0*j, 0*j, -0.5+0.86603*j/), (/2, 2/) )
195  Irr(3)%mat(:,:,12) =  RESHAPE( (/-0.5-0.86603*j, 0*j, 0*j, -0.5+0.86603*j/), (/2, 2/) )
196  Irr(3)%mat(:,:,13) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
197  Irr(3)%mat(:,:,14) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
198  Irr(3)%mat(:,:,15) =  RESHAPE( (/0*j, -0.5+0.86603*j, -0.5-0.86603*j, 0*j/), (/2, 2/) )
199  Irr(3)%mat(:,:,16) =  RESHAPE( (/0*j, -0.5+0.86603*j, -0.5-0.86603*j, 0*j/), (/2, 2/) )
200  Irr(3)%mat(:,:,17) =  RESHAPE( (/0*j, -0.5-0.86603*j, -0.5+0.86603*j, 0*j/), (/2, 2/) )
201  Irr(3)%mat(:,:,18) =  RESHAPE( (/0*j, -0.5-0.86603*j, -0.5+0.86603*j, 0*j/), (/2, 2/) )
202  Irr(3)%mat(:,:,19) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
203  Irr(3)%mat(:,:,20) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
204  Irr(3)%mat(:,:,21) =  RESHAPE( (/0*j, -0.5+0.86603*j, -0.5-0.86603*j, 0*j/), (/2, 2/) )
205  Irr(3)%mat(:,:,22) =  RESHAPE( (/0*j, -0.5+0.86603*j, -0.5-0.86603*j, 0*j/), (/2, 2/) )
206  Irr(3)%mat(:,:,23) =  RESHAPE( (/0*j, -0.5-0.86603*j, -0.5+0.86603*j, 0*j/), (/2, 2/) )
207  Irr(3)%mat(:,:,24) =  RESHAPE( (/0*j, -0.5-0.86603*j, -0.5+0.86603*j, 0*j/), (/2, 2/) )
208 
209  Irr(4)%name = "T1"
210  Irr(4)%dim = 3
211  Irr(4)%nsym = 24
212  ABI_MALLOC(Irr(4)%mat, (3,3,24))
213  Irr(4)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
214  Irr(4)%mat(:,:,2) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
215  Irr(4)%mat(:,:,3) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
216  Irr(4)%mat(:,:,4) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
217  Irr(4)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0/), (/3, 3/) )
218  Irr(4)%mat(:,:,6) =  RESHAPE( (/0.0, 1.0, 0.0, 0.0, 0.0, -1.0, -1.0, 0.0, 0.0/), (/3, 3/) )
219  Irr(4)%mat(:,:,7) =  RESHAPE( (/0.0, -1.0, 0.0, 0.0, 0.0, 1.0, -1.0, 0.0, 0.0/), (/3, 3/) )
220  Irr(4)%mat(:,:,8) =  RESHAPE( (/0.0, -1.0, 0.0, 0.0, 0.0, -1.0, 1.0, 0.0, 0.0/), (/3, 3/) )
221  Irr(4)%mat(:,:,9) =  RESHAPE( (/0.0, 0.0, 1.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0/), (/3, 3/) )
222  Irr(4)%mat(:,:,10) =  RESHAPE( (/0.0, 0.0, 1.0, -1.0, 0.0, 0.0, 0.0, -1.0, 0.0/), (/3, 3/) )
223  Irr(4)%mat(:,:,11) =  RESHAPE( (/0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0, 0.0/), (/3, 3/) )
224  Irr(4)%mat(:,:,12) =  RESHAPE( (/0.0, 0.0, -1.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0/), (/3, 3/) )
225  Irr(4)%mat(:,:,13) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 0.0, -1.0, 0.0, -1.0, 0.0/), (/3, 3/) )
226  Irr(4)%mat(:,:,14) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 1.0, 0.0/), (/3, 3/) )
227  Irr(4)%mat(:,:,15) =  RESHAPE( (/0.0, 0.0, -1.0, 0.0, -1.0, 0.0, -1.0, 0.0, 0.0/), (/3, 3/) )
228  Irr(4)%mat(:,:,16) =  RESHAPE( (/0.0, 0.0, 1.0, 0.0, -1.0, 0.0, 1.0, 0.0, 0.0/), (/3, 3/) )
229  Irr(4)%mat(:,:,17) =  RESHAPE( (/0.0, -1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
230  Irr(4)%mat(:,:,18) =  RESHAPE( (/0.0, 1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
231  Irr(4)%mat(:,:,19) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 0.0, -1.0, 0.0, 1.0, 0.0/), (/3, 3/) )
232  Irr(4)%mat(:,:,20) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, -1.0, 0.0/), (/3, 3/) )
233  Irr(4)%mat(:,:,21) =  RESHAPE( (/0.0, 0.0, -1.0, 0.0, 1.0, 0.0, 1.0, 0.0, 0.0/), (/3, 3/) )
234  Irr(4)%mat(:,:,22) =  RESHAPE( (/0.0, 0.0, 1.0, 0.0, 1.0, 0.0, -1.0, 0.0, 0.0/), (/3, 3/) )
235  Irr(4)%mat(:,:,23) =  RESHAPE( (/0.0, -1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
236  Irr(4)%mat(:,:,24) =  RESHAPE( (/0.0, 1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
237 
238  Irr(5)%name = "T2"
239  Irr(5)%dim = 3
240  Irr(5)%nsym = 24
241  ABI_MALLOC(Irr(5)%mat, (3,3,24))
242  Irr(5)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
243  Irr(5)%mat(:,:,2) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
244  Irr(5)%mat(:,:,3) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
245  Irr(5)%mat(:,:,4) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, -1.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
246  Irr(5)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 0.0, 0.0, 0.0, 1.0, 1.0, 0.0, 0.0/), (/3, 3/) )
247  Irr(5)%mat(:,:,6) =  RESHAPE( (/0.0, 1.0, 0.0, 0.0, 0.0, -1.0, -1.0, 0.0, 0.0/), (/3, 3/) )
248  Irr(5)%mat(:,:,7) =  RESHAPE( (/0.0, -1.0, 0.0, 0.0, 0.0, 1.0, -1.0, 0.0, 0.0/), (/3, 3/) )
249  Irr(5)%mat(:,:,8) =  RESHAPE( (/0.0, -1.0, 0.0, 0.0, 0.0, -1.0, 1.0, 0.0, 0.0/), (/3, 3/) )
250  Irr(5)%mat(:,:,9) =  RESHAPE( (/0.0, 0.0, 1.0, 1.0, 0.0, 0.0, 0.0, 1.0, 0.0/), (/3, 3/) )
251  Irr(5)%mat(:,:,10) =  RESHAPE( (/0.0, 0.0, 1.0, -1.0, 0.0, 0.0, 0.0, -1.0, 0.0/), (/3, 3/) )
252  Irr(5)%mat(:,:,11) =  RESHAPE( (/0.0, 0.0, -1.0, 1.0, 0.0, 0.0, 0.0, -1.0, 0.0/), (/3, 3/) )
253  Irr(5)%mat(:,:,12) =  RESHAPE( (/0.0, 0.0, -1.0, -1.0, 0.0, 0.0, 0.0, 1.0, 0.0/), (/3, 3/) )
254  Irr(5)%mat(:,:,13) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, 1.0, 0.0/), (/3, 3/) )
255  Irr(5)%mat(:,:,14) =  RESHAPE( (/1.0, 0.0, 0.0, 0.0, 0.0, -1.0, 0.0, -1.0, 0.0/), (/3, 3/) )
256  Irr(5)%mat(:,:,15) =  RESHAPE( (/0.0, 0.0, 1.0, 0.0, 1.0, 0.0, 1.0, 0.0, 0.0/), (/3, 3/) )
257  Irr(5)%mat(:,:,16) =  RESHAPE( (/0.0, 0.0, -1.0, 0.0, 1.0, 0.0, -1.0, 0.0, 0.0/), (/3, 3/) )
258  Irr(5)%mat(:,:,17) =  RESHAPE( (/0.0, 1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
259  Irr(5)%mat(:,:,18) =  RESHAPE( (/0.0, -1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, 1.0/), (/3, 3/) )
260  Irr(5)%mat(:,:,19) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0, -1.0, 0.0/), (/3, 3/) )
261  Irr(5)%mat(:,:,20) =  RESHAPE( (/-1.0, 0.0, 0.0, 0.0, 0.0, -1.0, 0.0, 1.0, 0.0/), (/3, 3/) )
262  Irr(5)%mat(:,:,21) =  RESHAPE( (/0.0, 0.0, 1.0, 0.0, -1.0, 0.0, -1.0, 0.0, 0.0/), (/3, 3/) )
263  Irr(5)%mat(:,:,22) =  RESHAPE( (/0.0, 0.0, -1.0, 0.0, -1.0, 0.0, 1.0, 0.0, 0.0/), (/3, 3/) )
264  Irr(5)%mat(:,:,23) =  RESHAPE( (/0.0, 1.0, 0.0, -1.0, 0.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
265  Irr(5)%mat(:,:,24) =  RESHAPE( (/0.0, -1.0, 0.0, 1.0, 0.0, 0.0, 0.0, 0.0, -1.0/), (/3, 3/) )
266 
267  RETURN
268   if (.FALSE.) write(std_out,*) j
269  end subroutine ptg_Td