QPMS
/
qpms
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

analysis: fix symmetry group info

This commit is contained in:
Marek Nečada 2022-01-03 23:51:13 +02:00
parent 63f7f216b0
commit 54936ee722
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
qpms/analysis.py
View File

@ -49,7 +49,8 @@ class FiniteRectLatAnalysis:
meta = data['meta'][()] meta = data['meta'][()]
scatter = 'symmetry_adapted' in meta.keys()# finiterectlat-constant-driving ; TODO unify scatter = 'symmetry_adapted' in meta.keys()# finiterectlat-constant-driving ; TODO unify
if scatter: if scatter:
thegroup = meta['symmetry_adapted'] #thegroup = meta['symmetry_adapted'] # always D2h, this option does something else!
thegroup = 'D2h'
else: # finiterectlat-modes else: # finiterectlat-modes
thegroup = 'D2h' if meta['D2'] else 'D4h' thegroup = 'D2h' if meta['D2'] else 'D4h'
Nx, Ny = meta['size'] Nx, Ny = meta['size']
@ -110,11 +111,13 @@ class FiniteRectLatAnalysis:
self.Nx, self.Ny, self.px, self.py = Nx, Ny, px, py self.Nx, self.Ny, self.px, self.py = Nx, Ny, px, py
self.posmap = posmap self.posmap = posmap
self.fullvec_poffsets = nelem * np.arange(Nx*Ny) self.fullvec_poffsets = nelem * np.arange(Nx*Ny)
self.group = thegroup
if scatter: if scatter:
self.typ = 'scatter' self.typ = 'scatter'
self.fvcs1_fromdata = data['fvcs1'] self.fvcs1_fromdata = data['fvcs1']
self.iris1_fromdata = data['iris1'] self.iris1_fromdata = data['iris1']
self.positions_fromdata = data['positions'] self.positions_fromdata = data['positions']
self.driving_symmetry_adapted = data['symmetry_adapted']
else: # modes else: # modes
self.typ = 'modes' self.typ = 'modes'
self.eigval = data['eigval'] self.eigval = data['eigval']