LCOV - code coverage report
Current view: top level - src/43_ptgroups - m_ptg_D4h.F90 (source / functions) Coverage Total Hit
Test: coverage.info Lines: 0.0 % 254 0
Test Date: 2026-09-19 15:24:51 Functions: 0.0 % 1 0

            Line data    Source code
       1              : !!****f* ABINIT/ptg_D4h
       2              : !!
       3              : !! NAME
       4              : !! ptg_D4h
       5              : !!
       6              : !! FUNCTION
       7              : !!
       8              : !! COPYRIGHT
       9              : !! Copyright (C) 2010-2026 ABINIT group (MG)
      10              : !! This file is distributed under the terms of the
      11              : !! GNU General Public License, see ~abinit/COPYING
      12              : !! or http://www.gnu.org/copyleft/gpl.txt .
      13              : !! For the initials of contributors, see ~abinit/doc/developers/contributors.txt .
      14              : !!
      15              : !! INPUTS
      16              : !!
      17              : !! OUTPUT
      18              : !!
      19              : !! SOURCE
      20              : !!
      21              : !********************************************************************************
      22              : ! This include file has been automatically generated by the script ptg.py
      23              : ! Do not edit! Change the script source instead.
      24              : !********************************************************************************
      25              : 
      26              : ! Point group name  D4h (4/mmm)
      27              : 
      28              : #if defined HAVE_CONFIG_H
      29              : #include "config.h"
      30              : #endif
      31              : 
      32              : #include "abi_common.h"
      33              : 
      34              : module m_ptg_D4h
      35              : 
      36              : #ifdef FC_INTEL
      37              : !DEC$ NOOPTIMIZE
      38              : #endif
      39              : 
      40              : contains
      41              : !!**
      42              : 
      43            0 :  subroutine ptg_D4h (nsym,nclass,sym,class_ids,class_names,Irr)
      44              :  use defs_basis
      45              :  use m_abicore
      46              :  use m_defs_ptgroups,  only : irrep_t
      47              :  implicit none
      48              : !Arguments ------------------------------------
      49              :  integer,intent(out) :: nclass,nsym
      50              :  !arrays
      51              :  integer,allocatable,intent(out) :: sym(:,:,:), class_ids(:,:)
      52              :  character(len=5),allocatable,intent(out) :: class_names(:)
      53              :  type(irrep_t),allocatable,intent(out) :: Irr(:)
      54              :  !Local variables-------------------------------
      55              :  complex(dp) :: j=(0.0_dp,1.0_dp)
      56              :  ! ********************************************************************************
      57              : ! List of symmetries packed in classes
      58            0 :  nsym = 16
      59            0 :  ABI_MALLOC(sym, (3,3,nsym))
      60            0 :  sym(:,:,1) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
      61            0 :  sym(:,:,2) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
      62            0 :  sym(:,:,3) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
      63            0 :  sym(:,:,4) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
      64            0 :  sym(:,:,5) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
      65            0 :  sym(:,:,6) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
      66            0 :  sym(:,:,7) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
      67            0 :  sym(:,:,8) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
      68            0 :  sym(:,:,9) = RESHAPE( (/-1, 0, 0, 0, -1, 0, 0, 0, -1/) ,(/3,3/) )
      69            0 :  sym(:,:,10) = RESHAPE( (/1, 0, 0, 0, 1, 0, 0, 0, -1/) ,(/3,3/) )
      70            0 :  sym(:,:,11) = RESHAPE( (/0, -1, 0, 1, 0, 0, 0, 0, -1/) ,(/3,3/) )
      71            0 :  sym(:,:,12) = RESHAPE( (/0, 1, 0, -1, 0, 0, 0, 0, -1/) ,(/3,3/) )
      72            0 :  sym(:,:,13) = RESHAPE( (/1, 0, 0, 0, -1, 0, 0, 0, 1/) ,(/3,3/) )
      73            0 :  sym(:,:,14) = RESHAPE( (/-1, 0, 0, 0, 1, 0, 0, 0, 1/) ,(/3,3/) )
      74            0 :  sym(:,:,15) = RESHAPE( (/0, -1, 0, -1, 0, 0, 0, 0, 1/) ,(/3,3/) )
      75            0 :  sym(:,:,16) = RESHAPE( (/0, 1, 0, 1, 0, 0, 0, 0, 1/) ,(/3,3/) )
      76              : 
      77              : ! Number of classes and corresponding indices
      78            0 :  nclass = 10
      79            0 :  ABI_MALLOC(class_ids, (2,nclass))
      80            0 :  class_ids(1,1) = 1
      81            0 :  class_ids(2,1) = 1
      82            0 :  class_ids(1,2) = 2
      83            0 :  class_ids(2,2) = 2
      84            0 :  class_ids(1,3) = 3
      85            0 :  class_ids(2,3) = 4
      86            0 :  class_ids(1,4) = 5
      87            0 :  class_ids(2,4) = 6
      88            0 :  class_ids(1,5) = 7
      89            0 :  class_ids(2,5) = 8
      90            0 :  class_ids(1,6) = 9
      91            0 :  class_ids(2,6) = 9
      92            0 :  class_ids(1,7) = 10
      93            0 :  class_ids(2,7) = 10
      94            0 :  class_ids(1,8) = 11
      95            0 :  class_ids(2,8) = 12
      96            0 :  class_ids(1,9) = 13
      97            0 :  class_ids(2,9) = 14
      98            0 :  class_ids(1,10) = 15
      99            0 :  class_ids(2,10) = 16
     100              : 
     101            0 : ABI_MALLOC(class_names,(10))
     102            0 :  class_names(1) = "1+"
     103            0 :  class_names(2) = "2+"
     104            0 :  class_names(3) = "4+"
     105            0 :  class_names(4) = "2+"
     106            0 :  class_names(5) = "2+"
     107            0 :  class_names(6) = "-2-"
     108            0 :  class_names(7) = "-2+"
     109            0 :  class_names(8) = "-4+"
     110            0 :  class_names(9) = "-2+"
     111            0 :  class_names(10) = "-2+"
     112              : 
     113              : ! List of irreducible representations.
     114            0 :  ABI_MALLOC(Irr, (10))
     115            0 :  Irr(1)%name = "A1g"
     116            0 :  Irr(1)%dim = 1
     117            0 :  Irr(1)%nsym = 16
     118            0 :  ABI_MALLOC(Irr(1)%mat, (1,1,16))
     119            0 :  Irr(1)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     120            0 :  Irr(1)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     121            0 :  Irr(1)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
     122            0 :  Irr(1)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
     123            0 :  Irr(1)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
     124            0 :  Irr(1)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
     125            0 :  Irr(1)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
     126            0 :  Irr(1)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
     127            0 :  Irr(1)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
     128            0 :  Irr(1)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
     129            0 :  Irr(1)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
     130            0 :  Irr(1)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
     131            0 :  Irr(1)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
     132            0 :  Irr(1)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
     133            0 :  Irr(1)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
     134            0 :  Irr(1)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
     135              : 
     136            0 :  Irr(2)%name = "A1u"
     137            0 :  Irr(2)%dim = 1
     138            0 :  Irr(2)%nsym = 16
     139            0 :  ABI_MALLOC(Irr(2)%mat, (1,1,16))
     140            0 :  Irr(2)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     141            0 :  Irr(2)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     142            0 :  Irr(2)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
     143            0 :  Irr(2)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
     144            0 :  Irr(2)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
     145            0 :  Irr(2)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
     146            0 :  Irr(2)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
     147            0 :  Irr(2)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
     148            0 :  Irr(2)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
     149            0 :  Irr(2)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
     150            0 :  Irr(2)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
     151            0 :  Irr(2)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
     152            0 :  Irr(2)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
     153            0 :  Irr(2)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
     154            0 :  Irr(2)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
     155            0 :  Irr(2)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
     156              : 
     157            0 :  Irr(3)%name = "A2g"
     158            0 :  Irr(3)%dim = 1
     159            0 :  Irr(3)%nsym = 16
     160            0 :  ABI_MALLOC(Irr(3)%mat, (1,1,16))
     161            0 :  Irr(3)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     162            0 :  Irr(3)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     163            0 :  Irr(3)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
     164            0 :  Irr(3)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
     165            0 :  Irr(3)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
     166            0 :  Irr(3)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
     167            0 :  Irr(3)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
     168            0 :  Irr(3)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
     169            0 :  Irr(3)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
     170            0 :  Irr(3)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
     171            0 :  Irr(3)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
     172            0 :  Irr(3)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
     173            0 :  Irr(3)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
     174            0 :  Irr(3)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
     175            0 :  Irr(3)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
     176            0 :  Irr(3)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
     177              : 
     178            0 :  Irr(4)%name = "A2u"
     179            0 :  Irr(4)%dim = 1
     180            0 :  Irr(4)%nsym = 16
     181            0 :  ABI_MALLOC(Irr(4)%mat, (1,1,16))
     182            0 :  Irr(4)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     183            0 :  Irr(4)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     184            0 :  Irr(4)%mat(:,:,3) =  RESHAPE( (/1.0/), (/1, 1/) )
     185            0 :  Irr(4)%mat(:,:,4) =  RESHAPE( (/1.0/), (/1, 1/) )
     186            0 :  Irr(4)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
     187            0 :  Irr(4)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
     188            0 :  Irr(4)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
     189            0 :  Irr(4)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
     190            0 :  Irr(4)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
     191            0 :  Irr(4)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
     192            0 :  Irr(4)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
     193            0 :  Irr(4)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
     194            0 :  Irr(4)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
     195            0 :  Irr(4)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
     196            0 :  Irr(4)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
     197            0 :  Irr(4)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
     198              : 
     199            0 :  Irr(5)%name = "B1g"
     200            0 :  Irr(5)%dim = 1
     201            0 :  Irr(5)%nsym = 16
     202            0 :  ABI_MALLOC(Irr(5)%mat, (1,1,16))
     203            0 :  Irr(5)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     204            0 :  Irr(5)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     205            0 :  Irr(5)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
     206            0 :  Irr(5)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
     207            0 :  Irr(5)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
     208            0 :  Irr(5)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
     209            0 :  Irr(5)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
     210            0 :  Irr(5)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
     211            0 :  Irr(5)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
     212            0 :  Irr(5)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
     213            0 :  Irr(5)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
     214            0 :  Irr(5)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
     215            0 :  Irr(5)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
     216            0 :  Irr(5)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
     217            0 :  Irr(5)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
     218            0 :  Irr(5)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
     219              : 
     220            0 :  Irr(6)%name = "B1u"
     221            0 :  Irr(6)%dim = 1
     222            0 :  Irr(6)%nsym = 16
     223            0 :  ABI_MALLOC(Irr(6)%mat, (1,1,16))
     224            0 :  Irr(6)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     225            0 :  Irr(6)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     226            0 :  Irr(6)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
     227            0 :  Irr(6)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
     228            0 :  Irr(6)%mat(:,:,5) =  RESHAPE( (/1.0/), (/1, 1/) )
     229            0 :  Irr(6)%mat(:,:,6) =  RESHAPE( (/1.0/), (/1, 1/) )
     230            0 :  Irr(6)%mat(:,:,7) =  RESHAPE( (/-1.0/), (/1, 1/) )
     231            0 :  Irr(6)%mat(:,:,8) =  RESHAPE( (/-1.0/), (/1, 1/) )
     232            0 :  Irr(6)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
     233            0 :  Irr(6)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
     234            0 :  Irr(6)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
     235            0 :  Irr(6)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
     236            0 :  Irr(6)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
     237            0 :  Irr(6)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
     238            0 :  Irr(6)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
     239            0 :  Irr(6)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
     240              : 
     241            0 :  Irr(7)%name = "B2g"
     242            0 :  Irr(7)%dim = 1
     243            0 :  Irr(7)%nsym = 16
     244            0 :  ABI_MALLOC(Irr(7)%mat, (1,1,16))
     245            0 :  Irr(7)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     246            0 :  Irr(7)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     247            0 :  Irr(7)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
     248            0 :  Irr(7)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
     249            0 :  Irr(7)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
     250            0 :  Irr(7)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
     251            0 :  Irr(7)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
     252            0 :  Irr(7)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
     253            0 :  Irr(7)%mat(:,:,9) =  RESHAPE( (/1.0/), (/1, 1/) )
     254            0 :  Irr(7)%mat(:,:,10) =  RESHAPE( (/1.0/), (/1, 1/) )
     255            0 :  Irr(7)%mat(:,:,11) =  RESHAPE( (/-1.0/), (/1, 1/) )
     256            0 :  Irr(7)%mat(:,:,12) =  RESHAPE( (/-1.0/), (/1, 1/) )
     257            0 :  Irr(7)%mat(:,:,13) =  RESHAPE( (/-1.0/), (/1, 1/) )
     258            0 :  Irr(7)%mat(:,:,14) =  RESHAPE( (/-1.0/), (/1, 1/) )
     259            0 :  Irr(7)%mat(:,:,15) =  RESHAPE( (/1.0/), (/1, 1/) )
     260            0 :  Irr(7)%mat(:,:,16) =  RESHAPE( (/1.0/), (/1, 1/) )
     261              : 
     262            0 :  Irr(8)%name = "B2u"
     263            0 :  Irr(8)%dim = 1
     264            0 :  Irr(8)%nsym = 16
     265            0 :  ABI_MALLOC(Irr(8)%mat, (1,1,16))
     266            0 :  Irr(8)%mat(:,:,1) =  RESHAPE( (/1.0/), (/1, 1/) )
     267            0 :  Irr(8)%mat(:,:,2) =  RESHAPE( (/1.0/), (/1, 1/) )
     268            0 :  Irr(8)%mat(:,:,3) =  RESHAPE( (/-1.0/), (/1, 1/) )
     269            0 :  Irr(8)%mat(:,:,4) =  RESHAPE( (/-1.0/), (/1, 1/) )
     270            0 :  Irr(8)%mat(:,:,5) =  RESHAPE( (/-1.0/), (/1, 1/) )
     271            0 :  Irr(8)%mat(:,:,6) =  RESHAPE( (/-1.0/), (/1, 1/) )
     272            0 :  Irr(8)%mat(:,:,7) =  RESHAPE( (/1.0/), (/1, 1/) )
     273            0 :  Irr(8)%mat(:,:,8) =  RESHAPE( (/1.0/), (/1, 1/) )
     274            0 :  Irr(8)%mat(:,:,9) =  RESHAPE( (/-1.0/), (/1, 1/) )
     275            0 :  Irr(8)%mat(:,:,10) =  RESHAPE( (/-1.0/), (/1, 1/) )
     276            0 :  Irr(8)%mat(:,:,11) =  RESHAPE( (/1.0/), (/1, 1/) )
     277            0 :  Irr(8)%mat(:,:,12) =  RESHAPE( (/1.0/), (/1, 1/) )
     278            0 :  Irr(8)%mat(:,:,13) =  RESHAPE( (/1.0/), (/1, 1/) )
     279            0 :  Irr(8)%mat(:,:,14) =  RESHAPE( (/1.0/), (/1, 1/) )
     280            0 :  Irr(8)%mat(:,:,15) =  RESHAPE( (/-1.0/), (/1, 1/) )
     281            0 :  Irr(8)%mat(:,:,16) =  RESHAPE( (/-1.0/), (/1, 1/) )
     282              : 
     283            0 :  Irr(9)%name = "Eu"
     284            0 :  Irr(9)%dim = 2
     285            0 :  Irr(9)%nsym = 16
     286            0 :  ABI_MALLOC(Irr(9)%mat, (2,2,16))
     287            0 :  Irr(9)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
     288            0 :  Irr(9)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
     289            0 :  Irr(9)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
     290            0 :  Irr(9)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
     291            0 :  Irr(9)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
     292            0 :  Irr(9)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
     293            0 :  Irr(9)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
     294            0 :  Irr(9)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
     295            0 :  Irr(9)%mat(:,:,9) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
     296            0 :  Irr(9)%mat(:,:,10) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
     297            0 :  Irr(9)%mat(:,:,11) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
     298            0 :  Irr(9)%mat(:,:,12) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
     299            0 :  Irr(9)%mat(:,:,13) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
     300            0 :  Irr(9)%mat(:,:,14) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
     301            0 :  Irr(9)%mat(:,:,15) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
     302            0 :  Irr(9)%mat(:,:,16) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
     303              : 
     304            0 :  Irr(10)%name = "Eg"
     305            0 :  Irr(10)%dim = 2
     306            0 :  Irr(10)%nsym = 16
     307            0 :  ABI_MALLOC(Irr(10)%mat, (2,2,16))
     308            0 :  Irr(10)%mat(:,:,1) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
     309            0 :  Irr(10)%mat(:,:,2) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
     310            0 :  Irr(10)%mat(:,:,3) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
     311            0 :  Irr(10)%mat(:,:,4) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
     312            0 :  Irr(10)%mat(:,:,5) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
     313            0 :  Irr(10)%mat(:,:,6) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
     314            0 :  Irr(10)%mat(:,:,7) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
     315            0 :  Irr(10)%mat(:,:,8) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
     316            0 :  Irr(10)%mat(:,:,9) =  RESHAPE( (/1.0, 0.0, 0.0, 1.0/), (/2, 2/) )
     317            0 :  Irr(10)%mat(:,:,10) =  RESHAPE( (/-1.0, 0.0, 0.0, -1.0/), (/2, 2/) )
     318            0 :  Irr(10)%mat(:,:,11) =  RESHAPE( (/1*j, 0*j, 0*j, -0-1*j/), (/2, 2/) )
     319            0 :  Irr(10)%mat(:,:,12) =  RESHAPE( (/-0-1*j, 0*j, 0*j, 1*j/), (/2, 2/) )
     320            0 :  Irr(10)%mat(:,:,13) =  RESHAPE( (/0.0, 1.0, 1.0, 0.0/), (/2, 2/) )
     321            0 :  Irr(10)%mat(:,:,14) =  RESHAPE( (/0.0, -1.0, -1.0, 0.0/), (/2, 2/) )
     322            0 :  Irr(10)%mat(:,:,15) =  RESHAPE( (/0*j, 1*j, -0-1*j, 0*j/), (/2, 2/) )
     323            0 :  Irr(10)%mat(:,:,16) =  RESHAPE( (/0*j, -0-1*j, 1*j, 0*j/), (/2, 2/) )
     324              : 
     325            0 :  RETURN
     326              :   if (.FALSE.) write(std_out,*) j
     327              :  end subroutine ptg_D4h
     328              : !!***
     329              : 
     330              : end module m_ptg_D4h
     331              : !!***
        

Generated by: LCOV version 2.3-1