← Back to branch summary

~ubuntu-branches/ubuntu/utopic/nwchem/utopic

« back to all changes in this revision

Viewing changes to src/symmetry/gensym_extra.f

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
http://bugs.debian.org/696361
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/symmetry/gensym_extra.f
 
1
C$Id: gensym_extra.f 22503 2012-05-20 06:58:57Z d3y133 $
 
2
 
 
3
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
4
c
 
5
c    The matrix SYMOPS contains the matrix reps. of all group operators
 
6
c   except the identity. The variable NOPS holds the number of operators
 
7
c   in SYMOPS.
 
8
c
 
9
c***********************************************************************
 
10
      subroutine gensym_extra(itype,numgrp,numset,symops,nops,oprint,
 
11
     >                        group_name)
 
12
      implicit none
 
13
      integer maxops
 
14
      parameter(maxops=192)
 
15
 
 
16
      integer itype,numgrp,numset
 
17
      real*8 symops(maxops*3,4)
 
18
      integer nops
 
19
      character*(*) group_name
 
20
      logical oprint
 
21
 
 
22
*     **** local varialbles ****
 
23
      integer iop
 
24
 
 
25
      nops = 0
 
26
 
 
27
******* Extra Triclinic groups ******
 
28
c 
 
29
c    A1 
 
30
c 
 
31
      if (numgrp.eq.231) then
 
32
      group_name = 'A1'
 
33
      nops = 1
 
34
      iop = 0
 
35
      symops(3*iop+1,1) = 1.0d0
 
36
      symops(3*iop+2,1) = 0.0d0
 
37
      symops(3*iop+3,1) = 0.0d0
 
38
      symops(3*iop+1,2) = 0.0d0
 
39
      symops(3*iop+2,2) = 1.0d0
 
40
      symops(3*iop+3,2) = 0.0d0
 
41
      symops(3*iop+1,3) = 0.0d0
 
42
      symops(3*iop+2,3) = 0.0d0
 
43
      symops(3*iop+3,3) = 1.0d0
 
44
      symops(3*iop+1,4) = 0.0d0
 
45
      symops(3*iop+2,4) = 0.5d0
 
46
      symops(3*iop+3,4) = 0.5d0
 
47
      
 
48
      end if
 
49
c 
 
50
c    B1 
 
51
c 
 
52
      if (numgrp.eq.232) then
 
53
      group_name = 'B1'
 
54
      nops = 1
 
55
      iop = 0
 
56
      symops(3*iop+1,1) = 1.0d0
 
57
      symops(3*iop+2,1) = 0.0d0
 
58
      symops(3*iop+3,1) = 0.0d0
 
59
      symops(3*iop+1,2) = 0.0d0
 
60
      symops(3*iop+2,2) = 1.0d0
 
61
      symops(3*iop+3,2) = 0.0d0
 
62
      symops(3*iop+1,3) = 0.0d0
 
63
      symops(3*iop+2,3) = 0.0d0
 
64
      symops(3*iop+3,3) = 1.0d0
 
65
      symops(3*iop+1,4) = 0.5d0
 
66
      symops(3*iop+2,4) = 0.0d0
 
67
      symops(3*iop+3,4) = 0.5d0
 
68
      end if
 
69
c 
 
70
c    C1 
 
71
c 
 
72
      if (numgrp.eq.233) then
 
73
      group_name = 'C1'
 
74
      nops = 1
 
75
      iop = 0
 
76
      symops(3*iop+1,1) = 1.0d0
 
77
      symops(3*iop+2,1) = 0.0d0
 
78
      symops(3*iop+3,1) = 0.0d0
 
79
      symops(3*iop+1,2) = 0.0d0
 
80
      symops(3*iop+2,2) = 1.0d0
 
81
      symops(3*iop+3,2) = 0.0d0
 
82
      symops(3*iop+1,3) = 0.0d0
 
83
      symops(3*iop+2,3) = 0.0d0
 
84
      symops(3*iop+3,3) = 1.0d0
 
85
      symops(3*iop+1,4) = 0.5d0
 
86
      symops(3*iop+2,4) = 0.5d0
 
87
      symops(3*iop+3,4) = 0.0d0
 
88
      end if
 
89
c 
 
90
c    F1 
 
91
c 
 
92
      if (numgrp.eq.234) then
 
93
      group_name = 'F1'
 
94
      nops = 3
 
95
      iop = 0
 
96
      symops(3*iop+1,1) = 1.0d0
 
97
      symops(3*iop+2,1) = 0.0d0
 
98
      symops(3*iop+3,1) = 0.0d0
 
99
      symops(3*iop+1,2) = 0.0d0
 
100
      symops(3*iop+2,2) = 1.0d0
 
101
      symops(3*iop+3,2) = 0.0d0
 
102
      symops(3*iop+1,3) = 0.0d0
 
103
      symops(3*iop+2,3) = 0.0d0
 
104
      symops(3*iop+3,3) = 1.0d0
 
105
      symops(3*iop+1,4) = 0.0d0
 
106
      symops(3*iop+2,4) = 0.5d0
 
107
      symops(3*iop+3,4) = 0.5d0
 
108
      iop = 1
 
109
      symops(3*iop+1,1) = 1.0d0
 
110
      symops(3*iop+2,1) = 0.0d0
 
111
      symops(3*iop+3,1) = 0.0d0
 
112
      symops(3*iop+1,2) = 0.0d0
 
113
      symops(3*iop+2,2) = 1.0d0
 
114
      symops(3*iop+3,2) = 0.0d0
 
115
      symops(3*iop+1,3) = 0.0d0
 
116
      symops(3*iop+2,3) = 0.0d0
 
117
      symops(3*iop+3,3) = 1.0d0
 
118
      symops(3*iop+1,4) = 0.5d0
 
119
      symops(3*iop+2,4) = 0.0d0
 
120
      symops(3*iop+3,4) = 0.5d0
 
121
      iop = 2
 
122
      symops(3*iop+1,1) = 1.0d0
 
123
      symops(3*iop+2,1) = 0.0d0
 
124
      symops(3*iop+3,1) = 0.0d0
 
125
      symops(3*iop+1,2) = 0.0d0
 
126
      symops(3*iop+2,2) = 1.0d0
 
127
      symops(3*iop+3,2) = 0.0d0
 
128
      symops(3*iop+1,3) = 0.0d0
 
129
      symops(3*iop+2,3) = 0.0d0
 
130
      symops(3*iop+3,3) = 1.0d0
 
131
      symops(3*iop+1,4) = 0.5d0
 
132
      symops(3*iop+2,4) = 0.5d0
 
133
      symops(3*iop+3,4) = 0.0d0
 
134
      end if
 
135
c 
 
136
c    I1 
 
137
c 
 
138
      if (numgrp.eq.235) then
 
139
      group_name = 'I1'
 
140
      nops = 1
 
141
      iop = 0
 
142
      symops(3*iop+1,1) = 1.0d0
 
143
      symops(3*iop+2,1) = 0.0d0
 
144
      symops(3*iop+3,1) = 0.0d0
 
145
      symops(3*iop+1,2) = 0.0d0
 
146
      symops(3*iop+2,2) = 1.0d0
 
147
      symops(3*iop+3,2) = 0.0d0
 
148
      symops(3*iop+1,3) = 0.0d0
 
149
      symops(3*iop+2,3) = 0.0d0
 
150
      symops(3*iop+3,3) = 1.0d0
 
151
      symops(3*iop+1,4) = 0.5d0
 
152
      symops(3*iop+2,4) = 0.5d0
 
153
      symops(3*iop+3,4) = 0.5d0
 
154
      end if
 
155
 
 
156
c 
 
157
c    A-1
 
158
c 
 
159
      if (numgrp.eq.236) then
 
160
      group_name = 'A-1'
 
161
      nops = 3
 
162
      iop = 0
 
163
      symops(3*iop+1,1) = -1.0d0
 
164
      symops(3*iop+2,1) = 0.0d0
 
165
      symops(3*iop+3,1) = 0.0d0
 
166
      symops(3*iop+1,2) = 0.0d0
 
167
      symops(3*iop+2,2) = -1.0d0
 
168
      symops(3*iop+3,2) = 0.0d0
 
169
      symops(3*iop+1,3) = 0.0d0
 
170
      symops(3*iop+2,3) = 0.0d0
 
171
      symops(3*iop+3,3) = -1.0d0
 
172
      symops(3*iop+1,4) = 0.0d0
 
173
      symops(3*iop+2,4) = 0.0d0
 
174
      symops(3*iop+3,4) = 0.0d0
 
175
      iop = 1
 
176
      symops(3*iop+1,1) = 1.0d0
 
177
      symops(3*iop+2,1) = 0.0d0
 
178
      symops(3*iop+3,1) = 0.0d0
 
179
      symops(3*iop+1,2) = 0.0d0
 
180
      symops(3*iop+2,2) = 1.0d0
 
181
      symops(3*iop+3,2) = 0.0d0
 
182
      symops(3*iop+1,3) = 0.0d0
 
183
      symops(3*iop+2,3) = 0.0d0
 
184
      symops(3*iop+3,3) = 1.0d0
 
185
      symops(3*iop+1,4) = 0.0d0
 
186
      symops(3*iop+2,4) = 0.5d0
 
187
      symops(3*iop+3,4) = 0.5d0
 
188
      iop = 2
 
189
      symops(3*iop+1,1) = -1.0d0
 
190
      symops(3*iop+2,1) = 0.0d0
 
191
      symops(3*iop+3,1) = 0.0d0
 
192
      symops(3*iop+1,2) = 0.0d0
 
193
      symops(3*iop+2,2) = -1.0d0
 
194
      symops(3*iop+3,2) = 0.0d0
 
195
      symops(3*iop+1,3) = 0.0d0
 
196
      symops(3*iop+2,3) = 0.0d0
 
197
      symops(3*iop+3,3) = -1.0d0
 
198
      symops(3*iop+1,4) = 0.0d0
 
199
      symops(3*iop+2,4) = 0.5d0
 
200
      symops(3*iop+3,4) = 0.5d0
 
201
      end if
 
202
c 
 
203
c    B-1
 
204
c 
 
205
      if (numgrp.eq.237) then
 
206
      group_name = 'B-1'
 
207
      nops = 3
 
208
      iop = 0
 
209
      symops(3*iop+1,1) = -1.0d0
 
210
      symops(3*iop+2,1) = 0.0d0
 
211
      symops(3*iop+3,1) = 0.0d0
 
212
      symops(3*iop+1,2) = 0.0d0
 
213
      symops(3*iop+2,2) = -1.0d0
 
214
      symops(3*iop+3,2) = 0.0d0
 
215
      symops(3*iop+1,3) = 0.0d0
 
216
      symops(3*iop+2,3) = 0.0d0
 
217
      symops(3*iop+3,3) = -1.0d0
 
218
      symops(3*iop+1,4) = 0.0d0
 
219
      symops(3*iop+2,4) = 0.0d0
 
220
      symops(3*iop+3,4) = 0.0d0
 
221
      iop = 1
 
222
      symops(3*iop+1,1) = 1.0d0
 
223
      symops(3*iop+2,1) = 0.0d0
 
224
      symops(3*iop+3,1) = 0.0d0
 
225
      symops(3*iop+1,2) = 0.0d0
 
226
      symops(3*iop+2,2) = 1.0d0
 
227
      symops(3*iop+3,2) = 0.0d0
 
228
      symops(3*iop+1,3) = 0.0d0
 
229
      symops(3*iop+2,3) = 0.0d0
 
230
      symops(3*iop+3,3) = 1.0d0
 
231
      symops(3*iop+1,4) = 0.5d0
 
232
      symops(3*iop+2,4) = 0.0d0
 
233
      symops(3*iop+3,4) = 0.5d0
 
234
      iop = 2
 
235
      symops(3*iop+1,1) = -1.0d0
 
236
      symops(3*iop+2,1) = 0.0d0
 
237
      symops(3*iop+3,1) = 0.0d0
 
238
      symops(3*iop+1,2) = 0.0d0
 
239
      symops(3*iop+2,2) = -1.0d0
 
240
      symops(3*iop+3,2) = 0.0d0
 
241
      symops(3*iop+1,3) = 0.0d0
 
242
      symops(3*iop+2,3) = 0.0d0
 
243
      symops(3*iop+3,3) = -1.0d0
 
244
      symops(3*iop+1,4) = 0.5d0
 
245
      symops(3*iop+2,4) = 0.0d0
 
246
      symops(3*iop+3,4) = 0.5d0
 
247
      end if
 
248
c 
 
249
c    C-1 
 
250
c 
 
251
      if (numgrp.eq.238) then
 
252
      group_name = 'C-1'
 
253
      nops = 3
 
254
      iop = 0
 
255
      symops(3*iop+1,1) = -1.0d0
 
256
      symops(3*iop+2,1) = 0.0d0
 
257
      symops(3*iop+3,1) = 0.0d0
 
258
      symops(3*iop+1,2) = 0.0d0
 
259
      symops(3*iop+2,2) = -1.0d0
 
260
      symops(3*iop+3,2) = 0.0d0
 
261
      symops(3*iop+1,3) = 0.0d0
 
262
      symops(3*iop+2,3) = 0.0d0
 
263
      symops(3*iop+3,3) = -1.0d0
 
264
      symops(3*iop+1,4) = 0.0d0
 
265
      symops(3*iop+2,4) = 0.0d0
 
266
      symops(3*iop+3,4) = 0.0d0
 
267
      iop = 1
 
268
      symops(3*iop+1,1) = 1.0d0
 
269
      symops(3*iop+2,1) = 0.0d0
 
270
      symops(3*iop+3,1) = 0.0d0
 
271
      symops(3*iop+1,2) = 0.0d0
 
272
      symops(3*iop+2,2) = 1.0d0
 
273
      symops(3*iop+3,2) = 0.0d0
 
274
      symops(3*iop+1,3) = 0.0d0
 
275
      symops(3*iop+2,3) = 0.0d0
 
276
      symops(3*iop+3,3) = 1.0d0
 
277
      symops(3*iop+1,4) = 0.5d0
 
278
      symops(3*iop+2,4) = 0.5d0
 
279
      symops(3*iop+3,4) = 0.0d0
 
280
      iop = 2
 
281
      symops(3*iop+1,1) = -1.0d0
 
282
      symops(3*iop+2,1) = 0.0d0
 
283
      symops(3*iop+3,1) = 0.0d0
 
284
      symops(3*iop+1,2) = 0.0d0
 
285
      symops(3*iop+2,2) = -1.0d0
 
286
      symops(3*iop+3,2) = 0.0d0
 
287
      symops(3*iop+1,3) = 0.0d0
 
288
      symops(3*iop+2,3) = 0.0d0
 
289
      symops(3*iop+3,3) = -1.0d0
 
290
      symops(3*iop+1,4) = 0.5d0
 
291
      symops(3*iop+2,4) = 0.5d0
 
292
      symops(3*iop+3,4) = 0.0d0
 
293
      end if
 
294
c 
 
295
c    F-1 
 
296
c 
 
297
      if (numgrp.eq.239) then
 
298
      group_name = 'F-1'
 
299
      nops = 7
 
300
      iop = 0
 
301
      symops(3*iop+1,1) = -1.0d0
 
302
      symops(3*iop+2,1) = 0.0d0
 
303
      symops(3*iop+3,1) = 0.0d0
 
304
      symops(3*iop+1,2) = 0.0d0
 
305
      symops(3*iop+2,2) = -1.0d0
 
306
      symops(3*iop+3,2) = 0.0d0
 
307
      symops(3*iop+1,3) = 0.0d0
 
308
      symops(3*iop+2,3) = 0.0d0
 
309
      symops(3*iop+3,3) = -1.0d0
 
310
      symops(3*iop+1,4) = 0.0d0
 
311
      symops(3*iop+2,4) = 0.0d0
 
312
      symops(3*iop+3,4) = 0.0d0
 
313
      iop = 1
 
314
      symops(3*iop+1,1) = 1.0d0
 
315
      symops(3*iop+2,1) = 0.0d0
 
316
      symops(3*iop+3,1) = 0.0d0
 
317
      symops(3*iop+1,2) = 0.0d0
 
318
      symops(3*iop+2,2) = 1.0d0
 
319
      symops(3*iop+3,2) = 0.0d0
 
320
      symops(3*iop+1,3) = 0.0d0
 
321
      symops(3*iop+2,3) = 0.0d0
 
322
      symops(3*iop+3,3) = 1.0d0
 
323
      symops(3*iop+1,4) = 0.0d0
 
324
      symops(3*iop+2,4) = 0.5d0
 
325
      symops(3*iop+3,4) = 0.5d0
 
326
      iop = 2
 
327
      symops(3*iop+1,1) = -1.0d0
 
328
      symops(3*iop+2,1) = 0.0d0
 
329
      symops(3*iop+3,1) = 0.0d0
 
330
      symops(3*iop+1,2) = 0.0d0
 
331
      symops(3*iop+2,2) = -1.0d0
 
332
      symops(3*iop+3,2) = 0.0d0
 
333
      symops(3*iop+1,3) = 0.0d0
 
334
      symops(3*iop+2,3) = 0.0d0
 
335
      symops(3*iop+3,3) = -1.0d0
 
336
      symops(3*iop+1,4) = 0.0d0
 
337
      symops(3*iop+2,4) = 0.5d0
 
338
      symops(3*iop+3,4) = 0.5d0
 
339
      iop = 3
 
340
      symops(3*iop+1,1) = 1.0d0
 
341
      symops(3*iop+2,1) = 0.0d0
 
342
      symops(3*iop+3,1) = 0.0d0
 
343
      symops(3*iop+1,2) = 0.0d0
 
344
      symops(3*iop+2,2) = 1.0d0
 
345
      symops(3*iop+3,2) = 0.0d0
 
346
      symops(3*iop+1,3) = 0.0d0
 
347
      symops(3*iop+2,3) = 0.0d0
 
348
      symops(3*iop+3,3) = 1.0d0
 
349
      symops(3*iop+1,4) = 0.5d0
 
350
      symops(3*iop+2,4) = 0.0d0
 
351
      symops(3*iop+3,4) = 0.5d0
 
352
      iop = 4
 
353
      symops(3*iop+1,1) = -1.0d0
 
354
      symops(3*iop+2,1) = 0.0d0
 
355
      symops(3*iop+3,1) = 0.0d0
 
356
      symops(3*iop+1,2) = 0.0d0
 
357
      symops(3*iop+2,2) = -1.0d0
 
358
      symops(3*iop+3,2) = 0.0d0
 
359
      symops(3*iop+1,3) = 0.0d0
 
360
      symops(3*iop+2,3) = 0.0d0
 
361
      symops(3*iop+3,3) = -1.0d0
 
362
      symops(3*iop+1,4) = 0.5d0
 
363
      symops(3*iop+2,4) = 0.0d0
 
364
      symops(3*iop+3,4) = 0.5d0
 
365
      iop = 5
 
366
      symops(3*iop+1,1) = 1.0d0
 
367
      symops(3*iop+2,1) = 0.0d0
 
368
      symops(3*iop+3,1) = 0.0d0
 
369
      symops(3*iop+1,2) = 0.0d0
 
370
      symops(3*iop+2,2) = 1.0d0
 
371
      symops(3*iop+3,2) = 0.0d0
 
372
      symops(3*iop+1,3) = 0.0d0
 
373
      symops(3*iop+2,3) = 0.0d0
 
374
      symops(3*iop+3,3) = 1.0d0
 
375
      symops(3*iop+1,4) = 0.5d0
 
376
      symops(3*iop+2,4) = 0.5d0
 
377
      symops(3*iop+3,4) = 0.0d0
 
378
      iop = 6
 
379
      symops(3*iop+1,1) = -1.0d0
 
380
      symops(3*iop+2,1) = 0.0d0
 
381
      symops(3*iop+3,1) = 0.0d0
 
382
      symops(3*iop+1,2) = 0.0d0
 
383
      symops(3*iop+2,2) = -1.0d0
 
384
      symops(3*iop+3,2) = 0.0d0
 
385
      symops(3*iop+1,3) = 0.0d0
 
386
      symops(3*iop+2,3) = 0.0d0
 
387
      symops(3*iop+3,3) = -1.0d0
 
388
      symops(3*iop+1,4) = 0.5d0
 
389
      symops(3*iop+2,4) = 0.5d0
 
390
      symops(3*iop+3,4) = 0.0d0
 
391
      end if
 
392
c 
 
393
c    I-1 
 
394
c 
 
395
      if (numgrp.eq.240) then
 
396
      group_name = 'I-1'
 
397
      nops = 3
 
398
      iop = 0
 
399
      symops(3*iop+1,1) = -1.0d0
 
400
      symops(3*iop+2,1) = 0.0d0
 
401
      symops(3*iop+3,1) = 0.0d0
 
402
      symops(3*iop+1,2) = 0.0d0
 
403
      symops(3*iop+2,2) = -1.0d0
 
404
      symops(3*iop+3,2) = 0.0d0
 
405
      symops(3*iop+1,3) = 0.0d0
 
406
      symops(3*iop+2,3) = 0.0d0
 
407
      symops(3*iop+3,3) = -1.0d0
 
408
      symops(3*iop+1,4) = 0.0d0
 
409
      symops(3*iop+2,4) = 0.0d0
 
410
      symops(3*iop+3,4) = 0.0d0
 
411
      iop = 1
 
412
      symops(3*iop+1,1) = 1.0d0
 
413
      symops(3*iop+2,1) = 0.0d0
 
414
      symops(3*iop+3,1) = 0.0d0
 
415
      symops(3*iop+1,2) = 0.0d0
 
416
      symops(3*iop+2,2) = 1.0d0
 
417
      symops(3*iop+3,2) = 0.0d0
 
418
      symops(3*iop+1,3) = 0.0d0
 
419
      symops(3*iop+2,3) = 0.0d0
 
420
      symops(3*iop+3,3) = 1.0d0
 
421
      symops(3*iop+1,4) = 0.5d0
 
422
      symops(3*iop+2,4) = 0.5d0
 
423
      symops(3*iop+3,4) = 0.5d0
 
424
      iop = 2
 
425
      symops(3*iop+1,1) = -1.0d0
 
426
      symops(3*iop+2,1) = 0.0d0
 
427
      symops(3*iop+3,1) = 0.0d0
 
428
      symops(3*iop+1,2) = 0.0d0
 
429
      symops(3*iop+2,2) = -1.0d0
 
430
      symops(3*iop+3,2) = 0.0d0
 
431
      symops(3*iop+1,3) = 0.0d0
 
432
      symops(3*iop+2,3) = 0.0d0
 
433
      symops(3*iop+3,3) = -1.0d0
 
434
      symops(3*iop+1,4) = 0.5d0
 
435
      symops(3*iop+2,4) = 0.5d0
 
436
      symops(3*iop+3,4) = 0.5d0
 
437
      end if
 
438
 
 
439
 
 
440
      return
 
441
      end
 
442