Total energy of G0W0 calculation

Dear all,
I am using PySCF to do G0W0 calculations.
I ran the example provided on PySCF forum.
But even though I set the ‘mol.verbose = 9’, there are only quasiparticle orbital energies given.
I wonder if it is possible to output the total energy and if possible, how to do that?
I post the input file below:
######################################
#!/usr/bin/env python
'''
A simple example to run a GW calculation
'''
from pyscf import gto, dft, gw
mol = gto.M(
atom = 'H 0 0 0; F 0 0 1.1',
basis = 'ccpvdz')
mol.output = '/home/scratch/zhishuo/PySCF/TEST_examples/GW'
mol.verbose = 9
mf = dft.RKS(mol)
mf.xc = 'pbe'
mf.kernel()
nocc = mol.nelectron//2
By default, GW is done with analytic continuation
gw = gw.GW(mf)
same as gw = gw.GW(mf, freq_int='ac')
gw.kernel()
print(gw.mo_energy)
######################################
with the output file shown as following:
######################################
******** <class 'pyscf.dft.rks.RKS'> ********
method = RKS-RHF
initial guess = minao
damping factor = 0
level_shift factor = 0
DIIS = <class 'pyscf.scf.diis.CDIIS'>
diis_start_cycle = 1
diis_space = 8
SCF conv_tol = 1e-09
SCF conv_tol_grad = None
SCF max_cycles = 50
direct_scf = True
direct_scf_tol = 1e-13
chkfile to save SCF result = /home/scratch/zhishuo/PySCF/TEST_examples/GW/tmpbawer1v3
max_memory 4000 MB (current use 75 MB)
XC library pyscf.dft.libxc version 5.1.7
S. Lehtola, C. Steigemann, M. J. Oliveira, and M. A. Marques, SoftwareX 7, 1 (2018)
XC functionals = pbe
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 78, 1396 (1997)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 78, 1396 (1997)
radial grids:
Treutler-Ahlrichs [JCP 102, 346 (1995); DOI:10.1063/1.469408] (M4) radial grids
becke partition: Becke, JCP 88, 2547 (1988); DOI:10.1063/1.454033
pruning grids: <function nwchem_prune at 0x7f4096f4e040>
grids dens level: 3
symmetrized grids: False
atomic radii adjust function: <function treutler_atomic_radii_adjust at 0x7f4096f3ff70>
atomic_radii : [3.77945036 0.66140414 2.64561657 2.74010288 1.98421243 1.60626721
1.32280829 1.22832198 1.13383567 0.94486306 2.83458919 3.40150702
2.83458919 2.36215766 2.07869874 1.88972612 1.88972612 1.88972612
3.40150702 4.15739747 3.40150702 3.0235618  2.64561657 2.55113027
2.64561657 2.64561657 2.64561657 2.55113027 2.55113027 2.55113027
2.55113027 2.45664396 2.36215766 2.17318504 2.17318504 2.17318504
3.59047964 4.44085639 3.77945225 3.40150702 2.92907549 2.74010288
2.74010288 2.55113027 2.45664396 2.55113027 2.64561657 3.0235618
2.92907549 2.92907549 2.74010288 2.74010288 2.64561657 2.64561657
3.96842486 4.91328792 4.06291117 3.68496594 3.49599333 3.49599333
3.49599333 3.49599333 3.49599333 3.49599333 3.40150702 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
2.92907549 2.74010288 2.55113027 2.55113027 2.45664396 2.55113027
2.55113027 2.55113027 2.83458919 3.59047964 3.40150702 3.0235618
3.59047964 2.74010288 3.96842486 3.40150702 4.06291117 3.68496594
3.40150702 3.40150702 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072]
small_rho_cutoff = 1e-07
Set gradient conv threshold to 3.16228e-05
Nelec from initial guess = 9.977410093935312
atom H rad-grids = 50, ang-grids = [ 50  50  50  50  50  50  50  50  50  50  50  50  50  50  86  86  86  86
266 266 266 266 302 302 302 302 302 302 302 302 302 302 302 302 302 266
266 266 266 266 266 266 266 266 266 266 266 266 266 266]
atom F rad-grids = 75, ang-grids = [ 50  50  50  50  50  50  50  50  50  50  50  50  50  50  50  50  50  50
50  50  50  86  86  86  86  86  86  86  86  86 266 266 266 266 266 266
302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302
266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266
266 266 266]
tot grids = 24466
Drop grids 2965
CPU time for setting up grids      0.02 sec, wall time      0.02 sec
nelec by numeric integration = 9.977410058506097
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -147.79366184144808  Ecoul = 53.81545980555761  Exc = -10.56315078337966
init E= -100.211721093561
cond(S) = 36.03560971202316
CPU time for initialize scf      0.09 sec, wall time      0.10 sec
HOMO = -0.433574392955141  LUMO = -0.0442325133887142
mo_energy =
[-24.42316777  -1.14925327  -0.5181955   -0.43357439  -0.43357439
-0.04423251   0.43679971   0.94825256   0.94825256   1.04310565
1.30276609   1.30276609   1.58236173   1.82722031   3.35715328
3.35715413   3.48370399   3.48370399   4.12568286]
nelec by numeric integration = 9.999999954153106
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -151.0872629864317  Ecoul = 57.46381765256835  Exc = -10.95113537035863
cycle= 1 E= -100.244948978513  delta_E= -0.0332  |g|= 0.628  |ddm|= 0.912
CPU time for cycle= 1      0.04 sec, wall time      0.04 sec
diis-norm(errvec)=1.02923
diis-c [-1.05932303  1.        ]
HOMO = -0.123004851649509  LUMO = 0.0172140447828145
mo_energy =
[-2.38123704e+01 -8.61303203e-01 -3.14577195e-01 -1.23004852e-01
-1.23004852e-01  1.72140448e-02  4.33449184e-01  1.09442013e+00
1.09442013e+00  1.27256228e+00  1.37012124e+00  1.37012124e+00
1.69046388e+00  1.96509937e+00  3.68634249e+00  3.68634349e+00
3.78246750e+00  3.78246750e+00  4.38773121e+00]
nelec by numeric integration = 9.999999942292893
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -144.99518711059403  Ecoul = 50.830645460069775  Exc = -10.20529865697312
cycle= 2 E= -100.040208581788  delta_E= 0.205  |g|= 1.25  |ddm|= 1.05
CPU time for cycle= 2      0.03 sec, wall time      0.03 sec
diis-norm(errvec)=2.02419
diis-c [-0.03198082  0.66546961  0.33453039]
HOMO = -0.296437558761834  LUMO = -0.0318248464538084
mo_energy =
[-24.23667377  -1.02673171  -0.42729085  -0.29643756  -0.29643756
-0.03182485   0.43198499   1.01987653   1.01987653   1.14175934
1.31753693   1.31753693   1.62424315   1.87830502   3.507597
3.50759796   3.60746965   3.60746965   4.22947123]
nelec by numeric integration = 9.999999949793388
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.21611744191034  Ecoul = 55.26765518445286  Exc = -10.698087527570502
cycle= 3 E= -100.316918059319  delta_E= -0.277  |g|= 0.0213  |ddm|= 0.723
CPU time for cycle= 3      0.03 sec, wall time      0.03 sec
diis-norm(errvec)=0.0437793
diis-c [-3.81030875e-04  3.78042053e-02  3.57987417e-02  9.26397053e-01]
HOMO = -0.292325571653969  LUMO = -0.0332430466973048
mo_energy =
[-24.23517775  -1.02185256  -0.42455774  -0.29232557  -0.29232557
-0.03324305   0.42891064   1.02131079   1.02131079   1.14461189
1.31701752   1.31701752   1.62592025   1.8799779    3.51228341
3.51228436   3.61268811   3.61268811   4.23368417]
nelec by numeric integration = 9.999999949142088
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.14948106689314  Ecoul = 55.19357036944616  Exc = -10.690773726178836
cycle= 4 E= -100.317052697917  delta_E= -0.000135  |g|= 0.00284  |ddm|= 0.0185
CPU time for cycle= 4      0.03 sec, wall time      0.03 sec
diis-norm(errvec)=0.00404271
diis-c [-4.86602792e-06  3.80053251e-04  3.08298350e-03  1.06646958e-01
8.89890005e-01]
HOMO = -0.292116387532438  LUMO = -0.0333014864608095
mo_energy =
[-24.23349534  -1.02179203  -0.42426063  -0.29211639  -0.29211639
-0.03330149   0.42895946   1.02156708   1.02156708   1.14475524
1.31706167   1.31706167   1.62592571   1.88001411   3.51240356
3.51240451   3.61272772   3.61272772   4.23375717]
nelec by numeric integration = 9.999999949190618
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.14500123850584  Ecoul = 55.18834477489713  Exc = -10.690029108801216
cycle= 5 E= -100.317053846701  delta_E= -1.15e-06  |g|= 0.000271  |ddm|= 0.00173
CPU time for cycle= 5      0.03 sec, wall time      0.03 sec
diis-norm(errvec)=0.000520753
diis-c [-2.58571251e-08 -1.09777162e-04  7.51218634e-06 -6.75273608e-03
8.64788948e-02  9.20376106e-01]
HOMO = -0.292130792516525  LUMO = -0.0333531699225601
mo_energy =
[-24.23351507  -1.0218054   -0.42432539  -0.29213079  -0.29213079
-0.03335317   0.42888589   1.02153383   1.02153383   1.14471416
1.31700999   1.31700999   1.62588855   1.87995637   3.51239195
3.51239291   3.61269377   3.61269377   4.23371024]
nelec by numeric integration = 9.999999949185455
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.14466195985864  Ecoul = 55.18797496026121  Exc = -10.689998587271836
cycle= 6 E= -100.31705386116  delta_E= -1.45e-08  |g|= 2.91e-05  |ddm|= 0.000155
CPU time for cycle= 6      0.03 sec, wall time      0.03 sec
diis-norm(errvec)=5.26053e-05
diis-c [-6.89065592e-10 -3.73588694e-05 -9.46190378e-06 -1.37443837e-03
7.88163419e-03  1.52532667e-01  8.41006958e-01]
HOMO = -0.292133728435749  LUMO = -0.0333585859092196
mo_energy =
[-24.23353017  -1.0218092   -0.42432717  -0.29213373  -0.29213373
-0.03335859   0.42887909   1.0215283    1.0215283    1.1447179
1.31700436   1.31700436   1.62587931   1.87995254   3.51238982
3.51239078   3.61269331   3.61269331   4.23370895]
nelec by numeric integration = 9.999999949184097
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.14465652319063  Ecoul = 55.18797201653816  Exc = -10.69000108041619
cycle= 7 E= -100.31705386136  delta_E= -1.99e-10  |g|= 8.54e-06  |ddm|= 3.32e-05
CPU time for cycle= 7      0.03 sec, wall time      0.03 sec
HOMO = -0.292130748240692  LUMO = -0.033355737045129
mo_energy =
[-24.23352062  -1.02180575  -0.42432416  -0.29213075  -0.29213075
-0.03335574   0.42888222   1.02153141   1.02153141   1.1447189
1.31700698   1.31700698   1.62588325   1.87995494   3.51239239
3.51239335   3.61269597   3.61269597   4.23371182]
nelec by numeric integration = 9.999999949184428
CPU time for vxc      0.03 sec, wall time      0.03 sec
E1 = -149.14461413263734  Ecoul = 55.187923145283825  Exc = -10.689994599717537
Extra cycle  E= -100.317053861362  delta_E= -2.4e-12  |g|= 8.47e-06  |ddm|= 9.75e-06
CPU time for scf_cycle      0.35 sec, wall time      0.35 sec
CPU time for SCF      0.35 sec, wall time      0.35 sec
converged SCF energy = -100.317053861362
Default auxbasis cc-pvdz-ri is used for F ccpvdz
Default auxbasis cc-pvdz-ri is used for H ccpvdz
******** <class 'pyscf.gw.gw_ac.GWAC'> ********
method = GWAC
GW nocc = 5, nvir = 14
use perturbative linearized QP eqn = False
analytic continuation method = pade
******** <class 'pyscf.df.df.DF'> ********
auxbasis = {'F': 'cc-pvdz-ri', 'H': 'cc-pvdz-ri'}
max_memory = 4000
_cderi_to_save = /home/scratch/zhishuo/PySCF/TEST_examples/GW/tmpgoo607ja
num shells = 24, num cGTOs = 70
size of aux basis 70
CPU time for 2c2e      0.00 sec, wall time      0.00 sec
shranges = [(0, 8, 190)]
int3c2e [1/1], AO [0:8], nrow = 190
CPU time for cholesky_eri      0.00 sec, wall time      0.00 sec
nelec by numeric integration = 9.999999949184426
CPU time for vxc      0.03 sec, wall time      0.03 sec
GW mo_energy =
[-26.20887482  -1.34556612  -0.63970147  -0.50371319  -0.50371319
0.13820756   0.61429425   1.29017698   1.2901769    1.38517064
1.56518126   1.56518087   1.91850741   2.15492359   3.86029428
3.86028626   4.02315672   4.02315594   4.3233935 ]
WARN: GW QP energies may not be sorted from min to max
CPU time for GW      0.19 sec, wall time      0.19 sec

[-26.20887482  -1.34556612  -0.63970147  -0.50371319  -0.50371319
0.13820756   0.61429425   1.29017698   1.2901769    1.38517064
1.56518126   1.56518087   1.91850741   2.15492359   3.86029428
3.86028626   4.02315672   4.02315594   4.3233935 ]
######################################
Best regards
Zhishuo Huang

The GW approximation is rarely used to calculate total energies. Most
similar is the RPA total energy, an efficient implementation of which is
not currently in PySCF but is a planned feature.
On Fri, Dec 31, 2021 at 6:48 AM ZhishuoHuang @.***>
wrote:

Dear all,
I am using PySCF to do G0W0 calculations.
I ran the example provided on PySCF forum.
But even though I set the ‘mol.verbose = 9’, there are only quasiparticle
orbital energies given.
I wonder if it is possible to output the total energy and if possible, how
to do that?
I post the input file below:
######################################
#!/usr/bin/env python
'''
A simple example to run a GW calculation
'''
from pyscf import gto, dft, gw
mol = gto.M(
atom = 'H 0 0 0; F 0 0 1.1',
basis = 'ccpvdz')
mol.output = '/home/scratch/zhishuo/PySCF/TEST_examples/GW'
mol.verbose = 9
mf = dft.RKS(mol)
mf.xc = 'pbe'
mf.kernel()
nocc = mol.nelectron//2
By default, GW is done with analytic continuation
gw = gw.GW(mf)
same as gw = gw.GW(mf, freq_int='ac')
gw.kernel()
print(gw.mo_energy)
######################################
with the output file shown as following:
######################################
******** <class 'pyscf.dft.rks.RKS'> ********
method = RKS-RHF
initial guess = minao
damping factor = 0
level_shift factor = 0
DIIS = <class 'pyscf.scf.diis.CDIIS'>
diis_start_cycle = 1
diis_space = 8
SCF conv_tol = 1e-09
SCF conv_tol_grad = None
SCF max_cycles = 50
direct_scf = True
direct_scf_tol = 1e-13
chkfile to save SCF result =
/home/scratch/zhishuo/PySCF/TEST_examples/GW/tmpbawer1v3
max_memory 4000 MB (current use 75 MB)
XC library pyscf.dft.libxc version 5.1.7
S. Lehtola, C. Steigemann, M. J. Oliveira, and M. A. Marques, SoftwareX 7,
1 (2018)
XC functionals = pbe
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 78, 1396 (1997)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996)
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 78, 1396 (1997)
radial grids:
Treutler-Ahlrichs [JCP 102, 346 (1995); DOI:10.1063/1.469408] (M4) radial
grids
becke partition: Becke, JCP 88, 2547 (1988); DOI:10.1063/1.454033
pruning grids: <function nwchem_prune at 0x7f4096f4e040>
grids dens level: 3
symmetrized grids: False
atomic radii adjust function: <function treutler_atomic_radii_adjust at
0x7f4096f3ff70>
atomic_radii : [3.77945036 0.66140414 2.64561657 2.74010288 1.98421243
1.60626721
1.32280829 1.22832198 1.13383567 0.94486306 2.83458919 3.40150702
2.83458919 2.36215766 2.07869874 1.88972612 1.88972612 1.88972612
3.40150702 4.15739747 3.40150702 3.0235618 2.64561657 2.55113027
2.64561657 2.64561657 2.64561657 2.55113027 2.55113027 2.55113027
2.55113027 2.45664396 2.36215766 2.17318504 2.17318504 2.17318504
3.59047964 4.44085639 3.77945225 3.40150702 2.92907549 2.74010288
2.74010288 2.55113027 2.45664396 2.55113027 2.64561657 3.0235618
2.92907549 2.92907549 2.74010288 2.74010288 2.64561657 2.64561657
3.96842486 4.91328792 4.06291117 3.68496594 3.49599333 3.49599333
3.49599333 3.49599333 3.49599333 3.49599333 3.40150702 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
2.92907549 2.74010288 2.55113027 2.55113027 2.45664396 2.55113027
2.55113027 2.55113027 2.83458919 3.59047964 3.40150702 3.0235618
3.59047964 2.74010288 3.96842486 3.40150702 4.06291117 3.68496594
3.40150702 3.40150702 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072 3.30702072
3.30702072 3.30702072 3.30702072 3.30702072 3.30702072]
small_rho_cutoff = 1e-07
Set gradient conv threshold to 3.16228e-05
Nelec from initial guess = 9.977410093935312
atom H rad-grids = 50, ang-grids = [ 50 50 50 50 50 50 50 50 50 50 50 50
50 50 86 86 86 86
266 266 266 266 302 302 302 302 302 302 302 302 302 302 302 302 302 266
266 266 266 266 266 266 266 266 266 266 266 266 266 266]
atom F rad-grids = 75, ang-grids = [ 50 50 50 50 50 50 50 50 50 50 50 50
50 50 50 50 50 50
50 50 50 86 86 86 86 86 86 86 86 86 266 266 266 266 266 266
302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302 302
266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266 266
266 266 266]
tot grids = 24466
Drop grids 2965
CPU time for setting up grids 0.02 sec, wall time 0.02 sec
nelec by numeric integration = 9.977410058506097
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -147.79366184144808 Ecoul = 53.81545980555761 Exc = -10.56315078337966
init E= -100.211721093561
cond(S) = 36.03560971202316
CPU time for initialize scf 0.09 sec, wall time 0.10 sec
HOMO = -0.433574392955141 LUMO = -0.0442325133887142
mo_energy =
[-24.42316777 -1.14925327 -0.5181955 -0.43357439 -0.43357439
-0.04423251 0.43679971 0.94825256 0.94825256 1.04310565
1.30276609 1.30276609 1.58236173 1.82722031 3.35715328
3.35715413 3.48370399 3.48370399 4.12568286]
nelec by numeric integration = 9.999999954153106
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -151.0872629864317 Ecoul = 57.46381765256835 Exc = -10.95113537035863
cycle= 1 E= -100.244948978513 delta_E= -0.0332 |g|= 0.628 |ddm|= 0.912
CPU time for cycle= 1 0.04 sec, wall time 0.04 sec
diis-norm(errvec)=1.02923
diis-c [-1.05932303 1. ]
HOMO = -0.123004851649509 LUMO = 0.0172140447828145
mo_energy =
[-2.38123704e+01 -8.61303203e-01 -3.14577195e-01 -1.23004852e-01
-1.23004852e-01 1.72140448e-02 4.33449184e-01 1.09442013e+00
1.09442013e+00 1.27256228e+00 1.37012124e+00 1.37012124e+00
1.69046388e+00 1.96509937e+00 3.68634249e+00 3.68634349e+00
3.78246750e+00 3.78246750e+00 4.38773121e+00]
nelec by numeric integration = 9.999999942292893
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -144.99518711059403 Ecoul = 50.830645460069775 Exc =
-10.20529865697312
cycle= 2 E= -100.040208581788 delta_E= 0.205 |g|= 1.25 |ddm|= 1.05
CPU time for cycle= 2 0.03 sec, wall time 0.03 sec
diis-norm(errvec)=2.02419
diis-c [-0.03198082 0.66546961 0.33453039]
HOMO = -0.296437558761834 LUMO = -0.0318248464538084
mo_energy =
[-24.23667377 -1.02673171 -0.42729085 -0.29643756 -0.29643756
-0.03182485 0.43198499 1.01987653 1.01987653 1.14175934
1.31753693 1.31753693 1.62424315 1.87830502 3.507597
3.50759796 3.60746965 3.60746965 4.22947123]
nelec by numeric integration = 9.999999949793388
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.21611744191034 Ecoul = 55.26765518445286 Exc =
-10.698087527570502
cycle= 3 E= -100.316918059319 delta_E= -0.277 |g|= 0.0213 |ddm|= 0.723
CPU time for cycle= 3 0.03 sec, wall time 0.03 sec
diis-norm(errvec)=0.0437793
diis-c [-3.81030875e-04 3.78042053e-02 3.57987417e-02 9.26397053e-01]
HOMO = -0.292325571653969 LUMO = -0.0332430466973048
mo_energy =
[-24.23517775 -1.02185256 -0.42455774 -0.29232557 -0.29232557
-0.03324305 0.42891064 1.02131079 1.02131079 1.14461189
1.31701752 1.31701752 1.62592025 1.8799779 3.51228341
3.51228436 3.61268811 3.61268811 4.23368417]
nelec by numeric integration = 9.999999949142088
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.14948106689314 Ecoul = 55.19357036944616 Exc =
-10.690773726178836
cycle= 4 E= -100.317052697917 delta_E= -0.000135 |g|= 0.00284 |ddm|= 0.0185
CPU time for cycle= 4 0.03 sec, wall time 0.03 sec
diis-norm(errvec)=0.00404271
diis-c [-4.86602792e-06 3.80053251e-04 3.08298350e-03 1.06646958e-01
8.89890005e-01]
HOMO = -0.292116387532438 LUMO = -0.0333014864608095
mo_energy =
[-24.23349534 -1.02179203 -0.42426063 -0.29211639 -0.29211639
-0.03330149 0.42895946 1.02156708 1.02156708 1.14475524
1.31706167 1.31706167 1.62592571 1.88001411 3.51240356
3.51240451 3.61272772 3.61272772 4.23375717]
nelec by numeric integration = 9.999999949190618
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.14500123850584 Ecoul = 55.18834477489713 Exc =
-10.690029108801216
cycle= 5 E= -100.317053846701 delta_E= -1.15e-06 |g|= 0.000271 |ddm|=
0.00173
CPU time for cycle= 5 0.03 sec, wall time 0.03 sec
diis-norm(errvec)=0.000520753
diis-c [-2.58571251e-08 -1.09777162e-04 7.51218634e-06 -6.75273608e-03
8.64788948e-02 9.20376106e-01]
HOMO = -0.292130792516525 LUMO = -0.0333531699225601
mo_energy =
[-24.23351507 -1.0218054 -0.42432539 -0.29213079 -0.29213079
-0.03335317 0.42888589 1.02153383 1.02153383 1.14471416
1.31700999 1.31700999 1.62588855 1.87995637 3.51239195
3.51239291 3.61269377 3.61269377 4.23371024]
nelec by numeric integration = 9.999999949185455
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.14466195985864 Ecoul = 55.18797496026121 Exc =
-10.689998587271836
cycle= 6 E= -100.31705386116 delta_E= -1.45e-08 |g|= 2.91e-05 |ddm|=
0.000155
CPU time for cycle= 6 0.03 sec, wall time 0.03 sec
diis-norm(errvec)=5.26053e-05
diis-c [-6.89065592e-10 -3.73588694e-05 -9.46190378e-06 -1.37443837e-03
7.88163419e-03 1.52532667e-01 8.41006958e-01]
HOMO = -0.292133728435749 LUMO = -0.0333585859092196
mo_energy =
[-24.23353017 -1.0218092 -0.42432717 -0.29213373 -0.29213373
-0.03335859 0.42887909 1.0215283 1.0215283 1.1447179
1.31700436 1.31700436 1.62587931 1.87995254 3.51238982
3.51239078 3.61269331 3.61269331 4.23370895]
nelec by numeric integration = 9.999999949184097
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.14465652319063 Ecoul = 55.18797201653816 Exc = -10.69000108041619
cycle= 7 E= -100.31705386136 delta_E= -1.99e-10 |g|= 8.54e-06 |ddm|=
3.32e-05
CPU time for cycle= 7 0.03 sec, wall time 0.03 sec
HOMO = -0.292130748240692 LUMO = -0.033355737045129
mo_energy =
[-24.23352062 -1.02180575 -0.42432416 -0.29213075 -0.29213075
-0.03335574 0.42888222 1.02153141 1.02153141 1.1447189
1.31700698 1.31700698 1.62588325 1.87995494 3.51239239
3.51239335 3.61269597 3.61269597 4.23371182]
nelec by numeric integration = 9.999999949184428
CPU time for vxc 0.03 sec, wall time 0.03 sec
E1 = -149.14461413263734 Ecoul = 55.187923145283825 Exc =
-10.689994599717537
Extra cycle E= -100.317053861362 delta_E= -2.4e-12 |g|= 8.47e-06 |ddm|=
9.75e-06
CPU time for scf_cycle 0.35 sec, wall time 0.35 sec
CPU time for SCF 0.35 sec, wall time 0.35 sec
converged SCF energy = -100.317053861362
Default auxbasis cc-pvdz-ri is used for F ccpvdz
Default auxbasis cc-pvdz-ri is used for H ccpvdz
******** <class 'pyscf.gw.gw_ac.GWAC'> ********
method = GWAC
GW nocc = 5, nvir = 14
use perturbative linearized QP eqn = False
analytic continuation method = pade
******** <class 'pyscf.df.df.DF'> ********
auxbasis = {'F': 'cc-pvdz-ri', 'H': 'cc-pvdz-ri'}
max_memory = 4000
_cderi_to_save = /home/scratch/zhishuo/PySCF/TEST_examples/GW/tmpgoo607ja
num shells = 24, num cGTOs = 70
size of aux basis 70
CPU time for 2c2e 0.00 sec, wall time 0.00 sec
shranges = [(0, 8, 190)]
int3c2e [1/1], AO [0:8], nrow = 190
CPU time for cholesky_eri 0.00 sec, wall time 0.00 sec
nelec by numeric integration = 9.999999949184426
CPU time for vxc 0.03 sec, wall time 0.03 sec
GW mo_energy =
[-26.20887482 -1.34556612 -0.63970147 -0.50371319 -0.50371319
0.13820756 0.61429425 1.29017698 1.2901769 1.38517064
1.56518126 1.56518087 1.91850741 2.15492359 3.86029428
3.86028626 4.02315672 4.02315594 4.3233935 ]
WARN: GW QP energies may not be sorted from min to max
CPU time for GW      0.19 sec, wall time      0.19 sec
[-26.20887482 -1.34556612 -0.63970147 -0.50371319 -0.50371319
0.13820756 0.61429425 1.29017698 1.2901769 1.38517064
1.56518126 1.56518087 1.91850741 2.15492359 3.86029428
3.86028626 4.02315672 4.02315594 4.3233935 ]
######################################
Best regards
Zhishuo Huang
—
Reply to this email directly, view it on GitHub
https://github.com/pyscf/pyscf/issues/1149, or unsubscribe
https://github.com/notifications/unsubscribe-auth/ADGAJCY7YU5V2AFGHTL45KDUTWKAZANCNFSM5LBGFWXA
.
Triage notifications on the go with GitHub Mobile for iOS
https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675
or Android
https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.
You are receiving this because you are subscribed to this thread.Message
ID: @.***>

Hello, thank you for your reply.
There is one approximate way of total energy calculation, of which the total energy is the summation of total energy at RPA level + the Galitskii-Migdal correction.
So the total energy at RPA level is not currently available in PySCF?

Indeed, given the full GW Green's function (and not just its quasiparticle approximation), one could in principle calculate the total energy by frequency integration, although controlling this integration error is numerically nontrivial.
An efficient RPA total energy calculation is not currently available in PySCF. Perhaps @zhutianyu1991, who wrote the latest GW implementation in PySCF, could tell us if/when he plans to contribute RPA total energies.
Two inefficient options are:

The plasmon formula for the RPA energy. PySCF does have implementations of RPA and TDA excitation energies, the summed difference of which gives the RPA total energy. But this requires enumerating all excitation energies and so will be O(N^6) cost.

The ring CCD version of the RPA energy. PySCF's CCSD codes can be modified to only include ring (or direct ring) diagrams, which also comes with O(N^6) or O(N^5) cost. My group is working a little on this, and has a few hacky codes that we could share. If you're interested in this, you can contact me offline and to discuss further.

Hi, current GW implementations only compute the quasiparticle energy, and the total energy is not computed. For non-self-consistent GW (G0W0), there're a few definitions of total energy (such as Galitskii-Migdal or Klein) that give different values, but I don't know if there's consensus on which one is the best.
Instead, I think RPA energy is what many literatures reported as "G0W0 energy". I have an O(N^4) implementation of RPA total energy (which re-uses many routines in GW-AC code) in my local repo, and should probably make a pull request to merge this code into PySCF's main branch now. In the meantime, if you want to use this code sooner, I'll email you about it.

Indeed, given the full GW Green's function (and not just its quasiparticle approximation), one could in principle calculate the total energy by frequency integration, although controlling this integration error is numerically nontrivial.
An efficient RPA total energy calculation is not currently available in PySCF. Perhaps @zhutianyu1991, who wrote the latest GW implementation in PySCF, could tell us if/when he plans to contribute RPA total energies.
Two inefficient options are:

The plasmon formula for the RPA energy. PySCF does have implementations of RPA and TDA excitation energies, the summed difference of which gives the RPA total energy. But this requires enumerating all excitation energies and so will be O(N^6) cost.
The ring CCD version of the RPA energy. PySCF's CCSD codes can be modified to only include ring (or direct ring) diagrams, which also comes with O(N^6) or O(N^5) cost. My group is working a little on this, and has a few hacky codes that we could share. If you're interested in this, you can contact me offline and to discuss further.

I see.
Thank you for your detailed explanation.

Indeed, given the full GW Green's function (and not just its quasiparticle approximation), one could in principle calculate the total energy by frequency integration, although controlling this integration error is numerically nontrivial.
An efficient RPA total energy calculation is not currently available in PySCF. Perhaps @zhutianyu1991, who wrote the latest GW implementation in PySCF, could tell us if/when he plans to contribute RPA total energies.
Two inefficient options are:

The plasmon formula for the RPA energy. PySCF does have implementations of RPA and TDA excitation energies, the summed difference of which gives the RPA total energy. But this requires enumerating all excitation energies and so will be O(N^6) cost.
The ring CCD version of the RPA energy. PySCF's CCSD codes can be modified to only include ring (or direct ring) diagrams, which also comes with O(N^6) or O(N^5) cost. My group is working a little on this, and has a few hacky codes that we could share. If you're interested in this, you can contact me offline and to discuss further.

I see.
Thank you for your detailed explanation.

Hi, current GW implementations only compute the quasiparticle energy, and the total energy is not computed. For non-self-consistent GW (G0W0), there're a few definitions of total energy (such as Galitskii-Migdal or Klein) that give different values, but I don't know if there's consensus on which one is the best. Instead, I think RPA energy is what many literatures reported as "G0W0 energy". I have an O(N^4) implementation of RPA total energy (which re-uses many routines in GW-AC code) in my local repo, and should probably make a pull request to merge this code into PySCF's main branch now. In the meantime, if you want to use this code sooner, I'll email you about it.

Hi, could you please email me the code, I want to use it sooner.
My email address is<EMAIL_ADDRESS>I will appreciate that.

Hi, current GW implementations only compute the quasiparticle energy, and the total energy is not computed. For non-self-consistent GW (G0W0), there're a few definitions of total energy (such as Galitskii-Migdal or Klein) that give different values, but I don't know if there's consensus on which one is the best. Instead, I think RPA energy is what many literatures reported as "G0W0 energy". I have an O(N^4) implementation of RPA total energy (which re-uses many routines in GW-AC code) in my local repo, and should probably make a pull request to merge this code into PySCF's main branch now. In the meantime, if you want to use this code sooner, I'll email you about it.

Hi, could you please email me the code, I want to use it sooner.
My email address is<EMAIL_ADDRESS>I will appreciate that.

Hi @ZhishuoHuang , has your question been answered?

Hi @ZhishuoHuang , has your question been answered?
Yes, to some extent I know that currently, total energy is only available at RPA level.

Indeed, given the full GW Green's function (and not just its quasiparticle approximation), one could in principle calculate the total energy by frequency integration, although controlling this integration error is numerically nontrivial.
An efficient RPA total energy calculation is not currently available in PySCF. Perhaps @zhutianyu1991, who wrote the latest GW implementation in PySCF, could tell us if/when he plans to contribute RPA total energies.
Two inefficient options are:

The plasmon formula for the RPA energy. PySCF does have implementations of RPA and TDA excitation energies, the summed difference of which gives the RPA total energy. But this requires enumerating all excitation energies and so will be O(N^6) cost.
The ring CCD version of the RPA energy. PySCF's CCSD codes can be modified to only include ring (or direct ring) diagrams, which also comes with O(N^6) or O(N^5) cost. My group is working a little on this, and has a few hacky codes that we could share. If you're interested in this, you can contact me offline and to discuss further.

@tberkel: I wonder if the ring CCD is available somewhere in PySCF or it is still an internal code?

No, it's not in PySCF. I don't know if it's of sufficiently broad interest that it needs to be in PySCF. But it's really only a few lines to modify the update_amps function. All the t1 stuff can be deleted, and only a few contractions need to be kept for ring CCD. If I have time, maybe I'll write something and place it in the examples/cc folder.

An example in pyscf/examples/cc would be helpful. I am interested in using it. I'll email you offline for further discussion then. Thanks for the prompt response.

Hi friends, I wrote a demo code by modifying the update_amps and energy in pyscf/cc/rccsd_slow.py according to the discussion above. I put the file here, maybe someone will find it helpful in future.
'''
Restricted direct ring-CCD

Ref: Scuseria et al., J. Chem. Phys. 129, 231101 (2008)
Ref: Masios et al., Phys. Rev. Lett. 131, 186401 (2023)

Original Source:
Modified from pyscf/cc/rccsd_slow.py

Contact:
<EMAIL_ADDRESS>'''

from functools import reduce
import numpy as np
from pyscf.cc import ccsd
from pyscf import lib
from pyscf.ao2mo import _ao2mo

einsum = lib.einsum

def update_amps(cc, t1, t2, eris):
    # Ref: Masios et al., Phys. Rev. Lett. 131, 186401 (2023) Eq.(3)
    nocc, nvir = t1.shape
    t1[:] = 0.0

    eris_ovvo = np.asarray(eris.ovvo)
    eris_ovov = np.asarray(eris.ovov)

    t2new = 2*einsum('iack, kjcb->iajb', eris_ovvo, t2)
    t2new += 2*einsum('ikac, jbck->iajb', t2, eris_ovvo)
    tmp = einsum('ikac, kcld->iald', t2, eris_ovov)
    t2new += 4*einsum('iald, ljdb->iajb', tmp, t2)
    t2new = t2new.transpose(0,2,1,3)

    t2new += np.asarray(eris.ovov).conj().transpose(0,2,1,3).copy()

    mo_e = eris.fock.diagonal().real
    eia = mo_e[:nocc,None] - mo_e[None,nocc:]
    eijab = lib.direct_sum('ia,jb->ijab',eia,eia)
    t2new /= eijab

    return t1, t2new

def energy(cc, t1, t2, eris):
    # Ref: Scuseria et al., J. Chem. Phys. 129, 231101 (2008) Eq.(11)
    nocc, nvir = t1.shape
    fock = eris.fock
    eris_ovov = np.asarray(eris.ovov)
    e = 2*einsum('ijab,iajb', t2, eris_ovov)
    return e.real

def init_amps(cc, eris):
    nocc = cc.nocc
    mo_e = eris.fock.diagonal().real
    eia = mo_e[:nocc,None] - mo_e[None,nocc:]
    eijab = lib.direct_sum('ia,jb->ijab', eia, eia)
    t1 = np.zeros_like(eris.fock[:nocc,nocc:])

    eris_ovov = np.asarray(eris.ovov)
    t2 = eris_ovov.transpose(0,2,1,3).conj() / eijab
    cc.emp2  = 2*einsum('ijab,iajb', t2, eris_ovov)
    cc.emp2 -=   einsum('ijab,ibja', t2, eris_ovov)
    lib.logger.info(cc, 'Init t2, MP2 energy = %.15g', cc.emp2)
    return cc.emp2, t1, t2

def ao2mo(cc, mo_coeff=None):
    return _ChemistsERIs(cc, mo_coeff)

class _ChemistsERIs:
    def __init__(self, cc, mo_coeff=None, method='incore'):
        if mo_coeff is None:
            self.mo_coeff = mo_coeff = ccsd._mo_without_core(cc, cc.mo_coeff)
        else:
            self.mo_coeff = mo_coeff = ccsd._mo_without_core(cc, mo_coeff)
        dm = cc._scf.make_rdm1(cc.mo_coeff, cc.mo_occ)
        fockao = cc._scf.get_hcore() + cc._scf.get_veff(cc.mol, dm)
        self.fock = reduce(np.dot, (mo_coeff.T, fockao, mo_coeff))

        nocc = cc.nocc
        nmo = cc.nmo
        nvir = nmo - nocc
        # eri1 = ao2mo.incore.full(cc._scf._eri, mo_coeff)
        # eri1 = ao2mo.restore(1, eri1, nmo)
        mo = cc._scf.mo_coeff
        Lpq = cc._scf._cderi
        ijslice = (0, nocc, nocc, nmo)
        Lov = _ao2mo.nr_e2(Lpq, mo, ijslice, aosym='s2', mosym='s1')
        eris_ovov = Lov.T @ Lov
        self.ovov = eris_ovov.reshape(nocc, nvir, nocc, nvir)
        self.ovvo = self.ovov.transpose(0,1,3,2).copy()
        self.Lov = Lov

class drCCD(ccsd.CCSD):
    init_amps = init_amps
    update_amps = update_amps
    energy = energy
    ao2mo = ao2mo

if __name__ == '__main__':
	from pyscf import gto, dft, scf
	mol = gto.Mole()
	mol.verbose = 4
	mol.atom = 'H 0 0 0; F 0 0 1.1'
	mol.basis = 'ccpvdz'

	mol.build()
	mf = scf.RHF(mol).density_fit()
	mf.kernel()

	mycc = drCCD(mf)
	e_drccd, _, t2 = mycc.kernel()
	from pyscf.gw import rpa
	myrpa = rpa.RPA(mf)
	e_rpa = myrpa.kernel()
	print(f"Difference between drccd and drpa: {e_drccd-e_rpa} Ha." )

Thanks a lot @changsuc for writing this. Since this question has come up multiple times, please feel free to turn your script into an example (to be placed into examples/cc) and submit a PR.


feat(tagging): add SessionService [EX-4875]

Pull request template
EX-4875
Motivation and context
This PR adds a new session id service, which allows us to get and remove the current session id.
Type of change

[ ] Bug fix (non-breaking change that fixes an issue)
[x] New feature (non-breaking change that adds functionality)
[ ] Breaking change (fix or feature that causes existing functionality to not work as expected)
[ ] Change requires a documentation update

What is the destination branch of this PR?

[x] Main
[ ] Other. Specify:

How has this been tested?
In order to see the new property in the localstorage you should register the tagging module in the Layout and add the wireServices in the tagging wiring.
Checklist:

[x] My code follows the style guidelines of this project.
[x] I have performed a self-review on my own code.
[x] I have commented my code, particularly in hard-to-understand areas.
[x] I have made corresponding changes to the documentation.
[x] My changes generate no new warnings.
[x] I have added tests that prove my fix is effective or that my feature works.
[x] New and existing unit tests pass locally with my changes.

amazing job, keep it up avilesino! 💪🏿


 Frequency filtering and gene list filtering work well. Very nice work!
But not consequence filtering, see below.
Also - to avoid adding lots of terms for filtering effects, in gemini we often use the LOW, MEDIUM, HIGH impact categories as summary terms (impact_severity). There is a detailed listing of impact_severity vs impacts vs SO_impacts here: http://gemini.readthedocs.org/en/latest/content/database_schema.html#details-of-the-impact-and-impact-severity-columns
[2016-02-04 09:15:58,769] ERROR in app: Exception on /variants/famall7 [GET]
Traceback (most recent call last):
  File "/Library/Python/2.7/site-packages/flask/app.py", line 1882, in wsgi_app
    response = self.full_dispatch_request()
  File "/Library/Python/2.7/site-packages/flask/app.py", line 1537, in full_dispatch_request
    rv = self.handle_user_exception(e)
  File "/Library/Python/2.7/site-packages/flask/app.py", line 1441, in handle_user_exception
    reraise(exc_type, exc_value, tb)
  File "/Library/Python/2.7/site-packages/flask/app.py", line 1535, in full_dispatch_request
    rv = self.dispatch_request()
  File "/Library/Python/2.7/site-packages/flask/app.py", line 1521, in dispatch_request
    return self.view_functions[rule.endpoint](**req.view_args)
  File "/Users/vaszac/proj/puzzle/puzzle/blueprints/variants/views.py", line 40, in variants
    return render_template('variants.html', **kwargs)
  File "/Library/Python/2.7/site-packages/flask/templating.py", line 123, in render_template
    context, ctx.app)
  File "/Library/Python/2.7/site-packages/flask/templating.py", line 105, in _render
    rv = template.render(context)
  File "/usr/local/lib/python2.7/site-packages/jinja2/environment.py", line 989, in render
    return self.environment.handle_exception(exc_info, True)
  File "/usr/local/lib/python2.7/site-packages/jinja2/environment.py", line 754, in handle_exception
    reraise(exc_type, exc_value, tb)
  File "/Users/vaszac/proj/puzzle/puzzle/blueprints/variants/templates/variants.html", line 3, in top-level template code
    {% from "macros/filters.html" import filters_form %}
  File "/Users/vaszac/proj/puzzle/puzzle/templates/layouts/base.html", line 30, in top-level template code
    {% block body %}
  File "/Users/vaszac/proj/puzzle/puzzle/templates/layouts/base.html", line 47, in block "body"
    {% block body_content %}{% endblock %}
  File "/Users/vaszac/proj/puzzle/puzzle/blueprints/variants/templates/variants.html", line 36, in block "body_content"
    {% for variant in variants %}
  File "/Users/vaszac/proj/puzzle/puzzle/plugins/gemini/variant_mixin.py", line 84, in variants
    if transcript['Consequence'] in consequences:
KeyError: 'Consequence'

Ok I know why this is happening, will fix asap

This is somehow related to #28 . I'm working on a branch to add this filter


"Invalid string length" error when querying states

Describe the bug
curl -Ssv "http://consensus:5052/eth/v2/debug/beacon/states/finalized"
*   Trying <IP_ADDRESS>:5052...
* Connected to consensus (<IP_ADDRESS>) port 5052
> GET /eth/v2/debug/beacon/states/finalized HTTP/1.1
> Host: consensus:5052
> User-Agent: curl/8.4.0
> Accept: */*
>
< HTTP/1.1 500 Internal Server Error
< vary: Origin
< access-control-allow-origin: *
< content-type: application/json; charset=utf-8
< content-length: 84
< Date: Thu, 23 Nov 2023 04:16:45 GMT
< Connection: keep-alive
< Keep-Alive: timeout=72
<
{ [84 bytes data]
* Connection #0 to host consensus left intact
{"statusCode":500,"error":"Internal Server Error","message":"Invalid string length"}

Expected behavior
State json or an error that is more descriptive
Steps to reproduce
(see above curl)
Additional context
No response
Operating system
Linux
Lodestar version or commit hash
v1.12.0

NB: This was on holesky, which I suppose has quite large states

We are hitting several limits now using json, there might be possible workarounds but those are potentially really slow and still require a lot of memory.
The error is thrown during serialization but the message it definitely not great, should be more descriptive but proper solution is to support ssz for all apis (https://github.com/ChainSafe/lodestar/pull/6080) and follow spec (e.g. set all metadata headers).

just tried again, this is still an issue

It's kinda good that it errors before converting the whole state to a stringified json as it would crash the node with default memory limit.
We support ssz on almost all apis now, just need to get alignment on the state validators api https://github.com/ethereum/beacon-APIs/pull/449 to be able to return a ssz response and this should be resolved.
I don't think it's worth to try and optimize the json response, if that's even technically possible.


Support field pages without version specified in URL

This URL works:
<EMAIL_ADDRESS>
This one doesn't:
https://registry.field33.com/field/@foaf_mirror/foaf

It should be possible to link to a field without having to specify an explicit version. This will make it a lot easier to allow for external sources to link to a field (e.g. also to support #20).

This has been fixed recently with https://registry.field33.com/field/@foaf_mirror/foaf redirecting to https://registry.field33.com/field/@foaf_mirror/foaf/latest, which displays the latest version.


Fix merge issue: remove incorrect and stale duplicates of targets

@raines Somehow the merge in #1105 brought in the fix from #1104 as a separate copy of the target and surrounding code, but above the old code that is still broken!
cc @eerhardt @srivatsn

@mmitche Is jenkins still having trouble getting through the queue?

@nguerrera Just a bit backed up.  I requested another core increase so that should help.  But it looks like you're stuff just started

Is the Mac build stuck?

@nguerrera looks like it. I'll be optimistic.


Make CARO ODK

Things to do:

[ ] caro-bridge.owl - put as component
[ ] other bridges keep as reference but don't fold in
[ ] Makefile - try with generic, the custom stuff probably not the most important
[ ] IMPORTS = caro uberon chebi pr go pco ro cl <- check for dangling class after


WSL hangs repeatedly when vscode indexing files.

Version
Microsoft Windows [Version 10.0.22622.598]
WSL Version

[X] WSL 2
[ ] WSL 1

Kernel Version
<IP_ADDRESS>
Distro Version
Ubuntu 20.04
Other Software
No response
Repro Steps
I have vscode and intellisense installed.  After a while, usually while indexing, vmmemWSl maxes out the CPU and the wsl instance freezes.  only killing the windows subsystem for linux process in task manager lets me restart it again.
Expected Behavior
it to not crash every hour or so.
Actual Behavior
Crashing, unresponsive, unusable.
Error message sometimes
A connection attempt failed because the connected party did not properly respond after a period of time, or established connection failed because connected host has failed to respond.
Error code: Wsl/Service/0x8007274c

[process exited with code<PHONE_NUMBER> (0xffffffff)]

Diagnostic Logs
WslLogs-2022-10-03_10-41-28.zip
WslLogs-2022-10-03_13-32-42.zip

my vs code editor terminal frequently terminated showing some exit code
any solution for this annoying kind of behaviour of vs code

Closing this and following #8824 as that is exactly the issue I am facing.


only try to use XDomainRequest for POST & GET

I'd actually prefer if XDomainRequest was removed entirely as there are too many caveats http://blogs.msdn.com/b/ieinternals/archive/2010/05/13/xdomainrequest-restrictions-limitations-and-workarounds.aspx and I'm working on a more general iframe + post message solution for ie9 https://github.com/Financial-Times/ft-xdr

:+1:


ScrollView: do scroll bounds checking early

The ScrollView will panic if scroll.y is given a number that is outside the content size.

Thanks! Not sure why that test randomly failed :/ But unrelated to this PR l


Make Key a clone type

This would be helpful, especially considering Key is a relatively small type.

Unfortunately this is not possible due to io::Error's lack of clone.

Have you considered taking the Error variant out of Key, and just iterating over io::Result<Key>?

That's a good idea, actually.


Huge memory consumption and/or crash (pythran 0.8.4.post0)

I am not able anymore to compile this code with Pythran 0.8.4.post0. It compiles fine with Pythran 0.8.3.
import numpy as np

# pythran export compute_subpix_2d_gaussian2(float32[][], int, int)
def compute_subpix_2d_gaussian2(correl, ix, iy):
    correl_crop = correl[iy-1:iy+2, ix-1:ix+2]
    # hoops, pythran crashes because of this line
    # correl_crop[correl_crop < 0] = 1e-6

    # we write it like this to please pythran
    tmp = np.where(correl_crop < 0)
    for i0, i1 in zip(tmp[0], tmp[1]):
        correl_crop[i0, i1] = 1e-6

    c10 = 0
    c01 = 0
    c11 = 0
    c20 = 0
    c02 = 0
    for i in range(3):
        for j in range(3):
            c10 += (i-1)*np.log(correl_crop[j, i])
            c01 += (j-1)*np.log(correl_crop[j, i])
            c11 += (i-1)*(j-1)*np.log(correl_crop[j, i])
            c20 += (3*(i-1)**2-2)*np.log(correl_crop[j, i])
            c02 += (3*(j-1)**2-2)*np.log(correl_crop[j, i])
            c00 = (5-3*(i-1)**2-3*(j-1)**2)*np.log(correl_crop[j, i])

    c00, c10, c01, c11, c20, c02 = \
        c00/9, c10/6, c01/6, c11/4, c20/6, c02/6
    deplx = (c11*c01-2*c10*c02)/(4*c20*c02-c11**2)
    deply = (c11*c10-2*c01*c20)/(4*c20*c02-c11**2)
    return deplx, deply, correl_crop

With Pythran 0.8.4.post0, g++ consumes too much memory (more than 12 Go) and finally crashes...
Same behavior with GCC 4.9 and 5.4...

Funny enough: we added getter computation of range values, so we can now infer that a[(i,j)] only involves positive index, so fast indexing can be used, but this was not implemented for generalized expression, so it created a big error, so big that it made... the compiler swap. Should be ok with #797


Suggestion: choose custom download name

Hello i use this great tool almost everyday and i think that it's the best gui for yt-dlp today available, but sometimes i need to give a different name to the output file before downloading it, do you think that it's possible to add the opportunity to choose a custom file name before clicking the download button?

Hey, glad to hear that you like this tool and thanks for the suggestion. I'll try to consider this in an upcoming update.

so the new version is out, and custom files names before downloading is not possible. thanks however for fixing downloading stats

https://twitter.com/i/status/1210978856628191233
https://twitter.com/i/status/1199121778058113026
can't download these because I can't select a new file name for the save....but also the errors show and then disappear.  And vividl says "No error log available."  so I can't see a whole history of my download/save errors, at least for the current session.

This has been implemented and the upcoming release will allow customizing download names. Thanks for your patience :)


Typo in readme: seamlessly not seemlessly

typo in readme: "seamlessly" not "seemlessly"

thank you. edited it.


Can not change the statistic values.

For dead by Daylight i try to change some values under the statistics tab and it says it stored it but then loads back the same as it was. Does anyone know why?

Are they flagged as increment only or protected?

true, not possible to change


Optimize images

It's nothing drastic, but I was able to losslessly optimize a number of files, which results in slightly smaller file sizes.

Thx!


[496] Features: Push notification, Incoming Transaction, Websockets

https://github.com/aeternity/superhero-wallet/blob/2518e011400cb3558216757ae8ae8f9856474dff/src/composables/websocket.ts#L165-L167
@ifaouibadi  Listening to each and every transaction in the Aeternity chain is not scalable and not correct. Hardly transactions are relevant to the account. Please make use of the Object payload type to listen to events for only required local accounts.
From Middleware documentation
Object, which takes a further field, target - can be any æternity entity. So you may subscribe to any æternity
 object type, and be sent all transactions which reference the object. For instance, if you have an oracle
ok_JcUaMCu9FzTwonCZkFE5BXsgxueagub9fVzywuQRDiCogTzse you may subscribe to this object and
be notified of any events which relate to it - presumable you would be interested in queries, to which you
would respond. Of course you can also subscribe to accounts, contracts, names, whatever you like.

https://github.com/aeternity/ae_mdw#supported-payloads
Perhaps we can discuss this in a meeting.

@ifaouibadi Listening to each and every transaction in the Aeternity chain is not scalable and not correct. Hardly transactions are relevant to the account. Please make use of the Object payload type to listen to events for only required local accounts.

fixed

I remember @thepiwo once said we will introduce a web socket for multisig backend too to detect new multisig vaults.
@ifaouibadi Please take this into consideration we should be able to use WebsocketClient.ts file even for other web socket connections.

I started work on websockets for the ga-multisig backend, but we never followed up with an actual plan on integration, so I abandoned the approach for the time https://github.com/aeternity/ga-multisig-backend/commit/515b598c47c541c2ef4739ae5f21acf625f8aaed
let's have a meeting to pick it up again

I started work on websockets for the ga-multisig backend, but we never followed up with an actual plan on integration, so I abandoned the approach for the time aeternity/ga-multisig-backend@515b598
let's have a meeting to pick it up again

This is important to optimize the performance of the SH wallet, especially to eliminate polling to the multisig backend. Let's talk soon ✌

Although the WebSocket methods are static it looks like there are many parallel connections that keep on getting added as we switch networks. If you observe the notifications are still received from the first connection made. Here is the screenshot

Improved network unsubscribe from previous when it's been changed

Notifications are received only from the active network. I believe notifications should be network agnostic with each notification indicating the network.

it's not a good idea to listen for all networks been added in the app as it can cause notifications duplicate or receiving unexpected notification from non active network

it's not a good idea to listen for all networks been added in the app as it can cause notifications duplicate or receiving unexpected notification from non active network

Fine.

I suggest to add a redirect from transaction notification to open this transaction details.


Feature/BoseFS2C

New address type BoseFS2C for two-component bosonic systems
Implementation of the two-component bosonic address type and Hamiltonians for studying 1D bosonic polaron and/or bosonic mixtures problems.
Exports in BitStringAddresses.jl

BoseFS2C{NA,NB,M,AA,AB} <: BosonicFockStateAddress <: BitStringAddressType: new address type that holds two BoseFS.

Exports in Hamiltonians.jl

TwoComponentBosonicHamiltonian{T} <: AbstractHamiltonian{T}: new subtype in AbstractHamiltonian;
BoseHubbardReal1D2C{T, BoseFS2C} <: TwoComponentBosonicHamiltonian{T}:  two-component 1D Bose-Hubbard Hamiltonian in real space;
BoseHubbardMom1D2C{T, BoseFS2C} <: TwoComponentBosonicHamiltonian{T}: two-component 1D Bose-Hubbard Hamiltonian in momentum space.

New testsets are added too for briefly testing the functionalities of BoseFS2C and the new Hamiltonians.

No idea why GitHub shows old commits from Joachim. They have been merged to develop branch already and are not part of this pull request.

No idea why GitHub shows old commits from Joachim. They have been merged to develop branch already and are not part of this pull request.

Have pulled recent changes from develop branch into this pull request.

Have pulled recent changes from develop branch into this pull request.

It is quite irritating to see these old commits. They don't seem to affect the actual code changes, because I see only three files changed. Nevertheless, would you mind opening a new (separate) pull request from the latest version of the BoseFS2C branch to develop? I would like to see whether this will look tidier.
Cheers!

It is quite irritating to see these old commits. They don't seem to affect the actual code changes, because I see only three files changed. Nevertheless, would you mind opening a new (separate) pull request from the latest version of the BoseFS2C branch to develop? I would like to see whether this will look tidier.
Cheers!

Close this pull request to submit a new one.

Close this pull request to submit a new one.

Reopen because the old commits will still be there even if I create a new pull request...

Reopen because the old commits will still be there even if I create a new pull request...

Ready for the final review before merge.

Ready for the final review before merge.


Thoughts on Changing the Way Project Magnitude is Calculated

 Issue by 3ullShark
Tuesday May 30, 2017 at 08:39 GMT
Originally opened as https://github.com/Erkan-Yilmaz/Gridcoin-tasks/issues/130

Here's an excerpt from a conversation we had on IRC about the decay of RAC being different on each project.
<BullShark> is it just me, or do some project RACS decay faster than others?
<jamezz> its not you
<jamezz> like dude i aint done collatz in months
<jamezz> aint done GPUgrid in 4 months
<pepperooo> nfs decays fast, gpugrid decays faster :(
<BullShark> Shouldn't we be calculating our own RAC then to make it fairer
<jamezz> its because every project has a different credit system

Would it be possible to change the way RAC is calculated so that we do the calculation in the NN code when we calculate the project magnitudes. We could do this by calculating the difference of the users total credit when they first registered their CPID in the NN to what they have now.
Some potential issue

They would not be rewarded for work that was done before they registered their CPID in the blockchain.
Would deleting beacons cause mag to drop to 0?

Is it fair the way it is now? Or should we look at changing it?

 Comment by grctest
Tuesday May 30, 2017 at 17:00 GMT

I like the current method of calculating magnitude on an individual project basis, proportional computation rewards (your_rac/team_rac).
If individual projects seem to have issues with their credit mechanisms it would be a good idea to contact their development team regarding your concerns.
Reason I think we should stick to using the project's RAC is that its a simple value that they already have, making it more abstract could be an additional hurdle for new users to understand how much they're capable of earning.
Aiding in development of the 4th gen credit mechanism could help prevent projects having to reinvent the wheel: https://boinc.berkeley.edu/dev/forum_thread.php?id=10953

 Comment by Foggyx420
Tuesday May 30, 2017 at 18:42 GMT

I agree about keeping it the same and each project admin tinkers with there own projects. A 4th credit system could be great however we may have similar situation as projects don't all just update to new system.

 Comment by tomasbrod
Tuesday May 30, 2017 at 22:12 GMT

I thing we should calculate the magnitude ourself. Calculate magnitude difference in total credits since last superblock. Why?  When projects go offline or project stats not updating, the mag calculation ceases tho be fair.

 Comment by Foggyx420
Tuesday May 30, 2017 at 22:58 GMT

so get the neural network to do all the calculations then? if projects are down that information would not be accurate either as boinc gathers that information from the projects and not in client. also would have to add code so the wallets would even keep an updated rac amount as it doesn't. and even then it can be manipulated as well still.

 Comment by 3ullShark
Tuesday May 30, 2017 at 22:59 GMT

Yeah, it doesn't seem fair that users can be rewarded for doing nothing just because they picked a project that takes ages for RAC to decay.
I agree with @tomasbrod if it's possible we should calculate mags ourselves and make it maybe make it decay down to 0 after 30 days. Then it does not matter which project you choose they will all be the same.
I don't think magnitude is hard to understand. We already use magnitude. It's actually harder to learn RAC and Magnitude. Especially if people are new to BOINC.

 Comment by Foggyx420
Tuesday May 30, 2017 at 23:53 GMT

@3ullShark
I hear you but I do have concerns.

How are you proposing exactly that the magnitude be calculated?
Who exactly is ourselves? the neural network? the person who owns the cpid?
Where is the RAC or MAG from the stats actually going to come from and where are you going to actually get these stats?
** We can't use the boinc client for the stats as it wouldn't be accurate either when project is down.
** Also if we did we would be adding more cpu usage as the wallet would have to keep either in sync or updated from the client_state.xml.
** At boot up the wallet checks client_state.xml for cpid information and if rac is over 100, etc but is not regularly updated least ive never seen as my list cpids been the same over days.
If the neural network is going to calculate all this where is it going to gather any kind of RAC or MAG information if we not using projects stats that are from the projects themselves? Remember the stats are not always perfect but they come from the project not a gridcoiner.
** Surely we can't allow a local client to calculate there own RAC or MAG as it can be manipulated unlike the neural network comes onto a consensus between the NN clients.
** Where will there be consensus on whether or not that information for RAC or MAG is correct?
If the neural network is going to do all these calculations it would be a over work of all the code as the network would have to also essentially come to a consensus which is like voting for popular stats.

Also there is a reason its 115000 * ((Your RAC / Team RAC) / WhitelistProject count)
Whitelist project count makes an even amount of magnitude spread across all projects meaning generally there is so many coins per day per project thou ur not neccessarily paid that day for that project. One project should not be favoured more then another and there not. Also Its not the NN fault that a project server is down and thus a greylist would help here however if the project is down then the NN should not allow rewarding for them during that time frame as there is no stats to fairly reward users to begin with so using the 1 sb stats before hand would be a disaster as well. yes it affects your mag when the projects are down and NN can't get that information.
Your RAC / Team RAC is because your in competition for the magnitude in each project. You bring more hardware to the table then your spending more money on hardware and electricity or if you make wise project choices for your hardware as well as other factors. It is a competition for the MAG and coins.
115000 Neural Multiplier keeps the Network Magnitude below 115000 Magnitude and that along with mag unit keeps the network steady in coin production for DPoR.
I'm greatly concerned about the what "ifs" in this situation. I understand people want to work towards change and thats one thing but you talking about redesigning gridcoin completely. Also there would have to be a consensus vote on this of course once you've laid out all the outcomes to this and complete direction to this.
I can't even begin to imagine the biggest hardest fork of gridcoin to ever happen. This kind sounding like a direction of other crypto coins leaning toward splitting into two different coins.
Not throwing your ideas under the bus at all but there is lack of direction and information to support the idea. Brainstorming is one thing but showing everything in complete documentation is more realistic then "what ifs"
To be honest its almost like a whole new coin being developed here and gridcoin is the base coin its derived from.


Example using ipv6?

Summary
Ability to use IPv6 public address.
Motivation
Because services like Fly.io mostly work well on ipv6 UDP.

+1 to this because I'm also trying to use fly.io.
In particular, I'd love to be able to have TURN listening on ipv4 and ipv6 UDP and TCP and have it relay for an SFU running on ipv6 UDP.

Unfortunately it requires the implementation of RFC 6156 to support IPv6.

@Sean-Der Is IPv6 support planned? Thanks in advance.

Seems like the TURN client currently does IPv4 lookups only:
https://github.com/pion/turn/blob/6880381c77d370eb8121fe7eb237b838ea8e4429/client.go#L95

See also #265


Fix utterances color scheme setting in the preferred mode

Reviously we set utterances' theme to preferred-color-scheme
under eureka theme. It means using the browser default color
scheme. It's fine if the default was light mode as we have
logic to switch dark mode. However, if the deafult was the dark
mode, utterances won't change its theme even if light mode is
explicity turned on for the website. This commit fixes it.

Do you mean that the browser default color scheme is dark and user switch to light mode using the button on the website? I've tested this use case and it seems normal.

Do you mean that the browser default color scheme is dark and user switch to light mode using the button on the website?

Yes. Probably due to different OS/browser settings I've on my side. The switching issue seems to be fixed by ccd5695125f1516db25424dc723eaa9fc6b3e87a; but with a defult browser light theme, I still see the following when first opening https://www.wangchucheng.com/en/posts/flutter-state-management/ (assuming you've incorporated the change there):

I still see the following when first opening

Wow, my browser works well. Could you provide the log in console which can show this issue? I will see what we can do.

I only see warnings about "Could not find the language 'dart', did you forget to load/include a language module?" in the console log. Nothing pointing to this issue.
Actually playing it a bit more, I think the issue appears when you first manually switches the website's theme to light via the toggle, and close the webpage. And then reopen it. It remembers the previous theme settings and will try to reuse that IIUC. Then you'll likely see the issue.

Thanks! I can reproduce this issue, I will check this later.

发件人: Lei Zhang @.>
发送时间: Wednesday, July 7, 2021 11:39:17 PM
收件人: wangchucheng/hugo-eureka @.>
抄送: C. Wang @.>; Comment @.>
主题: Re: [wangchucheng/hugo-eureka] Fix utterances color scheme setting in the preferred mode (#122)
I only see warnings about "Could not find the language 'dart', did you forget to load/include a language module?" in the console log. Nothing pointing to this issue.
Actually playing it a bit more, I think the issue appears when you first manually switches the website's theme to light via the toggle, and close the webpage. And then reopen it. It remembers the previous theme settings and will try to reuse that IIUC. Then you'll likely see the issue.
―
You are receiving this because you commented.
Reply to this email directly, view it on GitHubhttps://github.com/wangchucheng/hugo-eureka/pull/122#issuecomment-875710643, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AJEXGNFRPFTVPCHVZCTBHPTTWRYKLANCNFSM47I5YLLQ.

Awesome, thanks!


Remove version: apankrat.byenow version <IP_ADDRESS>

[Automated] It returns code over 400 in all urls
Microsoft Reviewers: Open in CodeFlow

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 797)

Hi @Exorcism0666,
This manifest has the highest version number for this package. Is it available from another location? (This might be in error if the version is switching from semantic to string, or string to semantic.)
(Automated response - build 795.)

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 807)

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 807)

Hi @Exorcism0666,
This manifest has the highest version number for this package. Is it available from another location? (This might be in error if the version is switching from semantic to string, or string to semantic.)
(Automated response - build 808.)

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 812)

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 813)

URL: https://github.com/apankrat/byenow/releases/download/byenow-<IP_ADDRESS>/byenow-0.11.zip
Response status code does not indicate success: 404 (Not Found).
(Automated message - build 813)

For moderators
This account is automated by Github Actions and the source code was created by CoolPlayLin. If you have any questions about any pull request, don't hesitate to ping Exorcism0666.

[!warning]
With the new version of Komac, numerous Pull Requests may contain bugs or unwanted elements. If such issues arise, please tag russellbanks so that these bugs can be addressed in future updates.

[!important]
Please carefully review these Pull Request before merging. If it is a Pull Request for removing incorrect content and the URLs are issue, free upon manual checking, please close this Pull Request directly. (It is best to inform Exorcism0666 of the closure.)

良い一日をお過ごしください！


Added spawn.datastore.type: MYSQL2 to use mysql connectorJ.

Adding mysql connectorJ (in addition to current drizzle driver) option to be compatible with oci mysql version.

Have to use a mysql driver that is compatible with oci mysql8.0. Current drizzle driver does not have newer versions supporting mysql8.0, therefore switched to connectorJ. Have not tried compatibility with non-oci mysql version, since drizzle works perfectly with older mysql for now.


Help pages on skpm.io broken

Just letting you know that all links on https://skpm.io/help/ are currently broken and return a 404.

I saw that yes, haven't have time to look into it yet
duplicate #157


Auto-scrolling in web-UI should only be active if user hasn't scrolled up.

Every time a new message appears in one of the feeds/logs of the web UI, we scroll to the bottom. However, if a user wants to read a certain part during the run and has scrolled up, we shouldn't scroll down.

[!tip]
General tips for frontend contributors:

Please open a separate PR for every issue that you address. Try to keep PRs small and concise (if you're thinking of bigger refactorings, please let us know)
In the web UI under "Extra Settings" there's a switch to run without needing API credits for testing purposes
Please install pre-commit (pip install pre-commit && pre-commit install) to ensure that your submission follows the styling guidelines

+1

@klieret we can close this, right?

Thanks for the reminder, yes

ty


Uasp force for the dev branch

An implementation to address:
#2488

Was able to merge the other PR into the dev branch - closing this one


Siri (Italian) - homebridge - iOs10

I installed public beta of iOs 10.
All work fine.
The problem is only Siri.
When ask in Italian "Siri tun on light Kitchen" "Siri accendi luce cucina", Siri don't is able to recognize message.
With iOs 9 worled fine.
There are  ohters with this problem?
Siri in English works fine?
Bye
Sandro

Same problem.

same here (ita)

Hello guys! To you it works iOS 10 - Homebridge? When there were beta Homebridge working, as soon as I installed the final or even RC funzioan no nothing .. :(
Ciao ragazzi! A voi funziona iOS 10 - Homebridge? Quando c'erano le beta homebridge funzionava, appena ho installato la definitiva o anche RC non funzioan più nulla.. :(

Personally running on iOS 10 Without any issues
Sent from my iPhone

On Sep 14, 2016, at 7:57 AM, steppige<EMAIL_ADDRESS>wrote:
Hello guys! To you it works iOS 10 - Homebridge? When there were beta Homebridge working, as soon as I installed the final or even RC funzioan no nothing .. :(
Ciao ragazzi! A voi funziona iOS 10 - Homebridge? Quando c'erano le beta homebridge funzionava, appena ho installato la definitiva o anche RC non funzioan più nulla.. :(
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub, or mute the thread.

Hello!! thanks for the reply! To you work with the beta version of iOS 10? Or with the final one? I have installed eDomoticz ...

I am also running iOS 10.0.1 (14A403) without any homebridge issues.

@steppige i'm running the final GM.
I have only ran the GM, So i can't give you any personal experience with other beta versions
Sent from my iPhone

On Sep 14, 2016, at 8:01 AM, gijoecool<EMAIL_ADDRESS>wrote:
I am also running iOS 10.0.1 (14A403) without any homebridge issues.
―
You are receiving this because you commented.
Reply to this email directly, view it on GitHub, or mute the thread.


Add Matchit support for Java statements

Adds matchit support for if/else, do/while, and try/catch/finally statements in Java.

Great, thank you!

FYI the new unit tests don't test jumping from try to catch to finally, might be worth testing too. I can add the test if you want, but for you it might be faster to just add it than deal with another PR.

Thank you, added this one as well


Add guacscanner

🗣 Description
This pull request modifies this Ansible role to allow for the addition of cisagov/guacscanner-docker to cisagov/guacamole-composition.
💭 Motivation and context
The functionality of cisagov/guacscanner-docker is necessary for the RPT UAT taking place in the COOL.
See also:

cisagov/guacscanner/pull/1
cisagov/guacscanner-docker#2
cisagov/guacamole-composition#30
cisagov/guacamole-packer#53
cisagov/cool-assessment-terraform#149

🧪 Testing
I have applied these changes to env0 in our staging COOL environment and verified that they function as expected.
✅ Checklist

[x] This PR has an informative and human-readable title.
[x] Changes are limited to a single goal - eschew scope creep!
[ ] All future TODOs are captured in issues, which are referenced in code comments.
[x] All relevant type-of-change labels have been added.
[x] I have read the CONTRIBUTING document.
[x] These code changes follow cisagov code standards.
[x] All relevant repo and/or project documentation has been updated to reflect the changes in this PR.
[x] Tests have been added and/or modified to cover the changes in this PR.
[x] All new and existing tests pass.
[ ] Revert back to the main branch of cisagov/guacamole-composition once cisagov/guacamole-composition#30 is merged.

Excellent work. I had a question about Debian 9-fu since it is not a supported platform for this role. I also went a bit too far down the rabbit hole and found that the Dockerfile for guacamole/guacamole is installing a PostgreSQL JDBC driver that's over six years old. We should probably make an issue in the Apache JIRA to either fix it ourselves with a PR or at least get some attention on it I guess. The closest issue I could find was someone else having a problem in GUACAMOLE-1433 and another user mentions the JDBC driver needing an update.

Good idea.  The old driver has already bitten us.  See cisagov/guacamole-composition@8197d0e804a145a8dac67078b1e0a14c68e8e99d for more details.


Add build pieces for Google Test

copied from https://groups.google.com/forum/?utm_medium=email&utm_source=footer#!topic/tesseract-dev/VLwwGFKtdPA

Ray | Jul 21
I have the OK to "throw the tests over the wall" already. ie provide them in a non-working form.
There are actually very few copyrighted images that would need to be replaced. Most of the tests run on synthetic data, existing test data, or don't require images.

If someone would put together the build pieces necessary to build and run an empty test (using Google test), then I will port at least one example, and then push the rest out.

Ray | Jul 13
There are the order of 50 tests, some of which complete in <1s. Some take ~10mins to run, but very few. You could probably run the whole lot on a single machine in about an hour.

Jul 13 Jeff Breidenbach
First step is to make a working but empty test for Tesseract under this
framework. Once that exemplar is in place, we can migrate the existing
tests that currently run inside Google.

Some tests for tesseract using java as part of Quan's tess4j project

https://github.com/nguyenq/tess4j/tree/master/src/test/java/net/sourceforge/tess4j

https://stackoverflow.com/questions/35998856/how-can-i-use-google-test-with-my-project-that-builds-via-autotools

@theraysmith, pull request #1051 adds the required pieces for GoogleTest.

Thanks!
I will port one test, and then determine how easy it is to write portable
tests that work both inside Google and outside.
On Mon, Jul 24, 2017 at 10:51 AM, Stefan Weil<EMAIL_ADDRESS>wrote:

@theraysmith https://github.com/theraysmith, pull request #1051
https://github.com/tesseract-ocr/tesseract/pull/1051 adds the required
pieces for GoogleTest.
—
You are receiving this because you were mentioned.
Reply to this email directly, view it on GitHub
https://github.com/tesseract-ocr/tesseract/issues/1050#issuecomment-317501971,
or mute the thread
https://github.com/notifications/unsubscribe-auth/AL056a53GU5T3Bo5ETJ3KJ6XOFSsIQ5Fks5sRNmPgaJpZM4Og8V1
.

--
Ray.

PR with Makefile.am for unittest
https://github.com/tesseract-ocr/tesseract/pull/1088


Event Hub components should expose initialOffset in manifest to control starting position (like Kafka)

Describe the feature

When an PubSub or Binding component starts up with an empty/new checkpoint store, it starts at the latest position in the backing store. There should be an option to allow starting at the earliest.
Kafka supports this with initialOffset

See https://github.com/Azure/azure-event-hubs-go/blob/3dc3aeb6144e5d79a70d32f023f9d5990f81569d/persist/checkpoint.go#L29 in the azure SDK go SDK for the mechanism to support this.
Release Note

RELEASE NOTE:

Bump.

Bump


total-functions/no-array-destructuring false positive with string literals

The code that reproduces the error:
interface Response {
  'error-code': number
}

function getError(
  { 'error-code': errorCode }: Response // <-- total-functions/no-array-destructuring
): string {
  return errorCode.toString()
}

There should be no total-functions/no-array-destructuring in this code.
For normal properties (not quoted), everything works fine.
I am using the library version:

<EMAIL_ADDRESS><EMAIL_ADDRESS>
Eslint config:
extends: [
  ...,
  'plugin:total-functions/recommended',
],
plugins: [
  ...,
  'total-functions'
]

Without rewriting the rules for eslint-plugin-total-functions.

Thanks for reporting @awwit . I'll fix this soon. In the meantime, PRs encouraged :P

Released as v1.39.0. Please let me know how you go @awwit

@danielnixon I confirm that the error has been fixed.
Thanks! =)

Yikes. Hold tight.

@awwit let me know how v1.40.0 goes. Sorry for the churn.


Visualization

Hi,
Very interesting work
I wanted to know if there was a way to produce the graphs given a video like shown in the README
Thank you

Thanks for your interest, we use this tool to visualize the learned graphs.

Thank you!
Is there any guidelines on how to use the tool on your data?
Best

No such guideline :), and you need to format the learned adj matrices to adapt to the tool. Basically, you need a JS file like this:
var nodes = [
    {
        "id": 0,
        "label": "V0",
        "group": 3
    },
    {
        "id": 1,
        "label": "C1",
        "group": 2
    },
}
var edges = [
    {
        "from": 0,
        "to": 1,
        "value": "0.2441"
    },
    {
        "from": 1,
        "to": 8,
        "value": "1.0000"
    },
}


Query.removed() is always empty

Bevy version
rev 72b2fc9
Operating system & version
Windows 10
What you did
Here is a simple example:
use bevy::prelude::*;

struct MyComponentOne;
struct MyComponentTwo;

fn main() {
    App::build()
        .add_plugins(DefaultPlugins)
        .add_system(spawn.system())
        .add_system(remove.system())
        .add_system(print.system())
        .run();
}

fn spawn(
    commands: &mut Commands,
) {
    commands.spawn((MyComponentOne, MyComponentTwo));
}

fn remove(
    commands: &mut Commands,
    query: Query<Entity>,
) {
    for entity in query.iter() {
        commands.remove_one::<MyComponentOne>(entity);
    }
}

fn print(
    query: Query<Entity>,
) {
    for thing in query.removed::<MyComponentOne>() {
        println!("{:?}", thing);
    }
}

What you expected to happen
To see something printed in the console.
What actually happened
Nothing is printed.
Additional information
It's been suggested in the Discord that perhaps Query.removed() is not meant to be used in this way. If so, how is it meant to be used? We should add this information to the documentation of the function.
If it's not meant to be used in this way, is there some best practice yet for how to check for removed components in Bevy? This would be great knowledge to have in the community as it's been very hard to find any information about this.

Your example above doesn't work for a couple of reasons:

all three systems were added to the same (default) stage, so they all run in parallel. the "remove" system should always run after "spawn" and the "print" system must run after components are actually removed.
remove_one and spawn are Commands, which are always applied after a stage has finished executing. this means that even if the parallel execution order happens to line up, query.removed::() will return nothing because the remove commands haven't been applied yet.

You can fix it by adding your systems to different stages to force execution order. This prints the expected output:
.add_system_to_stage(stage::PRE_UPDATE, spawn.system())
.add_system_to_stage(stage::UPDATE, remove.system())
.add_system_to_stage(stage::POST_UPDATE, print.system())

(note that you should create your own stages after stage::UPDATE for scenarios like this instead of adding them to PRE/POST update)

I do want to give optional, more fine grained control over execution order within a stage, but for now using stages is the recommended approach.
I'm closing this issue as it isn't a bug / the system is behaving as intended.

I'd say it's at least something that could be documented better,
possibly by fixing this example and adding it to the repo. I hit this
issue too, and have been all over the repo with grep looking for
instances of removed. It wasn't clear to me that the stage checking
for removal needed to be after the stage where the removal happened. It
also wasn't clear whether removal needed Changed or some other query
operator, or whether it worked on any component regardless of its
presence in the query. This example seems to clarify both questions for
me. Might be nice to submit it as a PR that closes this issue once
merged. It makes sense once you have a deeper understanding of how the
pieces interact, but the learning curve to get there is a bit steep. The
examples are great at shallowing it out a bit. :)
Thanks.

Yeah an example illustrating this issue is a good idea. Feel free to pr it :smile:

I'll get on a PR for the example, and then the documentation, unless someone else gets to it first.
Just announcing here to avoid any duplicated work.


[Feature] Plugin load priority

Is your feature request related to a problem? Please describe.
In some cases, the order plugins are loaded in breaks them.
Describe the feature you'd like
A good way of solving this would be an option in the plugin settings to choose the order. For example, grouping the plugins into high, medium, and low priority. By default, every plugin would be on medium priority and everything would load like always, but if the load order messes up something, you would be able to put the conflicting plugins into high or low priority, which would be loaded first or last respectively. It would probably be best to put a dropdown list next to the load/unload toggle switch in the plugin settings that allows you to choose the priority.

IIRC, the plugins are loaded in the order of their filenames.
What you could do in the meantime is rename the files, but I think it would break plugin updates.

oh, well, i did not know that. gonna try that out to see if that fixes the problem i'm personally having. still, quite a dirty workaround.

What plugins are having issues like that? This should not be the case.

for example, if "ServerFolders" by DevilBro loads after "Game Activity Toggle" by Egrodo, the button that GAT adds gets removed and only re-appears after manually reloading the plugin. Sure, the plugin devs could fix that themselves, but having configurable priority would be a catch-all solution for everything.

There is a bug with Game Activity Toggle where it doesn't add the button. The dev is working on a fix, it is not a bug exclusive to ServerFolders. GAT should load before ServerFolders anyways, just depends on when the code written by the plugin dev decides to execute, that is not something that BD can control.
I don't really intend to add loading priority, that's an unnecessarily complicated system to add for a very very rare issue.

Plus something like load priority sounds like a band-aid solution.


The positional embedding is fixed during pretraining but learnable during finetuning?

For instance, PIXELForSequenceClassification is used for finetuning GLUE. Its underlining ViT implementation is ViTModel instead of PIXELModel. It seems that PIXELModel follows MAE and uses fixed positional embeddings, while ViTModel follows the original ViT paper and uses learnable positional embeddings.

Hi, thanks for pointing this out! You're right the PIXEL finetuning models don't have their position embeddings fixed. I think behavior is the same in the original ViT-MAE (https://github.com/facebookresearch/mae).
I don't think this affects performance (noticeably), considering that the position embeddings are still loaded from the pretrained model, so they will be finetuned from sinusoids rather than from scratch. Since they are in the lowest layer of the model, the gradients for the position embeddings are also relatively small, so they don't change too much from their original sinusoidal form.
You could nevertheless try setting the requires_grad attribute of the position embeddigns to False after loading a finetuning model like PIXELForSequenceClassification and check whether that makes a difference downstream, or alternatively use a PIXELModel instead of ViTModel as the base model for finetuning. Note that in the latter case you would have to make sure that the mask_ratio is set to 0.0 as otherwise masking will be applied on the embeddings (this is the main reason why we use ViTModel actually).

Thanks for your prompt reply! That is really helpful! I will close this issue.


Mask Eopatches according to polygon

I have multiple eopatches created based on a river shapefile. Currently I am trying to mask these patches according to the shapefile so that if: pixels of the eopatch are located within the river they still contain the information and ifelse: 0 or no information. Is this possible?

Yes, this is possible using the VectorToRaster task, where for each eopatch you can burn your river vector data into a binary mask to be used to filter data.


When Discord bot token and channels are set, the server wont start up

The console is just left with press a key to continue

Sounds like the app crashes on startup
can you send me the manager log file?

closing as it is really old


Release sdk/resourcemanager/extendedlocation/armextendedlocation v0.1.0

…1.0 generation from spec commit: 466db813b82d7623b6882850636df09e0f98d831

[ ] The purpose of this PR is explained in this or a referenced issue.
[ ] The PR does not update generated files.

These files are managed by the codegen framework at Azure/autorest.go.

[ ] Tests are included and/or updated for code changes.
[ ] Updates to CHANGELOG.md are included.
[ ] MIT license headers are included in each file.

/check-enforcer override


Enable mergify for bots

This mergify configuration enables automatic merging of PRs, if:

Build passes
Target branch is master
The PR is tagged with "dependencies"

I modified the renovate configuration to assign the correct labels to PRs. This should help reduce your workload, and we let the bots to the job, until they all rise up and claim humanity :)


Fix _GET_X_LPARAM and _GET_Y_LPARAM from api_windows.go

Negative values should not disappear. For example, if _HIWORD returns 63000, then _GET_Y_LPARAM should return -2536, but now it returns 63000. This fixes that.
Here's how it's defined in windowsx.h:
#define GET_X_LPARAM(lp)                        ((int)(short)LOWORD(lp))
#define GET_Y_LPARAM(lp)                        ((int)(short)HIWORD(lp))

Did the original code cause an actual issue? I'd like to know how serious the issue is. If this is serious, I'll back-port the fix to the stable branch. Thanks,

These functions (macros) are used at WM_MOUSEMOVE (a mouse move) and WM_MOVE (a window move). For a mouse move, the cursor position is a relative position in a window so these cannot be a negative value. For a window move, I have confirmed a too-big value can be passed to inputWindowPos. However, Ebitengine doesn't use any callbacks for window position changes. Therefore, I think there is no serious issue, which means we don't have to backport the fix to the stable branch.

Did the original code cause an actual issue?

Yes, but not while using Ebiten engine. I'm using yours goglfw (together with OpenGL (no CGO!) with my own ImGui Go bindings, and encountered an issue where cursorPosCallback is called with wrong cursor position parameters (huge values) when moving or resizing secondary viewports up or to the left.


Arguments not identified at all (cc3200/toniebox)

I am currently try to use your lib on one of my projects.
First of all:
I am using it on a ti cc3200 (toniebox) with energia (arduino fork).
But it doesn't seem to be able to read my input. (I am sending it as string over the parse interface).
Definition
    cmdI2C = cli.addCmd("i2c");
    cmdI2C.setDescription("Access I2C");
    cmdI2C.addFlagArg("r/ead");
    cmdI2C.addFlagArg("w/rite");
    cmdI2C.addArg("a/ddress");
    cmdI2C.addArg("r/egister");
    cmdI2C.addArg("v/alue", "");
    cmdI2C.addArg("l/ength", "");
    cmdI2C.addArg("o/utput", "b");

    cmdHelp = cli.addCmd("h/elp");
    cmdHelp.setDescription("Get help!");

Input
i2c read -address 0x10 -register 0b10000000
Output
ERROR: Unknown argument at command 'i2c' at 'read'
Did you mean "i2c [-r/ead] [-w/rite] -a/ddress <value> -r/egister <value> [-v/alue <value>] [-l/ength <value>] [-o/utput <b>]
Access I2C"?

Full code:
https://github.com/toniebox-reverse-engineering/hackiebox_cfw/blob/develop/BoxCLI.cpp

Hi, shouldn't it be i2c -read -address 0x10 -register 0b10000000? There is  a - missing before read.

facepalm
Thanks for the hint.


Bug: rename "login" to "logIn"

Function names should be in the imperative form, which means they should be verbs. Since the verb for the action is a compound verb, it's two words: log and in. In camelCase, this would be logIn (just as it is with its complementary function logOut), not login (which is the noun form).

 "Logout" is not a verb. Thanks for pointing this out :slightly_smiling_face:
Agreed we should fix this (with backwards compatibility), so any PR is welcome.

Sure, I'll get going on a PR for this.

@getsnoopy Did you start working on a little PR for this?

@sebastianvitterso I haven't yet; I'll probably get to it next week.


[PrintAsObjC] Handle typealiases to non-type Clang declarations.

For historic reasons, Clang's representation of an Objective-C class
declaration ObjCInterfaceDecl) and compatibility alias
(ObjCCompatibleAliasDecl) are not actually Clang TypeDecl nodes. Cope
with this in Objective-C printing, fixing rdar://problem/32308192.

@swift-ci please smoke test and merge


Feature request: private key for encryption and public key for decryption

Hi @z-hao-wang , I'm very interesting in and appreciate your project.
But my scenario is:
I need encrypt some data on server, and decrypt them on client.
So, the private key should be kept on server for encryption, and the public key will be distribute to any client for decryption.
But I found that, This project can only use public key for encryption and private key for decryption. It will not fit this scenario. Could you help to support this case.

Hi @xing-zheng,
Looks like you need to store the private key on client.
RSA works like this, you keep private key to yourself, and open you public key to everyone. So any one can send you encrypted message using public key, however this message can only be decrypted using your private key.
My solution is every client has it's own private key and public key. I store each client's public key in server database. Then if I need to send an encrypted message to a client, I need to look up its public key, then encrypt using its public key.
Another note is RSA can only encrypt certain size data, if you need to encrypt larger data, you must use something like AES. I use RSA to encrypt the AES cipher(password), then use AES to encrypt/decrypt the actual data blob.

Thank you @z-hao-wang .
Thank you for your explanation.
May be I haven't declare my requirement clearly.
I should keep the private key on server, and use server to broadcast encrypted message to any client which hold the public key.

@xing-zheng
From my understanding, that is not possible with RSA.
Encryption using private key is only for signing purpose. The message is still not encrypted but client can verify using public key that this message is truly from the server who holds the private key.
see http://crypto.stackexchange.com/questions/2123/rsa-encryption-with-private-key-and-decryption-with-a-public-key
If you want to broadcast the same encrypted message, you would have to go with symmetric encryption method.
Or using the method I purposed, you just encrypt the same message for each client (each client would get a special encrypted message using its public key). Then client use its private key to decrypt the message.

Thank you @z-hao-wang .
I've change the code on my side, Then it works for my scenario.
xing-zheng/react-native-rsa@6ebeaaee2379c51aaf119d9dac1180c4e1166175

@xing-zheng: How can you use this library with .pem file?


Wait for powerdown

+1

+1

(Note, still need to test this change more thoroughly)

test failed on DC5 want_sles12=1 :
onadmin_rebootcloud
Warning: Permanently added '<IP_ADDRESS>' (ECDSA) to the list of known hosts.
Waiting for: node d52-54-01-77-77-01.g15.cloud.suse.de to go down
.......................
Waiting for: unknown process
....................................................................................................
Error: Waiting for 'unknown process' timed out.
This check was used: crowbar node_state status | grep d52-54-01-77-77-01.g15.cloud.suse.de | grep -q 'Power'

return 11

which is probably the bug that is tested for...

@bmwiedemann yeah, it hangs for 2 minutes on shutdown due to the cinder/tgt changes .. I've extended the timeout.

@bmwiedemann Can trigger the test again pls?


weather: move Env and SenseEnv out of physic

physic is now just a bucket for units.

The next step is to rename physic to metric or measurement.

Another option is to inverse the names. Make the package named environment, and change the struct Env to Weather. Then it'd make it more sensible to have a struct Position for geospatial position sensing.

Thinking about it more; I'll drop this change, as all in all, it's not really worth the churn.


Unable to watch Discord Streaming

There seems to be a problem with the unlocker intervening with discord's streaming. After enabling the unlocker, I am unable to watch my friend's share screen smoothly. Discord will only respond if I click on it, but otherwise it will hang. This issue only happens with discord and when the FPS unlock is running.
Has anyone encountered similar issue

Sounds like spec issues

If you mean the PC spec, I'm running a Ryzen 9 7900, RX 6700XT with 64GB of ram. The game is loaded on a Samsung 980 Pro.
I do not think the PC is the issue.

I've also tried with other games like Assetto Corsa, PUGB and Planet Coaster and those game do not have issues with the stream

I could not replicate the issue. Can you try recoding it?


new attributes for ogit/ServiceManagement/Ticket

In Ticket we currently have optional attributes:
ogit/ServiceManagement/externalTicketSource
ogit/ServiceManagement/externalTicketId
which could correspond to an external source of the tool (e.g. BMC) we are importing the ticket from (e.g. Ticket coming to BMC from ATOS)
we also need optional attributes:
ogit/ServiceManagement/ticketSource
ogit/ServiceManagement/ticketId
which would be identifying the tool (e.g. BMC) and the id of the ticket in it (in BMC)
Tks.

Descriptions:
ticketSource:
Source of the ticket, e.g. BMC
If a ticket was imported/duplicated from another system this attribute will define the ticket ID from the source system.
ticketId:
ID of the ticket in the source system, e.g. the ID of the ticket in BMC

Should the connection be "ogit/Software/Application -> ogit/creates ogit/ServiceManagement/Ticket" or
"ogit/Software/Application -> ogit/creates ogit/ServiceManagement/Incident" because of talking of entity Incident?

postponed


Array stats tests failing on main

I've seen this on several PRs now:
FAILED dask/array/tests/test_stats.py::test_two[power_divergence-kwargs7] - V...
FAILED dask/array/tests/test_stats.py::test_two[power_divergence-kwargs8] - V...
FAILED dask/array/tests/test_stats.py::test_two[chisquare-kwargs4] - ValueErr...
FAILED dask/array/tests/test_stats.py::test_two[power_divergence-kwargs5] - V...
FAILED dask/array/tests/test_stats.py::test_two[power_divergence-kwargs6] - V...

ref
My guess is that a new version of scipy was released, but I'm checking now.

I can reproduce locally and it seems like this is an intentional change https://github.com/scipy/scipy/issues/14298


NameError Exception: uninitialized constant Aws::CF

I'm getting this error when trying to use the aws_signer option according to the documentation in the README.
Here's my carrierwave config file:
CarrierWave.configure do |config|
  config.aws_credentials = {
    access_key_id:      'my-key',
    secret_access_key:  'my-secret',
    region:             'us-east-1'
  }

  config.asset_host = 'my-cloudfront-endpoint'

  config.aws_bucket = 'my-bucket'

  config.aws_signer = -> (unsigned_url, options) { Aws::CF::Signer.sign_url(unsigned_url, options) }
end

I tried switching to use
config.aws_signer = -> (unsigned_url, options) { Aws::CloudFront::UrlSigner.new.sign_url(unsigned_url, options) }

as I was troubleshooting but it gave me an error about a missing key_pair_id (documentation here )
Any pointers are appreciated 😄

As noted in #70, you'll need to use the aws_cf_signer gem. This has come up enough that I think I need to  make a note in the README.

Yep that would be a good idea. Thanks!

Documentation has been updated!