mirror of
				https://github.com/paboyle/Grid.git
				synced 2025-10-29 19:14:33 +00:00 
			
		
		
		
	
		
			
				
	
	
		
			70 lines
		
	
	
		
			2.1 KiB
		
	
	
	
		
			Plaintext
		
	
	
	
	
	
			
		
		
	
	
			70 lines
		
	
	
		
			2.1 KiB
		
	
	
	
		
			Plaintext
		
	
	
	
	
	
| * Sign off 2+1f HMC with Hasenbush and strange RHMC
 | |
| 
 | |
| - Wilson plaquette cross checked against CPS and literature GwilsonFnone
 | |
| - Timesteps matched
 | |
| 
 | |
| - Use 16^3x32 
 | |
| 
 | |
| ********************************************************************
 | |
| * From previous CPS runs:
 | |
| ********************************************************************
 | |
| 
 | |
| Strange (m=0.04)  has eigenspan 
 | |
| **** 
 | |
| 16^3 done as 1+1+1 with separate PV's. 
 | |
| /dirac1/archive/QCDOC/host/QCDDWF/DWF/2+1f/16nt32/IWASAKI/b2.13/ls16/M1_8/ms0.04/mu0.01/rhmc_multitimescale/evol5/work
 | |
| ****
 | |
| 2+1f 16^3  - [ 4e^-4, 2.42 ]    for strange
 | |
| 
 | |
| ****
 | |
| 24^3 done as 1+1+1 at strange, and single quotient https://arxiv.org/pdf/0804.0473.pdf Eq 83,
 | |
| ****
 | |
| double lambda_low =   4.0000000000000002e-04 <- strange
 | |
| double lambda_low =   1.0000000000000000e-02 <- pauli villars
 | |
| And high = 2.5
 | |
| 
 | |
| Array bsn_mass[3] = { 
 | |
| double bsn_mass[0] =   1.0000000000000000e+00
 | |
| double bsn_mass[1] =   1.0000000000000000e+00
 | |
| double bsn_mass[2] =   1.0000000000000000e+00
 | |
| }
 | |
| Array frm_mass[3] = { 
 | |
| double frm_mass[0] =   4.0000000000000001e-02
 | |
| double frm_mass[1] =   4.0000000000000001e-02
 | |
| double frm_mass[2] =   4.0000000000000001e-02
 | |
| }
 | |
| 
 | |
| ***
 | |
| 32^3 
 | |
| /dirac1/archive/QCDOC/host/QCDDWF/DWF/2+1f/32nt64/IWASAKI/b2.25/ls16/M1_8/ms0.03/mu0.004/evol6/work
 | |
| ***
 | |
| Similar det scheme
 | |
| double lambda_low =   4.0000000000000002e-04
 | |
| double lambda_low =   1.0000000000000000e-02
 | |
| 
 | |
| Array bsn_mass[3] = { 
 | |
| double bsn_mass[0] =   1.0000000000000000e+00
 | |
| double bsn_mass[1] =   1.0000000000000000e+00
 | |
| double bsn_mass[2] =   1.0000000000000000e+00
 | |
| }
 | |
| Array frm_mass[3] = { 
 | |
| double frm_mass[0] =   3.0000000000000002e-02
 | |
| double frm_mass[1] =   3.0000000000000002e-02
 | |
| double frm_mass[2] =   3.0000000000000002e-02
 | |
| }
 | |
| 
 | |
| ********************************************************************
 | |
| * Grid: Power method bounds check
 | |
| ********************************************************************
 | |
| - Finding largest eigenvalue approx 25 not 2.5
 | |
| - Conventions:
 | |
| 
 | |
| Grid MpcDagMpc based on:
 | |
| 
 | |
|    (Moo-Moe Mee^-1 Meo)^dag(Moo-Moe Mee^-1 Meo)
 | |
| 
 | |
| - with  Moo = 5-M5 = 3.2
 | |
| - CPS use(d) Moo = 1
 | |
| - Eigenrange in Grid is 3.2^2 rescaled so factor of 10 accounted for
 | |
| 
 |