======================================================
Command used: wombat w_hinvMF.par
======================================================

wombat w_hinvMF.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_hinvMF.par"
**RUNOP --hinv -v**
**COM Ecample20/J: Calculate H^{-1} for A{Gamma} with 2 metafounders**
**PED ../InputFiles/pedigree2**
  Pedigree file : "../InputFiles/pedigree2"
**MRK ../InputFiles/MarkerCsMF.dat**
  Marker counts file : "../InputFiles/MarkerCsMF.dat"
**SPECIAL**
**HINVERSE SNP 500**
  HINV: nsnp, msnp	    500 	500
 HACK_NUMBER NSNPS
**HINVERSE METAF 2 SCALE**
  HINV: mf_gamma1    2.00000000000000
 HACK_NUMBER METAF
 no. of meta founders		2
 Expect to read UPPER triangle of GAMMA matrix from "MFGamma.dat"
  HINV: set up A-gamma inverse scaled  T glsqu	       200
**HINVERSE CENTER half**
  HINV: how_center HALF  5.000000000000000E-002
**HINVERSE out AGAMMA**
**HINVERSE EPSILON 0.d0**
  HINV: epsilon   0.000000000000000E+000
**END**

 End of parameter file reached :  11 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					   = ../InputFiles/pedigree2
 Total no. of animal IDs found			   =	       14
 no. of NRM matrices		1	    1
**************************** Count of "animal" types ****************************
 No. of non-genotyped animals			   =	       12
 No. of genotyped animals in total		   =		2
 Total no. of "animals" 			   =	       14

 Begin: read M2 from ../InputFiles/MarkerCsMF. sys-time 	    0.0 	    0.0
 tvec4	       500
 chevck marks	      500	  500
 No. of rows in ../InputFiles/MarkerCsMF.dat read  =		2
 No. of marker allele counts per row read	   =	      500
 End: read M2 from .dat 		       sys-time 	    0.0 	    0.0
 Fitting meta-founder(s) - forcing re-build of NRM Inverse
 MF	      1 	  1	      0 	  0
 MF	      2 	  2	      0 	  0
 start get_mfgamma22
 Matrix GAMMA for metafounders has been read from "MFGamma.dat"
    1  0.5500  0.4800
    2  0.4800  0.7700

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

 File "X.gin"	 with A^{-gamma} written out
 File "X.codes"  written out
 X.codes: All animals given "type" =1
 NRM Inverse no. :	     1	 "X"
 No. of elements =	    42	 log determinant =  0.000000000000000E+000
 A22: No. genotyped animals + ancestors 	   =		4
 End: calculating A22 (Colleau) 	       sys-time 	    0.0 	    0.0
 End: centering M2			       sys-time 	    0.0 	    0.0
 End: calculating GRM			       sys-time 	    0.0 	    0.0
 end calc g22 F
 F
 End: inverting g22			       sys-time 	    0.0 	    0.0
 End: inverting a22			       sys-time 	    0.0 	    0.0
 Start writing out "Hinverse.gin"
 Successfully allocated arrays for in-core NRM
 NRM Inverse no. :	     1	 "X"
 No. of elements =	    42	 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   =	   43.000000
	 0.00000		  alpha 	  0.00000
 Total time used (secs): CPU =		 0  System =	       0
 Machine used = "didgeridoo"
 "WOMBAT" has finished !  26/05/2025


0.002u 0.010s 0:00.00 0.0%	0+0k 0+96io 0pf+0w

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