print("Start")
  print(date())

res <- spfabia(paste("dataSim",run,"fabia",sep=""),p=SS,alpha=0.03,cyc=50,spl=0,non_negative=1,write_file=0,norm=0,lap=100.0,samples=(1:haploN),iter=A,quant=0.01,lowerB=1.1,upperB=1000,eps=1e-5)


  print("End Fabia")
  print(date())


sPF <- samplesPerFeature(X=paste("dataSim",run,"fabia",sep=""),samples=(1:haploN),lowerB=1.1,upperB=1000)


# save(res,sPF,annot,file=paste("resA",run,".Rda",sep=""))



off1 <- 0

resIBD1 <- getIBDs(res=res,sPF=sPF,annot=annot,chrom=1,labelsA=indi,ps=ps,psZ=psZ,inteA=inteA,thresA=thresA,minSNPs=minSNPs,off=off1)


if ( length(resIBD1) > 1) {
  comp <- compareIBDres(resIBD1,simv="minD")

  mergedIBD1 <- mergeIBDres(resIBD1=resIBD1,resIBD2=NULL,comp=comp,cut=0.8) 
} else {
  
 mergedIBD1 <- resIBD1 
}




off2=inteA%/%2

resIBD2 <- getIBDs(res=res,sPF=sPF,annot=annot,chrom=1,labelsA=indi,ps=ps,psZ=psZ,inteA=inteA,thresA=thresA,minSNPs=minSNPs,off=off2)


if ( length(resIBD2) > 1) {
  comp <- compareIBDres(resIBD2,simv="minD")

  mergedIBD2 <- mergeIBDres(resIBD1=resIBD2,resIBD2=NULL,comp=comp,cut=0.8) 
} else {
  mergedIBD2 <- resIBD2 
}

if ( length(mergedIBD1) > 0) {

  if ( length(mergedIBD2) > 0) {

    comp12 <- compareIBDres(mergedIBD1,mergedIBD2,simv="minD")

    mergedIBD <- mergeIBDres(resIBD1=mergedIBD1,resIBD2=mergedIBD2,comp=comp12,cut=0.8) 

} else {

mergedIBD <- mergedIBD1
}

} else {

  if ( length(mergedIBD2) > 0) {
    mergedIBD <- mergedIBD2
  } else {
    mergedIBD <- mergedIBD1
  }
  
  
}
  
  print("End HapFabia")

  print(date())


# save(mergedIBD,resIBD1,resIBD2,mergedIBD1,mergedIBD2,file=paste("resIBDA",run,".Rda",sep=""))