10
0
mirror of https://github.com/LCPQ/QUESTDB_website.git synced 2024-12-30 16:15:50 +01:00
QUESTDB_website/tools/lib/formats/default/doubleColumnHandler.py

62 lines
2.7 KiB
Python
Raw Normal View History

2020-07-03 14:51:12 +02:00
from ..formatHandlerBase import formatHandlerBase
from ..formatName import formatName
2020-07-08 13:24:29 +02:00
from ...data import dataFileBase,DataType,method,excitationValue,datafileSelector,AbsDataFile,getSubtablesRange,state
from ...LaTeX import newCommand,extractMath
2020-07-03 14:51:12 +02:00
import re
2020-07-08 13:24:29 +02:00
from TexSoup import TexSoup
import numpy as np
2020-07-03 14:51:12 +02:00
from ...utils import getValFromCell
@formatName("doubleColumn")
class doubleColumnHandler(formatHandlerBase):
def readFromTable(self,table):
2020-07-08 13:24:29 +02:00
datalist=list()
datacls=dict()
subtablesMol=getSubtablesRange(table)
for rangeMol in subtablesMol:
mymolecule=str(table[rangeMol[0],0])
moltable=table[rangeMol,:]
subtablestrans=getSubtablesRange(moltable,firstindex=0,column=1,count=2)
for rangeTrans in subtablestrans:
mytrans=moltable[rangeTrans,:]
mytransdesc=mytrans[0:2,1]
for i in range(2):
2020-10-08 18:41:13 +02:00
mathsoup=extractMath(mytransdesc[i],Soup=True,commands=self.Commands)
2020-07-08 13:24:29 +02:00
mytransdesc[i]=str(mathsoup)
for colindex in range(3,np.size(table,1)):
col=mytrans[:,colindex]
mybasis=str(table[1,colindex])
for index,cell in enumerate(col):
methodnameAT1=str(mytrans[index,2])
PTString=r"($\%T_1$)"
HasT1=methodnameAT1.endswith(PTString)
if HasT1:
methodname=methodnameAT1[:-len(PTString)]
else:
methodname=str(methodnameAT1)
mymethod=method(methodname,mybasis)
strcell=str(cell)
if strcell!="":
2020-07-03 14:51:12 +02:00
if HasT1:
2020-07-08 13:24:29 +02:00
m=re.match(r"^(?P<value>[-+]?\d+\.?\d*)\s*(?:\((?P<T1>\d+\.?\d*)\\\%\))?",strcell)
val,unsafe=getValFromCell(TexSoup(m.group("value")))
T1=m.group("T1")
2020-07-03 14:51:12 +02:00
else:
2020-07-08 13:24:29 +02:00
m=re.match(r"^[-+]?\d+\.?\d*",strcell)
val,unsafe=getValFromCell(TexSoup(m.group(0)))
T1=None
if (mymolecule,mymethod.name,mymethod.basis) in datacls:
data=datacls[(mymolecule,mymethod.name,mymethod.basis)]
else:
data=AbsDataFile()
data.molecule=mymolecule
data.method=mymethod
datacls[(mymolecule,mymethod.name,mymethod.basis)]=data
infin=mytransdesc[0].split(r"\rightarrow")
for i,item in enumerate(infin):
m=re.match(r"^(?P<number>\d)\\[,:;\s]\s*\^(?P<multiplicity>\d)(?P<sym>\S*)",item.strip())
infin[i]=state(m.group("number"),m.group("multiplicity"),m.group("sym"))
data.excitations.append(excitationValue(infin[0],infin[1],val,type=mytransdesc[1],isUnsafe=unsafe,T1=T1))
for value in datacls.values():
datalist.append(value)
return datalist