diff --git a/dpdata/abacus/scf.py b/dpdata/abacus/scf.py index 3a42b853c..531837ef8 100644 --- a/dpdata/abacus/scf.py +++ b/dpdata/abacus/scf.py @@ -1,11 +1,10 @@ import os,sys import numpy as np +from ..unit import EnergyConversion, PressureConversion, LengthConversion -bohr2ang = 0.5291770 -ry2ev = 13.605698 -kbar2evperang3 = 1e3 / 1.6021892e6 -# The consts are cited from $ABACUS_ROOT/source/src_global/constant.h - +bohr2ang = LengthConversion("bohr", "angstrom").value() +ry2ev = EnergyConversion("rydberg", "eV").value() +kbar2evperang3 = PressureConversion("kbar", "eV/angstrom^3").value() def get_block (lines, keyword, skip = 0, nlines = None): ret = [] diff --git a/dpdata/amber/md.py b/dpdata/amber/md.py index eb9d21b9f..64227a6a6 100644 --- a/dpdata/amber/md.py +++ b/dpdata/amber/md.py @@ -3,9 +3,11 @@ from scipy.io import netcdf import numpy as np from dpdata.amber.mask import pick_by_amber_mask +from dpdata.unit import EnergyConversion +from ..periodic_table import ELEMENTS -kcalmol2eV= 0.04336410390059322 -symbols = ['X', 'H', 'He', 'Li', 'Be', 'B', 'C', 'N', 'O', 'F', 'Ne', 'Na', 'Mg', 'Al', 'Si', 'P', 'S', 'Cl', 'Ar', 'K', 'Ca', 'Sc', 'Ti', 'V', 'Cr', 'Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', 'Ga', 'Ge', 'As', 'Se', 'Br', 'Kr', 'Rb', 'Sr', 'Y', 'Zr', 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', 'In', 'Sn', 'Sb', 'Te', 'I', 'Xe', 'Cs', 'Ba', 'La', 'Ce', 'Pr', 'Nd', 'Pm', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Ho', 'Er', 'Tm', 'Yb', 'Lu', 'Hf', 'Ta', 'W', 'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg', 'Tl', 'Pb', 'Bi', 'Po', 'At', 'Rn', 'Fr', 'Ra', 'Ac', 'Th', 'Pa', 'U', 'Np', 'Pu', 'Am', 'Cm', 'Bk', 'Cf', 'Es', 'Fm', 'Md', 'No', 'Lr', 'Rf', 'Db', 'Sg', 'Bh', 'Hs', 'Mt', 'Ds', 'Rg', 'Cn', 'Nh', 'Fl', 'Mc', 'Lv', 'Ts', 'Og'] +kcalmol2eV = EnergyConversion("kcal_mol", "eV").value() +symbols = ['X'] + ELEMENTS energy_convert = kcalmol2eV force_convert = energy_convert diff --git a/dpdata/amber/sqm.py b/dpdata/amber/sqm.py index 840f4e139..f569ebe4b 100644 --- a/dpdata/amber/sqm.py +++ b/dpdata/amber/sqm.py @@ -1,6 +1,8 @@ import numpy as np -from ..periodic_table import ELEMENTS +from dpdata.periodic_table import ELEMENTS +from dpdata.unit import EnergyConversion +kcal2ev = EnergyConversion("kcal_mol", "eV").value() START = 0 READ_ENERGY = 1 @@ -47,7 +49,8 @@ def parse_sqm_out(fname): if len(coords) == len(charges): flag = START elif flag == READ_FORCES: - forces.append([float(x) for x in line.strip().split()[-3:]]) + ll = line.strip() + forces.append([float(ll[-60:-40]), float(ll[-40:-20]), float(ll[-20:])]) if len(forces) == len(charges): flag = START @@ -62,7 +65,7 @@ def parse_sqm_out(fname): data['coords'] = np.array([coords]) energies = np.array(energies) - forces = np.array([forces], dtype=np.float32) + forces = np.array([forces], dtype=np.float32) * kcal2ev if len(forces) > 0: data['energies'] = energies data['forces'] = forces @@ -99,4 +102,4 @@ def make_sqm_in(data, fname=None, frame_idx=0, **kwargs): if fname is not None: with open(fname, 'w') as fp: fp.write(ret) - return ret \ No newline at end of file + return ret diff --git a/dpdata/cp2k/output.py b/dpdata/cp2k/output.py index d295d7ef6..9d13f7c99 100644 --- a/dpdata/cp2k/output.py +++ b/dpdata/cp2k/output.py @@ -3,11 +3,12 @@ import re from collections import OrderedDict from .cell import cell_to_low_triangle +from ..unit import EnergyConversion, LengthConversion, ForceConversion, PressureConversion #%% -AU_TO_ANG = 5.29177208590000E-01 -AU_TO_EV = 2.72113838565563E+01 -AU_TO_EV_EVERY_ANG = AU_TO_EV/AU_TO_ANG +AU_TO_ANG = LengthConversion("bohr", "angstrom").value() +AU_TO_EV = EnergyConversion("hartree", "eV").value() +AU_TO_EV_EVERY_ANG = ForceConversion("hartree/bohr", "eV/angstrom").value() delimiter_patterns=[] delimiter_p1 = re.compile(r'^ \* GO CP2K GO! \*+') delimiter_p2 = re.compile(r'^ \*+') @@ -221,9 +222,9 @@ def get_frames (fname) : coord_flag = False force_flag = False stress_flag = False - eV = 2.72113838565563E+01 # hatree to eV - angstrom = 5.29177208590000E-01 # Bohr to Angstrom - GPa = 160.21766208 # 1 eV/(Angstrom^3) = 160.21 GPa + eV = EnergyConversion("hartree", "eV").value() + angstrom = LengthConversion("bohr", "angstrom").value() + GPa = PressureConversion("eV/angstrom^3", "GPa").value() atom_symbol_list = [] cell = [] coord = [] diff --git a/dpdata/gaussian/log.py b/dpdata/gaussian/log.py index 9c20977db..5ed976fea 100644 --- a/dpdata/gaussian/log.py +++ b/dpdata/gaussian/log.py @@ -1,14 +1,12 @@ import numpy as np +from ..unit import LengthConversion, EnergyConversion, ForceConversion +from ..periodic_table import ELEMENTS +length_convert = LengthConversion("bohr", "angstrom").value() +energy_convert = EnergyConversion("hartree", "eV").value() +force_convert = ForceConversion("hartree/bohr", "eV/angstrom").value() -hartree2ev = 27.211386018 -bohr2ang = 0.52917721067 - -length_convert = bohr2ang -energy_convert = hartree2ev -force_convert = energy_convert / length_convert - -symbols = ['X', 'H', 'He', 'Li', 'Be', 'B', 'C', 'N', 'O', 'F', 'Ne', 'Na', 'Mg', 'Al', 'Si', 'P', 'S', 'Cl', 'Ar', 'K', 'Ca', 'Sc', 'Ti', 'V', 'Cr', 'Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', 'Ga', 'Ge', 'As', 'Se', 'Br', 'Kr', 'Rb', 'Sr', 'Y', 'Zr', 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', 'In', 'Sn', 'Sb', 'Te', 'I', 'Xe', 'Cs', 'Ba', 'La', 'Ce', 'Pr', 'Nd', 'Pm', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Ho', 'Er', 'Tm', 'Yb', 'Lu', 'Hf', 'Ta', 'W', 'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg', 'Tl', 'Pb', 'Bi', 'Po', 'At', 'Rn', 'Fr', 'Ra', 'Ac', 'Th', 'Pa', 'U', 'Np', 'Pu', 'Am', 'Cm', 'Bk', 'Cf', 'Es', 'Fm', 'Md', 'No', 'Lr', 'Rf', 'Db', 'Sg', 'Bh', 'Hs', 'Mt', 'Ds', 'Rg', 'Cn', 'Nh', 'Fl', 'Mc', 'Lv', 'Ts', 'Og'] +symbols = ["X"] + ELEMENTS def to_system_data(file_name, md=False): data = {} diff --git a/dpdata/gromacs/gro.py b/dpdata/gromacs/gro.py index d2349bb14..2114810ef 100644 --- a/dpdata/gromacs/gro.py +++ b/dpdata/gromacs/gro.py @@ -1,11 +1,10 @@ #!/usr/bin/env python3 import re import numpy as np -import json -import os +from ..unit import LengthConversion -nm2ang = 10. -ang2nm = 1. / nm2ang +nm2ang = LengthConversion("nm", "angstrom").value() +ang2nm = LengthConversion("angstrom", "nm").value() cell_idx_gmx2dp = [0, 4, 8, 1, 2, 3, 5, 6, 7] def _format_atom_name(atom_name): diff --git a/dpdata/pwmat/atomconfig.py b/dpdata/pwmat/atomconfig.py index 13034e70d..6535edc64 100644 --- a/dpdata/pwmat/atomconfig.py +++ b/dpdata/pwmat/atomconfig.py @@ -1,5 +1,5 @@ #!/usr/bin/python3 - +from ..periodic_table import ELEMENTS import numpy as np def _to_system_data_lower(lines) : @@ -39,13 +39,6 @@ def _to_system_data_lower(lines) : for jj in range(ii) : atom_types.append(idx) system['atom_types'] = np.array(atom_types, dtype = int) - ELEMENTS=['H', 'He', 'Li', 'Be', 'B', 'C', 'N', 'O', 'F', 'Ne', 'Na', 'Mg', 'Al', 'Si', 'P', 'S', 'Cl', 'Ar', 'K', 'Ca', \ - 'Sc', 'Ti', 'V', 'Cr','Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', 'Ga', 'Ge', 'As', 'Se', 'Br', 'Kr', 'Rb', 'Sr', 'Y', 'Zr', \ - 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', 'In', 'Sn', 'Sb', 'Te', 'I', 'Xe', 'Cs', 'Ba', 'La', 'Ce', 'Pr', 'Nd',\ - 'Pm', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Ho', 'Er', 'Tm', 'Yb', 'Lu', 'Hf', 'Ta', 'W', 'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg', \ - 'Tl', 'Pb', 'Bi', 'Po', 'At', 'Rn', 'Fr', 'Ra', 'Ac', 'Th', 'Pa', 'U', 'Np', 'Pu', 'Am', 'Cm', 'Bk', 'Cf', 'Es', 'Fm', \ - 'Md', 'No', 'Lr'] - system['atom_names'] = [ELEMENTS[ii-1] for ii in np.unique(sorted(atomic_number))] return system @@ -80,12 +73,6 @@ def from_system_data(system, f_idx = 0, skip_zeros = True) : for ii, jj in zip(atom_numbs, atom_names): for kk in range(ii): symbal.append(jj) - ELEMENTS=['H', 'He', 'Li', 'Be', 'B', 'C', 'N', 'O', 'F', 'Ne', 'Na', 'Mg', 'Al', 'Si', 'P', 'S', 'Cl', 'Ar', 'K', 'Ca', \ - 'Sc', 'Ti', 'V', 'Cr','Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', 'Ga', 'Ge', 'As', 'Se', 'Br', 'Kr', 'Rb', 'Sr', 'Y', 'Zr', \ - 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', 'In', 'Sn', 'Sb', 'Te', 'I', 'Xe', 'Cs', 'Ba', 'La', 'Ce', 'Pr', 'Nd',\ - 'Pm', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Ho', 'Er', 'Tm', 'Yb', 'Lu', 'Hf', 'Ta', 'W', 'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg', \ - 'Tl', 'Pb', 'Bi', 'Po', 'At', 'Rn', 'Fr', 'Ra', 'Ac', 'Th', 'Pa', 'U', 'Np', 'Pu', 'Am', 'Cm', 'Bk', 'Cf', 'Es', 'Fm', \ - 'Md', 'No', 'Lr'] atomic_numbers = [] for ii in symbal: atomic_numbers.append(ELEMENTS.index(ii)+1) diff --git a/dpdata/pwmat/movement.py b/dpdata/pwmat/movement.py index ef79906a7..121f38855 100644 --- a/dpdata/pwmat/movement.py +++ b/dpdata/pwmat/movement.py @@ -1,4 +1,5 @@ import numpy as np +from ..periodic_table import ELEMENTS def system_info (lines, type_idx_zero = False) : atom_names = [] @@ -32,12 +33,6 @@ def system_info (lines, type_idx_zero = False) : atom_types.append(idx) else : atom_types.append(idx+1) - ELEMENTS=['H', 'He', 'Li', 'Be', 'B', 'C', 'N', 'O', 'F', 'Ne', 'Na', 'Mg', 'Al', 'Si', 'P', 'S', 'Cl', 'Ar', 'K', 'Ca', \ - 'Sc', 'Ti', 'V', 'Cr','Mn', 'Fe', 'Co', 'Ni', 'Cu', 'Zn', 'Ga', 'Ge', 'As', 'Se', 'Br', 'Kr', 'Rb', 'Sr', 'Y', 'Zr', \ - 'Nb', 'Mo', 'Tc', 'Ru', 'Rh', 'Pd', 'Ag', 'Cd', 'In', 'Sn', 'Sb', 'Te', 'I', 'Xe', 'Cs', 'Ba', 'La', 'Ce', 'Pr', 'Nd',\ - 'Pm', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Ho', 'Er', 'Tm', 'Yb', 'Lu', 'Hf', 'Ta', 'W', 'Re', 'Os', 'Ir', 'Pt', 'Au', 'Hg', \ - 'Tl', 'Pb', 'Bi', 'Po', 'At', 'Rn', 'Fr', 'Ra', 'Ac', 'Th', 'Pa', 'U', 'Np', 'Pu', 'Am', 'Cm', 'Bk', 'Cf', 'Es', 'Fm', \ - 'Md', 'No', 'Lr'] for ii in np.unique(sorted(atomic_number)): atom_names.append(ELEMENTS[ii-1]) return atom_names, atom_numbs, np.array(atom_types, dtype = int), nelm diff --git a/dpdata/qe/traj.py b/dpdata/qe/traj.py index e98003307..49d3338dd 100644 --- a/dpdata/qe/traj.py +++ b/dpdata/qe/traj.py @@ -1,16 +1,14 @@ #!/usr/bin/python3 - import numpy as np -import dpdata,warnings +import dpdata, warnings +from ..unit import EnergyConversion, LengthConversion, ForceConversion, PressureConversion -ry2ev = 13.605693009 -hartree2ev = 27.211386018 -bohr2ang = 0.52917721067 -kbar2evperang3 = 1e3 / 1.602176621e6 +ry2ev = EnergyConversion("rydberg", "eV").value() +kbar2evperang3 = PressureConversion("kbar", "eV/angstrom^3").value() -length_convert = bohr2ang -energy_convert = hartree2ev -force_convert = energy_convert / length_convert +length_convert = LengthConversion("bohr", "angstrom").value() +energy_convert = EnergyConversion("hartree", "eV").value() +force_convert = ForceConversion("hartree/bohr", "eV/angstrom").value() def load_key (lines, key) : for ii in lines : diff --git a/dpdata/unit.py b/dpdata/unit.py new file mode 100644 index 000000000..7b5aeffc4 --- /dev/null +++ b/dpdata/unit.py @@ -0,0 +1,169 @@ +from abc import ABC +from scipy import constants + +AVOGADRO = constants.Avogadro # Avagadro constant +ELE_CHG = constants.elementary_charge # Elementary Charge, in C +BOHR = constants.value("atomic unit of length") # Bohr, in m +HARTREE = constants.value("atomic unit of energy") # Hartree, in Jole +RYDBERG = constants.Rydberg * constants.h * constants.c # Rydberg, in Jole + +# energy conversions +econvs = { + "eV": 1.0, + "hartree": HARTREE / ELE_CHG, + "kJ_mol": 1 / (ELE_CHG * AVOGADRO / 1000), + "kcal_mol": 1 / (ELE_CHG * AVOGADRO / 1000 / 4.184), + "rydberg": RYDBERG / ELE_CHG, + "J": 1 / ELE_CHG, + "kJ": 1000 / ELE_CHG +} + +# length conversions +lconvs = { + "angstrom": 1.0, + "bohr": BOHR * 1E10, + "nm": 10.0, + "m": 1E10, +} + +def check_unit(unit): + if unit not in econvs.keys() and unit not in lconvs.keys(): + try: + eunit = unit.split("/")[0] + lunit = unit.split("/")[1] + if eunit not in econvs.keys(): + raise RuntimeError(f"Invaild unit: {unit}") + if lunit not in lconvs.keys(): + raise RuntimeError(f"Invalid unit: {unit}") + except: + raise RuntimeError(f"Invalid unit: {unit}") + +class Conversion(ABC): + def __init__(self, unitA, unitB, check=True): + """ + Parent class for unit conversion + + Parameters + ---------- + unitA : str, unit to be converted + unitB : str, unit which unitA is converted to, i.e. `1 unitA = self._value unitB` + check : bool, whether to check unit validity + + Examples + -------- + >>> conv = Conversion("foo", "bar", check=False) + >>> conv.set_value("10.0") + >>> print(conv) + 1 foo = 10.0 bar + >>> conv.value() + 10.0 + """ + if check: + check_unit(unitA) + check_unit(unitB) + self.unitA = unitA + self.unitB = unitB + self._value = 0.0 + + def value(self): + return self._value + + def set_value(self, value): + self._value = value + + def __repr__(self): + return f"1 {self.unitA} = {self._value} {self.unitB}" + + def __str__(self): + return self.__repr__() + +class EnergyConversion(Conversion): + def __init__(self, unitA, unitB): + """ + Class for energy conversion + + Examples + -------- + >>> conv = LengthConversion("eV", "kcal_mol") + >>> conv.value() + 23.06054783061903 + """ + super().__init__(unitA, unitB) + self.set_value(econvs[unitA] / econvs[unitB]) + +class LengthConversion(Conversion): + def __init__(self, unitA, unitB): + """ + Class for length conversion + + Examples + -------- + >>> conv = LengthConversion("angstrom", "nm") + >>> conv.value() + 0.1 + """ + super().__init__(unitA, unitB) + self.set_value(lconvs[unitA] / lconvs[unitB]) + +class ForceConversion(Conversion): + def __init__(self, unitA, unitB): + """ + Class for force conversion + + Parameters + ---------- + unitA, unitB : str, in format of "energy_unit/length_unit" + + Examples + -------- + >>> conv = ForceConversion("kJ_mol/nm", "eV/angstrom") + >>> conv.value() + 0.0010364269656262175 + """ + super().__init__(unitA, unitB) + econv = EnergyConversion(unitA.split("/")[0], unitB.split("/")[0]).value() + lconv = LengthConversion(unitA.split("/")[1], unitB.split("/")[1]).value() + self.set_value(econv / lconv) + +class PressureConversion(Conversion): + def __init__(self, unitA, unitB): + """ + Class for pressure conversion + + Parameters + ---------- + unitA, unitB : str, in format of "energy_unit/length_unit^3", or in `["Pa", "pa", "kPa", "kpa", "bar", "kbar"]` + + Examples + -------- + >>> conv = PressureConversion("kbar", "eV/angstrom^3") + >>> conv.value() + 0.0006241509074460763 + """ + super().__init__(unitA, unitB, check=False) + unitA, factorA = self._convert_unit(unitA) + unitB, factorB = self._convert_unit(unitB) + eunitA, lunitA = self._split_unit(unitA) + eunitB, lunitB = self._split_unit(unitB) + econv = EnergyConversion(eunitA, eunitB).value() * factorA / factorB + lconv = LengthConversion(lunitA, lunitB).value() + self.set_value(econv / lconv**3) + + def _convert_unit(self, unit): + if unit == "Pa" or unit == "pa": + return "J/m^3", 1 + elif unit == "kPa" or unit == "kpa": + return "kJ/m^3", 1 + elif unit == "GPa" or unit == "Gpa": + return "kJ/m^3", 1E6 + elif unit == "bar": + return "J/m^3", 1E5 + elif unit == "kbar": + return "kJ/m^3", 1E5 + else: + return unit, 1 + + def _split_unit(self, unit): + eunit = unit.split("/")[0] + lunit = unit.split("/")[1][:-2] + return eunit, lunit \ No newline at end of file diff --git a/tests/cp2k/aimd/deepmd/energy.raw b/tests/cp2k/aimd/deepmd/energy.raw index 5944c7889..ed4325d62 100644 --- a/tests/cp2k/aimd/deepmd/energy.raw +++ b/tests/cp2k/aimd/deepmd/energy.raw @@ -1,3 +1,3 @@ -7.597057812500000000e+04 --7.597049218750000000e+04 --7.597026562500000000e+04 +-7.597050000000000000e+04 +-7.597027343750000000e+04 diff --git a/tests/cp2k/aimd/deepmd/force.raw b/tests/cp2k/aimd/deepmd/force.raw index d17e071a8..e3dd01762 100644 --- a/tests/cp2k/aimd/deepmd/force.raw +++ b/tests/cp2k/aimd/deepmd/force.raw @@ -1,3 +1,3 @@ --1.181987579911947250e-02 -6.981573533266782761e-03 -2.535107731819152832e-02 2.352045150473713875e-03 6.612877477891743183e-04 -1.095135696232318878e-02 9.170096367597579956e-03 3.188682254403829575e-03 -1.015020068734884262e-02 4.770424682646989822e-03 7.456199382431805134e-04 -8.743807673454284668e-03 -5.867771804332733154e-03 1.279020961374044418e-02 -8.277409709990024567e-03 1.680318824946880341e-02 8.171994239091873169e-03 -1.793961599469184875e-02 -7.178519736044108868e-04 -1.655224710702896118e-02 -8.576685562729835510e-03 2.648287639021873474e-02 2.245652861893177032e-02 -2.479314804077148438e-02 -2.599848061800003052e-02 1.871763146482408047e-03 -1.798075251281261444e-02 -4.216609231662005186e-04 3.547093831002712250e-03 -4.467548802495002747e-03 1.272541843354701996e-02 -2.251772210001945496e-03 -1.833504997193813324e-02 4.828531818930059671e-04 -8.190506137907505035e-03 -1.827283017337322235e-02 1.298047136515378952e-02 -8.585427887737751007e-03 -1.525641232728958130e-02 1.295270305126905441e-02 1.116835791617631912e-02 -6.951748859137296677e-03 -2.566835097968578339e-02 2.900975756347179413e-02 -1.581999845802783966e-02 -3.503899322822690010e-03 2.792218001559376717e-04 -7.236626930534839630e-03 -1.866620848886668682e-03 2.224518451839685440e-03 5.501646548509597778e-03 -3.064292110502719879e-02 4.510228987783193588e-03 -9.799502789974212646e-03 4.271630663424730301e-03 -1.547238416969776154e-02 -2.558453381061553955e-02 -1.630285009741783142e-02 -1.694099791347980499e-02 7.224800065159797668e-03 9.183466434478759766e-03 -7.137896493077278137e-03 -1.933829486370086670e-02 6.450743228197097778e-02 -3.525239601731300354e-02 -9.753017127513885498e-02 -3.070719912648200989e-02 -5.610661115497350693e-03 -3.235116228461265564e-02 -1.719913817942142487e-02 -5.414743209257721901e-04 8.389509283006191254e-03 6.940178573131561279e-02 3.834697604179382324e-02 -5.660694837570190430e-02 -8.162738755345344543e-03 -1.202093530446290970e-02 -3.740955144166946411e-03 1.456787064671516418e-02 5.491876509040594101e-04 8.945382200181484222e-03 3.462864458560943604e-02 7.256064563989639282e-02 8.536062203347682953e-03 7.866033352911472321e-03 -5.486220121383666992e-03 1.132776588201522827e-02 1.698522083461284637e-02 -1.576600410044193268e-02 -1.265805494040250778e-02 -2.978571504354476929e-02 9.144899435341358185e-03 2.029114635661244392e-03 -6.899709254503250122e-02 1.742369830608367920e-01 -2.683743238449096680e-01 -4.429239407181739807e-02 6.901354994624853134e-03 2.609206922352313995e-02 2.234185859560966492e-02 2.111544087529182434e-02 -5.753100384026765823e-03 -8.917613886296749115e-03 -4.039614275097846985e-02 4.068822041153907776e-02 5.109296180307865143e-03 -1.294293347746133804e-02 3.259593248367309570e-02 -1.423830687999725342e-01 7.929538935422897339e-02 -4.696274176239967346e-02 -2.426092978566884995e-03 -1.106515377759933472e-01 -9.850924462080001831e-02 1.005635559558868408e-01 -3.671586886048316956e-02 -9.566508978605270386e-02 1.548637151718139648e-01 2.369528636336326599e-03 2.690098881721496582e-01 -9.338349848985671997e-02 1.726326048374176025e-01 -1.000776235014200211e-02 -5.573894456028938293e-02 -2.007779628038406372e-01 -3.753501921892166138e-02 -2.343843281269073486e-01 1.200124099850654602e-01 4.204010590910911560e-02 -3.067891672253608704e-02 -5.148376803845167160e-03 1.371575593948364258e-01 -3.283195942640304565e-02 -1.383109539747238159e-01 1.536984890699386597e-01 2.477211654186248779e-01 2.673551440238952637e-01 -5.387232303619384766e-01 -3.399512544274330139e-02 -3.363141715526580811e-01 -6.386671960353851318e-02 -2.391789257526397705e-01 -3.341148421168327332e-02 -5.645525455474853516e-01 -2.335435897111892700e-02 2.483330816030502319e-01 1.459018737077713013e-01 9.233356118202209473e-01 4.329773783683776855e-01 -4.229253828525543213e-01 -2.839886210858821869e-02 -3.265270292758941650e-01 6.476197391748428345e-02 -6.324636340141296387e-01 -1.054913327097892761e-01 2.155288904905319214e-01 -2.199295908212661743e-01 -4.232226014137268066e-01 2.678050696849822998e-01 3.901438117027282715e-01 5.183776021003723145e-01 2.177307903766632080e-01 -3.374685943126678467e-01 5.977681279182434082e-01 -7.905463576316833496e-01 -1.668028905987739563e-02 -1.058150410652160645e+00 -3.875943124294281006e-01 -3.925251588225364685e-02 4.368921369314193726e-02 1.205425709486007690e-01 -7.931236177682876587e-02 1.952598616480827332e-02 2.281134203076362610e-02 -1.511140167713165283e-02 -2.594140172004699707e-02 -2.723168432712554932e-01 -3.352718485984951258e-04 6.987435370683670044e-02 -2.582673169672489166e-02 5.285776779055595398e-02 6.644758861511945724e-03 2.587023377418518066e-01 -1.801109127700328827e-02 -3.912601992487907410e-02 2.104566246271133423e-01 -4.200668260455131531e-02 1.420020312070846558e-02 2.380419820547103882e-01 2.954865992069244385e-02 -2.349371276795864105e-02 2.397044599056243896e-01 -5.117009393870830536e-03 -2.002117969095706940e-02 1.887647360563278198e-01 4.477936029434204102e-02 2.717604674398899078e-02 -9.436257183551788330e-02 -8.993718773126602173e-03 -3.621655981987714767e-03 -2.380774617195129395e-01 1.019843518733978271e-01 4.974570591002702713e-03 2.735705114901065826e-02 6.167408078908920288e-02 1.060477178543806076e-02 -3.488719463348388672e-01 8.042410947382450104e-03 -8.855650573968887329e-02 1.035048961639404297e-01 -1.011960506439208984e-01 -6.915342062711715698e-02 2.238952666521072388e-01 7.932367734611034393e-03 -1.187361180782318115e-01 -9.918287396430969238e-02 6.672424077987670898e-02 -6.066209450364112854e-02 7.064928859472274780e-02 -5.082299560308456421e-02 -1.655636169016361237e-02 2.118203341960906982e-01 4.289011284708976746e-02 3.366689980030059814e-01 -7.251539081335067749e-02 9.767106920480728149e-02 2.942092716693878174e-01 -1.480986326932907104e-01 6.034224852919578552e-02 -4.308999180793762207e-01 1.146930575370788574e+00 2.093196809291839600e-01 -8.582702279090881348e-02 -2.207106947898864746e-01 -1.496844831854104996e-02 1.843213587999343872e-01 -3.934044837951660156e-01 2.001809477806091309e-01 -8.990695476531982422e-01 8.304765820503234863e-01 -7.793455719947814941e-01 -3.349324762821197510e-01 2.378815561532974243e-01 3.620493710041046143e-01 1.195638060569763184e+00 4.119195342063903809e-01 2.902950346469879150e-01 -6.464575976133346558e-02 -3.813989758491516113e-01 -1.106751933693885803e-01 1.081159189343452454e-01 -3.471554815769195557e-02 9.297623299062252045e-03 -1.349829137325286865e-03 -2.145842649042606354e-03 1.335945213213562965e-03 4.220723174512386322e-03 -1.248779669404029846e-01 -3.870024578645825386e-03 4.961200524121522903e-03 -7.309851795434951782e-02 -4.484003875404596329e-04 -4.894866142421960831e-03 1.382533553987741470e-02 3.704959759488701820e-03 -4.498402122408151627e-03 -1.058152914047241211e-01 -1.672759768553078175e-03 -4.989482928067445755e-03 -8.488651365041732788e-02 -1.313319546170532703e-03 -7.889172993600368500e-03 -9.988170117139816284e-02 -9.862751699984073639e-03 3.771294141188263893e-03 -5.090115591883659363e-02 6.931694224476814270e-04 -1.896959962323307991e-03 -8.603682368993759155e-02 5.891426000744104385e-03 1.052301097661256790e-02 -6.807098537683486938e-02 -1.224050763994455338e-02 -8.844595140544697642e-05 -1.171086076647043228e-02 9.057481773197650909e-03 -6.648872513324022293e-04 -7.828289270401000977e-02 -4.207353107631206512e-03 1.285037375055253506e-03 -5.484163016080856323e-02 4.715403076261281967e-03 -1.758480258285999298e-02 -8.468031138181686401e-02 -9.094506502151489258e-03 -5.433769430965185165e-03 -1.080428957939147949e-01 -8.356085163541138172e-04 -8.645591326057910919e-03 -8.232723921537399292e-02 -2.507082745432853699e-02 4.364704713225364685e-02 1.537128835916519165e-01 5.729600414633750916e-02 5.963776633143424988e-02 1.158852726221084595e-01 -6.959976255893707275e-02 1.238217577338218689e-01 1.664686501026153564e-01 2.036262303590774536e-02 9.312792867422103882e-02 1.595683246850967407e-01 8.057528734207153320e-02 5.345786362886428833e-02 2.299466133117675781e-01 5.173162370920181274e-02 -4.056172445416450500e-03 2.901870310306549072e-01 1.976664178073406219e-02 -3.087277896702289581e-02 1.369883716106414795e-01 -1.363353151828050613e-02 4.292251169681549072e-02 1.564999669790267944e-01 1.284152865409851074e-01 -6.114494800567626953e-02 3.792459368705749512e-01 -4.886947199702262878e-02 6.220578402280807495e-02 2.492766827344894409e-01 4.413247108459472656e-02 1.222508102655410767e-01 2.082742750644683838e-01 -1.421974319964647293e-02 -4.172386135905981064e-03 1.484585851430892944e-01 2.711948193609714508e-02 -2.002323791384696960e-02 2.530520856380462646e-01 -3.040535189211368561e-02 -8.523207157850265503e-03 2.680878937244415283e-01 -9.880235418677330017e-03 3.497625887393951416e-02 1.259254366159439087e-01 3.049739636480808258e-02 6.191216409206390381e-03 1.380322426557540894e-01 -3.762089461088180542e-02 4.634927585721015930e-02 1.610230468213558197e-02 2.170730940997600555e-02 -7.616019248962402344e-02 -1.117401407100260258e-03 2.092877984978258610e-04 5.715562496334314346e-03 -1.721147894859313965e-01 -1.635478809475898743e-02 2.458540350198745728e-02 8.387092500925064087e-02 -1.793704368174076080e-02 -6.610923260450363159e-02 4.725636169314384460e-02 -4.781737457960844040e-03 4.136905074119567871e-03 -1.320595741271972656e-01 9.640093892812728882e-03 1.156122237443923950e-02 2.084830403327941895e-01 -3.595944959670305252e-03 2.512173540890216827e-02 -1.716694682836532593e-01 5.623002536594867706e-03 1.370912231504917145e-02 -1.162035763263702393e-01 3.378943726420402527e-02 -3.141476586461067200e-02 4.943099990487098694e-02 9.966830164194107056e-02 3.148675709962844849e-02 3.660736605525016785e-03 6.501960009336471558e-02 -1.878396607935428619e-02 1.421562954783439636e-02 1.706852577626705170e-02 -9.174210019409656525e-03 -1.557013839483261108e-01 -1.136684697121381760e-02 4.417669493705034256e-03 7.971448451280593872e-02 4.406356718391180038e-03 -2.304119803011417389e-02 1.649928279221057892e-02 -4.582734312862157822e-03 -3.686447627842426300e-03 -7.609025388956069946e-02 --1.104494463652372360e-02 3.476131334900856018e-02 2.383937090635299683e-01 -2.983302436769008636e-02 6.788226217031478882e-02 -7.082360796630382538e-03 -6.491459608078002930e-01 3.486107289791107178e-02 1.303477287292480469e-01 -7.131726294755935669e-02 1.783661693334579468e-01 -2.583156526088714600e-01 -1.776951104402542114e-01 1.444620639085769653e-01 -1.592680066823959351e-01 8.508963137865066528e-02 8.697990328073501587e-02 -1.015986874699592590e-01 7.329495251178741455e-02 3.711572587490081787e-01 1.321593374013900757e-01 5.504114925861358643e-02 -6.570916622877120972e-02 2.412003278732299805e-02 2.140499949455261230e-01 -7.536160200834274292e-02 1.978813558816909790e-01 -9.274689108133316040e-02 -5.994388461112976074e-01 -1.937819868326187134e-01 -1.141302436590194702e-01 3.028306961059570312e-01 -1.305004507303237915e-01 1.491460949182510376e-01 5.856458842754364014e-02 1.309930831193923950e-01 3.258333206176757812e-01 -8.520995825529098511e-02 -2.395059764385223389e-01 -2.042067795991897583e-01 -5.118007063865661621e-01 -4.316213726997375488e-02 -7.233747094869613647e-02 1.305323392152786255e-01 -1.028775498270988464e-01 -1.012448966503143311e-02 3.805541247129440308e-02 -1.114979460835456848e-01 9.408849477767944336e-02 -1.270032376050949097e-01 -9.719438105821609497e-02 -1.077755019068717957e-01 -1.132776588201522827e-02 -1.396108977496623993e-02 1.897062808275222778e-01 -1.100956648588180542e-01 -1.400613635778427124e-01 -2.287099137902259827e-02 -2.160446532070636749e-02 -5.284748226404190063e-02 3.776950389146804810e-02 -1.654561460018157959e-01 1.634496599435806274e-01 -1.084990128874778748e-01 8.678758889436721802e-02 2.098215520381927490e-01 6.759327650070190430e-01 -4.112494885921478271e-01 -3.361912667751312256e-01 1.149864196777343750e-01 -1.030446738004684448e-01 -1.138047426939010620e-01 -6.231747269630432129e-01 -3.622010648250579834e-01 -1.078145802021026611e-01 1.205538883805274963e-01 1.539226919412612915e-01 2.042021602392196655e-02 3.630397608503699303e-03 -1.659852825105190277e-02 1.282276064157485962e-01 3.265661001205444336e-01 2.516091763973236084e-01 6.139280274510383606e-02 3.997751772403717041e-01 -3.187895417213439941e-01 6.818771362304687500e-02 6.308767199516296387e-02 -3.385834395885467529e-02 -2.106453478336334229e-02 -4.315339494496583939e-03 -2.343010306358337402e-01 -1.357208192348480225e-01 -7.493480294942855835e-02 1.520822942256927490e-01 -2.871083915233612061e-01 -1.007563918828964233e-01 3.371898829936981201e-02 1.746123582124710083e-01 -4.844523966312408447e-02 2.442136593163013458e-02 1.908154487609863281e-01 5.124671384692192078e-02 1.412152685225009918e-02 3.879537433385848999e-02 2.635992765426635742e-01 -4.273878037929534912e-01 -4.049178957939147949e-02 -1.563117653131484985e-01 3.645001351833343506e-02 -9.467676281929016113e-02 -7.148797810077667236e-02 -2.167239487171173096e-01 -1.213452741503715515e-01 2.505318820476531982e-01 -1.560865342617034912e-02 -3.080011904239654541e-01 2.920737266540527344e-01 -3.169656032696366310e-03 6.070307493209838867e-01 -2.779676318168640137e-01 4.885317087173461914e-01 2.117822766304016113e-01 8.048170059919357300e-02 -1.967896670103073120e-01 7.976076006889343262e-02 -3.200874328613281250e-01 2.314280867576599121e-01 5.967787653207778931e-02 2.285798192024230957e-01 5.649356245994567871e-01 -5.458348989486694336e-02 -1.947559267282485962e-01 6.182063370943069458e-02 2.359593957662582397e-01 -1.382492482662200928e-01 -6.764315068721771240e-02 -2.699308693408966064e-01 -7.123858481645584106e-02 -6.276844143867492676e-01 -2.283124178647994995e-01 5.317709967494010925e-02 6.410788744688034058e-02 -3.124867379665374756e-01 -3.441524505615234375e-02 2.965572178363800049e-01 1.769356131553649902e-01 6.501157879829406738e-01 3.025216460227966309e-01 -3.078834414482116699e-01 1.825997489504516125e-03 7.202688604593276978e-02 -2.511818706989288330e-01 -9.704058170318603516e-01 -2.338984012603759766e-01 1.342753469944000244e-01 -2.212151437997817993e-01 -3.379781842231750488e-01 2.369970828294754028e-01 7.971648573875427246e-01 8.250258564949035645e-01 4.743896126747131348e-01 -3.288477063179016113e-01 5.826305150985717773e-01 -6.866187453269958496e-01 3.443632647395133972e-02 -9.923686981201171875e-01 -3.228811919689178467e-01 7.178880274295806885e-02 -1.244876757264137268e-01 2.564778178930282593e-02 -1.374753415584564209e-01 -9.046375006437301636e-02 1.591620892286300659e-01 3.975151777267456055e-01 2.328185364603996277e-02 -6.311548948287963867e-01 3.031587786972522736e-02 -5.235480666160583496e-01 1.079878732562065125e-01 -1.716468483209609985e-01 7.066404223442077637e-01 6.992629170417785645e-01 9.035884588956832886e-02 -2.569920420646667480e-01 2.911028563976287842e-01 -2.731956541538238525e-01 1.626562178134918213e-01 -1.209534332156181335e-01 2.627821639180183411e-02 -3.783429786562919617e-02 4.602382779121398926e-01 3.054902553558349609e-01 5.480573773384094238e-01 1.684962213039398193e-01 -1.160549670457839966e-01 1.016665622591972351e-01 1.942725502885878086e-03 -2.473349869251251221e-02 -5.631847307085990906e-02 -2.918592989444732666e-01 -5.663929581642150879e-01 3.910899758338928223e-01 1.680426746606826782e-01 -1.523800343275070190e-01 -1.604018658399581909e-01 -3.901119530200958252e-01 1.903649866580963135e-01 -1.977188587188720703e-01 -1.694372445344924927e-01 4.718992412090301514e-01 4.436777830123901367e-01 3.941043615341186523e-01 -6.879305243492126465e-01 2.727374807000160217e-02 -1.745609343051910400e-01 9.348119795322418213e-02 1.866610646247863770e-01 2.327002584934234619e-01 -4.432057440280914307e-01 3.967453837394714355e-01 5.185770988464355469e-01 9.205937385559082031e-02 4.654653072357177734e-01 -5.490436404943466187e-02 7.019934058189392090e-02 3.512594997882843018e-01 -4.940852820873260498e-01 1.627919673919677734e-01 -1.851636469364166260e-01 6.728849411010742188e-01 3.807649388909339905e-02 -7.737894654273986816e-01 4.971485119313001633e-03 -1.257660239934921265e-01 4.114834666252136230e-01 -5.589094758033752441e-01 3.938091993331909180e-01 -1.208497166633605957e+00 5.434005856513977051e-01 -3.033346235752105713e-01 -4.138216376304626465e-01 4.525141412159428000e-05 6.551947593688964844e-01 1.291819453239440918e+00 4.318990409374237061e-01 -1.062909439206123352e-01 -4.860830008983612061e-01 -5.503559708595275879e-01 -1.151345148682594299e-01 -3.370392322540283203e-01 2.496016621589660645e-01 7.138925138860940933e-03 4.087025765329599380e-03 -5.033191759139299393e-03 -2.961396705359220505e-03 2.248378284275531769e-02 -1.294812709093093872e-01 7.782215252518653870e-03 -3.985723946243524551e-03 -7.532303780317306519e-02 1.993221975862979889e-02 8.204390294849872589e-03 1.980006508529186249e-02 -5.281046032905578613e-04 6.842733826488256454e-03 -1.057103946805000305e-01 -9.847325272858142853e-03 1.671731239184737206e-03 -8.024258911609649658e-02 4.962229286320507526e-04 -7.765245623886585236e-03 -9.957728534936904907e-02 -1.425522472709417343e-02 8.225987665355205536e-03 -4.023519158363342285e-02 -5.335553083568811417e-03 8.743293583393096924e-03 -8.677421510219573975e-02 -6.039007101207971573e-03 1.079194806516170502e-02 -5.816143751144409180e-02 -1.022990513592958450e-02 9.163411450572311878e-04 -1.089119259268045425e-02 1.388190034776926041e-02 1.676359213888645172e-03 -8.718302100896835327e-02 -2.216290915384888649e-03 8.119029924273490906e-03 -5.274360999464988708e-02 1.222816645167768002e-03 -1.122080814093351364e-02 -8.391463756561279297e-02 -1.298869866877794266e-02 -1.625605672597885132e-02 -1.077976152300834656e-01 -1.185638550668954849e-02 -1.393897831439971924e-02 -8.446948230266571045e-02 -1.552997715771198273e-02 4.289371520280838013e-02 1.527399867773056030e-01 6.501240283250808716e-02 5.235177278518676758e-02 1.159608662128448486e-01 -5.968558788299560547e-02 1.221289411187171936e-01 1.672075837850570679e-01 3.620216250419616699e-02 1.012747213244438171e-01 1.585223972797393799e-01 6.760561466217041016e-02 5.228698253631591797e-02 2.168802618980407715e-01 5.048926547169685364e-02 -2.149133756756782532e-02 2.762387990951538086e-01 1.242562755942344666e-02 -2.042330056428909302e-02 1.278681606054306030e-01 -1.921539567410945892e-02 4.353494569659233093e-02 1.464042663574218750e-01 1.288631707429885864e-01 -6.012164801359176636e-02 3.738620579242706299e-01 -3.939598426222801208e-02 5.242067947983741760e-02 2.547222673892974854e-01 4.888849705457687378e-02 1.314548552036285400e-01 2.049518823623657227e-01 -2.355078980326652527e-02 -1.129537075757980347e-02 1.436465084552764893e-01 2.126508019864559174e-02 -2.505334280431270599e-02 2.450934946537017822e-01 -4.344187304377555847e-02 -5.161746870726346970e-03 2.797349989414215088e-01 -6.929637398570775986e-03 3.583140671253204346e-02 1.322745084762573242e-01 3.419772908091545105e-02 2.383412560448050499e-03 1.350003927946090698e-01 -5.564895644783973694e-02 2.226575277745723724e-02 1.167949289083480835e-02 2.564058266580104828e-02 -9.204909205436706543e-02 -1.253567077219486237e-02 -2.450312674045562744e-02 1.618663780391216278e-02 -1.755765229463577271e-01 -2.139569260179996490e-02 3.943043947219848633e-02 7.631908357143402100e-02 -1.542301941663026810e-02 -5.633389949798583984e-02 4.985009133815765381e-02 -4.008349962532520294e-03 1.968128047883510590e-02 -1.294668763875961304e-01 1.023813267238438129e-03 -3.601652756333351135e-02 2.150717526674270630e-01 1.322781108319759369e-02 1.532223261892795563e-02 -1.604856848716735840e-01 -1.355588436126708984e-02 2.606070227921009064e-02 -1.114218384027481079e-01 5.919913575053215027e-02 1.487126108258962631e-03 3.223289176821708679e-02 9.704525768756866455e-02 4.147806391119956970e-02 7.003684877417981625e-04 5.199490487575531006e-02 1.038211397826671600e-03 1.214074902236461639e-02 2.444193512201309204e-02 9.088849648833274841e-03 -1.592788100242614746e-01 1.247602142393589020e-02 -2.815512195229530334e-02 7.731667160987854004e-02 5.528386216610670090e-03 -1.595163904130458832e-02 8.866705931723117828e-03 -3.608954697847366333e-02 -1.909198425710201263e-02 -7.481499016284942627e-02 --9.893090464174747467e-03 7.110025733709335327e-02 4.683017432689666748e-01 -6.107398495078086853e-02 1.347078084945678711e-01 -1.578657305799424648e-03 -1.261842370033264160e+00 4.915383458137512207e-02 2.711197435855865479e-01 -1.311180293560028076e-01 3.344485759735107422e-01 -4.808214902877807617e-01 -3.216167092323303223e-01 2.555764019489288330e-01 -2.926167249679565430e-01 1.411504745483398438e-01 1.599719822406768799e-01 -1.761498749256134033e-01 1.470336765050888062e-01 6.908518075942993164e-01 2.576049864292144775e-01 7.998599112033843994e-02 -1.543093770742416382e-01 7.660036534070968628e-02 4.400098621845245361e-01 -1.362982988357543945e-01 3.931432664394378662e-01 -1.867906451225280762e-01 -1.155184268951416016e+00 -3.807911872863769531e-01 -2.209132909774780273e-01 5.606989860534667969e-01 -2.357346713542938232e-01 2.812612056732177734e-01 1.144310608506202698e-01 2.691909074783325195e-01 5.828824639320373535e-01 -1.458211392164230347e-01 -4.240993559360504150e-01 -4.254646003246307373e-01 -1.020523309707641602e+00 -8.228095620870590210e-02 -1.227393224835395813e-01 2.231588959693908691e-01 -1.860954165458679199e-01 -9.655521251261234283e-03 7.481139153242111206e-02 -2.035866379737854004e-01 1.776267141103744507e-01 -2.375889569520950317e-01 -1.883476972579956055e-01 -1.761766225099563599e-01 -2.894804999232292175e-02 -1.338002085685729980e-02 3.621239364147186279e-01 -1.977126896381378174e-01 -2.472341954708099365e-01 -2.230174839496612549e-02 -1.701401732861995697e-02 -1.129578202962875366e-01 5.937705561518669128e-02 -3.134961426258087158e-01 3.318846821784973145e-01 -2.587897479534149170e-01 1.865679919719696045e-01 4.676409959793090820e-01 1.349116921424865723e+00 -8.273280262947082520e-01 -6.393083930015563965e-01 2.348702698945999146e-01 -1.994476765394210815e-01 -2.319345921277999878e-01 -1.305436015129089355e+00 -7.539456486701965332e-01 -1.434228122234344482e-01 2.326097637414932251e-01 3.093150258064270020e-01 3.671895340085029602e-02 -3.112577600404620171e-03 -3.046294488012790680e-02 2.379859387874603271e-01 5.817290544509887695e-01 3.898522555828094482e-01 1.008556336164474487e-01 7.762391567230224609e-01 -6.269609332084655762e-01 1.192472502589225769e-01 9.827578812837600708e-02 -4.342541843652725220e-02 -3.569976612925529480e-02 1.863998360931873322e-02 -4.682919979095458984e-01 -2.610754668712615967e-01 -7.949799299240112305e-02 1.178094893693923950e-01 -2.877820432186126709e-01 -1.527363806962966919e-01 5.917342379689216614e-02 3.124430179595947266e-01 -1.266787648200988770e-01 2.466356381773948669e-02 3.761935234069824219e-01 1.092194318771362305e-01 6.830958276987075806e-02 3.851101174950599670e-02 5.131099224090576172e-01 -8.361957669258117676e-01 -1.178449690341949463e-01 -1.629236042499542236e-01 -8.984463289380073547e-03 -1.315494626760482788e-01 -1.241508573293685913e-01 -3.026718199253082275e-01 -1.452256739139556885e-01 3.789070844650268555e-01 -2.524823416024446487e-03 -4.800285696983337402e-01 3.944936096668243408e-01 -1.162138651125133038e-03 8.724945783615112305e-01 -4.213492870330810547e-01 7.462729811668395996e-01 3.899874985218048096e-01 2.097485363483428955e-01 -1.942535340785980225e-01 1.942555904388427734e-01 -3.750139176845550537e-01 3.149148821830749512e-01 6.796608120203018188e-02 4.504814445972442627e-01 1.019416213035583496e+00 -2.055576443672180176e-01 -3.281714916229248047e-01 2.426272928714752197e-01 3.063994050025939941e-01 -5.537286996841430664e-01 -4.357289671897888184e-01 1.795658469200134277e-02 -1.048259362578392029e-01 -8.470062017440795898e-01 -3.549505472183227539e-01 3.670609593391418457e-01 1.665252149105072021e-01 -2.213976904749870300e-02 -4.017965868115425110e-02 3.127880692481994629e-01 1.877275556325912476e-01 2.954963743686676025e-01 1.175163835287094116e-01 -1.510980725288391113e-01 4.575740918517112732e-02 4.829621911048889160e-01 -5.672902464866638184e-01 -1.197457313537597656e+00 -3.361994922161102295e-01 6.976174563169479370e-02 -2.216095477342605591e-01 -2.158086299896240234e-01 1.900338381528854370e-01 1.096209406852722168e+00 1.025184154510498047e+00 6.810111999511718750e-01 -3.208042681217193604e-01 5.260307192802429199e-01 -5.059360265731811523e-01 9.677066653966903687e-02 -8.423340320587158203e-01 -1.967562437057495117e-01 1.529091596603393555e-01 -2.647413611412048340e-01 -5.476912483572959900e-02 -1.835166066884994507e-01 -1.947086155414581299e-01 2.880592048168182373e-01 7.779361009597778320e-01 7.493274658918380737e-02 -9.699820876121520996e-01 4.461738094687461853e-02 -1.026048064231872559e+00 2.282728254795074463e-01 -3.384255766868591309e-01 1.343114972114562988e+00 1.096434593200683594e+00 1.817188858985900879e-01 -4.267115890979766846e-01 3.620889782905578613e-01 -4.821353256702423096e-01 3.120321631431579590e-01 -4.564875364303588867e-01 3.491455316543579102e-02 -5.989950522780418396e-02 6.529316306114196777e-01 6.110997796058654785e-01 1.099410414695739746e+00 1.527127325534820557e-01 -2.757677733898162842e-01 1.618555784225463867e-01 8.699533343315124512e-02 -5.625213682651519775e-02 -1.119247451424598694e-01 -3.362915515899658203e-01 -1.203253149986267090e+00 7.881099581718444824e-01 3.196734786033630371e-01 -3.434947431087493896e-01 -3.287900984287261963e-01 -3.869170844554901123e-01 3.485536575317382812e-01 -2.719147205352783203e-01 -4.036462306976318359e-01 1.044531226158142090e+00 9.418208003044128418e-01 5.448548197746276855e-01 -1.327193737030029297e+00 2.001449614763259888e-01 -2.329656034708023071e-01 1.114362403750419617e-01 4.299717545509338379e-01 3.690345585346221924e-01 -8.036394119262695312e-01 8.094382882118225098e-01 7.861651778221130371e-01 1.207611188292503357e-01 5.581638813018798828e-01 -3.813305869698524475e-02 5.217076838016510010e-02 3.965011239051818848e-01 -8.052031993865966797e-01 2.708004117012023926e-01 9.015933424234390259e-02 1.360761523246765137e-01 -1.289413422346115112e-01 -1.314676046371459961e+00 2.050768435001373291e-01 -2.040915936231613159e-01 5.657923221588134766e-01 -6.838779449462890625e-01 5.440202355384826660e-01 -1.415313124656677246e+00 2.235332578420639038e-01 2.245771139860153198e-01 -4.179127514362335205e-01 -2.605195939540863037e-01 8.231499791145324707e-01 1.273505926132202148e+00 3.739458620548248291e-01 -3.938071429729461670e-01 -8.364529013633728027e-01 -6.571554541587829590e-01 -1.396422684192657471e-01 -7.621161341667175293e-01 5.516409873962402344e-01 4.955544136464595795e-03 9.621581993997097015e-03 -6.986201275140047073e-03 -7.330729160457849503e-03 4.049744457006454468e-02 -1.321984082460403442e-01 1.939485967159271240e-02 -1.333425566554069519e-02 -7.554672658443450928e-02 3.975645452737808228e-02 2.119051851332187653e-02 2.583701536059379578e-02 -4.761168733239173889e-03 1.824557594954967499e-02 -1.068581342697143555e-01 -1.780694536864757538e-02 7.972991093993186951e-03 -7.572824507951736450e-02 2.534079365432262421e-03 -7.397064007818698883e-03 -1.013287156820297241e-01 -1.824506185948848724e-02 1.250841654837131500e-02 -2.939645014703273773e-02 -1.128045842051506042e-02 1.917631551623344421e-02 -8.776049315929412842e-02 -1.817410066723823547e-02 1.099866535514593124e-02 -4.806677252054214478e-02 -8.188963867723941803e-03 1.888732425868511200e-03 -1.017076987773180008e-02 1.872277259826660156e-02 3.620113246142864227e-03 -9.292635321617126465e-02 -3.630397550296038389e-04 1.480698306113481522e-02 -4.907670244574546814e-02 -2.207034965977072716e-03 -4.811048042029142380e-03 -8.566195517778396606e-02 -1.708395220339298248e-02 -2.702949382364749908e-02 -1.082830354571342468e-01 -2.264319173991680145e-02 -1.930178515613079071e-02 -8.770547062158584595e-02 -5.845660343766212463e-03 4.200719669461250305e-02 1.513104438781738281e-01 7.280849665403366089e-02 4.511514678597450256e-02 1.158076301217079163e-01 -4.996527731418609619e-02 1.204417869448661804e-01 1.678513884544372559e-01 5.110633000731468201e-02 1.099300831556320190e-01 1.553280502557754517e-01 5.640126019716262817e-02 5.075923353433609009e-02 2.010067850351333618e-01 4.866635426878929138e-02 -3.678374364972114563e-02 2.593600451946258545e-01 5.267161875963211060e-03 -9.717226959764957428e-03 1.177632063627243042e-01 -2.470110170543193817e-02 4.407642036676406860e-02 1.362098455429077148e-01 1.288472414016723633e-01 -5.881449952721595764e-02 3.665066361427307129e-01 -2.970858290791511536e-02 4.231470078229904175e-02 2.592504918575286865e-01 5.263562500476837158e-02 1.400428414344787598e-01 2.000724524259567261e-01 -3.160451352596282959e-02 -1.877728104591369629e-02 1.365667134523391724e-01 1.508774794638156891e-02 -2.999140322208404541e-02 2.362375855445861816e-01 -5.841957777738571167e-02 -2.730511478148400784e-04 2.864455878734588623e-01 -3.346033627167344093e-03 3.664439171552658081e-02 1.374090015888214111e-01 3.791400045156478882e-02 -1.554489019326865673e-03 1.316636204719543457e-01 -7.310777902603149414e-02 -2.165897283703088760e-03 7.529732771217823029e-03 2.937588281929492950e-02 -1.078212708234786987e-01 -2.507854066789150238e-02 -4.747644811868667603e-02 2.562515623867511749e-02 -1.774122864007949829e-01 -2.635997720062732697e-02 5.408058315515518188e-02 6.827718764543533325e-02 -1.295218896120786667e-02 -4.677093774080276489e-02 5.198718979954719543e-02 -3.239589976146817207e-03 3.646955639123916626e-02 -1.268859952688217163e-01 -7.840836420655250549e-03 -8.373003453016281128e-02 2.224749922752380371e-01 3.069074451923370361e-02 5.268704611808061600e-03 -1.480256021022796631e-01 -3.176289424300193787e-02 3.841330856084823608e-02 -1.046834960579872131e-01 8.282603323459625244e-02 3.443067148327827454e-02 1.815867237746715546e-02 9.441193938255310059e-02 5.159072577953338623e-02 -1.701041823253035545e-03 3.985929861664772034e-02 1.937891915440559387e-02 9.308421984314918518e-03 3.155926242470741272e-02 2.736630849540233612e-02 -1.607952415943145752e-01 3.609520196914672852e-02 -6.125190481543540955e-02 7.475945353507995605e-02 6.751202512532472610e-03 -8.983434177935123444e-03 1.465528854168951511e-03 -6.945423781871795654e-02 -3.552904725074768066e-02 -7.128743827342987061e-02 +-1.181987673044204712e-02 -6.981573998928070068e-03 -2.535107918083667755e-02 2.352045383304357529e-03 6.612878059968352318e-04 -1.095135789364576340e-02 9.170097298920154572e-03 3.188682487234473228e-03 -1.015020161867141724e-02 4.770425148308277130e-03 7.456199964508414268e-04 -8.743808604776859283e-03 -5.867772269994020462e-03 1.279021054506301880e-02 -8.277410641312599182e-03 1.680318824946880341e-02 8.171995170414447784e-03 -1.793961599469184875e-02 -7.178520900197327137e-04 -1.655224896967411041e-02 -8.576686494052410126e-03 2.648287825286388397e-02 2.245653048157691956e-02 -2.479314990341663361e-02 -2.599848248064517975e-02 1.871763262897729874e-03 -1.798075437545776367e-02 -4.216609522700309753e-04 3.547094296663999557e-03 -4.467549268156290054e-03 1.272541936486959457e-02 -2.251772442832589149e-03 -1.833505183458328247e-02 4.828532109968364239e-04 -8.190507069230079651e-03 -1.827283203601837158e-02 1.298047229647636414e-02 -8.585428819060325623e-03 -1.525641325861215591e-02 1.295270491391420364e-02 1.116835884749889374e-02 -6.951749324798583984e-03 -2.566835284233093262e-02 2.900975942611694336e-02 -1.581999845802783966e-02 -3.503899788483977318e-03 2.792218292597681284e-04 -7.236627396196126938e-03 -1.866621081717312336e-03 2.224518684670329094e-03 5.501647014170885086e-03 -3.064292483031749725e-02 4.510229453444480896e-03 -9.799503721296787262e-03 4.271631129086017609e-03 -1.547238603234291077e-02 -2.558453567326068878e-02 -1.630285196006298065e-02 -1.694099977612495422e-02 7.224800530821084976e-03 9.183467365801334381e-03 -7.137896958738565445e-03 -1.933829672634601593e-02 6.450743973255157471e-02 -3.525239974260330200e-02 -9.753017872571945190e-02 -3.070720098912715912e-02 -5.610661581158638000e-03 -3.235116600990295410e-02 -1.719913817942142487e-02 -5.414743791334331036e-04 8.389510214328765869e-03 6.940179318189620972e-02 3.834697976708412170e-02 -5.660695582628250122e-02 -8.162738755345344543e-03 -1.202093716710805893e-02 -3.740955376997590065e-03 1.456787157803773880e-02 5.491877091117203236e-04 8.945383131504058838e-03 3.462864831089973450e-02 7.256065309047698975e-02 8.536063134670257568e-03 7.866033352911472321e-03 -5.486220587044954300e-03 1.132776681333780289e-02 1.698522269725799561e-02 -1.576600596308708191e-02 -1.265805587172508240e-02 -2.978571876883506775e-02 9.144900366663932800e-03 2.029114868491888046e-03 -6.899709999561309814e-02 1.742369979619979858e-01 -2.683743536472320557e-01 -4.429239779710769653e-02 6.901355460286140442e-03 2.609207108616828918e-02 2.234186045825481415e-02 2.111544273793697357e-02 -5.753100849688053131e-03 -8.917614817619323730e-03 -4.039614647626876831e-02 4.068822413682937622e-02 5.109296645969152451e-03 -1.294293440878391266e-02 3.259593620896339417e-02 -1.423830837011337280e-01 7.929539680480957031e-02 -4.696274548768997192e-02 -2.426093211397528648e-03 -1.106515452265739441e-01 -9.850925207138061523e-02 1.005635634064674377e-01 -3.671586886048316956e-02 -9.566509723663330078e-02 1.548637300729751587e-01 2.369528869166970253e-03 2.690099179744720459e-01 -9.338350594043731689e-02 1.726326197385787964e-01 -1.000776235014200211e-02 -5.573894828557968140e-02 -2.007779777050018311e-01 -3.753502294421195984e-02 -2.343843579292297363e-01 1.200124174356460571e-01 4.204010963439941406e-02 -3.067892044782638550e-02 -5.148377269506454468e-03 1.371575593948364258e-01 -3.283196315169334412e-02 -1.383109688758850098e-01 1.536985039710998535e-01 2.477211803197860718e-01 2.673551440238952637e-01 -5.387232899665832520e-01 -3.399512916803359985e-02 -3.363142013549804688e-01 -6.386671960353851318e-02 -2.391789555549621582e-01 -3.341148793697357178e-02 -5.645526051521301270e-01 -2.335436083376407623e-02 2.483331114053726196e-01 1.459018886089324951e-01 9.233356714248657227e-01 4.329774081707000732e-01 -4.229254126548767090e-01 -2.839886583387851715e-02 -3.265270292758941650e-01 6.476198136806488037e-02 -6.324636340141296387e-01 -1.054913401603698730e-01 2.155289053916931152e-01 -2.199296057224273682e-01 -4.232226312160491943e-01 2.678050994873046875e-01 3.901438415050506592e-01 5.183776617050170898e-01 2.177308052778244019e-01 -3.374686241149902344e-01 5.977681875228881836e-01 -7.905464172363281250e-01 -1.668029092252254486e-02 -1.058150410652160645e+00 -3.875943422317504883e-01 -3.925251960754394531e-02 4.368921741843223572e-02 1.205425858497619629e-01 -7.931236922740936279e-02 1.952598802745342255e-02 2.281134389340877533e-02 -1.511140260845422745e-02 -2.594140544533729553e-02 -2.723168730735778809e-01 -3.352718777023255825e-04 6.987436115741729736e-02 -2.582673355937004089e-02 5.285777151584625244e-02 6.644759327173233032e-03 2.587023675441741943e-01 -1.801109313964843750e-02 -3.912602365016937256e-02 2.104566395282745361e-01 -4.200668632984161377e-02 1.420020405203104019e-02 2.380420118570327759e-01 2.954866178333759308e-02 -2.349371463060379028e-02 2.397044748067855835e-01 -5.117009859532117844e-03 -2.002118155360221863e-02 1.887647509574890137e-01 4.477936401963233948e-02 2.717604860663414001e-02 -9.436257928609848022e-02 -8.993719704449176788e-03 -3.621656214818358421e-03 -2.380774915218353271e-01 1.019843593239784241e-01 4.974570591002702713e-03 2.735705487430095673e-02 6.167408451437950134e-02 1.060477271676063538e-02 -3.488719761371612549e-01 8.042410947382450104e-03 -8.855651319026947021e-02 1.035049110651016235e-01 -1.011960580945014954e-01 -6.915342807769775391e-02 2.238952815532684326e-01 7.932367734611034393e-03 -1.187361255288124084e-01 -9.918288141489028931e-02 6.672424823045730591e-02 -6.066209822893142700e-02 7.064929604530334473e-02 -5.082299932837486267e-02 -1.655636355280876160e-02 2.118203490972518921e-01 4.289011657238006592e-02 3.366690278053283691e-01 -7.251539826393127441e-02 9.767107665538787842e-02 2.942093014717102051e-01 -1.480986326932907104e-01 6.034225225448608398e-02 -4.308999478816986084e-01 1.146930694580078125e+00 2.093196958303451538e-01 -8.582703024148941040e-02 -2.207107096910476685e-01 -1.496844924986362457e-02 1.843213737010955811e-01 -3.934045135974884033e-01 2.001809626817703247e-01 -8.990696072578430176e-01 8.304767012596130371e-01 -7.793456315994262695e-01 -3.349325060844421387e-01 2.378815710544586182e-01 3.620494008064270020e-01 1.195638179779052734e+00 4.119195640087127686e-01 2.902950644493103027e-01 -6.464576721191406250e-02 -3.813990056514739990e-01 -1.106752008199691772e-01 1.081159263849258423e-01 -3.471555188298225403e-02 9.297624230384826660e-03 -1.349829253740608692e-03 -2.145842881873250008e-03 1.335945329628884792e-03 4.220723174512386322e-03 -1.248779743909835815e-01 -3.870024811476469040e-03 4.961200989782810211e-03 -7.309852540493011475e-02 -4.484004166442900896e-04 -4.894866608083248138e-03 1.382533740252256393e-02 3.704959992319345474e-03 -4.498402588069438934e-03 -1.058152988553047180e-01 -1.672759884968400002e-03 -4.989483393728733063e-03 -8.488652110099792480e-02 -1.313319546170532703e-03 -7.889173924922943115e-03 -9.988170862197875977e-02 -9.862752631306648254e-03 3.771294374018907547e-03 -5.090116336941719055e-02 6.931694806553423405e-04 -1.896960078738629818e-03 -8.603683114051818848e-02 5.891426466405391693e-03 1.052301190793514252e-02 -6.807099282741546631e-02 -1.224050857126712799e-02 -8.844595868140459061e-05 -1.171086169779300690e-02 9.057482704520225525e-03 -6.648873095400631428e-04 -7.828290015459060669e-02 -4.207353573292493820e-03 1.285037491470575333e-03 -5.484163388609886169e-02 4.715403541922569275e-03 -1.758480444550514221e-02 -8.468031883239746094e-02 -9.094506502151489258e-03 -5.433769896626472473e-03 -1.080429032444953918e-01 -8.356085745617747307e-04 -8.645592257380485535e-03 -8.232724666595458984e-02 -2.507082931697368622e-02 4.364705085754394531e-02 1.537128984928131104e-01 5.729601159691810608e-02 5.963777005672454834e-02 1.158852875232696533e-01 -6.959977000951766968e-02 1.238217651844024658e-01 1.664686650037765503e-01 2.036262489855289459e-02 9.312793612480163574e-02 1.595683246850967407e-01 8.057529479265213013e-02 5.345786735415458679e-02 2.299466282129287720e-01 5.173162743449211121e-02 -4.056172911077737808e-03 2.901870608329772949e-01 1.976664364337921143e-02 -3.087278082966804504e-02 1.369883865118026733e-01 -1.363353244960308075e-02 4.292251542210578918e-02 1.564999818801879883e-01 1.284153014421463013e-01 -6.114495173096656799e-02 3.792459666728973389e-01 -4.886947572231292725e-02 6.220578774809837341e-02 2.492766976356506348e-01 4.413247480988502502e-02 1.222508251667022705e-01 2.082742899656295776e-01 -1.421974413096904755e-02 -4.172386601567268372e-03 1.484586000442504883e-01 2.711948379874229431e-02 -2.002323977649211884e-02 2.530521154403686523e-01 -3.040535375475883484e-02 -8.523208089172840118e-03 2.680879235267639160e-01 -9.880236349999904633e-03 3.497626259922981262e-02 1.259254515171051025e-01 3.049740009009838104e-02 6.191216874867677689e-03 1.380322575569152832e-01 -3.762089833617210388e-02 4.634927958250045776e-02 1.610230654478073120e-02 2.170731127262115479e-02 -7.616019248962402344e-02 -1.117401523515582085e-03 2.092878130497410893e-04 5.715562961995601654e-03 -1.721148043870925903e-01 -1.635478809475898743e-02 2.458540536463260651e-02 8.387093245983123779e-02 -1.793704554438591003e-02 -6.610924005508422852e-02 4.725636541843414307e-02 -4.781737923622131348e-03 4.136905539780855179e-03 -1.320595890283584595e-01 9.640094824135303497e-03 1.156122330576181412e-02 2.084830552339553833e-01 -3.595945192500948906e-03 2.512173727154731750e-02 -1.716694831848144531e-01 5.623003002256155014e-03 1.370912324637174606e-02 -1.162035912275314331e-01 3.378944098949432373e-02 -3.141476958990097046e-02 4.943100363016128540e-02 9.966830909252166748e-02 3.148676082491874695e-02 3.660737071186304092e-03 6.501960754394531250e-02 -1.878396794199943542e-02 1.421563047915697098e-02 1.706852763891220093e-02 -9.174210950732231140e-03 -1.557013988494873047e-01 -1.136684790253639221e-02 4.417669959366321564e-03 7.971449196338653564e-02 4.406357184052467346e-03 -2.304119989275932312e-02 1.649928465485572815e-02 -4.582734778523445129e-03 -3.686448093503713608e-03 -7.609026134014129639e-02 +-1.104494556784629822e-02 3.476131707429885864e-02 2.383937388658523560e-01 -2.983302623033523560e-02 6.788226962089538574e-02 -7.082361262291669846e-03 -6.491460204124450684e-01 3.486107662320137024e-02 1.303477436304092407e-01 -7.131726294755935669e-02 1.783661842346191406e-01 -2.583156824111938477e-01 -1.776951253414154053e-01 1.444620639085769653e-01 -1.592680215835571289e-01 8.508963882923126221e-02 8.697991073131561279e-02 -1.015986949205398560e-01 7.329495996236801147e-02 3.711572885513305664e-01 1.321593374013900757e-01 5.504115298390388489e-02 -6.570917367935180664e-02 2.412003464996814728e-02 2.140500098466873169e-01 -7.536160945892333984e-02 1.978813707828521729e-01 -9.274689853191375732e-02 -5.994388461112976074e-01 -1.937820017337799072e-01 -1.141302511096000671e-01 3.028307259082794189e-01 -1.305004656314849854e-01 1.491461098194122314e-01 5.856459215283393860e-02 1.309930831193923950e-01 3.258333504199981689e-01 -8.520996570587158203e-02 -2.395059913396835327e-01 -2.042067945003509521e-01 -5.118007659912109375e-01 -4.316214099526405334e-02 -7.233747839927673340e-02 1.305323541164398193e-01 -1.028775572776794434e-01 -1.012449059635400772e-02 3.805541619658470154e-02 -1.114979535341262817e-01 9.408850222826004028e-02 -1.270032525062561035e-01 -9.719439595937728882e-02 -1.077755093574523926e-01 -1.132776681333780289e-02 -1.396109163761138916e-02 1.897062957286834717e-01 -1.100956723093986511e-01 -1.400613635778427124e-01 -2.287099324166774750e-02 -2.160446718335151672e-02 -5.284748598933219910e-02 3.776950761675834656e-02 -1.654561609029769897e-01 1.634496748447418213e-01 -1.084990203380584717e-01 8.678759634494781494e-02 2.098215818405151367e-01 6.759327650070190430e-01 -4.112495183944702148e-01 -3.361912965774536133e-01 1.149864271283149719e-01 -1.030446812510490417e-01 -1.138047501444816589e-01 -6.231747865676879883e-01 -3.622010946273803711e-01 -1.078145951032638550e-01 1.205538958311080933e-01 1.539227068424224854e-01 2.042021788656711578e-02 3.630398074164986610e-03 -1.659852825105190277e-02 1.282276064157485962e-01 3.265661299228668213e-01 2.516092061996459961e-01 6.139280647039413452e-02 3.997752070426940918e-01 -3.187895715236663818e-01 6.818772107362747192e-02 6.308767944574356079e-02 -3.385834768414497375e-02 -2.106453664600849152e-02 -4.315339960157871246e-03 -2.343010455369949341e-01 -1.357208341360092163e-01 -7.493481040000915527e-02 1.520823091268539429e-01 -2.871084213256835938e-01 -1.007563993334770203e-01 3.371899202466011047e-02 1.746123731136322021e-01 -4.844524338841438293e-02 2.442136779427528381e-02 1.908154636621475220e-01 5.124671757221221924e-02 1.412152778357267380e-02 3.879537805914878845e-02 2.635992765426635742e-01 -4.273878335952758789e-01 -4.049179330468177795e-02 -1.563117653131484985e-01 3.645001724362373352e-02 -9.467677026987075806e-02 -7.148798555135726929e-02 -2.167239636182785034e-01 -1.213452816009521484e-01 2.505319118499755859e-01 -1.560865435749292374e-02 -3.080012202262878418e-01 2.920737564563751221e-01 -3.169656265527009964e-03 6.070308089256286621e-01 -2.779676318168640137e-01 4.885317385196685791e-01 2.117823064327239990e-01 8.048170804977416992e-02 -1.967896819114685059e-01 7.976076751947402954e-02 -3.200874626636505127e-01 2.314281016588211060e-01 5.967788025736808777e-02 2.285798341035842896e-01 5.649356842041015625e-01 -5.458349734544754028e-02 -1.947559416294097900e-01 6.182063743472099304e-02 2.359594106674194336e-01 -1.382492631673812866e-01 -6.764315813779830933e-02 -2.699308991432189941e-01 -7.123859226703643799e-02 -6.276844739913940430e-01 -2.283124327659606934e-01 5.317710340023040771e-02 6.410789489746093750e-02 -3.124867677688598633e-01 -3.441524878144264221e-02 2.965572476387023926e-01 1.769356280565261841e-01 6.501158475875854492e-01 3.025216758251190186e-01 -3.078834712505340576e-01 1.825997605919837952e-03 7.202689349651336670e-02 -2.511819005012512207e-01 -9.704058766365051270e-01 -2.338984161615371704e-01 1.342753618955612183e-01 -2.212151587009429932e-01 -3.379782140254974365e-01 2.369971126317977905e-01 7.971649169921875000e-01 8.250259160995483398e-01 4.743896424770355225e-01 -3.288477063179016113e-01 5.826305150985717773e-01 -6.866188049316406250e-01 3.443633019924163818e-02 -9.923687577247619629e-01 -3.228812217712402344e-01 7.178880274295806885e-02 -1.244876831769943237e-01 2.564778551459312439e-02 -1.374753564596176147e-01 -9.046375751495361328e-02 1.591621041297912598e-01 3.975152075290679932e-01 2.328185550868511200e-02 -6.311549544334411621e-01 3.031587973237037659e-02 -5.235481262207031250e-01 1.079878807067871094e-01 -1.716468632221221924e-01 7.066404819488525391e-01 6.992629766464233398e-01 9.035885334014892578e-02 -2.569920718669891357e-01 2.911028861999511719e-01 -2.731956839561462402e-01 1.626562327146530151e-01 -1.209534481167793274e-01 2.627821825444698334e-02 -3.783430159091949463e-02 4.602383077144622803e-01 3.054902851581573486e-01 5.480574369430541992e-01 1.684962362051010132e-01 -1.160549819469451904e-01 1.016665697097778320e-01 1.942725735716521740e-03 -2.473350055515766144e-02 -5.631847679615020752e-02 -2.918592989444732666e-01 -5.663930177688598633e-01 3.910900056362152100e-01 1.680426895618438721e-01 -1.523800343275070190e-01 -1.604018807411193848e-01 -3.901119828224182129e-01 1.903650015592575073e-01 -1.977188736200332642e-01 -1.694372594356536865e-01 4.718992710113525391e-01 4.436778128147125244e-01 3.941043913364410400e-01 -6.879305839538574219e-01 2.727374993264675140e-02 -1.745609492063522339e-01 9.348120540380477905e-02 1.866610795259475708e-01 2.327002882957458496e-01 -4.432057738304138184e-01 3.967454135417938232e-01 5.185771584510803223e-01 9.205938130617141724e-02 4.654653668403625488e-01 -5.490437150001525879e-02 7.019934803247451782e-02 3.512595295906066895e-01 -4.940853416919708252e-01 1.627919822931289673e-01 -1.851636618375778198e-01 6.728850007057189941e-01 3.807649761438369751e-02 -7.737895250320434570e-01 4.971485584974288940e-03 -1.257660388946533203e-01 4.114834964275360107e-01 -5.589095354080200195e-01 3.938092291355133057e-01 -1.208497285842895508e+00 5.434006452560424805e-01 -3.033346533775329590e-01 -4.138216674327850342e-01 4.525141775957308710e-05 6.551947593688964844e-01 1.291819572448730469e+00 4.318991005420684814e-01 -1.062909588217735291e-01 -4.860830307006835938e-01 -5.503559708595275879e-01 -1.151345223188400269e-01 -3.370392620563507080e-01 2.496016919612884521e-01 7.138925604522228241e-03 4.087025765329599380e-03 -5.033191759139299393e-03 -2.961396938189864159e-03 2.248378470540046692e-02 -1.294812858104705811e-01 7.782215718179941177e-03 -3.985724411904811859e-03 -7.532304525375366211e-02 1.993222162127494812e-02 8.204391226172447205e-03 1.980006694793701172e-02 -5.281046614982187748e-04 6.842734292149543762e-03 -1.057104021310806274e-01 -9.847326204180717468e-03 1.671731355600059032e-03 -8.024259656667709351e-02 4.962229286320507526e-04 -7.765246555209159851e-03 -9.957729279994964600e-02 -1.425522565841674805e-02 8.225988596677780151e-03 -4.023519530892372131e-02 -5.335553549230098724e-03 8.743294514715671539e-03 -8.677422255277633667e-02 -6.039007566869258881e-03 1.079194899648427963e-02 -5.816144496202468872e-02 -1.022990606725215912e-02 9.163412614725530148e-04 -1.089119352400302887e-02 1.388190127909183502e-02 1.676359446719288826e-03 -8.718302845954895020e-02 -2.216291148215532303e-03 8.119029924273490906e-03 -5.274361371994018555e-02 1.222816761583089828e-03 -1.122080907225608826e-02 -8.391464501619338989e-02 -1.298869960010051727e-02 -1.625605858862400055e-02 -1.077976226806640625e-01 -1.185638643801212311e-02 -1.393898017704486847e-02 -8.446948975324630737e-02 -1.552997902035713196e-02 4.289371892809867859e-02 1.527400016784667969e-01 6.501240283250808716e-02 5.235178023576736450e-02 1.159608736634254456e-01 -5.968559533357620239e-02 1.221289560198783875e-01 1.672075986862182617e-01 3.620216250419616699e-02 1.012747362256050110e-01 1.585224121809005737e-01 6.760562211275100708e-02 5.228698626160621643e-02 2.168802767992019653e-01 5.048927292227745056e-02 -2.149133943021297455e-02 2.762388288974761963e-01 1.242562849074602127e-02 -2.042330242693424225e-02 1.278681755065917969e-01 -1.921539753675460815e-02 4.353494942188262939e-02 1.464042812585830688e-01 1.288631856441497803e-01 -6.012165173888206482e-02 3.738620877265930176e-01 -3.939598798751831055e-02 5.242068320512771606e-02 2.547222971916198730e-01 4.888850077986717224e-02 1.314548552036285400e-01 2.049519121646881104e-01 -2.355079352855682373e-02 -1.129537168890237808e-02 1.436465233564376831e-01 2.126508206129074097e-02 -2.505334466695785522e-02 2.450935095548629761e-01 -4.344187676906585693e-02 -5.161747336387634277e-03 2.797350287437438965e-01 -6.929637864232063293e-03 3.583141043782234192e-02 1.322745233774185181e-01 3.419773280620574951e-02 2.383412793278694153e-03 1.350004076957702637e-01 -5.564896017313003540e-02 2.226575464010238647e-02 1.167949382215738297e-02 2.564058639109134674e-02 -9.204909950494766235e-02 -1.253567170351743698e-02 -2.450312860310077667e-02 1.618663780391216278e-02 -1.755765378475189209e-01 -2.139569446444511414e-02 3.943044319748878479e-02 7.631909102201461792e-02 -1.542302034795284271e-02 -5.633390322327613831e-02 4.985009506344795227e-02 -4.008349962532520294e-03 1.968128234148025513e-02 -1.294668763875961304e-01 1.023813383653759956e-03 -3.601653128862380981e-02 2.150717675685882568e-01 1.322781294584274292e-02 1.532223355025053024e-02 -1.604856997728347778e-01 -1.355588529258966446e-02 2.606070414185523987e-02 -1.114218533039093018e-01 5.919913947582244873e-02 1.487126224674284458e-03 3.223289549350738525e-02 9.704526513814926147e-02 4.147806763648986816e-02 7.003685459494590759e-04 5.199490860104560852e-02 1.038211514241993427e-03 1.214074995368719101e-02 2.444193698465824127e-02 9.088850580155849457e-03 -1.592788249254226685e-01 1.247602235525846481e-02 -2.815512381494045258e-02 7.731667906045913696e-02 5.528386682271957397e-03 -1.595163904130458832e-02 8.866706863045692444e-03 -3.608955070376396179e-02 -1.909198611974716187e-02 -7.481499761343002319e-02 +-9.893091395497322083e-03 7.110026478767395020e-02 4.683018028736114502e-01 -6.107398867607116699e-02 1.347078233957290649e-01 -1.578657422214746475e-03 -1.261842489242553711e+00 4.915383830666542053e-02 2.711197733879089355e-01 -1.311180442571640015e-01 3.344486057758331299e-01 -4.808215200901031494e-01 -3.216167390346527100e-01 2.555764317512512207e-01 -2.926167547702789307e-01 1.411504894495010376e-01 1.599719971418380737e-01 -1.761498898267745972e-01 1.470336914062500000e-01 6.908518671989440918e-01 2.576050162315368652e-01 7.998599857091903687e-02 -1.543093919754028320e-01 7.660036534070968628e-02 4.400098919868469238e-01 -1.362982988357543945e-01 3.931432962417602539e-01 -1.867906600236892700e-01 -1.155184388160705566e+00 -3.807912170886993408e-01 -2.209133207798004150e-01 5.606990456581115723e-01 -2.357347011566162109e-01 2.812612354755401611e-01 1.144310683012008667e-01 2.691909372806549072e-01 5.828825235366821289e-01 -1.458211541175842285e-01 -4.240993857383728027e-01 -4.254646301269531250e-01 -1.020523428916931152e+00 -8.228096365928649902e-02 -1.227393299341201782e-01 2.231589108705520630e-01 -1.860954314470291138e-01 -9.655521251261234283e-03 7.481139153242111206e-02 -2.035866528749465942e-01 1.776267290115356445e-01 -2.375889718532562256e-01 -1.883477270603179932e-01 -1.761766374111175537e-01 -2.894805371761322021e-02 -1.338002178817987442e-02 3.621239662170410156e-01 -1.977127045392990112e-01 -2.472342103719711304e-01 -2.230175025761127472e-02 -1.701401919126510620e-02 -1.129578277468681335e-01 5.937706306576728821e-02 -3.134961724281311035e-01 3.318847119808197021e-01 -2.587897777557373047e-01 1.865680068731307983e-01 4.676410257816314697e-01 1.349117040634155273e+00 -8.273280858993530273e-01 -6.393084526062011719e-01 2.348702996969223022e-01 -1.994476914405822754e-01 -2.319346070289611816e-01 -1.305436134338378906e+00 -7.539457082748413086e-01 -1.434228271245956421e-01 2.326097786426544189e-01 3.093150556087493896e-01 3.671895712614059448e-02 -3.112577833235263824e-03 -3.046294674277305603e-02 2.379859536886215210e-01 5.817291140556335449e-01 3.898522853851318359e-01 1.008556410670280457e-01 7.762392759323120117e-01 -6.269609928131103516e-01 1.192472577095031738e-01 9.827579557895660400e-02 -4.342542216181755066e-02 -3.569976985454559326e-02 1.863998547196388245e-02 -4.682920277118682861e-01 -2.610754966735839844e-01 -7.949800044298171997e-02 1.178094968199729919e-01 -2.877820432186126709e-01 -1.527363955974578857e-01 5.917343124747276306e-02 3.124430477619171143e-01 -1.266787797212600708e-01 2.466356568038463593e-02 3.761935532093048096e-01 1.092194393277168274e-01 6.830959022045135498e-02 3.851101547479629517e-02 5.131099820137023926e-01 -8.361958265304565430e-01 -1.178449839353561401e-01 -1.629236191511154175e-01 -8.984463289380073547e-03 -1.315494775772094727e-01 -1.241508722305297852e-01 -3.026718199253082275e-01 -1.452256888151168823e-01 3.789071142673492432e-01 -2.524823416024446487e-03 -4.800285995006561279e-01 3.944936394691467285e-01 -1.162138767540454865e-03 8.724946975708007812e-01 -4.213493466377258301e-01 7.462730407714843750e-01 3.899875283241271973e-01 2.097485512495040894e-01 -1.942535489797592163e-01 1.942556053400039673e-01 -3.750139474868774414e-01 3.149149119853973389e-01 6.796608865261077881e-02 4.504814743995666504e-01 1.019416213035583496e+00 -2.055576592683792114e-01 -3.281715214252471924e-01 2.426273077726364136e-01 3.063994050025939941e-01 -5.537287592887878418e-01 -4.357289969921112061e-01 1.795658655464649200e-02 -1.048259437084197998e-01 -8.470063209533691406e-01 -3.549505770206451416e-01 3.670609891414642334e-01 1.665252298116683960e-01 -2.213977091014385223e-02 -4.017966240644454956e-02 3.127880990505218506e-01 1.877275705337524414e-01 2.954964041709899902e-01 1.175163909792900085e-01 -1.510980874300003052e-01 4.575741291046142578e-02 4.829622209072113037e-01 -5.672903060913085938e-01 -1.197457432746887207e+00 -3.361995220184326172e-01 6.976174563169479370e-02 -2.216095775365829468e-01 -2.158086448907852173e-01 1.900338530540466309e-01 1.096209406852722168e+00 1.025184273719787598e+00 6.810112595558166504e-01 -3.208042681217193604e-01 5.260307788848876953e-01 -5.059360861778259277e-01 9.677067399024963379e-02 -8.423340916633605957e-01 -1.967562586069107056e-01 1.529091745615005493e-01 -2.647413611412048340e-01 -5.476912856101989746e-02 -1.835166215896606445e-01 -1.947086304426193237e-01 2.880592346191406250e-01 7.779361605644226074e-01 7.493275403976440430e-02 -9.699821472167968750e-01 4.461738467216491699e-02 -1.026048183441162109e+00 2.282728403806686401e-01 -3.384256064891815186e-01 1.343115091323852539e+00 1.096434712409973145e+00 1.817189007997512817e-01 -4.267116189002990723e-01 3.620890080928802490e-01 -4.821353554725646973e-01 3.120321929454803467e-01 -4.564875662326812744e-01 3.491455689072608948e-02 -5.989950895309448242e-02 6.529316902160644531e-01 6.110998392105102539e-01 1.099410533905029297e+00 1.527127474546432495e-01 -2.757678031921386719e-01 1.618555784225463867e-01 8.699534088373184204e-02 -5.625214055180549622e-02 -1.119247600436210632e-01 -3.362915813922882080e-01 -1.203253269195556641e+00 7.881100177764892578e-01 3.196735084056854248e-01 -3.434947729110717773e-01 -3.287901282310485840e-01 -3.869171142578125000e-01 3.485536873340606689e-01 -2.719147503376007080e-01 -4.036462604999542236e-01 1.044531345367431641e+00 9.418208599090576172e-01 5.448548793792724609e-01 -1.327193737030029297e+00 2.001449763774871826e-01 -2.329656183719635010e-01 1.114362478256225586e-01 4.299717843532562256e-01 3.690345883369445801e-01 -8.036394715309143066e-01 8.094383478164672852e-01 7.861652374267578125e-01 1.207611262798309326e-01 5.581639409065246582e-01 -3.813306242227554321e-02 5.217077210545539856e-02 3.965011537075042725e-01 -8.052032589912414551e-01 2.708004415035247803e-01 9.015934169292449951e-02 1.360761672258377075e-01 -1.289413422346115112e-01 -1.314676165580749512e+00 2.050768584012985229e-01 -2.040916085243225098e-01 5.657923817634582520e-01 -6.838780045509338379e-01 5.440202951431274414e-01 -1.415313243865966797e+00 2.235332727432250977e-01 2.245771437883377075e-01 -4.179128110408782959e-01 -2.605196237564086914e-01 8.231500387191772461e-01 1.273506045341491699e+00 3.739458918571472168e-01 -3.938071727752685547e-01 -8.364529609680175781e-01 -6.571555137634277344e-01 -1.396422833204269409e-01 -7.621161937713623047e-01 5.516410470008850098e-01 4.955544602125883102e-03 9.621582925319671631e-03 -6.986202206462621689e-03 -7.330730091780424118e-03 4.049744829535484314e-02 -1.321984231472015381e-01 1.939486153423786163e-02 -1.333425659686326981e-02 -7.554673403501510620e-02 3.975645825266838074e-02 2.119052037596702576e-02 2.583701722323894501e-02 -4.761169198900461197e-03 1.824557781219482422e-02 -1.068581417202949524e-01 -1.780694723129272461e-02 7.972992025315761566e-03 -7.572825253009796143e-02 2.534079598262906075e-03 -7.397064473479986191e-03 -1.013287231326103210e-01 -1.824506372213363647e-02 1.250841747969388962e-02 -2.939645387232303619e-02 -1.128045935183763504e-02 1.917631737887859344e-02 -8.776050060987472534e-02 -1.817410066723823547e-02 1.099866628646850586e-02 -4.806677624583244324e-02 -8.188963867723941803e-03 1.888732542283833027e-03 -1.017077080905437469e-02 1.872277446091175079e-02 3.620113478973507881e-03 -9.292636066675186157e-02 -3.630397841334342957e-04 1.480698399245738983e-02 -4.907670617103576660e-02 -2.207035198807716370e-03 -4.811048507690429688e-03 -8.566196262836456299e-02 -1.708395406603813171e-02 -2.702949568629264832e-02 -1.082830503582954407e-01 -2.264319360256195068e-02 -1.930178701877593994e-02 -8.770547807216644287e-02 -5.845660809427499771e-03 4.200720041990280151e-02 1.513104587793350220e-01 7.280850410461425781e-02 4.511515051126480103e-02 1.158076375722885132e-01 -4.996528103947639465e-02 1.204417943954467773e-01 1.678514033555984497e-01 5.110633745789527893e-02 1.099300980567932129e-01 1.553280651569366455e-01 5.640126764774322510e-02 5.075923725962638855e-02 2.010067999362945557e-01 4.866635799407958984e-02 -3.678374737501144409e-02 2.593600451946258545e-01 5.267162341624498367e-03 -9.717227891087532043e-03 1.177632212638854980e-01 -2.470110356807708740e-02 4.407642409205436707e-02 1.362098604440689087e-01 1.288472414016723633e-01 -5.881450325250625610e-02 3.665066659450531006e-01 -2.970858477056026459e-02 4.231470450758934021e-02 2.592505216598510742e-01 5.263562873005867004e-02 1.400428563356399536e-01 2.000724673271179199e-01 -3.160451725125312805e-02 -1.877728290855884552e-02 1.365667283535003662e-01 1.508774887770414352e-02 -2.999140694737434387e-02 2.362376004457473755e-01 -5.841958150267601013e-02 -2.730511769186705351e-04 2.864455878734588623e-01 -3.346033859997987747e-03 3.664439544081687927e-02 1.374090164899826050e-01 3.791400417685508728e-02 -1.554489135742187500e-03 1.316636353731155396e-01 -7.310777902603149414e-02 -2.165897516533732414e-03 7.529733236879110336e-03 2.937588468194007874e-02 -1.078212782740592957e-01 -2.507854253053665161e-02 -4.747645184397697449e-02 2.562515810132026672e-02 -1.774123013019561768e-01 -2.635998092591762543e-02 5.408058688044548035e-02 6.827719509601593018e-02 -1.295218989253044128e-02 -4.677094146609306335e-02 5.198719725012779236e-02 -3.239590208977460861e-03 3.646956011652946472e-02 -1.268860101699829102e-01 -7.840836420655250549e-03 -8.373003453016281128e-02 2.224750071763992310e-01 3.069074638187885284e-02 5.268705077469348907e-03 -1.480256170034408569e-01 -3.176289796829223633e-02 3.841331228613853455e-02 -1.046835035085678101e-01 8.282604068517684937e-02 3.443067520856857300e-02 1.815867424011230469e-02 9.441194683313369751e-02 5.159073323011398315e-02 -1.701041939668357372e-03 3.985930234193801880e-02 1.937892101705074310e-02 9.308422915637493134e-03 3.155926614999771118e-02 2.736631035804748535e-02 -1.607952564954757690e-01 3.609520569443702698e-02 -6.125190854072570801e-02 7.475946098566055298e-02 6.751203443855047226e-03 -8.983435109257698059e-03 1.465528970584273338e-03 -6.945424526929855347e-02 -3.552904725074768066e-02 -7.128743827342987061e-02 diff --git a/tests/cp2k/restart_aimd/deepmd/energy.raw b/tests/cp2k/restart_aimd/deepmd/energy.raw index 852484ea9..e97cf5960 100644 --- a/tests/cp2k/restart_aimd/deepmd/energy.raw +++ b/tests/cp2k/restart_aimd/deepmd/energy.raw @@ -1,5 +1,5 @@ --2.196387023925781250e+02 --2.196259307861328125e+02 --2.196165771484375000e+02 --2.196100006103515625e+02 --2.196051788330078125e+02 +-2.196387329101562500e+02 +-2.196259460449218750e+02 +-2.196165924072265625e+02 +-2.196100158691406250e+02 +-2.196051940917968750e+02 diff --git a/tests/cp2k/restart_aimd/deepmd/force.raw b/tests/cp2k/restart_aimd/deepmd/force.raw index 2d03f2c87..d08bd63a2 100644 --- a/tests/cp2k/restart_aimd/deepmd/force.raw +++ b/tests/cp2k/restart_aimd/deepmd/force.raw @@ -1,5 +1,5 @@ --2.169162631034851074e-01 3.279730081558227539e-01 6.035397052764892578e-01 7.185621261596679688e-01 1.594181656837463379e-01 4.608764052391052246e-01 -3.765185177326202393e-01 3.857482671737670898e-01 -1.248831078410148621e-01 -6.494076848030090332e-01 1.203697919845581055e-01 -4.328375160694122314e-01 5.217282772064208984e-01 -9.955028891563415527e-01 -5.069978833198547363e-01 --1.989514529705047607e-01 2.921981513500213623e-01 5.636984109878540039e-01 6.776738762855529785e-01 1.945625692605972290e-01 5.051451325416564941e-01 -4.636886715888977051e-01 3.652833104133605957e-01 -1.360401511192321777e-01 -7.025374770164489746e-01 9.427926689386367798e-02 -3.964312076568603516e-01 6.848884224891662598e-01 -9.483894705772399902e-01 -5.365054011344909668e-01 --1.648406237363815308e-01 2.307343930006027222e-01 4.815881848335266113e-01 5.865107774734497070e-01 2.308490723371505737e-01 5.361675620079040527e-01 -5.500597357749938965e-01 3.308701217174530029e-01 -1.434295028448104858e-01 -7.340555787086486816e-01 5.354630947113037109e-02 -3.259690999984741211e-01 8.598077297210693359e-01 -8.481215834617614746e-01 -5.483376383781433105e-01 --1.165429651737213135e-01 1.469107717275619507e-01 3.615300059318542480e-01 4.501595497131347656e-01 2.673597633838653564e-01 5.535713434219360352e-01 -6.303640604019165039e-01 2.890907227993011475e-01 -1.467755436897277832e-01 -7.438977956771850586e-01 -1.830625406000763178e-04 -2.249499261379241943e-01 1.038020730018615723e+00 -7.053260803222656250e-01 -5.432190895080566406e-01 --5.708466097712516785e-02 4.539745301008224487e-02 2.101059257984161377e-01 2.773819267749786377e-01 3.029201626777648926e-01 5.577098131179809570e-01 -6.990006566047668457e-01 2.466197013854980469e-01 -1.458370834589004517e-01 -7.330569624900817871e-01 -6.438967585563659668e-02 -9.865168482065200806e-02 1.209161996841430664e+00 -5.326966047286987305e-01 -5.230426192283630371e-01 +-2.169162780046463013e-01 3.279730379581451416e-01 6.035397648811340332e-01 7.185621857643127441e-01 1.594181805849075317e-01 4.608764648437500000e-01 -3.765185475349426270e-01 3.857482969760894775e-01 -1.248831152915954590e-01 -6.494077444076538086e-01 1.203698068857192993e-01 -4.328375458717346191e-01 5.217282772064208984e-01 -9.955029487609863281e-01 -5.069979429244995117e-01 +-1.989514678716659546e-01 2.921981811523437500e-01 5.636984705924987793e-01 6.776739358901977539e-01 1.945625841617584229e-01 5.051451921463012695e-01 -4.636887311935424805e-01 3.652833402156829834e-01 -1.360401660203933716e-01 -7.025375366210937500e-01 9.427927434444427490e-02 -3.964312374591827393e-01 6.848884820938110352e-01 -9.483895301818847656e-01 -5.365054607391357422e-01 +-1.648406386375427246e-01 2.307344228029251099e-01 4.815882444381713867e-01 5.865108370780944824e-01 2.308490872383117676e-01 5.361676216125488281e-01 -5.500597953796386719e-01 3.308701515197753906e-01 -1.434295177459716797e-01 -7.340556979179382324e-01 5.354631319642066956e-02 -3.259690999984741211e-01 8.598078489303588867e-01 -8.481216430664062500e-01 -5.483376979827880859e-01 +-1.165429726243019104e-01 1.469107866287231445e-01 3.615300357341766357e-01 4.501595795154571533e-01 2.673597931861877441e-01 5.535714030265808105e-01 -6.303641200065612793e-01 2.890907526016235352e-01 -1.467755436897277832e-01 -7.438978552818298340e-01 -1.830625551519915462e-04 -2.249499410390853882e-01 1.038020730018615723e+00 -7.053261399269104004e-01 -5.432191491127014160e-01 +-5.708466470241546631e-02 4.539745673537254333e-02 2.101059406995773315e-01 2.773819565773010254e-01 3.029201924800872803e-01 5.577098727226257324e-01 -6.990007162094116211e-01 2.466197162866592407e-01 -1.458370983600616455e-01 -7.330570816993713379e-01 -6.438968330621719360e-02 -9.865169227123260498e-02 1.209162116050720215e+00 -5.326966047286987305e-01 -5.230426192283630371e-01 diff --git a/tests/test_abacus_pw_scf.py b/tests/test_abacus_pw_scf.py index 8c1dfd29b..3813943f2 100644 --- a/tests/test_abacus_pw_scf.py +++ b/tests/test_abacus_pw_scf.py @@ -2,6 +2,9 @@ import numpy as np import unittest from context import dpdata +from dpdata.unit import LengthConversion + +bohr2ang = LengthConversion("bohr", "angstrom").value() class TestABACUSSinglePointEnergy: @@ -25,7 +28,8 @@ def test_atom_types(self) : # self.assertEqual(self.system_h2o.data['atom_types'][ii], ref_type[ii]) def test_cell(self) : - cell = 5.29177 * np.eye(3) + # cell = 5.29177 * np.eye(3) + cell = bohr2ang * 10 * np.eye(3) for ii in range(cell.shape[0]) : for jj in range(cell.shape[1]) : self.assertAlmostEqual(self.system_ch4.data['cells'][0][ii][jj], cell[ii][jj]) @@ -49,7 +53,7 @@ def test_coord(self) : coord = np.array(coord) for ii in range(coord.shape[0]) : for jj in range(coord.shape[1]) : - self.assertAlmostEqual(self.system_ch4.data['coords'][0][ii][jj], coord[ii][jj]) + self.assertAlmostEqual(self.system_ch4.data['coords'][0][ii][jj], coord[ii][jj], places=5) fp.close() # fp = open('qe.scf/h2o_coord') diff --git a/tests/test_cp2k_output.py b/tests/test_cp2k_output.py index d19750e0e..30dae7885 100644 --- a/tests/test_cp2k_output.py +++ b/tests/test_cp2k_output.py @@ -16,14 +16,14 @@ def test_cell(self): cell = np.loadtxt('cp2k/ref_cell') for ii in range(cell.shape[0]) : for jj in range(cell.shape[1]) : - self.assertEqual(self.system.data['cells'][0][ii][jj], cell[ii][jj]) + self.assertAlmostEqual(self.system.data['cells'][0][ii][jj], cell[ii][jj]) def test_coord(self): coord = np.loadtxt('cp2k/ref_coord') for ii in range(coord.shape[0]) : for jj in range(coord.shape[1]) : - self.assertEqual(self.system.data['coords'][0][ii][jj], coord[ii][jj]) + self.assertAlmostEqual(self.system.data['coords'][0][ii][jj], coord[ii][jj]) def test_force(self): #eV = 2.72113838565563E+01 @@ -31,18 +31,18 @@ def test_force(self): force = np.loadtxt('cp2k/ref_force') for ii in range(force.shape[0]) : for jj in range(force.shape[1]) : - self.assertEqual(self.system.data['forces'][0][ii][jj], force[ii][jj]) + self.assertAlmostEqual(self.system.data['forces'][0][ii][jj], force[ii][jj], places=6) def test_energy(self): #eV = 2.72113838565563E+01 - ref_energy = -48061.44424374075 + ref_energy = -48061.44846401638 self.assertEqual(self.system.data['energies'][0], ref_energy) def test_virial(self): virial = np.loadtxt("cp2k/ref_virial") for ii in range(virial.shape[0]) : for jj in range(virial.shape[1]) : - self.assertEqual(self.system.data['virials'][0][ii][jj], virial[ii][jj]) + self.assertAlmostEqual(self.system.data['virials'][0][ii][jj], virial[ii][jj], places=6) diff --git a/tests/test_qe_cp_traj.py b/tests/test_qe_cp_traj.py index b98a47a4f..cad93fa6c 100644 --- a/tests/test_qe_cp_traj.py +++ b/tests/test_qe_cp_traj.py @@ -3,6 +3,7 @@ import unittest from context import dpdata +bohr2ang = dpdata.unit.LengthConversion("bohr", "angstrom").value() class TestCPTRAJProps : def test_atom_names(self) : @@ -18,7 +19,7 @@ def test_atom_types(self) : self.assertEqual(self.system.data['atom_types'][ii], 1) def test_cell(self) : - ref = 0.52917721067 * 23.5170 * np.eye(3) + ref = bohr2ang * 23.5170 * np.eye(3) self.assertEqual(self.system.get_nframes(), 2) for ff in range(self.system.get_nframes()) : for ii in range(3) : @@ -32,7 +33,6 @@ def test_coord(self) : coords = [] for ii in lines : coords.append([float(jj) for jj in ii.split()]) - bohr2ang = 0.52917721067 coords = bohr2ang * np.array(coords) celll = bohr2ang * 23.5170 for ii in range(coords.shape[0]) : diff --git a/tests/test_qe_cp_traj_skipload.py b/tests/test_qe_cp_traj_skipload.py index c87baa314..b0d92e9d0 100644 --- a/tests/test_qe_cp_traj_skipload.py +++ b/tests/test_qe_cp_traj_skipload.py @@ -6,6 +6,8 @@ from comp_sys import CompLabeledSys from comp_sys import IsPBC +bohr2ang = dpdata.unit.LengthConversion("bohr", "angstrom").value() + class TestPWSCFTrajSkip(unittest.TestCase, CompSys, IsPBC): def setUp(self): self.system_1 = dpdata.System(os.path.join('qe.traj', 'traj6'), @@ -42,7 +44,7 @@ def test_cell(self): ref_cell = [5.359985500701728967e+00, 0, 0, 3.585941820098031974e-01, 5.317218997480877896e+00, 0, 7.606780476053129902e-01, 7.811107228901693622e-01, 5.715864930517207121e+00 ] - ref_cell = 0.52917721067 * np.array(ref_cell).reshape(3,3) + ref_cell = bohr2ang * np.array(ref_cell).reshape(3,3) for ii in range(3) : for jj in range(3) : @@ -51,7 +53,7 @@ def test_cell(self): ref_cell = [5.308510801020571712e+00, 0, 0, 3.076052782312116429e-01, 5.279388982187173340e+00, 0, 4.321921336152507731e-01, 8.121110815096156399e-01, 5.301664983741235737e+00] - ref_cell = 0.52917721067 * np.array(ref_cell).reshape(3,3) + ref_cell = bohr2ang * np.array(ref_cell).reshape(3,3) for ii in range(3) : for jj in range(3) :