Comments (5)
I observe this issue for the fourth-order symmetric tensor for hexagonal (6/mmm).
Lines 113 to 122 in a928f36
This issue seems to happen when input representation matrices have noise. The noise is amplified for higher ranks. Currently, I have no idea to properly handle it.
from spgrep.
I tested the above example and found that len(sym_tensors)
is 19, so the test fails.
But to be honest, I am still not very clear about the physical (mathematical) meaning of the following test steps performed in this test:
rep = get_representation_on_symmetric_matrix(rotations)
tensors = get_symmetry_adapted_tensors(rep, rotations, rank=4, real=True)
sym_tensors = apply_intrinsic_symmetry(tensors)
assert len(sym_tensors) == 18
As a user of GAP, I try to check the representation of these rotations, aka, the point group, as follows:
gap> rotations:=[[[ 1, 0, 0],
> [ 0, 1, 0],
> [ 0, 0, 1]],
>
> [[-1, 0, 0],
> [ 0, -1, 0],
> [ 0, 0, -1]],
>
> [[ 1, -1, 0],
> [ 1, 0, 0],
> [ 0, 0, 1]],
>
> [[-1, 1, 0],
> [-1, 0, 0],
> [ 0, 0, -1]],
>
> [[ 0, -1, 0],
> [ 1, -1, 0],
> [ 0, 0, 1]],
>
> [[ 0, 1, 0],
> [-1, 1, 0],
> [ 0, 0, -1]],
>
> [[-1, 0, 0],
> [ 0, -1, 0],
> [ 0, 0, 1]],
>
> [[ 1, 0, 0],
> [ 0, 1, 0],
> [ 0, 0, -1]],
>
> [[-1, 1, 0],
> [-1, 0, 0],
> [ 0, 0, 1]],
>
> [[ 1, -1, 0],
> [ 1, 0, 0],
> [ 0, 0, -1]],
>
> [[ 0, 1, 0],
> [-1, 1, 0],
> [ 0, 0, 1]],
>
> [[ 0, -1, 0],
> [ 1, -1, 0],
> [ 0, 0, -1]],
>
> [[ 0, -1, 0],
> [-1, 0, 0],
> [ 0, 0, -1]],
>
> [[ 0, 1, 0],
> [ 1, 0, 0],
> [ 0, 0, 1]],
>
> [[-1, 0, 0],
> [-1, 1, 0],
> [ 0, 0, -1]],
>
> [[ 1, 0, 0],
> [ 1, -1, 0],
> [ 0, 0, 1]],
>
> [[-1, 1, 0],
> [ 0, 1, 0],
> [ 0, 0, -1]],
>
> [[ 1, -1, 0],
> [ 0, -1, 0],
> [ 0, 0, 1]],
>
> [[ 0, 1, 0],
> [ 1, 0, 0],
> [ 0, 0, -1]],
>
> [[ 0, -1, 0],
> [-1, 0, 0],
> [ 0, 0, 1]],
>
> [[ 1, 0, 0],
> [ 1, -1, 0],
> [ 0, 0, -1]],
>
> [[-1, 0, 0],
> [-1, 1, 0],
> [ 0, 0, 1]],
>
> [[ 1, -1, 0],
> [ 0, -1, 0],
> [ 0, 0, -1]],
>
> [[-1, 1, 0],
> [ 0, 1, 0],
> [ 0, 0, 1]]];;
gap>
gap> P:=Group(rotations);;
gap> char:=Irr(P);;
gap> rep:=IrreducibleRepresentationsDixon(P, char:unitary);;
gap> List(rep, r -> List(P, g -> g^r));
[ [ [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ],
[ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ] ],
[ [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ],
[ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ] ],
[ [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ],
[ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ] ],
[ [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ],
[ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ] ],
[ [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ],
[ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ] ],
[ [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ],
[ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ] ],
[ [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ],
[ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ -1 ] ] ],
[ [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ],
[ [ -1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ], [ [ 1 ] ], [ [ -1 ] ], [ [ 1 ] ] ],
[ [ [ 1, 0 ], [ 0, 1 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1, 0 ], [ 0, -1 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ -1, 0 ], [ 0, -1 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ 1, 0 ], [ 0, 1 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ -1, 0 ], [ 0, 1 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ 1, 0 ], [ 0, -1 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ -1, 0 ], [ 0, 1 ] ] ],
[ [ [ 1, 0 ], [ 0, 1 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1, 0 ], [ 0, -1 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ 1, 0 ], [ 0, 1 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1, 0 ], [ 0, -1 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1, 0 ], [ 0, 1 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1, 0 ], [ 0, 1 ] ] ],
[ [ [ 1, 0 ], [ 0, 1 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ 1, 0 ], [ 0, 1 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ -1, 0 ], [ 0, -1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ -1, 0 ], [ 0, -1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ -1, 0 ], [ 0, 1 ] ], [ [ 1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ -1, 0 ], [ 0, 1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ 1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ 1, 0 ], [ 0, -1 ] ] ],
[ [ [ 1, 0 ], [ 0, 1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ 1, 0 ], [ 0, 1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ 1, 0 ], [ 0, 1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ], [ [ 1, 0 ], [ 0, 1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, -1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, -1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ 1, 0 ], [ 0, -1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ],
[ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ], [ [ 1, 0 ], [ 0, -1 ] ],
[ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ], [ [ -1/2, -1/2*E(12)^7+1/2*E(12)^11 ], [ -1/2*E(12)^7+1/2*E(12)^11, 1/2 ] ],
[ [ 1, 0 ], [ 0, -1 ] ], [ [ -1/2, 1/2*E(12)^7-1/2*E(12)^11 ], [ 1/2*E(12)^7-1/2*E(12)^11, 1/2 ] ] ] ]
But this is very different to the restults given by the rep variable of your code, as shown below:
In [31]: rep
Out[31]:
array([[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 0. , 0. , 0. ,
1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , -1. ,
1. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 0. , 0. , 0. , -1. ,
0. , 0. ],
[ 0. , 0. , 0. , 0. ,
0. , 1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 0. , 0. , 0. ,
0. , 0. ],
[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 0. , 0. , 0. ,
-1. , 0. ],
[ 1.41421356, 0. , 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]],
[[ 1. , 1. , 0. , 0. ,
0. , -1.41421356],
[ 0. , 1. , 0. , 0. ,
0. , 0. ],
[ 0. , 0. , 1. , 0. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
0. , 0. ],
[ 0. , 0. , 0. , 1. ,
-1. , 0. ],
[ 0. , 1.41421356, 0. , 0. ,
0. , -1. ]]])
If you are willing to provide a more intuitive explanation of these results given by your code, I would be willing to conduct more GAP based comparative studies in order to pinpoint the cause of the problem you are referring to.
Anyway, the point group here is very simple, as shown below:
gap> IdGroup(P);
[ 24, 14 ]
gap> StructureDescription(P);
"C2 x C2 x S3"
So, I don't understand why such a simple finite unimodular matrtix group can trigger such a strange problem.
Regards,
Zhao
from spgrep.
In this test, I consider the action of rotations on 3x3 symmetric matrices (for example, elastic constants). The procedure to generate representation matrices is described in an example page: https://spglib.github.io/spgrep/examples/symmetry_adapted_tensor.html.
from spgrep.
#100 (comment)
This issue seems to happen when input representation matrices have noise. The noise is amplified for higher ranks. Currently, I have no idea to properly handle it.
If I understand you correctly, you are referring to the cumulative error caused by decimal number approximations. As far as I know, it's a difficult thing to do float point computation in group theory. Strictly speaking, it is impossible to conduct exact computational group theory calculation over the real number field, so GAP only has extensive supports on the complete Cyclotomic field.
See here for the related discussion.
#100 (comment)
In this test, I consider the action of rotations on 3x3 symmetric matrices (for example, elastic constants). The procedure to generate representation matrices is described in an example page: https://spglib.github.io/spgrep/examples/symmetry_adapted_tensor.html.
I'll try to see if I can understand and reproduce it in GAP.
One thing that puzzles me is that your example page doesn't use any imprecise numbers, such as approximate floating point representations of irrational numbers, so this example cannot reflect the noise problem you've mentioned above.
from spgrep.
#100 (comment)
This issue seems to happen when input representation matrices have noise. The noise is amplified for higher ranks. Currently, I have no idea to properly handle it.
I'm not sure if using higher precision floating point numbers would help avoid this problem.
from spgrep.
Related Issues (13)
- Add co-representation of space group HOT 5
- Add double representation HOT 2
- Decomposition of Kronecker product of small representations
- Construct the little group based on Wyckoff positions. HOT 2
- Interface with DFT codes to calculate/analyze the irreducible representations of electronic states. HOT 2
- Calculate topological invariants to facilitate topological material analysis. HOT 2
- About the k-point used in the document example. HOT 7
- Identify a possible isomorphism between two space groups. HOT 3
- Missing data file in lattice vibration example HOT 2
- Some confusion on the description in the joss paper of this project. HOT 3
- Get the irreducible representations for a space-group represented on non-standard basis. HOT 5
- Irreducible Force Constant Examples HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from spgrep.