!{\src2tex{textfont=tt}}
!!****f* ABINIT/ptg_D4h
!!
!! NAME
!! ptg_D4h
!!
!! FUNCTION
!!
!! COPYRIGHT
!! Copyright (C) 2010-2014 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
!!
!********************************************************************************
! This include file has been automatically generated by the script ptg.py
! Do not edit! Change the script source instead.
!********************************************************************************

! Point group name  D4h (4/mmm)

#if defined HAVE_CONFIG_H
#include "config.h"
#endif

#include "abi_common.h"

 subroutine ptg_D4h (nsym,nclass,sym,class_ids,class_names,Irr)
 use defs_basis
 use m_profiling
 use m_defs_ptgroups,  only : irrep_t

!This section has been created automatically by the script Abilint (TD).
!Do not modify the following lines by hand.
#undef ABI_FUNC
#define ABI_FUNC 'ptg_D4h'
!End of the abilint section

 implicit none
!Arguments ------------------------------------
 integer,intent(out) :: nclass,nsym 
 !arrays
 integer,allocatable,intent(out) :: sym(:,:,:), class_ids(:,:)
 character(len=5),allocatable,intent(out) :: class_names(:)
 type(irrep_t),allocatable,intent(out) :: Irr(:)
 !Local variables-------------------------------
 complex(dpc) :: j=(0.0_dp,1.0_dp) 
 ! ********************************************************************************
! List of symmetries packed in classes
 nsym = 16
 ABI_MALLOC(sym, (3,3,nsym))
 sym(:,:,1) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,2) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,3) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,4) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,5) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,6) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,7) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,8) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,9) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,10) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,11) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,12) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
 sym(:,:,13) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,14) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,15) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
 sym(:,:,16) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )

! Number of classes and corresponding indeces
 nclass = 10
 ABI_MALLOC(class_ids, (2,nclass))
 class_ids(1,1) = 1
 class_ids(2,1) = 1
 class_ids(1,2) = 2
 class_ids(2,2) = 2
 class_ids(1,3) = 3
 class_ids(2,3) = 4
 class_ids(1,4) = 5
 class_ids(2,4) = 6
 class_ids(1,5) = 7
 class_ids(2,5) = 8
 class_ids(1,6) = 9
 class_ids(2,6) = 9
 class_ids(1,7) = 10
 class_ids(2,7) = 10
 class_ids(1,8) = 11
 class_ids(2,8) = 12
 class_ids(1,9) = 13
 class_ids(2,9) = 14
 class_ids(1,10) = 15
 class_ids(2,10) = 16

ABI_MALLOC(class_names,(10))
 class_names(1) = "1+" 
 class_names(2) = "2+" 
 class_names(3) = "4+" 
 class_names(4) = "2+" 
 class_names(5) = "2+" 
 class_names(6) = "-2-" 
 class_names(7) = "-2+" 
 class_names(8) = "-4+" 
 class_names(9) = "-2+" 
 class_names(10) = "-2+" 

! List of irreducible representations.
 ABI_DT_MALLOC(Irr, (10))
 Irr(1)%name = "A1g"
 Irr(1)%dim = 1
 Irr(1)%nsym = 16
 ABI_MALLOC(Irr(1)%mat, (1,1,16))
 Irr(1)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(1)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )

 Irr(2)%name = "A1u"
 Irr(2)%dim = 1
 Irr(2)%nsym = 16
 ABI_MALLOC(Irr(2)%mat, (1,1,16))
 Irr(2)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(2)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )

 Irr(3)%name = "A2g"
 Irr(3)%dim = 1
 Irr(3)%nsym = 16
 ABI_MALLOC(Irr(3)%mat, (1,1,16))
 Irr(3)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(3)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )

 Irr(4)%name = "A2u"
 Irr(4)%dim = 1
 Irr(4)%nsym = 16
 ABI_MALLOC(Irr(4)%mat, (1,1,16))
 Irr(4)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(4)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )

 Irr(5)%name = "B1g"
 Irr(5)%dim = 1
 Irr(5)%nsym = 16
 ABI_MALLOC(Irr(5)%mat, (1,1,16))
 Irr(5)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(5)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )

 Irr(6)%name = "B1u"
 Irr(6)%dim = 1
 Irr(6)%nsym = 16
 ABI_MALLOC(Irr(6)%mat, (1,1,16))
 Irr(6)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(6)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )

 Irr(7)%name = "B2g"
 Irr(7)%dim = 1
 Irr(7)%nsym = 16
 ABI_MALLOC(Irr(7)%mat, (1,1,16))
 Irr(7)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(7)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )

 Irr(8)%name = "B2u"
 Irr(8)%dim = 1
 Irr(8)%nsym = 16
 ABI_MALLOC(Irr(8)%mat, (1,1,16))
 Irr(8)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
 Irr(8)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )

 Irr(9)%name = "Eu"
 Irr(9)%dim = 2
 Irr(9)%nsym = 16
 ABI_MALLOC(Irr(9)%mat, (2,2,16))
 Irr(9)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
 Irr(9)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
 Irr(9)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
 Irr(9)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
 Irr(9)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
 Irr(9)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
 Irr(9)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
 Irr(9)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
 Irr(9)%mat(:,:,9) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
 Irr(9)%mat(:,:,10) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
 Irr(9)%mat(:,:,11) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
 Irr(9)%mat(:,:,12) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
 Irr(9)%mat(:,:,13) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
 Irr(9)%mat(:,:,14) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
 Irr(9)%mat(:,:,15) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
 Irr(9)%mat(:,:,16) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )

 Irr(10)%name = "Eg"
 Irr(10)%dim = 2
 Irr(10)%nsym = 16
 ABI_MALLOC(Irr(10)%mat, (2,2,16))
 Irr(10)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
 Irr(10)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
 Irr(10)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
 Irr(10)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
 Irr(10)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
 Irr(10)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
 Irr(10)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
 Irr(10)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
 Irr(10)%mat(:,:,9) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
 Irr(10)%mat(:,:,10) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
 Irr(10)%mat(:,:,11) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
 Irr(10)%mat(:,:,12) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
 Irr(10)%mat(:,:,13) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
 Irr(10)%mat(:,:,14) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
 Irr(10)%mat(:,:,15) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
 Irr(10)%mat(:,:,16) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )

 RETURN
  if (.FALSE.) write(std_out,*) j
 end subroutine ptg_D4h 
!!***
