TABLE OF CONTENTS


ABINIT/ptg_D4h [ Functions ]

[ Top ] [ Functions ]

NAME

 ptg_D4h

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  D4h (4/mmm)
 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_D4h
 42 contains
 43 !!**
 44 
 45 
 46 
 47  subroutine ptg_D4h (nsym,nclass,sym,class_ids,class_names,Irr)
 48  use defs_basis
 49  use m_abicore
 50  use m_defs_ptgroups,  only : irrep_t
 51 
 52 !This section has been created automatically by the script Abilint (TD).
 53 !Do not modify the following lines by hand.
 54 #undef ABI_FUNC
 55 #define ABI_FUNC 'ptg_D4h'
 56 !End of the abilint section
 57 
 58  implicit none
 59 !Arguments ------------------------------------
 60  integer,intent(out) :: nclass,nsym 
 61  !arrays
 62  integer,allocatable,intent(out) :: sym(:,:,:), class_ids(:,:)
 63  character(len=5),allocatable,intent(out) :: class_names(:)
 64  type(irrep_t),allocatable,intent(out) :: Irr(:)
 65  !Local variables-------------------------------
 66  complex(dpc) :: j=(0.0_dp,1.0_dp) 
 67  ! ********************************************************************************
 68 ! List of symmetries packed in classes
 69  nsym = 16
 70  ABI_MALLOC(sym, (3,3,nsym))
 71  sym(:,:,1) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
 72  sym(:,:,2) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
 73  sym(:,:,3) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 74  sym(:,:,4) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 75  sym(:,:,5) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
 76  sym(:,:,6) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
 77  sym(:,:,7) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 78  sym(:,:,8) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 79  sym(:,:,9) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
 80  sym(:,:,10) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
 81  sym(:,:,11) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 82  sym(:,:,12) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 83  sym(:,:,13) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
 84  sym(:,:,14) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
 85  sym(:,:,15) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 86  sym(:,:,16) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 87 
 88 ! Number of classes and corresponding indeces
 89  nclass = 10
 90  ABI_MALLOC(class_ids, (2,nclass))
 91  class_ids(1,1) = 1
 92  class_ids(2,1) = 1
 93  class_ids(1,2) = 2
 94  class_ids(2,2) = 2
 95  class_ids(1,3) = 3
 96  class_ids(2,3) = 4
 97  class_ids(1,4) = 5
 98  class_ids(2,4) = 6
 99  class_ids(1,5) = 7
100  class_ids(2,5) = 8
101  class_ids(1,6) = 9
102  class_ids(2,6) = 9
103  class_ids(1,7) = 10
104  class_ids(2,7) = 10
105  class_ids(1,8) = 11
106  class_ids(2,8) = 12
107  class_ids(1,9) = 13
108  class_ids(2,9) = 14
109  class_ids(1,10) = 15
110  class_ids(2,10) = 16
111 
112 ABI_MALLOC(class_names,(10))
113  class_names(1) = "1+" 
114  class_names(2) = "2+" 
115  class_names(3) = "4+" 
116  class_names(4) = "2+" 
117  class_names(5) = "2+" 
118  class_names(6) = "-2-" 
119  class_names(7) = "-2+" 
120  class_names(8) = "-4+" 
121  class_names(9) = "-2+" 
122  class_names(10) = "-2+" 
123 
124 ! List of irreducible representations.
125  ABI_DT_MALLOC(Irr, (10))
126  Irr(1)%name = "A1g"
127  Irr(1)%dim = 1
128  Irr(1)%nsym = 16
129  ABI_MALLOC(Irr(1)%mat, (1,1,16))
130  Irr(1)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
131  Irr(1)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
132  Irr(1)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
133  Irr(1)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
134  Irr(1)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
135  Irr(1)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
136  Irr(1)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
137  Irr(1)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
138  Irr(1)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
139  Irr(1)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
140  Irr(1)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
141  Irr(1)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
142  Irr(1)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
143  Irr(1)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
144  Irr(1)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
145  Irr(1)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
146 
147  Irr(2)%name = "A1u"
148  Irr(2)%dim = 1
149  Irr(2)%nsym = 16
150  ABI_MALLOC(Irr(2)%mat, (1,1,16))
151  Irr(2)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
152  Irr(2)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
153  Irr(2)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
154  Irr(2)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
155  Irr(2)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
156  Irr(2)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
157  Irr(2)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
158  Irr(2)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
159  Irr(2)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
160  Irr(2)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
161  Irr(2)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
162  Irr(2)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
163  Irr(2)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
164  Irr(2)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
165  Irr(2)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
166  Irr(2)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
167 
168  Irr(3)%name = "A2g"
169  Irr(3)%dim = 1
170  Irr(3)%nsym = 16
171  ABI_MALLOC(Irr(3)%mat, (1,1,16))
172  Irr(3)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
173  Irr(3)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
174  Irr(3)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
175  Irr(3)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
176  Irr(3)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
177  Irr(3)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
178  Irr(3)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
179  Irr(3)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
180  Irr(3)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
181  Irr(3)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
182  Irr(3)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
183  Irr(3)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
184  Irr(3)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
185  Irr(3)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
186  Irr(3)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
187  Irr(3)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
188 
189  Irr(4)%name = "A2u"
190  Irr(4)%dim = 1
191  Irr(4)%nsym = 16
192  ABI_MALLOC(Irr(4)%mat, (1,1,16))
193  Irr(4)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
194  Irr(4)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
195  Irr(4)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
196  Irr(4)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
197  Irr(4)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
198  Irr(4)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
199  Irr(4)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
200  Irr(4)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
201  Irr(4)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
202  Irr(4)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
203  Irr(4)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
204  Irr(4)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
205  Irr(4)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
206  Irr(4)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
207  Irr(4)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
208  Irr(4)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
209 
210  Irr(5)%name = "B1g"
211  Irr(5)%dim = 1
212  Irr(5)%nsym = 16
213  ABI_MALLOC(Irr(5)%mat, (1,1,16))
214  Irr(5)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
215  Irr(5)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
216  Irr(5)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
217  Irr(5)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
218  Irr(5)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
219  Irr(5)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
220  Irr(5)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
221  Irr(5)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
222  Irr(5)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
223  Irr(5)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
224  Irr(5)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
225  Irr(5)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
226  Irr(5)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
227  Irr(5)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
228  Irr(5)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
229  Irr(5)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
230 
231  Irr(6)%name = "B1u"
232  Irr(6)%dim = 1
233  Irr(6)%nsym = 16
234  ABI_MALLOC(Irr(6)%mat, (1,1,16))
235  Irr(6)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
236  Irr(6)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
237  Irr(6)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
238  Irr(6)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
239  Irr(6)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
240  Irr(6)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
241  Irr(6)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
242  Irr(6)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
243  Irr(6)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
244  Irr(6)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
245  Irr(6)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
246  Irr(6)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
247  Irr(6)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
248  Irr(6)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
249  Irr(6)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
250  Irr(6)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
251 
252  Irr(7)%name = "B2g"
253  Irr(7)%dim = 1
254  Irr(7)%nsym = 16
255  ABI_MALLOC(Irr(7)%mat, (1,1,16))
256  Irr(7)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
257  Irr(7)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
258  Irr(7)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
259  Irr(7)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
260  Irr(7)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
261  Irr(7)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
262  Irr(7)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
263  Irr(7)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
264  Irr(7)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
265  Irr(7)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
266  Irr(7)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
267  Irr(7)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
268  Irr(7)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
269  Irr(7)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
270  Irr(7)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
271  Irr(7)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
272 
273  Irr(8)%name = "B2u"
274  Irr(8)%dim = 1
275  Irr(8)%nsym = 16
276  ABI_MALLOC(Irr(8)%mat, (1,1,16))
277  Irr(8)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
278  Irr(8)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
279  Irr(8)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
280  Irr(8)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
281  Irr(8)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
282  Irr(8)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
283  Irr(8)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
284  Irr(8)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
285  Irr(8)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
286  Irr(8)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
287  Irr(8)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
288  Irr(8)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
289  Irr(8)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
290  Irr(8)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
291  Irr(8)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
292  Irr(8)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
293 
294  Irr(9)%name = "Eu"
295  Irr(9)%dim = 2
296  Irr(9)%nsym = 16
297  ABI_MALLOC(Irr(9)%mat, (2,2,16))
298  Irr(9)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
299  Irr(9)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
300  Irr(9)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
301  Irr(9)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
302  Irr(9)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
303  Irr(9)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
304  Irr(9)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
305  Irr(9)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
306  Irr(9)%mat(:,:,9) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
307  Irr(9)%mat(:,:,10) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
308  Irr(9)%mat(:,:,11) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
309  Irr(9)%mat(:,:,12) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
310  Irr(9)%mat(:,:,13) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
311  Irr(9)%mat(:,:,14) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
312  Irr(9)%mat(:,:,15) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
313  Irr(9)%mat(:,:,16) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
314 
315  Irr(10)%name = "Eg"
316  Irr(10)%dim = 2
317  Irr(10)%nsym = 16
318  ABI_MALLOC(Irr(10)%mat, (2,2,16))
319  Irr(10)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
320  Irr(10)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
321  Irr(10)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
322  Irr(10)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
323  Irr(10)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
324  Irr(10)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
325  Irr(10)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
326  Irr(10)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
327  Irr(10)%mat(:,:,9) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
328  Irr(10)%mat(:,:,10) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
329  Irr(10)%mat(:,:,11) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
330  Irr(10)%mat(:,:,12) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
331  Irr(10)%mat(:,:,13) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
332  Irr(10)%mat(:,:,14) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
333  Irr(10)%mat(:,:,15) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
334  Irr(10)%mat(:,:,16) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
335 
336  RETURN
337   if (.FALSE.) write(std_out,*) j
338  end subroutine ptg_D4h