forked from spglib/spglib
-
Notifications
You must be signed in to change notification settings - Fork 0
/
ChangeLog
866 lines (571 loc) · 27.5 KB
/
ChangeLog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
2014-02-17 Atsushi Togo <atz.togo@gmail.com>
* Version 1.6.0:
* src/spglib.c (spg_get_spacegroup_type):
A function is made for accessing space-group-type database.
2014-01-07 Atsushi Togo <atz.togo@gmail.com>
* src/tetrahedron_method.*:
Tetrahedron method is implemented. Example is found in the example
directory.
2013-09-10 Atsushi Togo <atz.togo@gmail.com>
* Version 1.5.2:
* kpoint.c (kpt_relocate_BZ_grid_address):
New function kpt_relocate_BZ_grid_address is implemented in
kpoint.c. k-points on the mesh of the reciprocal primitive lattice
are relocated to fit in the Brillouin zone. To treat
translationally equivalent k-points, the memory space of
bz_grid_address[prod(mesh + 1)][3] is required. bz_map[prod(mesh
* 2 - 1)] is used to map k-points on the expanded k-space where
the shape is same as the reciprocal primitive lattice but the size
along each lattice vector is multiplied by mesh * 2 - 1. BZ zone
is included in this expanded k-space, so the bz_map can hold the
location of relocated k-points. Irreducible k-points are stored
from 0 to index of prod(mesh) - 1 in bz_grid_address in the same
order of grid_address, and the translationally equivalent k-points
are stored after that. Fast indexing of the location of relocated
k-points can be executed using bz_map with similar way to usual
grid_address scheme.
* kpoint.c (kpt_get_BZ_triplets_at_q):
New function kpt_get_BZ_triplets_at_q is implemented using
relocate_BZ_grid_address. This generates a set of triplets that
cleanly fit in BZ.
2013-09-10 Atsushi Togo <atz.togo@gmail.com>
* kpoint.c:
Remove get_triplets_reciprocal_mesh and
extract_triplets_reciprocal_mesh_at_q. Add
get_grid_triplets_at_q. This function treats q-points on boundary
correctly.
* spglib.c (spg_get_symmetry):
Previously symmetry operations were obtained by sym_get_operation.
Now it is changed so that they are extracted from spglib_dataset.
2013-04-16 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c (get_operation):
In some case, get_space_group_operation returns more than 48
symmetry operations for a primitive cell. To avoid this, an
iterative approach, which is just reducing tolerance gradually to
achieve number of symmetry operations <= 48, is implemented.
* symmetry.c:
In the symmetry operation search, rot[] and trans[] were used. Now
they are replaced by Symmetry *symmetry in symmetry.c.
2013-04-14 Atsushi Togo <atz.togo@gmail.com>
* Fortran interface given by Dimitar Pashov:
The new fortran interface is given by Dimitar Pashov, which is
much safer than the previous one. The new fortran interface, its
example, and makefile are found in the example directory as
spglib_f08.f90, example_f08.f90 and Makefile_f08.
* spglib.h, spglib.c, kpoint.h, kpoint,c:
Symprec parameter was removed from spg_get_stabilized_reciprocal_mesh.
2013-01-28 Atsushi Togo <atz.togo@gmail.com>
* spglib.c, spglib_f.c:
Pointerd out by Michael Rutter that the length of charasters of
space-group type symbol (International table) was wrong in
spglib_f.c. That of Schoenflies symbol was also found wrong. The
former and latter were fixed as from 21 to 11 and from 10 to 7,
respectively.
2012-10-10 Atsushi Togo <atz.togo@gmail.com>
* Version 1.3
Rethink algorithms to search symmetry operations for speed-up and
stability. Many small modifications are applied. Some parts of the
codes that look unnecessary were removed but this change might
make the code unstable.
* symmetry.c (get_translation):
Improve OpenMP efficiency. Number of atoms to enable OpenMP
concurrency is set as 1000 (NUM_ATOMS_CRITERION_FOR_OPENMP).
2012-10-09 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c (sym_get_pure_translation, get_translation):
In 'sym_get_pure_translation', iterative approach to search
pure_translation was employed. But 'sym_get_pure_translation' was
only used from two functions in primitive.c. In these two
functions, tolerance is iteratively reduced when the search
fails. Because there is no need to do iteration twice, the
iteration part in 'sym_get_pure_translation' was removed.
* primitive.{c,h}:
'get_primitive_with_pure_translation' was replaced by
'get_primitive_and_pure_translation'. The former was only used
from 'get_operation' in symmetry.c and the interative approach was
implemented in 'get_operation'. The iterative part was moved into
'get_primitive_and_pure_translation'.
The Primitive structure was defined in primitive.h. This just
packs Cell * primitive_cell and VecDBL * pure_trans. The Primitive
structure is only used in 'get_operation' in symmetry.c and
'(prm_)get_primitive_and_pure_translation' in primitive.c.
2012-10-08 Atsushi Togo <atz.togo@gmail.com>
* spacegroup.c (get_hall_number_local_iteration):
REDUCE_RATE to reduce number of symmetry operations is changed
from 0.2 to 0.95. Symmetry operations depend on tolerance
delicately.
* Remove many lines of debug_print:
Now to follow the flow is more important than to see the
detailed behavior of implementation.
2012-10-03 Atsushi Togo <atz.togo@gmail.com>
* Version 1.2.5
* primitive.c:
Separate out get_primitive_with_mapping_table from
prm_get_primitive_with_mapping_table. Rename
prm_get_primitive_with_all to
get_primitive_with_pure_translation.
In get_primitive_with_pure_translation, it was
if (pure_trans->size > 0)
but not it was fixed to
if (pure_trans->size > 1)
not to search primitive cell when pure_trans->size == 1.
2012-10-01 Atsushi Togo <atz.togo@gmail.com>
* Version 1.2.4
* kpoint.c (get_ir_reciprocal_mesh):
In spglib, k-point grid is indexed as, for example,
0 1 2 3 4 5 6 -5 -4 -3 -2 -1
for 12 points mesh. Now the different index style such as
0 1 2 3 4 5 -6 -5 -4 -3 -2 -1
is activated by CFLAGS='-DGRID_BOUNDARY_AS_NEGATIVE' before
running configure script at the compilation time.
* symmetry.c (get_operation):
In the iterative reduction of tolerance of primitive cell search,
pure_trans was searched by a fixed tolerance rather than the
reduced one. This is fixed to use the reduced tolerance at each
iteration.
2012-08-13 Atsushi Togo <atz.togo@gmail.com>
* Version 1.2.3
* symmetry.c (is_identity_metric):
To check identity of metric tensors, too small angle induces
numerical size error, so if the dtheta is calculated smaller than
1e-6, automatically set the angles be identical.
* spacegroup.c (get_hall_number_local_iteration):
Implementation of iterative reduction of symmetry tolerance to
search symmetry operation to get Hall symbol is separated as the
function get_hall_number_local_iteration and the tolerance
reduction rate of REDUCE_RATE is changed from 0.2 from 0.95.
2012-07-17 Atsushi Togo <atz.togo@gmail.com>
* Version 1.2.2
* lattice.c (get_centering):
Tranformation A-centring monoclinic lattice to C-centring
monoclinic lattice preserving b-axis is not well handled and now
the matrix 'monocli_a2c' is used for it.
2012-07-14 Atsushi Togo <atz.togo@gmail.com>
* Version 1.2.1
* pointgroup.c (laue_one_axis):
Transformation matrix is expected to be always positive, but it
was not correctly implemented.
2012-05-28 Atsushi Togo <atz.togo@gmail.com>
* refinement.c (reduce_symmetry_in_frame):
Fix memory leak by freeing lattice_trans.
2012-04-26 Atsushi Togo <atz.togo@gmail.com>
The numbering of the hall symbol is included in the 'struct
SpglibDataset'. This starts from 1 to 530.
2012-04-25 Atsushi Togo <atz.togo@gmail.com>
* spglib.c, spglib.h:
To find symmetry with angle tolerance, for spg_* functions,
respective spgat_* functions are prepared. In the spgat_*
functions, the last variable is the angle tolerance (in degrees).
* symmetry.c:
The way to find lattice symmetry is changed. Now a new static
variable, 'angle_tolerance', can be used for the tolerance of
angle. Therefore the non-diagonal elements of metric tensor are
directry compared with this angle. 'angle_tolerance' can be
accessed via 'sym_set_angle_tolerance' inside spglib. However the
default behavior does not use the angle tolerance. In the default
behavior, averaged_lattice_vector_length * sin(delta_theta) is
compared with usual 'symprec'. The default behavior is invoked
when 'angle_tolerance' has negative value. These are implemented
in 'is_identity_metric' in symmetry.c.
2012-04-12 Atsushi Togo <atz.togo@gmail.com>
* lattice.c (get_centering):
Correction matrix for I-centring monoclinic was implemented.
2012-03-03 Atsushi Togo <atz.togo@gmail.com>
* primitive.c:
The functions of 'prm_get_current_tolerance',
'set_current_tolerance' and the value of static double
current_tolerance were introduced in 'primitive.c'. These are used
to return reduced tolerance value that is iteratively determined
in 'prm_get_primitive_with_mapping_table'.
The interative technique is used in several parts to obtain
possible primitive cell parameters by reducing symmetry
tolerance. Only in 'get_overlap_table' in 'primitive.c', the
interative techinique is used with reducing and increasing the
tolerance, i.e., in the other parts, only reducing.
2012-03-02 Atsushi Togo <atz.togo@gmail.com>
* Version 1.1.2
* refinement.c, primitive.c, symmetry.c:
Modifications at many small parts of the code for improve
robustness with large symmetry tolerance values.
2012-02-28 Atsushi Togo <atz.togo@gmail.com>
* primitive.c (get_overlap_table):
Simplify the algorithm for search overlapping table. The new
algorithm is as follows: Symmetry tolerance monotonically
decreases by REDUCE_RATE = 0.95 when the search is failed. When
number of overlapping atoms is found smaller than the ratio of
number of atoms between original and primitive cells, the symmetry
tolerance is doubled. For this algorithm, 'check_overlap_table'
function was implemented.
2012-02-09 Atsushi Togo <atz.togo@gmail.com>
* spin.{c,h}:
'spin.{c,h}' are added. Colliear spins on atoms are employed as
additional constraints to search symmetry. The API is
'spg_get_symmetry_with_collinear_spin'. (reduce_symmetry_in_frame):
Bug fix to get refined symmetry operations. When the original cell
is not primitive and the lattice symmetry is lower than that
primitive lattice, the symmetry of primitive has to be
reduced to those belongs to the lower symmetry. But this was not
correctly done.
2011-10-05 Atsushi Togo <atz.togo@gmail.com>
* spglib.c (spg_get_pointgroup):
'spg_get_pointgroup' is implemented. This is used to obtain point
group symbols.
2011-09-30 Atsushi Togo <atz.togo@gmail.com>
* site_symmetry.c (get_Wyckoff_notation):
Bug fix of assignment of Wyckoff letters
'sitesym_database.c' is compressed. Now only the first element of
'coordinates' is used to aasign a Wyckoff letter.
2011-08-09 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c (sym_get_pure_translation):
Iterative reduction of symmetry tolerance to find pure
translations is implemented as well as different parts.
2011-07-29 Atsushi Togo <atz.togo@gmail.com>
* kpoint.c, kpoint.h:
symmetry_kpoint.c and symmetry_kpoint.h are renamed as kpoint.c
and kpoint.h, respectively.
2011-07-06 Atsushi Togo <atz.togo@gmail.com>
* site_symmetry.c:
A set of equivalent atoms is returned. It is accessible
through spg_get_dataset. (spglib-1.0.7)
2011-07-05 Atsushi Togo <atz.togo@gmail.com>
* site_symmetry.c, sitesym_database.c, spglib-1.0.6:
Newly added two source files (and also *.h). The part related to
site symmetry in refinement.c is moved to site_symmetry.c. Wyckoff
position determination is implemented in site_symmetry.c and the
database is in sitesym_database.c. The pointers in the structure
SpglibDataset are made free by calling spg_free_dataset. It was
difficult to implement Wyckoff position determination without
Yusuke Seko's formated crystallographic database. I appreciate his
kind offer. (spglib-1.0.6)
* primitive.c (prm_get_primitive):
The function is modified to return mapping table from input cell
to primitive cell as the first argument.
2011-06-18 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c:
OpenMP is included to search translations with each rotation.
Compilation with gcc is achived by '-lgomp' and '-fopenmp'.
* refinement.c:
Performance tuning was made for 'get_refined_symmetry_operations'.
2011-06-18 Atsushi Togo <atz.togo@gmail.com>
* refinement.c:
In 'ref_get_refined_symmetry_operations', it was forgotten to
consider origin shift to obtain symmetry operations from
database. This is fixed by modifying translations in
'set_translation_with_origin_shift'.
2011-06-15 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.4:
The new version of spglib is released. To the next version,
implementation of Wyckoff positions is planed.
2011-06-15 Atsushi Togo <atz.togo@gmail.com>
* pyspglib:
pyspglib didn't work correctly. It was fixed.
2011-06-15 Atsushi Togo <atz.togo@gmail.com>
* spglib.c:
A new function 'spg_get_dataset is implmented. This returns a
'SpglibDataset' type structure that is defined in 'spglib.h'. This
includes refined informations of crystal symmetry. For example,
the symmetry operations are derived based on the operations in
database and so they are different from those obtained from
'spg_get_symmetry'. Wyckoff positions are to be impmemented.
2011-06-13 Atsushi Togo <atz.togo@gmail.com>
* spg_database.c:
A static array of 'symmetry_operations' is compressed because it
was very large. It contains symmetry operations of 530 hall
symbols. Each space group operation requires a rotation 3x3 matrix
and a translation vector, therefore 12 elements were implemented
for that. Now it is compressed using ternary numerical system for
rotation and duodecimal system for translation. This is achieved because
each element of rotation matrix can have only one of {-1,0,1}, and
the translation can have one of {0,2,3,4,6,8,9,10} divided by
12. Therefore 3^9 * 12^3 = 34012224 different values can map space
group operations. In principle, octal numerical system can be used
for translation, but duodecimal system is used for the convenient.
The fact that determinant of rotation matrix is -1 or 1 could be
used as a constraint to reduce the number of mapping space. These
may be implemented if more compression is required.
2011-06-10 Atsushi Togo <atz.togo@gmail.com>
* spglib.c:
spg_show_symmetry was removed.
2011-05-17 Atsushi Togo <atz.togo@gmail.com>
* primitive.c, etc
A tolerance to find primitive lattice vectors sometimes fails
because found pure translations do not give possible lattice
vectors in which the volume ratio of the input cell and primitive
cell does not match to the number of pure translations. So
iterative tolerance reduction is implemented. The reduced
tolerance should be used to find spacegroup and refined cell
because to hold consistency with symmetry operations to be found
for primitive cell. The reduced tolerance is not used inside
symmetry opration search. The tolerance used to find primitive
cell is stored in a static variable 'tolerance' in primitive.c and
this can be accessed via prm_get_tolerance().
Old 'prm_get_primitive' was removed and new one is that of
'prm_get_primitive_with_pure_trans' renamed. So now pure
translations are always required to call to obtain primitive cell.
When calling 'prm_get_primitive' and if the input cell is already
primitive, primitive->size=0 was returned in the previous
implementation, but in the new implementation, primitive cell that
has smallest lattice vectors is returned.
2011-04-21 Atsushi Togo <atz.togo@gmail.com>
* refinement.c (refine_cell):
When an input structure is a primitive cell but the lattice is not
smallest, 'refine_cell' didn't work correctly. This is fixed by
using new function of 'prm_get_cell_with_smallest_lattice' in
primitive.c. This function was originally in spacegroup.c as
'get_cell_with_smallest_lattice'. Following this, spacegroup.c,
primitive.c, and primitive.h were modified.
2011-03-18 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.3.1:
A critical bug in the structure refinment in refinement.c was
removed. A tolerance parameter in hall_symbol.c came back to
is_hall_symbol.
2011-03-15 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.3:
Crystal structure refinement by spg_refine_cell is
implemented. Code around tolerance is tunded up, hopefully
becoming more robust with large tolerance.
* refinement.c:
bravais.{c,h} are renamed to refinement.{c,h}.
2011-03-14 Atsushi Togo <atz.togo@gmail.com>
* debug.h (SPGWARNING):
Debug option of SPGWARNING has been defined. This may be used for
beta-tests.
* bravais.c, spacegrouop.c, hall_symbol.c, lattice.c, primitive.c:
Tolerance parameter has been reconsidered to make finding
spacegroup robust with a loose tolerance value. Tolerance is used
for Cartesian distance, norm squared, or volume. It wold be
avoided for ratio, reduced position or distance.
* bravais.c:
Exact locations of atoms are obtained by get_exact_location in
bravais.c.
* sapcegroup.c:
Iterative search for hall symbol is implemented in
get_hall_number. The tolerance is sequentially reduced until
finding a hall symbol. With the reduced tolerance, symmetry
operations are researched by sym_reduce_operation. The cost of the
reduction of symmetry operations is expeted much smaller than that
of the increase of symmetry operations with increasing tolerance.
* hall_symbol.c:
All tolerance parameters were kicked out.
* lattice.c:
get_base_center is modified to remove the tolerance paremeter. The
centring check is now based on
integers. lat_smallest_lattice_vector returns 1 or 0 to show if
the lattice has non-zero volume or not.
* primitive.c:
Primitive lattice is made of pure translations. Previously the
primitive lattice vectors are constructed from pure
translations. Currently, it is set such that the primitive lattice
parameters multiplied with 'intergers' become the original
latttice parameters.
* mathfunc.c (mat_is_int_matrix):
A macro of INT_PREC (currently set 0.1 ) is used for checking a
value is integer or not in mat_is_int_matrix.
2011-02-18 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.2:
Bug fix in finding lattice point symmetry (symmetry.c).
spg_get_bravais is added and spg_get_international_with_bravais
was removed.
2010-12-26 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.1:
By the help of Pawel T. Jochym, the setup.py in the ASE-interface,
pysgpblib, has been updated for matching the latest spglib.
2010-11-08 Atsushi Togo <atz.togo@gmail.com>
* spglib-1.0.0:
A lot of source code was re-written for this one month. This
release contains the algorithms of symmetry opertation finding and
space-group type determination which are different from the
versions before spglib 0.7.4.
2010-11-08 Atsushi Togo <atz.togo@gmail.com>
Algorithm to search basis vectors is changed to that of
[Grosse-Kunstleve, Acta Cryst. A55, 383-395 (1999)], which is
implemented in lattice.c. bravais.{c,h}, bravais_virtual.{c,h} are
removed. spacegroup_type.{c,h} are moved to spacegroup.{c,h}.
2010-11-05 Atsushi Togo <atz.togo@gmail.com>
* lattice.c:
lattice.c and lattice.h are added. Delaunay reduction is moved
there from bravais.c.
2010-11-02 Atsushi Togo <atz.togo@gmail.com>
* spglib 1.0-beta:
New release under the modified BSD license.
* pointgroup.c:
Sometimes symmetry operations that don't match point-group type
can be found. In the current implementation, some of non-necessary
operations are abandoned. For this, get_cast_symmetry,
check_pointgroup_condition, ptg_get_symmetry_pointgroup are added.
However still this implementation is not perfect. Maybe another
strategy is required.
2010-10-27 Atsushi Togo <atz.togo@gmail.com>
* primitive.c:
trim_cell is divided into trim_cell, set_primitive_positions,
get_positions_primitive, get_overlap_table. is_overlap is
modified to check Cartesian distance, and moved to cel_is_overlap
in cell.c.
* symmetry.c (get_translation):
Atomic position overlap check part is modified to use
cel_is_overlap.
* cell.c:
cel_frac_to_cart and cel_cart_to_frac are removed.
2010-10-25 Atsushi Togo <atz.togo@gmail.com>
* bravais_arc.* is renamed to bravais_virtual.*
The function is_holohedry is removed.
* bravais.c:
brv_set_relative_lattice is added. brv_get_brv_lattice_in_loop is
modified for handling RHOMB case flexibly. get_metric was removed
and moved to mathfunc.c as mat_get_metric.
* pointgroup.c:
The first argument of ptg_get_holohedry is removed, which was used
just for debugging.
* symmetry.c:
get_space_group_operation is separated from
get_operation. sym_get_operation_direct, set_axes and
get_lattice_symmetry are newly added. sym_get_operation_direct is
the function for finding symmetry operations without relying on
Bravais lattice. The lattice symmetry is exhaustively searched
through the matrices whose elements are -1, 0, or 1.
* mathfunc.c:
mat_get_metric and mat_check_identity_matrix_d3 are added.
2010-10-22 Atsushi Togo <atz.togo@gmail.com>
* bravais.c:
The names of local functions starting with brv_... is changed to
get_brv_.... Add functions get_brv_cubic_primitive,
get_brv_tetra_primitive, get_brv_ortho_primitive to find more
variety of lattice transformation matrices.
2010-10-22 Atsushi Togo <atz.togo@gmail.com>
The space group type finder is replaced completely. spacegroup.*,
spacegroup_data.*, spacegroup_database.* are
removed. spacegroup_type.*, spg_database.*, and hall_symbol.* are
added. The helper tools to create hall symbol related database are
in database directories.
2010-10-09 Atsushi Togo <atz.togo@gmail.com>
spglib-0.7.4 is released.
cel_new_cell is renamed to cel_alloc_cell, and cel_alloc_cell
returns a poiter of the Cell structure. cel_delete_cell is renamed
to cel_free_cell.
sym_new_symmetry is renamed to sym_alloc_symmetry, and
sym_alloc_symmetry returns a poiter of the Symmetry
structure. sym_delete_symmetry is renamed to sym_free_symmetry.
The other functions returns Cell and Symmetry structures are
modified to return those pointers.
2010-10-05 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c (get_index_with_least_atoms):
To search translation, it is efficient if an atom in least atoms
with an atomic type is used for the search origine. This function
looks for it. This is expected to accelerate especially when the
cell is huge and the lattice has high symmetry, but the number of
atoms of some atomic type is few.
2010-10-04 Atsushi Togo <atz.togo@gmail.com>
* primitive.c (prm_get_primitive_with_pure_trans):
Create new function 'prm_get_primitive_with_pure_trans' to avoid
searching same pure translations many times, which leads some
speed up.
2010-09-26 Atsushi Togo <atz.togo@gmail.com>
'const' of multidimensional arrays are removed. Just to notify it
should be 'const', 'SPGCONST' replaces those 'const', however it
is defined as nothing in mathfunc.h. This is done because of
avoiding lots of waring in C-compiler of 'gcc: incompatible
pointer type'.
* spglib_f.c:
spglib_f.c is updated and some examples are found in example.f90.
2010-04-12 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c, mathfunc.c:
To make the translations clean, sym_get_fractional_translation in
symmetry.c and mat_Dmod1 in mathfunc.c are
added. sym_get_fractional_translation is supposed to look for the
translation that is most closest to one of 0, 1/6,
1/4, 1/3, 1/2, 2/3, 3/4, and 5/6. The first implementation
contained a bug and Jorg Meyer solved this.
2010-03-25 Atsushi Togo <atz.togo@gmail.com>
* bravais_art.c:
Space group of the virtual structures in the test directory
failed. Fixed points are as follows:
P-Cubic, P-Tetra, P-Ortho to Base center monoclinic, which was
collect in the old code, and I-Cubic to F-Orthorhombic.
Cubic to Rhombohedral. To obtain correct principal axes,
'is_rhombo' was implemented.
2010-03-22 Atsushi Togo <atz.togo@gmail.com>
* spglib.h, spglib.c:
spg_get_bravais_lattice is removed.
spg_get_conventional_symmetry is implemented. This is used to get
symmetry operations for the conventional unit cell. The
conventional unit cell is searched from the primitive cell found
from the input cell, i.e., even if the input cell is not
primitive, the primitive cell is automatically searched and the
conventional cell is determined. The Bravais lattice of the
conventinal cell is returned as bravais_lattice. This
bravais_lattice is found including the internal atomic
positions. Therefore artificial cell is expected to be handled
correctly.
2010-02-11 Atsushi Togo <atz.togo@gmail.com>
* primitive.c:
There was a problem that primitive cell is not well found in naive
cases (inconsistency between tolerances in finding pure
translations and in finding overlapping atoms). David Lonie
implemented an iterative way to adjust symmetry finder tolerance
on the fly.
2009-07-15 Atsushi Togo <atz.togo@gmail.com>
* bravais.c:
Bug to find bravais lattice in FCC and rhombohedral, when the cell
is a supercell. This is probably ad-hoc fix.
2009-05-19 Atsushi Togo <atz.togo@gmail.com>
* symmetry_kpoint.c:
Time reversal symmetry is included.
* bravais.c (brv_smallest_lattice_vector):
Delaunay reduction is implemented instead of ABINIT
implementation.
2009-03-09 Atsushi Togo <atz.togo@gmail.com>
* symmetry.c:
Critical bug for finding symmetry for primitive cell which is
reduced from convensional cell.
2009-03-05 Atsushi Togo <atz.togo@gmail.com>
* symmetry_kpoint.c:
Functions to look for the irreducible points in reciprocal space
are implemented.
2009-02-25 Atsushi Togo <atz.togo@gmail.com>
* bravais_art.c:
Hexagonal to Orthorhombic symmetry lowering case is implemented.
* python module:
Python module is implemented. To use them, numpy is
required. Interface for ASE is implemented.
2008-08-22 Atsushi Togo <atz.togo@gmail.com>
* bravais_art.c:
This file is used for the cases that the lattice has higher
symmetry than the symmetry given by atomic positions. This does
not affect to the symmetry operations that is firstly obtained by
assuming the bravais lattice having high symmetry. Thus this file
works only for changing its bravais lattice. The information about
bravais lattice is used for only determing space group. Even if
the code in this file fails, or let's say, space group
determination fails), we can not say the symmetry operations are
wrong.
2008-07-11 Atsushi Togo <atz.togo@gmail.com>
* Most of the code
The procedure for searching a Bravais lattice from a primitive
cell is rewritten. The algorithm is also totally changed.
* test directory
For the test of space group search, 'test' directory is
created. Many POSCAR's are prepared in this directory. For running
tests, a small ruby code is written and put into this directory.
2008-05-28 Atsushi Togo <atz.togo@gmail.com>
* spacegroup_data.c (get_spacegroup_data_special_case):
Bug fix for space group determinations for special cases (Pcc2,
Pbam, Pccn).
2008-05-25 Atsushi Togo <atz.togo@gmail.com>
* *.c, *.h :
Many codes are re-built. Many functions are changed to
"static". The first three characters of public functions are
re-named using those file names.
* bravais.c :
Monoclinic is possibly improved. Monoclinic is difficult. It may
contain further bugs...
2008-05-01 Atsushi Togo <atz.togo@gmail.com>
* spacegroup.c (get_proper_rotation_class_6axis):
exit(1) is inserted.
* bravais.c/h:
"lattice_symmetry" is removed from the "Bravais" structure.
* spglib.h:
Misspelled function prototypes of spg_get_symmetry and
spg_check_symmetry are corrected.
2008-04-23 Atsushi Togo <atz.togo@gmail.com>
* Many patches by Xavier Andrade
For fixing many warning in many codes and gcc-4.3 related problem
in bravais.c.