======================================================
Command used: wombat hinv_gin.par
======================================================

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

 Parameter file opened : "hinv_gin.par"
**RUNOP  --hinv -v**
**COM Example12/A1: Step 1 - generate H^{-1}**
**PED ../pedigree.dat**
  Pedigree file : "../pedigree.dat"
**MRK ../MarkerCounts.dat**
  Marker counts file : "../MarkerCounts.dat"
**SPECIAL**
**HINVERSE SNP 1000**
  HINV: nsnp, msnp	   1000        1000
 HACK_NUMBER NSNPS
**HINVERSE lambda 1.0**
  HINV: lambda	  1.00000000000000
**HINVERSE epsilon 0.05d0**
  HINV: epsilon   5.000000000000000E-002
**HINVERSE OUT GIN animal**
  HINV: how_houtput GIN
 animal.gin
**END**

 End of parameter file reached :  10 lines read
 Colleau's method specified for calculation of A22 - but not needed
 "In core" storage for NRM Inverse switched on
 ***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					   = ../pedigree.dat
 Total no. of animal IDs found			   =	      600
 no. of NRM matrices		1	    1
**************************** Count of "animal" types ****************************
 No. of non-genotyped animals			   =	      313
 No. of genotyped animals in total		   =	      287
 Total no. of "animals" 			   =	      600

 Begin: read M2 from ../MarkerCounts.dat       sys-time 	    0.0 	    0.0
 tvec4	      1000
 chevck marks	     1000	 1000
 No. of rows in ../MarkerCounts.dat read	   =	      287
 No. of marker allele counts per row read	   =	     1000
 End: read M2 from .dat 		       sys-time 	    0.1 	    0.1
 Successfully allocated arrays for in-core NRM

 Inbreeding coefficients for random effect	   =		1 computed
 ... TIER's algorithm
 Covariance structure				   =	NRM
 No. of inbred animals				   =	       24
 Average inbreeding coefficient 		   = 0.500000	  (in %)
 ... amongst inbred animals			   =  12.5000	  (in %)

 NRM Inverse no. :	     1	 "X"
 No. of elements =	  1780	 log determinant =  -277.258872223979
 A22: No. genotyped animals + ancestors 	   =	      368
 End: calculating A22 (Colleau) 	       sys-time 	    0.2 	    0.2
 End: centering M2			       sys-time 	    0.2 	    0.2
 End: calculating GRM			       sys-time 	    0.2 	    0.2
 end calc g22 F
 F
 End: inverting g22			       sys-time 	    0.2 	    0.2
 End: inverting a22			       sys-time 	    0.2 	    0.2
 Start writing out "animal.gin"
 H^{-1} : No. of non-zero elements in A^11, A^21   =	   42505.000
	 0.00000		  alpha 	  0.00000
 Total time used (secs): CPU =		 1  System =	       0
 Machine used = "didgeridoo"
 "WOMBAT" has finished !  26/05/2025


0.978u 0.129s 0:00.30 363.3%	0+0k 0+4184io 0pf+0w

Script done on 2025-05-26 17:31:14+1000
======================================================
Command used: wombat --hdiag animal.gin
======================================================

wombat --hdiag animal.gin
Script started on 2025-05-26 17:31:14+1000
 WOMBAT 26-05-2025 Maximum no. of threads to be used set to	      8

 WOMBAT: Run to invert a matrix only
 INVERSION: Input file	"			   = animal.gin"
    No. of elements read			   =	    42505
    Max. row no. found				   =	      600
    Max. column no. found			   =	      600
    Matrix storage mode 			   = SPAR
    Skip (single) "header" line 		   =  T

 Start of "AMD": no. of rows			   =	      600
 End of "AMD": no of subscripts 		   =	    85636
 End of "SMBFCT": maxsub			   =	     1612
		  maxlnz			   =	    42966
		  nnops 			   = 0.410327E+07
 "Super-nodes" :  No. of nodes			   =	      284
		  Max. size			   =	      289
 Output file: Diagonal elements of inverse	   = "animal.hdiags"
 Output file: Log determinant of input matrix	   = "animal.gin.det"
 Rank of matrix found				   =	      600
 Log determinant of input matrix		   =  667.102
 Operational zero used				   = 0.100000E-11
 Total time used (secs): CPU =		 0  System =	       0
 Machine used = "didgeridoo"
 "WOMBAT" has finished !  26/05/2025


0.195u 0.011s 0:00.17 117.6%	0+0k 0+88io 0pf+0w

Script done on 2025-05-26 17:31:14+1000
======================================================
Command used: wombat s1blup.par
======================================================

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

	   1 xone
	   2 xtwo
	  50	      50 0   1.00000
 File						   = "animal.gin"
 No. of records 				   =	    42505
 Denseness of matrix				   =  23.5746	 %
 No. of diagonal elements found 		   =	      600
 No. of zero elements found			   =		0
 No. of elements close to zero			   =		0    0.00%
 Log determinant given				   =  0.00000
 Have read					   =	      600 codes from file "animal.codes"
 Max. number of levels for "animal" set to	   =	      600
 No. of "subjects" in data			   =	      600
 No. of combinations of records 		   =		1
 No. of effect codes/ind. (maxnr)		   =		4
 Random effect no.	     0
  1   1.19670
  2   1.16150	    12.5249
 Eigenvalues   12.6428	     1.07884
 in %	       92.1377	     7.86232
 Random effect no.	     1
  1   1.19670
  2  0.666500	    3.13120
 Eigenvalues   3.33860	    0.989303
 in %	       77.1413	     22.8587
 adjact ivcol	  3911544     3911544
 Start of "GENMMD": no. of rows 		   =	     1814
 End of "GENMMD": no of subscripts		   =	    29150
 adjact ivcol	  3911544     3911544
 End of "SMBFCT": maxsub			   =	     8870	  1815
		  maxlnz			   =	   181101
		  nnops 			   = 0.331167E+08 0.100000E+16
 "Super-nodes" :  No. of nodes			   =	      591
		  Max. size			   =	      575
   1  animal		      1200  zero soln.s       0
 SS residuals	420.862       4661.99
 Total time used (secs): CPU =		 0  System =	       0
 Machine used = "didgeridoo"
 "WOMBAT" has finished !  26/05/2025


0.429u 0.140s 0:00.38 147.3%	0+0k 0+3632io 0pf+0w

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