Last updated on 2023-09-29 07:55:15 CEST.
| Flavor | Version | Tinstall | Tcheck | Ttotal | Status | Flags | 
|---|---|---|---|---|---|---|
| r-devel-linux-x86_64-debian-clang | 2.0.3 | 32.45 | 306.84 | 339.29 | NOTE | |
| r-devel-linux-x86_64-debian-gcc | 2.0.3 | 21.20 | 226.15 | 247.35 | NOTE | |
| r-devel-linux-x86_64-fedora-clang | 2.0.3 | 403.14 | OK | |||
| r-devel-linux-x86_64-fedora-gcc | 2.0.3 | 999.44 | ERROR | |||
| r-devel-windows-x86_64 | 2.0.3 | 21.00 | 208.00 | 229.00 | NOTE | |
| r-patched-linux-x86_64 | 2.0.3 | 26.93 | 290.63 | 317.56 | NOTE | |
| r-release-linux-x86_64 | 2.0.3 | 24.24 | 298.60 | 322.84 | NOTE | |
| r-release-macos-arm64 | 2.0.3 | 105.00 | NOTE | |||
| r-release-macos-x86_64 | 2.0.3 | 351.00 | NOTE | |||
| r-release-windows-x86_64 | 2.0.3 | 30.00 | 283.00 | 313.00 | NOTE | |
| r-oldrel-macos-arm64 | 2.0.3 | 107.00 | NOTE | |||
| r-oldrel-macos-x86_64 | 2.0.3 | 134.00 | NOTE | |||
| r-oldrel-windows-x86_64 | 2.0.3 | 29.00 | 241.00 | 270.00 | NOTE | 
Version: 2.0.3
Check: package dependencies
Result: NOTE
    Package which this enhances but not available for checking: ‘INLA’
Flavors: r-devel-linux-x86_64-debian-clang, r-devel-linux-x86_64-debian-gcc, r-devel-windows-x86_64, r-patched-linux-x86_64, r-release-linux-x86_64, r-release-macos-arm64, r-release-macos-x86_64, r-release-windows-x86_64, r-oldrel-macos-arm64, r-oldrel-macos-x86_64, r-oldrel-windows-x86_64
Version: 2.0.3
Check: re-building of vignette outputs
Result: ERROR
    Error(s) in re-building vignettes:
    --- re-building ‘bym.Rnw’ using knitr
    
    
    *************************************************************************** 
    CONTAINS Runtime Modules of Parallel Sparse Linear Solver PARDISO Vers. 7.2
    CUSTOMIZED FOR THE R-INLA PACKAGE WHICH SOLVES A LARGE CLASS OF STATISTICAL
    MODELS USING THE INLA APPROACH.                                            
    Copyright Universita della Svizzera italiana 2000-2020 All Rights Reserved.
    
    No PARDISO license file found.  Please see `http://www.pardiso-project.org/r-inla
    where to place the license file pardiso.lic 
    *************************************************************************** 
    	Read ntt 2 1 with max.threads 24
    	Found num.threads = 2:1 max_threads = 2
    
    	882c43d68fb015d72b699bcc458cf277dfe110c2 - Mon Apr 24 16:19:16 2023 +0100
    Report bugs to <help@r-inla.org>
    Set reordering to id=[0] and name=[default]
    Process file[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/Model.ini] threads[2] max.threads[24] blas_threads[2] nested[2:1]
    inla_build...
    	number of sections=[131]
    	parse section=[0] name=[INLA.libR] type=[LIBR]
    	inla_parse_libR...
    		section[INLA.libR]
    			R_HOME=[/data/gannet/ripley/R/R-devel]
    	parse section=[8] name=[INLA.Expert] type=[EXPERT]
    	inla_parse_expert...
    		section[INLA.Expert]
    			disable.gaussian.check=[0]
    			Measure dot.product.gain=[No]
    			cpo.manual=[0]
    			jp.file=[(null)]
    			jp.model=[(null)]
    	parse section=[1] name=[INLA.Model] type=[PROBLEM]
    	inla_parse_problem...
    		name=[INLA.Model]
    		R-INLA version = [23.04.24]
    		R-INLA build date = [Mon Apr 24 10:06:18 PM +03 2023 (Version_23.04.24)]
    		Build tag = [Version_23.04.24]
    		System memory = [62.8Gb]
    		Cores = (Physical= 24, Logical= 24)
    		'char' is signed
    		openmp.strategy=[default]
    		pardiso-library installed and working? = [no]
    		smtp = [taucs]
    		strategy = [default]
    		store results in directory=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/results.files]
    		output:
    			gcpo=[0]
    				num.level.sets=[-1]
    				size.max=[32]
    				strategy=[Posterior]
    				correct.hyperpar=[1]
    				epsilon=[0.005]
    				prior.diagonal=[0.0001]
    				keep=[]
    				remove.fixed=[1]
    				remove=[]
    			cpo=[0]
    			po=[0]
    			dic=[0]
    			kld=[1]
    			mlik=[1]
    			q=[0]
    			graph=[0]
    			hyperparameters=[1]
    			config=[0]
    			likelihood.info=[0]
    			internal.opt=[1]
    			summary=[1]
    			return.marginals=[1]
    			return.marginals.predictor=[0]
    			nquantiles=[3]  [ 0.025 0.5 0.975 ]
    			ncdf=[0]  [ ]
    	parse section=[3] name=[Predictor] type=[PREDICTOR]
    	inla_parse_predictor ...
    		section=[Predictor]
    		dir=[predictor]
    		PRIOR->name=[loggamma]
    		hyperid=[53001|Predictor]
    		PRIOR->from_theta=[function (x) <<NEWLINE>>exp(x)]
    		PRIOR->to_theta = [function (x) <<NEWLINE>>log(x)]
    		PRIOR->PARAMETERS=[1, 1e-05]
    		initialise log_precision[13.8155]
    		fixed=[1]
    		user.scale=[1]
    		n=[120]
    		m=[0]
    		ndata=[120]
    		compute=[1]
    		read offsets from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c]
    		read n=[240] entries from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c]
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 0/120  (idx,y) = (0, -0.0877389)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 1/120  (idx,y) = (1, -0.0270927)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 2/120  (idx,y) = (2, -0.0783941)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 3/120  (idx,y) = (3, -0.801757)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 4/120  (idx,y) = (4, 0.702145)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 5/120  (idx,y) = (5, -0.555352)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 6/120  (idx,y) = (6, 0.363164)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 7/120  (idx,y) = (7, 1.49259)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 8/120  (idx,y) = (8, -0.0179155)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804a52645c] 9/120  (idx,y) = (9, 0.945142)
    		A=[(null)]
    		Aext=[(null)]
    		AextPrecision=[1e+08]
    		output:
    			summary=[1]
    			return.marginals=[1]
    			return.marginals.predictor=[0]
    			nquantiles=[3]  [ 0.025 0.5 0.975 ]
    			ncdf=[0]  [ ]
    	parse section=[2] name=[INLA.Data1] type=[DATA]
    	inla_parse_data [section 1]...
    		tag=[INLA.Data1]
    		family=[POISSON]
    		likelihood=[POISSON]
    		file->name=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828025575eb9]
    		file->name=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582804266e41d]
    		file->name=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582801dbd3bed]
    		file->name=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828056082521]
    		read n=[360] entries from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828025575eb9]
    		mdata.nattributes = 0
    			0/120  (idx,a,y,d) = (0, 1, 4, 1)
    			1/120  (idx,a,y,d) = (1, 1, 0, 1)
    			2/120  (idx,a,y,d) = (2, 1, 3, 1)
    			3/120  (idx,a,y,d) = (3, 1, 0, 1)
    			4/120  (idx,a,y,d) = (4, 1, 5, 1)
    			5/120  (idx,a,y,d) = (5, 1, 1, 1)
    			6/120  (idx,a,y,d) = (6, 1, 0, 1)
    			7/120  (idx,a,y,d) = (7, 1, 3, 1)
    			8/120  (idx,a,y,d) = (8, 1, 1, 1)
    			9/120  (idx,a,y,d) = (9, 1, 4, 1)
    		likelihood.variant=[0]
    		Link model   [LOG]
    		Link order   [-1]
    		Link variant [-1]
    		Link a       [1]
    		Link ntheta  [0]
    		mix.use[0]
    	parse section=[6] name=[region.indexS] type=[FFIELD]
    	inla_parse_ffield...
    		section=[region.indexS]
    		dir=[random.effect00000001]
    		model=[bym2]
    		PRIOR0->name=[pcprec]
    		hyperid=[11001|region.indexS]
    		PRIOR0->from_theta=[function (x) <<NEWLINE>>exp(x)]
    		PRIOR0->to_theta = [function (x) <<NEWLINE>>log(x)]
    		PRIOR0->PARAMETERS0=[1 0.05]
    		PRIOR1->name=[table]
    		hyperid=[11002|region.indexS]
    		PRIOR1->from_theta=[function (x) <<NEWLINE>>exp(x)/(1 + exp(x))]
    		PRIOR1->to_theta = [function (x) <<NEWLINE>>log(x/(1 - x))]
    		PRIOR1->table=[table: /tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file58280398e6fa6]
    		vb.correct n[1]  -1
    		correct=[-1]
    		constr=[0]
    		diagonal=[0.0001]
    		id.names=<not present>
    		compute=[1]
    		nrep=[1]
    		ngroup=[1]
    		read covariates from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6]
    		read n=[240] entries from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6]
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 0/120  (idx,y) = (0, 0)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 1/120  (idx,y) = (1, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 2/120  (idx,y) = (2, 2)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 3/120  (idx,y) = (3, 3)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 4/120  (idx,y) = (4, 4)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 5/120  (idx,y) = (5, 5)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 6/120  (idx,y) = (6, 6)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 7/120  (idx,y) = (7, 7)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 8/120  (idx,y) = (8, 8)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582802bf750e6] 9/120  (idx,y) = (9, 9)
    		read graph from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file582801d4aeb26]
    		file for locations=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file58280534fc4fa]
    			nlocations=[120]
    			locations[0]=[1]
    			locations[1]=[2]
    			locations[2]=[3]
    			locations[3]=[4]
    			locations[4]=[5]
    			locations[5]=[6]
    			locations[6]=[7]
    			locations[7]=[8]
    			locations[8]=[9]
    			locations[9]=[10]
    		initialise log_precision [4]
    		fixed=[0]
    		initialise phi_intern [-3]
    		fixed=[0]
    		adjust.for.con.comp[1]
    		scale.model[1]
    		connected component[0] size[120] scale[0.507989]
    		scale.model: prec_scale[0.507989]
    		read extra constraint from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828063538b12]
    		Constraint[0]
    			A[120] = 1.000000
    			A[121] = 1.000000
    			A[122] = 1.000000
    			A[123] = 1.000000
    			A[124] = 1.000000
    			A[125] = 1.000000
    			A[126] = 1.000000
    			A[127] = 1.000000
    			A[128] = 1.000000
    			A[129] = 1.000000
    			A[130] = 1.000000
    			e[0] = 0.000000
    		rank-deficiency is *defined* [1]
    		output:
    			summary=[1]
    			return.marginals=[1]
    			return.marginals.predictor=[0]
    			nquantiles=[3]  [ 0.025 0.5 0.975 ]
    			ncdf=[0]  [ ]
    	section=[4] name=[(Intercept)] type=[LINEAR]
    	inla_parse_linear...
    		section[(Intercept)]
    		dir=[fixed.effect00000001]
    		file for covariates=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198]
    		read n=[240] entries from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198]
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 0/120  (idx,y) = (0, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 1/120  (idx,y) = (1, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 2/120  (idx,y) = (2, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 3/120  (idx,y) = (3, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 4/120  (idx,y) = (4, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 5/120  (idx,y) = (5, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 6/120  (idx,y) = (6, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 7/120  (idx,y) = (7, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 8/120  (idx,y) = (8, 1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828024a44198] 9/120  (idx,y) = (9, 1)
    		prior mean=[0]
    		prior precision=[0]
    		compute=[1]
    		output:
    			summary=[1]
    			return.marginals=[1]
    			return.marginals.predictor=[0]
    			nquantiles=[3]  [ 0.025 0.5 0.975 ]
    			ncdf=[0]  [ ]
    	section=[5] name=[poverty] type=[LINEAR]
    	inla_parse_linear...
    		section[poverty]
    		dir=[fixed.effect00000002]
    		file for covariates=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6]
    		read n=[240] entries from file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6]
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 0/120  (idx,y) = (0, 23.1)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 1/120  (idx,y) = (1, 20.5)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 2/120  (idx,y) = (2, 10.8)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 3/120  (idx,y) = (3, 16.5)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 4/120  (idx,y) = (4, 21.7)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 5/120  (idx,y) = (5, 25.6)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 6/120  (idx,y) = (6, 41.6)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 7/120  (idx,y) = (7, 8.7)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 8/120  (idx,y) = (8, 18.2)
    		file=[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/data.files/file5828057c227c6] 9/120  (idx,y) = (9, 20.5)
    		prior mean=[0]
    		prior precision=[0.001]
    		compute=[1]
    		output:
    			summary=[1]
    			return.marginals=[1]
    			return.marginals.predictor=[0]
    			nquantiles=[3]  [ 0.025 0.5 0.975 ]
    			ncdf=[0]  [ ]
    	parse section=[130] name=[INLA.pardiso] type=[PARDISO]
    	inla_parse_pardiso...
    		section[INLA.pardiso]
    		verbose[0]
    		debug[0]
    		parallel.reordering[1]
    		nrhs[-1]
    	parse section=[129] name=[INLA.lp.scale] type=[LP.SCALE]
    	inla_parse_lp_scale...
    		section[INLA.lp.scale]
    
    	Index table: number of entries[4], total length[362]
    		tag                            start-index     length
    		Predictor                               0        120
    		region.indexS                         120        240
    		(Intercept)                           360          1
    		poverty                               361          1
    	List of hyperparameters: 
    		theta[0] = [Log precision for region.indexS]
    		theta[1] = [Logit phi for region.indexS]
    
    	parse section=[7] name=[INLA.Parameters] type=[INLA]
    	inla_parse_INLA...
    		section[INLA.Parameters]
    			lincomb.derived.correlation.matrix = [No]
    		global_node.factor = 2.000
    		global_node.degree = 2147483647
    		reordering = -1
    		constr.marginal.diagonal = 1.49e-08
    Contents of ai_param 0x7fa01c2801c0
    	Optimiser: DEFAULT METHOD
    		Option for GSL-BFGS2: tol  = 0.1
    		Option for GSL-BFGS2: step_size = 1
    		Option for GSL-BFGS2: epsx = 0.005
    		Option for GSL-BFGS2: epsf = 0.01
    		Option for GSL-BFGS2: epsg = 0.005
    		Restart: 0
    		Optimise: try to be smart: Yes
    		Optimise: use directions: Yes
    		Mode known: No
    		parallel linesearch [0]
    	Gaussian approximation:
    		tolerance_func = 0.0005
    		tolerance_step = 0.001
    		optpar_fp = 0
    		optpar_nr_step_factor = -0.1
    	Gaussian data: No
    	Strategy: 	Use a mean-skew corrected Gaussian by fitting a Skew-Normal
    	Fast mode: 	On
    	Use linear approximation to log(|Q +c|)? Yes
    		Method:	 Compute the derivative exact
    	Parameters for improved approximations
    		Number of points evaluate:	 9
    		Step length to compute derivatives numerically:	 0.0001
    		Stencil to compute derivatives numerically:	 5
    		Cutoff value to construct local neigborhood:	 0.0001
    	Log calculations:	 On
    	Log calculated marginal for the hyperparameters:	 On
    	Integration strategy:	 Automatic (GRID for dim(theta)=1 and 2 and otherwise CCD)
    		f0 (CCD only):	 1.100
    		dz (GRID only):	 0.750
    		Adjust weights (GRID only):	 On
    		Difference in log-density limit (GRID only):	 6.000
    		Skip configurations with (presumed) small density (GRID only):	 On
    	Gradient is computed using Central difference with step-length 0.005000
    	Hessian is computed using Central difference with step-length 0.070711
    	Hessian matrix is forced to be a diagonal matrix? [No]
    	Compute effective number of parameters? [Yes]
    	Perform a Monte Carlo error-test? [No]
    	Interpolator [Auto]
    	CPO required diff in log-density [3]
    	Stupid search mode:
    		Status     [On]
    		Max iter   [1000]
    		Factor     [1.05]
    	Numerical integration of hyperparameters:
    		Maximum number of function evaluations [100000]
    		Relative error ....................... [1e-05]
    		Absolute error ....................... [1e-06]
    	To stabilise the numerical optimisation:
    		Minimum value of the -Hessian [-inf]
    		Strategy for the linear term [Skip]
    	CPO manual calculation[No]
    	VB correction is [Enabled]
    		strategy            = [mean]
    		verbose             = [Yes]
    		f_enable_limit_mean = [30]
    		f_enable_limit_var  = [25]
    		iter_max            = [25]
    		emergency           = [25.00]
    		hessian_update      = [2]
    		hessian_strategy    = [full]
    	Read [120] sections with mode=[LINCOMB]
    inla_build: check for unused entries in[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/Model.ini]
    inla_INLA_preopt_experimental...
    	Strategy = [DEFAULT]
    	Mode....................... [Compact]
    	Setup...................... [0.11s]
    	Sparse-matrix library...... [taucs]
    	OpenMP strategy............ [small]
    	num.threads................ [2:1]
    	blas.num.threads........... [2]
    	Density-strategy........... [High]
    	Size of graph.............. [242]
    	Number of constraints...... [1]
    	Optimizing sort2 cut-off... [64]
    	Optimizing Qx-strategy..... serial[0.365] parallel [0.635] choose[serial]
    	Optimizing pred-strategy... plain [0.512] data-rich[0.488] choose[plain]
    	Optimizing dot-products.... plain....[0.333] group....[0.205]
    	                            plain.mkl[0.258] group.mkl[0.205]
    	                            ==> optimal.mix.strategy  [0.195]
    	                                plain....[ 0.0%] group....[40.0%]
    	                                plain.mkl[10.0%] group.mkl[50.0%]
    	Found optimal reordering=[amdc] nnz(L)=[1004] and use_global_nodes(user)=[no]
    
    	List of hyperparameters: 
    		theta[0] = [Log precision for region.indexS]
    		theta[1] = [Logit phi for region.indexS]
    
    
    Compute initial values...
    	Iter[0] RMS(err) = 1.000, update with step-size = 1.473
    	Iter[1] RMS(err) = 0.447, update with step-size = 0.441
    	Iter[2] RMS(err) = 0.201, update with step-size = 1.475
    	Initial values computed in 0.1026 seconds
    		x[0] = 0.2783
    		x[1] = -0.1126
    		x[2] = 0.2207
    		x[3] = 0.6318
    		x[4] = -0.3818
    		x[237] = 0.0000
    		x[238] = 0.0000
    		x[239] = 0.0000
    		x[240] = 0.0037
    		x[241] = 0.0047
    
    Optimise using DEFAULT METHOD
    Smart optimise part I: estimate gradient using forward differences
    maxld= -133.154 fn=  1 theta= 4.000 -3.000 [0.2, 0.02]
    maxld= -133.153 fn=  3 theta= 4.000 -2.995 [0.2, 0.02]
    maxld= -132.891 fn=  4 theta= 3.405 -2.196 [0.4, 0.02]
    maxld= -132.891 fn=  5 theta= 3.410 -2.196 [0.4, 0.02]
    maxld= -132.890 fn=  6 theta= 3.405 -2.196 [0.4, 0.01]
    maxld= -132.890 fn=  7 theta= 3.405 -2.191 [0.4, 0.01]
    maxld= -132.886 fn=  9 theta= 3.333 -2.098 [0.4, 0.01]
    maxld= -132.886 fn= 10 theta= 3.338 -2.098 [0.4, 0.01]
    maxld= -132.886 fn= 12 theta= 3.333 -2.093 [0.4, 0.01]
    New directions for numerical gradient
    	  dir01	  dir02
    	 -0.595	  0.804
    	  0.804	  0.595
    Iter=1 |grad|=0.0411 |x-x.old|=0.793 |f-f.old|=0.269 
    maxld= -132.884 fn= 14 theta= 3.349 -2.004 [0.5, 0.01]
    maxld= -132.884 fn= 16 theta= 3.349 -2.004 [0.5, 0.01]
    maxld= -132.883 fn= 17 theta= 3.353 -2.001 [0.5, 0.01]
    New directions for numerical gradient
    	  dir01	  dir02
    	   . 	  0.986
    	  0.986	   . 
    Iter=2 |grad| = 0.0227(pass) |x-x.old|=0.0679 |f-f.old|=0.00226(pass) 
    Smart optimise part II: estimate gradient using central differences
    Smart optimise part II: restart optimiser
    maxld= -132.883 fn= 18 theta= 3.354 -1.996 [0.5, 0.01]
    maxld= -132.883 fn= 21 theta= 3.358 -2.002 [0.5, 0.01]
    maxld= -132.883 fn= 23 theta= 3.409 -2.008 [0.4, 0.01]
    maxld= -132.883 fn= 25 theta= 3.404 -2.007 [0.4, 0.01]
    maxld= -132.883 fn= 31 theta= 3.398 -2.001 [0.5, 0.01]
    New directions for numerical gradient
    	  dir01	  dir02
    	  0.992	   . 
    	   . 	  0.992
    Iter=1 |grad| = 0.00493(pass) |x-x.old|=0.0315 |f-f.old|=0.000721(pass) 
    Optim: Number of function evaluations = 32
    Compute the Hessian using central differences and step_size[0.0707107]. Matrix-type [dense]
    Smart optimise part IV: estimate Hessian using central differences
    
    GMRFLib_opt_estimate_hessian: ensure spd for Hessian: mode seems fine
    GMRFLib_opt_estimate_hessian: set tol=[0.6732]. number of negative eigenvalues=[0]
    
          0.539     -0.060
         -0.060      0.400
    Eigenvectors of the Hessian
    	  0.938	  0.348
    	 -0.348	  0.938
    Eigenvalues of the Hessian
    	       0.561
    	       0.378
    StDev/Correlation matrix (scaled inverse Hessian)
       1.374   0.129
               1.595
    Compute corrected stdev for theta[0]: negative 0.542  positive 1.479
    Compute corrected stdev for theta[1]: negative 0.991  positive 1.071
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.002sec] hess/grad[4.004]
    		Number of nodes corrected for [32] max(dx/sd)[0.0253]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.002sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0003]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.003]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[240] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.025]
    		Node[241] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.024]
    		Implied correction for [330] nodes
    config  0/45=[  0.000  0.833 ] log(rel.dens)= -0.350, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.048]
    		Number of nodes corrected for [32] max(dx/sd)[0.0070]
    		Node[4] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[12] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[36] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[-0.001] dx/sd[-0.006] |x-mode|/sd[-0.006]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[84] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[108] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[116] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[124] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[132] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[156] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[172] delta[-0.010] dx/sd[-0.006] |x-mode|/sd[-0.006]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[188] delta[0.003] dx/sd[0.004] |x-mode|/sd[0.004]
    		Node[196] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[204] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[212] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[228] delta[0.003] dx/sd[0.004] |x-mode|/sd[0.004]
    		Node[236] delta[-0.003] dx/sd[-0.003] |x-mode|/sd[-0.003]
    		Node[240] delta[-0.002] dx/sd[-0.007] |x-mode|/sd[-0.007]
    		Node[241] delta[0.000] dx/sd[0.007] |x-mode|/sd[0.007]
    		Implied correction for [330] nodes
    config  1/45=[  0.000  2.083 ] log(rel.dens)= -1.865, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.016sec] hess/grad[7.604]
    		Number of nodes corrected for [32] max(dx/sd)[0.2843]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.016sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0368]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.016sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0080]
    		Node[4] delta[0.001] dx/sd[0.003] |x-mode|/sd[0.226]
    		Node[12] delta[0.004] dx/sd[0.008] |x-mode|/sd[0.190]
    		Node[20] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[-0.068]
    		Node[28] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.059]
    		Node[36] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.018]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.090]
    		Node[52] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[-0.059]
    		Node[60] delta[0.003] dx/sd[0.006] |x-mode|/sd[0.167]
    		Node[68] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.115]
    		Node[76] delta[0.002] dx/sd[0.003] |x-mode|/sd[0.050]
    		Node[84] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.002]
    		Node[92] delta[0.001] dx/sd[0.002] |x-mode|/sd[0.051]
    		Node[100] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.044]
    		Node[108] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.031]
    		Node[116] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.021]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.008]
    		Node[132] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.040]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.013]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.025]
    		Node[172] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.058]
    		Node[180] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.038]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.033]
    		Node[196] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.039]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.010]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.008]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[236] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.045]
    		Node[240] delta[-0.002] dx/sd[-0.005] |x-mode|/sd[-0.253]
    		Node[241] delta[0.000] dx/sd[0.004] |x-mode|/sd[0.189]
    		Implied correction for [330] nodes
    config  2/45=[ -1.896 -1.928 ] log(rel.dens)= -20.618, [0] accept, compute, 0.02s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.013sec] hess/grad[4.598]
    		Number of nodes corrected for [32] max(dx/sd)[0.0051]
    		Node[4] delta[0.000] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[28] delta[0.000] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[0.000] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[92] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[100] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[108] delta[0.000] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[0.002] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.002] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[156] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[-0.003] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[180] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[188] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.002] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[212] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[220] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[228] delta[0.002] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[0.001] dx/sd[0.005] |x-mode|/sd[0.005]
    		Node[241] delta[-0.000] dx/sd[-0.005] |x-mode|/sd[-0.005]
    		Implied correction for [330] nodes
    config  3/45=[  0.000  3.750 ] log(rel.dens)= -4.420, [1] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[7.404]
    		Number of nodes corrected for [32] max(dx/sd)[0.6773]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0363]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0011]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.093]
    		Node[12] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.063]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.025]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.018]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.030]
    		Node[52] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.059]
    		Node[60] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.070]
    		Node[68] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.040]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.014]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.012]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.016]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.034]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.037]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.006]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.012]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.019]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.066]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.039]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.032]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.034]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.014]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.007]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.048]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.235]
    		Node[241] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.167]
    		Implied correction for [330] nodes
    config  4/45=[ -1.896 -0.771 ] log(rel.dens)= -10.102, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[6.050]
    		Number of nodes corrected for [32] max(dx/sd)[0.0721]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0020]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.010]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[240] delta[0.000] dx/sd[0.002] |x-mode|/sd[-0.070]
    		Node[241] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[0.067]
    		Implied correction for [330] nodes
    config  5/45=[  1.150 -3.470 ] log(rel.dens)= -4.098, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.740]
    		Number of nodes corrected for [32] max(dx/sd)[0.5044]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0219]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.002sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0005]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.040]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.025]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.007]
    		Node[52] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.057]
    		Node[60] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.035]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.022]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.009]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.005]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[116] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.029]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.028]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.005]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.011]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.010]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.068]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.031]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.029]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.025]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.014]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.007]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.040]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.188]
    		Node[241] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.121]
    		Implied correction for [330] nodes
    config  6/45=[ -1.896  0.000 ] log(rel.dens)= -5.840, [0] accept, compute, 0.02s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.014sec] hess/grad[176.933]
    		Number of nodes corrected for [32] max(dx/sd)[0.0281]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.014sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0005]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[60] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[0.000] dx/sd[0.001] |x-mode|/sd[-0.028]
    		Node[241] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.026]
    		Implied correction for [330] nodes
    config  7/45=[  1.150 -1.928 ] log(rel.dens)= -1.843, [1] accept, compute, 0.02s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.013sec] hess/grad[185.333]
    		Number of nodes corrected for [32] max(dx/sd)[0.3284]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.013sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0096]
    		Node[4] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.013]
    		Node[12] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.012]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[36] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.010]
    		Node[44] delta[0.000] dx/sd[0.001] |x-mode|/sd[-0.008]
    		Node[52] delta[0.003] dx/sd[0.005] |x-mode|/sd[-0.060]
    		Node[60] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[0.014]
    		Node[68] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.018]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[116] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.022]
    		Node[124] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.006]
    		Node[132] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[0.017]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[148] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.006]
    		Node[156] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.015]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[172] delta[0.008] dx/sd[0.006] |x-mode|/sd[-0.068]
    		Node[180] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[0.016]
    		Node[188] delta[-0.001] dx/sd[-0.002] |x-mode|/sd[0.025]
    		Node[196] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.012]
    		Node[204] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.005]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[0.011]
    		Node[236] delta[0.001] dx/sd[0.002] |x-mode|/sd[-0.030]
    		Node[240] delta[0.001] dx/sd[0.002] |x-mode|/sd[-0.140]
    		Node[241] delta[0.000] dx/sd[0.005] |x-mode|/sd[-0.072]
    		Implied correction for [330] nodes
    config  8/45=[ -1.896  0.833 ] log(rel.dens)= -3.803, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.189]
    		Number of nodes corrected for [32] max(dx/sd)[0.0094]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.004] |x-mode|/sd[-0.004]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[132] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.001] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[240] delta[-0.002] dx/sd[-0.009] |x-mode|/sd[-0.009]
    		Node[241] delta[0.000] dx/sd[0.009] |x-mode|/sd[0.009]
    		Implied correction for [330] nodes
    config  9/45=[  1.150 -0.771 ] log(rel.dens)= -0.713, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.376]
    		Number of nodes corrected for [32] max(dx/sd)[0.0855]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0027]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.019]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[52] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.056]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.020]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.014]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.020]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[156] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.020]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[172] delta[0.002] dx/sd[0.001] |x-mode|/sd[-0.058]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.022]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.016]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.016]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.022]
    		Node[240] delta[0.001] dx/sd[0.003] |x-mode|/sd[-0.083]
    		Node[241] delta[-0.000] dx/sd[-0.003] |x-mode|/sd[0.077]
    		Implied correction for [330] nodes
    config 10/45=[ -1.896  2.083 ] log(rel.dens)= -3.807, [0] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.157]
    		Number of nodes corrected for [32] max(dx/sd)[0.0031]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.003] |x-mode|/sd[-0.003]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[132] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[156] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[180] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.001] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[212] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[220] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[240] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[241] delta[0.000] dx/sd[0.002] |x-mode|/sd[0.002]
    		Implied correction for [330] nodes
    config 11/45=[  1.150 -0.000 ] log(rel.dens)= -0.385, [1] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.770]
    		Number of nodes corrected for [32] max(dx/sd)[0.0039]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[92] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[124] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[132] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[156] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[180] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[204] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[212] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[220] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[228] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[240] delta[0.001] dx/sd[0.004] |x-mode|/sd[0.004]
    		Node[241] delta[-0.000] dx/sd[-0.004] |x-mode|/sd[-0.004]
    		Implied correction for [330] nodes
    config 12/45=[  1.150  0.833 ] log(rel.dens)= -0.525, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.088]
    		Number of nodes corrected for [32] max(dx/sd)[0.0090]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[28] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[36] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[92] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[100] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[108] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[116] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[124] delta[0.001] dx/sd[0.002] |x-mode|/sd[0.002]
    		Node[132] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[156] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.001] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[204] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[212] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[220] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[228] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[240] delta[0.002] dx/sd[0.009] |x-mode|/sd[0.009]
    		Node[241] delta[-0.000] dx/sd[-0.009] |x-mode|/sd[-0.009]
    		Implied correction for [330] nodes
    config 13/45=[  1.150  2.083 ] log(rel.dens)= -1.819, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.830]
    		Number of nodes corrected for [32] max(dx/sd)[0.2959]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0354]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0050]
    		Node[4] delta[0.001] dx/sd[0.002] |x-mode|/sd[0.187]
    		Node[12] delta[0.003] dx/sd[0.005] |x-mode|/sd[0.142]
    		Node[20] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.055]
    		Node[28] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.048]
    		Node[36] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.017]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.069]
    		Node[52] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.048]
    		Node[60] delta[0.002] dx/sd[0.004] |x-mode|/sd[0.131]
    		Node[68] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.085]
    		Node[76] delta[0.001] dx/sd[0.002] |x-mode|/sd[0.026]
    		Node[84] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.009]
    		Node[92] delta[0.001] dx/sd[0.002] |x-mode|/sd[0.041]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.034]
    		Node[108] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.029]
    		Node[116] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.022]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.003]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.010]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.006]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.014]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.010]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.008]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.009]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[240] delta[-0.002] dx/sd[-0.005] |x-mode|/sd[-0.265]
    		Node[241] delta[0.000] dx/sd[0.004] |x-mode|/sd[0.214]
    		Implied correction for [330] nodes
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.862]
    		Number of nodes corrected for [32] max(dx/sd)[0.0125]
    config 14/45=[ -1.053 -3.470 ] log(rel.dens)= -19.589, [0] accept, compute, 0.00s
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[60] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Implied correction for [330] nodes
    config 15/45=[  1.150  3.750 ] log(rel.dens)= -4.488, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.012]
    		Number of nodes corrected for [32] max(dx/sd)[0.0029]
    		Node[4] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.003] |x-mode|/sd[-0.003]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[241] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Implied correction for [330] nodes
    config 16/45=[  2.875 -3.470 ] log(rel.dens)= -5.257, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.029]
    		Number of nodes corrected for [32] max(dx/sd)[0.0092]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[0.002] dx/sd[0.009] |x-mode|/sd[0.009]
    		Node[241] delta[-0.000] dx/sd[-0.009] |x-mode|/sd[-0.009]
    		Implied correction for [330] nodes
    config 17/45=[  2.875 -1.928 ] log(rel.dens)= -3.309, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.917]
    		Number of nodes corrected for [32] max(dx/sd)[0.0120]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.011]
    		Implied correction for [330] nodes
    config 18/45=[  2.875 -0.771 ] log(rel.dens)= -2.080, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.005sec] hess/grad[5.915]
    		Number of nodes corrected for [32] max(dx/sd)[0.2409]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.005sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0177]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.006sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0013]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.055]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.024]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.015]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.007]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.013]
    		Node[52] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.022]
    		Node[60] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.036]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.007]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.013]
    		Node[116] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.020]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.004]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.012]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.006]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.007]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.019]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.013]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.010]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.011]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.006]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.005]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.017]
    		Node[240] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.225]
    		Node[241] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.207]
    		Implied correction for [330] nodes
    config 19/45=[ -1.053 -1.928 ] log(rel.dens)= -7.036, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.683]
    		Number of nodes corrected for [32] max(dx/sd)[0.0131]
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[6.337]
    		Number of nodes corrected for [32] max(dx/sd)[0.1724]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.013]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Implied correction for [330] nodes
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0098]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.015]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.012]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.005]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.002]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.011]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[172] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.022]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.019]
    		Node[240] delta[0.003] dx/sd[0.010] |x-mode|/sd[-0.163]
    		Node[241] delta[-0.000] dx/sd[-0.009] |x-mode|/sd[0.153]
    		Implied correction for [330] nodes
    config 20/45=[  2.875  0.000 ] log(rel.dens)= -1.560, [1] accept, compute, 0.00s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.213]
    		Number of nodes corrected for [32] max(dx/sd)[0.0138]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.014]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Implied correction for [330] nodes
    config 21/45=[  2.875  0.833 ] log(rel.dens)= -1.406, [1] accept, compute, 0.00s
    config 22/45=[ -1.053 -0.771 ] log(rel.dens)= -2.280, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.868]
    		Number of nodes corrected for [32] max(dx/sd)[0.0144]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.014]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 23/45=[  2.875  2.083 ] log(rel.dens)= -2.202, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.014sec] hess/grad[6.266]
    		Number of nodes corrected for [32] max(dx/sd)[0.1277]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.014sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0059]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.005]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.018]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.008]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.011]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[172] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.025]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[236] delta[0.001] dx/sd[0.001] |x-mode|/sd[-0.019]
    		Node[240] delta[0.002] dx/sd[0.006] |x-mode|/sd[-0.122]
    		Node[241] delta[-0.000] dx/sd[-0.006] |x-mode|/sd[0.115]
    		Implied correction for [330] nodes
    config 24/45=[ -1.053  0.000 ] log(rel.dens)= -1.055, [0] accept, compute, 0.02s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.480]
    		Number of nodes corrected for [32] max(dx/sd)[0.0148]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 25/45=[  2.875  3.750 ] log(rel.dens)= -4.755, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.399]
    		Number of nodes corrected for [32] max(dx/sd)[0.0852]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0028]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.025]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[172] delta[0.001] dx/sd[0.000] |x-mode|/sd[-0.029]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.013]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.017]
    		Node[240] delta[0.001] dx/sd[0.003] |x-mode|/sd[-0.082]
    		Node[241] delta[-0.000] dx/sd[-0.003] |x-mode|/sd[0.077]
    		Implied correction for [330] nodes
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.041]
    		Number of nodes corrected for [32] max(dx/sd)[0.0148]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 26/45=[  5.175 -1.928 ] log(rel.dens)= -5.656, [1] accept, compute, 0.02s
    config 27/45=[ -1.053  0.833 ] log(rel.dens)= -1.079, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.005sec] hess/grad[53.928]
    		Number of nodes corrected for [32] max(dx/sd)[0.0150]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.005sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0001]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 28/45=[  5.175 -0.771 ] log(rel.dens)= -4.308, [1] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.922]
    		Number of nodes corrected for [32] max(dx/sd)[0.0150]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0002]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.010sec] hess/grad[5.547]
    		Number of nodes corrected for [32] max(dx/sd)[0.0414]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.010sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0008]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.003]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.025]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.010]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.003]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.026]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.013]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.008]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.009]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.011]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Node[240] delta[0.000] dx/sd[0.001] |x-mode|/sd[-0.041]
    		Node[241] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.038]
    		Implied correction for [330] nodes
    config 29/45=[ -1.053  2.083 ] log(rel.dens)= -2.449, [0] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.719]
    		Number of nodes corrected for [32] max(dx/sd)[0.0136]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0002]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[240] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[241] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Implied correction for [330] nodes
    config 30/45=[ -1.053  3.750 ] log(rel.dens)= -4.590, [0] accept, compute, 0.00s
    config 31/45=[  5.175 -0.000 ] log(rel.dens)= -3.567, [1] accept, compute, 0.02s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.967]
    		Number of nodes corrected for [32] max(dx/sd)[0.0151]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0002]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 32/45=[  5.175  0.833 ] log(rel.dens)= -3.071, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.009sec] hess/grad[6.018]
    		Number of nodes corrected for [32] max(dx/sd)[0.2484]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.010sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0131]
    	[0]Iter [2/25] VB correct with strategy [MEAN] in total[0.010sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0007]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.070]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.039]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.018]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.015]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.008]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.014]
    		Node[52] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.021]
    		Node[60] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.039]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.018]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.013]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.009]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.016]
    		Node[116] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.019]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[172] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.006]
    		Node[180] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[188] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[220] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[236] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.005]
    		Node[240] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.236]
    		Node[241] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.216]
    		Implied correction for [330] nodes
    config 33/45=[ -0.421 -3.470 ] log(rel.dens)= -10.787, [0] accept, compute, 0.01s
    	[1]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.207]
    		Number of nodes corrected for [32] max(dx/sd)[0.0151]
    	[1]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0002]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[124] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[156] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[212] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[240] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.015]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.014]
    		Implied correction for [330] nodes
    config 34/45=[  5.175  2.083 ] log(rel.dens)= -3.175, [1] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[4.419]
    		Number of nodes corrected for [32] max(dx/sd)[0.1653]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0081]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.012]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.004]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.014]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[240] delta[0.002] dx/sd[0.008] |x-mode|/sd[-0.157]
    		Node[241] delta[-0.000] dx/sd[-0.008] |x-mode|/sd[0.149]
    		Implied correction for [330] nodes
    config 35/45=[ -0.421 -1.928 ] log(rel.dens)= -2.880, [0] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.355]
    		Number of nodes corrected for [32] max(dx/sd)[0.1047]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.002sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0039]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.004]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.011]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[240] delta[0.001] dx/sd[0.004] |x-mode|/sd[-0.101]
    		Node[241] delta[-0.000] dx/sd[-0.004] |x-mode|/sd[0.096]
    		Implied correction for [330] nodes
    config 36/45=[ -0.421 -0.771 ] log(rel.dens)= -0.563, [0] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.010sec] hess/grad[106.531]
    		Number of nodes corrected for [32] max(dx/sd)[0.0726]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.011sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0021]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.007]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.009]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[240] delta[0.001] dx/sd[0.002] |x-mode|/sd[-0.070]
    		Node[241] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[0.067]
    		Implied correction for [330] nodes
    config 37/45=[ -0.421 -0.000 ] log(rel.dens)= -0.115, [0] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.247]
    		Number of nodes corrected for [32] max(dx/sd)[0.0452]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0010]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.010]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.005]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.003]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[240] delta[0.000] dx/sd[0.001] |x-mode|/sd[-0.044]
    		Node[241] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.042]
    		Implied correction for [330] nodes
    config 38/45=[ -0.421  0.833 ] log(rel.dens)= -0.457, [0] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.008sec] hess/grad[4.727]
    		Number of nodes corrected for [32] max(dx/sd)[0.0181]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.008sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0003]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[28] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[44] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[68] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[84] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[108] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[148] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.004]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.012]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.007]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.005]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.003]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[240] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.018]
    		Node[241] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.017]
    		Implied correction for [330] nodes
    config 39/45=[ -0.421  2.083 ] log(rel.dens)= -2.009, [0] accept, compute, 0.01s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.176]
    		Number of nodes corrected for [32] max(dx/sd)[0.0044]
    		Node[4] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[12] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[20] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[28] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[36] delta[0.000] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[52] delta[-0.001] dx/sd[-0.004] |x-mode|/sd[-0.004]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[68] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[76] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[84] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[116] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[124] delta[0.002] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[132] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[140] delta[-0.001] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[148] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[156] delta[0.001] dx/sd[0.001] |x-mode|/sd[0.001]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[172] delta[-0.008] dx/sd[-0.004] |x-mode|/sd[-0.004]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[0.002] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[204] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.003] dx/sd[0.003] |x-mode|/sd[0.003]
    		Node[236] delta[-0.002] dx/sd[-0.002] |x-mode|/sd[-0.002]
    		Node[240] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.001]
    		Node[241] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Implied correction for [330] nodes
    config 40/45=[ -0.421  3.750 ] log(rel.dens)= -4.448, [0] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.017sec] hess/grad[192.927]
    		Number of nodes corrected for [32] max(dx/sd)[0.1982]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.018sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0095]
    		Node[4] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.026]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.006]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.007]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.008]
    		Node[60] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.015]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[76] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[-0.005]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[92] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.012]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.008]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.010]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[240] delta[0.003] dx/sd[0.010] |x-mode|/sd[-0.189]
    		Node[241] delta[-0.000] dx/sd[-0.009] |x-mode|/sd[0.177]
    		Implied correction for [330] nodes
    config 41/45=[  0.000 -3.470 ] log(rel.dens)= -6.834, [0] accept, compute, 0.02s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.254]
    		Number of nodes corrected for [32] max(dx/sd)[0.1160]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0043]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.004]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[20] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[36] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.003]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.013]
    		Node[100] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[240] delta[0.001] dx/sd[0.004] |x-mode|/sd[-0.112]
    		Node[241] delta[-0.000] dx/sd[-0.004] |x-mode|/sd[0.106]
    		Implied correction for [330] nodes
    config 42/45=[  0.000 -1.928 ] log(rel.dens)= -1.875, [0] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.367]
    		Number of nodes corrected for [32] max(dx/sd)[0.0671]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0019]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.002]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.009]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[164] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[196] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[240] delta[0.000] dx/sd[0.002] |x-mode|/sd[-0.065]
    		Node[241] delta[-0.000] dx/sd[-0.002] |x-mode|/sd[0.062]
    		Implied correction for [330] nodes
    config 43/45=[  0.000 -0.771 ] log(rel.dens)= -0.299, [0] accept, compute, 0.00s
    	[0]Iter [0/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[5.423]
    		Number of nodes corrected for [32] max(dx/sd)[0.0437]
    	[0]Iter [1/25] VB correct with strategy [MEAN] in total[0.001sec] hess/grad[nan]
    		Number of nodes corrected for [32] max(dx/sd)[0.0010]
    		Node[4] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[12] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[20] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[28] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[36] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[44] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[52] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[60] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[68] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[76] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.001]
    		Node[84] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[92] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.006]
    		Node[100] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[108] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.001]
    		Node[116] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.001]
    		Node[124] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[132] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[140] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.000]
    		Node[148] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[156] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.001]
    		Node[164] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.000]
    		Node[172] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.003]
    		Node[180] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[188] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[196] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.000]
    		Node[204] delta[0.000] dx/sd[0.000] |x-mode|/sd[0.002]
    		Node[212] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[220] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[-0.000]
    		Node[228] delta[-0.000] dx/sd[-0.000] |x-mode|/sd[0.002]
    		Node[236] delta[0.000] dx/sd[0.000] |x-mode|/sd[-0.002]
    		Node[240] delta[0.000] dx/sd[0.001] |x-mode|/sd[-0.043]
    		Node[241] delta[-0.000] dx/sd[-0.001] |x-mode|/sd[0.041]
    		Implied correction for [330] nodes
    config 44/45=[  0.000  0.000 ] log(rel.dens)= -0.001, [0] accept, compute, 0.02s
    
    Combine the densities with relative weights:
    config  0/45=[ -1.896 -1.928 ] weight=  0.000 adjusted weight=  0.000
    config  1/45=[ -1.896 -0.771 ] weight=  0.000 adjusted weight=  0.000
    config  2/45=[ -1.896  0.000 ] weight=  0.006 adjusted weight=  0.006
    config  3/45=[ -1.896  0.833 ] weight=  0.028 adjusted weight=  0.028
    config  4/45=[ -1.896  2.083 ] weight=  0.051 adjusted weight=  0.051
    config  5/45=[ -1.053 -3.470 ] weight=  0.000 adjusted weight=  0.000
    config  6/45=[ -1.053 -1.928 ] weight=  0.001 adjusted weight=  0.001
    config  7/45=[ -1.053 -0.771 ] weight=  0.078 adjusted weight=  0.078
    config  8/45=[ -1.053  0.000 ] weight=  0.411 adjusted weight=  0.411
    config  9/45=[ -1.053  0.833 ] weight=  0.259 adjusted weight=  0.259
    config 10/45=[ -1.053  2.083 ] weight=  0.120 adjusted weight=  0.120
    config 11/45=[ -1.053  3.750 ] weight=  0.023 adjusted weight=  0.023
    config 12/45=[ -0.421 -3.470 ] weight=  0.000 adjusted weight=  0.000
    config 13/45=[ -0.421 -1.928 ] weight=  0.043 adjusted weight=  0.043
    config 14/45=[ -0.421 -0.771 ] weight=  0.238 adjusted weight=  0.238
    config 15/45=[ -0.421 -0.000 ] weight=  0.577 adjusted weight=  0.577
    config 16/45=[ -0.421  0.833 ] weight=  0.265 adjusted weight=  0.265
    config 17/45=[ -0.421  2.083 ] weight=  0.102 adjusted weight=  0.102
    config 18/45=[ -0.421  3.750 ] weight=  0.015 adjusted weight=  0.015
    config 19/45=[  0.000 -3.470 ] weight=  0.002 adjusted weight=  0.002
    config 20/45=[  0.000 -1.928 ] weight=  0.181 adjusted weight=  0.181
    config 21/45=[  0.000 -0.771 ] weight=  0.480 adjusted weight=  0.480
    config 22/45=[  0.000  0.000 ] weight=  1.000 adjusted weight=  1.000
    config 23/45=[  0.000  0.833 ] weight=  0.456 adjusted weight=  0.456
    config 24/45=[  0.000  2.083 ] weight=  0.183 adjusted weight=  0.183
    config 25/45=[  0.000  3.750 ] weight=  0.023 adjusted weight=  0.023
    config 26/45=[  1.150 -3.470 ] weight=  0.021 adjusted weight=  0.021
    config 27/45=[  1.150 -1.928 ] weight=  0.121 adjusted weight=  0.121
    config 28/45=[  1.150 -0.771 ] weight=  0.205 adjusted weight=  0.205
    config 29/45=[  1.150 -0.000 ] weight=  0.440 adjusted weight=  0.440
    config 30/45=[  1.150  0.833 ] weight=  0.248 adjusted weight=  0.248
    config 31/45=[  1.150  2.083 ] weight=  0.124 adjusted weight=  0.124
    config 32/45=[  1.150  3.750 ] weight=  0.014 adjusted weight=  0.014
    config 33/45=[  2.875 -3.470 ] weight=  0.012 adjusted weight=  0.012
    config 34/45=[  2.875 -1.928 ] weight=  0.051 adjusted weight=  0.051
    config 35/45=[  2.875 -0.771 ] weight=  0.095 adjusted weight=  0.095
    config 36/45=[  2.875  0.000 ] weight=  0.248 adjusted weight=  0.248
    config 37/45=[  2.875  0.833 ] weight=  0.187 adjusted weight=  0.187
    config 38/45=[  2.875  2.083 ] weight=  0.154 adjusted weight=  0.154
    config 39/45=[  2.875  3.750 ] weight=  0.020 adjusted weight=  0.020
    config 40/45=[  5.175 -1.928 ] weight=  0.008 adjusted weight=  0.008
    config 41/45=[  5.175 -0.771 ] weight=  0.017 adjusted weight=  0.017
    config 42/45=[  5.175 -0.000 ] weight=  0.055 adjusted weight=  0.055
    config 43/45=[  5.175  0.833 ] weight=  0.058 adjusted weight=  0.058
    config 44/45=[  5.175  2.083 ] weight=  0.095 adjusted weight=  0.095
    
    Marginal likelihood: Integration -130.490 Gaussian-approx -130.269
    Compute the marginal for each of the 2 hyperparameters
    Interpolation method: Auto
    	Compute the marginal for theta[0] to theta[1] using numerical integration...
    	Compute the marginal for theta[0] to theta[1] using numerical integration... Done.
    Compute the marginal for the hyperparameters... done.
    Store results in directory[/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/results.files]
    
    Wall-clock time used on [/tmp/RtmpipCVyJ/working_dir/Rtmp7czjB8/file58280503d42d7/Model.ini]
    	Preparations             :   0.214 seconds
    	Approx inference (stage1):   1.327 seconds
    	Approx inference (stage2):   0.000 seconds
    	Approx inference (total) :   1.327 seconds
    	Output                   :   0.029 seconds
    	------------------------------------------
    	Total                    :   1.570 seconds
    
    Stage1:	Number of fn-calls= 90 with 0.0147 sec/fn-call
    Stage2:	Number of fn-calls= 0 with 1.3273 sec/fn-call
    Stage1:	Accumulated CPU-time is equivalent to 0.47 threads running at 100%
    	Efficiency using 2 threads = 23.42%
    Stage2:	Accumulated CPU-time is equivalent to 0.43 threads running at 100%
    	Efficiency using 2 threads = 21.66%
    Total:	Accumulated CPU-time is equivalent to 0.47 threads running at 100%
    	Efficiency using 2 threads = 23.42%
    
    --- finished re-building ‘bym.Rnw’
    
    make[1]: Entering directory '/data/gannet/ripley/R/packages/tests-devel/diseasemapping.Rcheck/vign_test/diseasemapping/vignettes'
    xelatex bym;
    This is XeTeX, Version 3.141592653-2.6-0.999995 (TeX Live 2023) (preloaded format=xelatex)
     restricted \write18 enabled.
    entering extended mode
    (./bym.tex
    LaTeX2e <2023-06-01> patch level 1
    L3 programming layer <2023-08-29>
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/base/article.cls
    Document Class: article 2023/05/17 v1.4n Standard LaTeX document class
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/base/size12.clo))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics/graphicx.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics/keyval.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics/graphics.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics/trig.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics-cfg/graphics.cfg)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics-def/xetex.def)))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/xcolor/xcolor.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics-cfg/color.cfg)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/graphics/mathcolor.ltx))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/framed/framed.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/base/alltt.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/geometry/geometry.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/generic/iftex/ifvtex.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/generic/iftex/iftex.sty)))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/amsmath/amsmath.sty
    For additional information on amsmath, use the `?' option.
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/amsmath/amstext.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/amsmath/amsgen.sty))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/amsmath/amsbsy.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/amsmath/amsopn.sty))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/caption/caption.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/caption/caption3.sty))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/caption/subcaption.sty)
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/upquote/upquote.sty
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/base/textcomp.sty))
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/l3backend/l3backend-xetex.def)
    No file bym.aux.
    (/data/gannet2/ripley/2023/texmf-dist/tex/latex/base/ts1cmr.fd)
    *geometry* driver: auto-detecting
    *geometry* detected driver: xetex
    
    Overfull \hbox (38.1448pt too wide) in paragraph at lines 99--99
    []    []\TU/lmtt/m/n/12 larynxRates[] []=[] []\TU/lmtt/b/n/12 structure[][]\TU/
    lmtt/m/n/12 ([][]\TU/lmtt/b/n/12 c[][]\TU/lmtt/m/n/12 ([][]0[][],[] []0[][],[] 
    []0[][],[] []0[][],[] []1e-06[][],[] []6e-06[][],[] []2.3e-05[][],[] []4.5e-05[
    ][],[] []9.9e-05[][],[][] 
    
    Overfull \hbox (24.59683pt too wide) in paragraph at lines 100--100
    []        []\TU/lmtt/m/n/12 0.000163[][],[] []0.000243[][],[] []0.000299[][],[]
     []0.000343[][],[] []0.000308[][],[] []0.000291[][],[] []0.000217[][],[] []0[][
    ],[][] 
    
    Overfull \hbox (61.66481pt too wide) in paragraph at lines 101--101
    []        []\TU/lmtt/m/n/12 0[][],[] []0[][],[] []1e-06[][],[] []1e-06[][],[] [
    ]3e-06[][],[] []8e-06[][],[] []1.3e-05[][],[] []2.3e-05[][],[] []3.5e-05[][],[]
     []5.8e-05[][],[] []6.8e-05[][],[][] 
    
    Overfull \hbox (55.66483pt too wide) in paragraph at lines 102--102
    []        []\TU/lmtt/m/n/12 7.5e-05[][],[] []5.5e-05[][],[] []4.1e-05[][],[] []
    3e-05[][]),[] [].Names[] []=[] []\TU/lmtt/b/n/12 c[][]\TU/lmtt/m/n/12 ([][]"M_1
    0"[][],[] []"M_15"[][],[] []"M_20"[][],[] []"M_25"[][],[][] 
    
    Overfull \hbox (67.88081pt too wide) in paragraph at lines 103--103
    []        []\TU/lmtt/m/n/12 "M_30"[][],[] []"M_35"[][],[] []"M_40"[][],[] []"M_
    45"[][],[] []"M_50"[][],[] []"M_55"[][],[] []"M_60"[][],[] []"M_65"[][],[] []"M
    _70"[][],[] []"M_75"[][],[][] 
    
    Overfull \hbox (67.9048pt too wide) in paragraph at lines 104--104
    []        []\TU/lmtt/m/n/12 "M_80"[][],[] []"M_85"[][],[] []"F_10"[][],[] []"F_
    15"[][],[] []"F_20"[][],[] []"F_25"[][],[] []"F_30"[][],[] []"F_35"[][],[] []"F
    _40"[][],[] []"F_45"[][],[][] 
    [1]
    Overfull \hbox (44.62482pt too wide) in paragraph at lines 162--162
    [][]\TU/lmtt/m/n/12 kBYM[] []=[] []\TU/lmtt/b/n/12 try[][]\TU/lmtt/m/n/12 ([][]
    \TU/lmtt/b/n/12 bym[][]\TU/lmtt/m/n/12 ([][]formula[] []= observed[] []~[] []\T
    U/lmtt/b/n/12 offset[][]\TU/lmtt/m/n/12 (logExpected)[] []+[] []poverty,[] []da
    ta[] []= kentucky,[][] 
    
    Overfull \hbox (87.16478pt too wide) in paragraph at lines 163--163
    []    []\TU/lmtt/m/n/12 prior[] []=[] []\TU/lmtt/b/n/12 list[][]\TU/lmtt/m/n/12
     ([][]sdSpatial[] []=[] []\TU/lmtt/b/n/12 c[][]\TU/lmtt/m/n/12 ([][]0.01[][],[]
     []0.2[][]),[] []sdIndep[] []=[] []\TU/lmtt/b/n/12 c[][]\TU/lmtt/m/n/12 ([][]0.
    01[][],[] []0.2[][])),[] []region.id[] []=[] []"County"[][]))[][] 
    ! Argument of \caption@@@withoptargs has an extra }.
    <inserted text> 
                    \par 
    l.210 ...dth=0.48\textwidth]{figure/priorPost-1} }
                                                      
    ? 
    ! Emergency stop.
    <inserted text> 
                    \par 
    l.210 ...dth=0.48\textwidth]{figure/priorPost-1} }
                                                      
    Output written on bym.pdf (1 page).
    Transcript written on bym.log.
    make[1]: *** [Makefile:7: bym.pdf] Error 1
    make[1]: Leaving directory '/data/gannet/ripley/R/packages/tests-devel/diseasemapping.Rcheck/vign_test/diseasemapping/vignettes'
    Error in tools::buildVignettes(dir = "/data/gannet/ripley/R/packages/tests-devel/diseasemapping.Rcheck/vign_test/diseasemapping",  : 
      running 'make' failed
    Execution halted
Flavor: r-devel-linux-x86_64-fedora-gcc