======================================================
Command used: wombat w_hinvmeta1.par
======================================================

wombat w_hinvmeta1.par
Script started on 2025-05-26 17:31:12+1000
 WOMBAT 26-05-2025 Maximum no. of threads to be used set to	      8

 Parameter file opened : "w_hinvmeta1.par"
**RUNOP --hinv -v**
**COM Example20/K: calculate H^{-1} with A-gamma for single MF; estimate gamma**
**PED ../InputFiles/pedsMF1.dat**
  Pedigree file : "../InputFiles/pedsMF1.dat"
**MRK ../InputFiles/MarkerCounts.dat**
  Marker counts file : "../InputFiles/MarkerCounts.dat"
**SPECIAL**
**HINVERSE SNP 3000**
  HINV: nsnp, msnp	   3000        3000
 HACK_NUMBER NSNPS
**HINVERSE META -1 SCALE**
  HINV: mf_gamma1   -1.00000000000000
 HACK_NUMBER METAF
 ONE meta founder  -1.00000000000000
  HINV: set up A-gamma inverse scaled  T glsqu	       200
**HINVERSE CENTER half**
  HINV: how_center HALF  5.000000000000000E-002
**END**

 End of parameter file reached :   9 lines read
 Colleau's method specified for calculation of A22 - but not needed
 ***Time: begin idvec			      0 	  0 Used since last	    0.0 	0.0
 ***Time: end idvec			      0 	  0 Used since last	    0.0 	0.0
 Pedigree file					   = tmpped.dat
 Total no. of animal IDs found			   =	      750
 no. of NRM matrices		1	    1
**************************** Count of "animal" types ****************************
 No. of non-genotyped animals			   =	      380
 No. of genotyped animals in total		   =	      370
 Total no. of "animals" 			   =	      750

 Begin: read M2 from ../InputFiles/MarkerCount sys-time 	    0.0 	    0.0
 tvec4	      3000
 chevck marks	     3000	 3000
 No. of rows in ../InputFiles/MarkerCounts.dat read=	      370
 No. of marker allele counts per row read	   =	     3000
 End: read M2 from .dat 		       sys-time 	    0.3 	    0.3

 Inbreeding coefficients for random effect	   =		1 computed
 ... MEUWISSEN & LIU algorithm
 Covariance structure				   =	NRM
 No. of inbred animals				   =	       73
 Average inbreeding coefficient 		   = 0.270152	  (in %)
 ... amongst inbred animals			   =  2.77183	  (in %)

 NRM Inverse no. :	     1	 "X"
 No. of elements =	  2545	 log determinant =  -417.061508887861
 A22: No. genotyped animals + ancestors 	   =	      624
 End: calculating A22 (Colleau) 	       sys-time 	    0.3 	    0.3
 ************* Estimating degree(s) of self-relationship for metafounders ************
 No. of metafounders				   =		1
 Method to deal with missing genotype info.	   = SETUNITY
 Diagonals of Gamma for "missing" MF set to	   = 0.500000E-01
 Minimum allele frequency allowed		   =  0.00000
 No. of markers used				   =	     3000	  3000
 Threshold to disregard MF with Diag( Q'A22^{-1}Q) = 0.100000E-07
 Method to estimate "Gamma"			   = glsqu
 Estimate					   = 0.256577

 ***Time: begin idvec			      0 	  0 Used since last	    0.5 	0.3
 ***Time: end idvec			      0 	  0 Used since last	    0.0 	0.0
 Pedigree file					   = ../InputFiles/pedsMF1.dat
 Total no. of animal IDs found			   =	      751
 no. of NRM matrices		1	    1
**************************** Count of "animal" types ****************************
 No. of non-genotyped animals			   =	      381
 No. of genotyped animals in total		   =	      370
 Total no. of "animals" 			   =	      751

 Fitting meta-founder(s) - forcing re-build of NRM Inverse
 MF	      1 	  1	      0 	  0

 Inbreeding coefficients for random effect	   =		1 computed
 ... with Meta-founder(s) - modified MEUWISSEN & LIU algorithm
 Covariance structure				   =	NRM
 No. of inbred animals				   =	      750
 Average inbreeding coefficient 		   =  13.0640	  (in %)
 ... amongst inbred animals			   =  13.0640	  (in %)

 NRM Inverse no. :	     1	 "X"
 No. of elements =	  2696	 log determinant =  0.000000000000000E+000
 A22: No. genotyped animals + ancestors 	   =	      625
 End: calculating A22 (Colleau) 	       sys-time 	    0.3 	    0.3
 End: centering M2			       sys-time 	    0.3 	    0.3
 End: calculating GRM			       sys-time 	    0.3 	    0.3
 end calc g22 F
 F
 End: inverting g22			       sys-time 	    0.3 	    0.3
 End: inverting a22			       sys-time 	    0.3 	    0.3
 Start writing out "Hinverse.gin"
 Successfully allocated arrays for in-core NRM
 NRM Inverse no. :	     1	 "X"
 No. of elements =	  2696	 log determinant =  0.000000000000000E+000
 Succesfully read from file for in-core storage
 H^{-1} : No. of non-zero elements in A^11, A^21   =	   70545.000
	 0.00000		  alpha 	  0.00000
 Total time used (secs): CPU =		 1  System =	       0
 Machine used = "didgeridoo"
 "WOMBAT" has finished !  26/05/2025


1.840u 0.125s 0:00.51 384.3%	0+0k 0+7488io 0pf+0w

Script done on 2025-05-26 17:31:12+1000
