From 0a67dc497244d9ef4d0947fb4cde429d6854f0ef Mon Sep 17 00:00:00 2001
From: PharmCat <v.s.arnautov@yandex.ru>
Date: Sun, 21 Feb 2021 23:06:42 +0300
Subject: [PATCH] Inf fix

---
 src/reml.jl | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/reml.jl b/src/reml.jl
index 07beb9f1..fd69828f 100644
--- a/src/reml.jl
+++ b/src/reml.jl
@@ -101,7 +101,7 @@ function reml_sweep_β(lmm, data::AbstractLMMDataBlocks, θ::Vector{T}) where T
             θ₁  += logdetv(V)
         catch
             lmmlog!(lmm, LMMLogMsg(:ERROR, "θ₁ not estimated during REML calculation, V isn't positive definite or |V| less zero."))
-            return (Inf, nothing, nothing, Inf)
+            return (1e100, nothing, nothing, 1e100)
         end
         #sweepb!(view(akk, 1:qswm), Vp, 1:q)
         sweep!(Vp, 1:q)
@@ -120,7 +120,7 @@ function reml_sweep_β(lmm, data::AbstractLMMDataBlocks, θ::Vector{T}) where T
         logdetθ₂ = logdet(θ₂)
     catch
         lmmlog!(lmm, LMMLogMsg(:ERROR, "logdet(θ₂) not estimated during REML calculation"))
-        return (Inf, nothing, nothing, Inf)
+        return (1e100, nothing, nothing, 1e100)
     end
     return   θ₁ + logdetθ₂ + θ₃ + c, β, θ₂, θ₃ #REML, β, iC, θ₃
 end
@@ -158,7 +158,7 @@ function reml_sweep_β(lmm, data::AbstractLMMDataBlocks, θ::Vector{T}, β::Vect
             θ₁  += logdetv(V)
         catch
             lmmlog!(lmm, LMMLogMsg(:ERROR, "θ₁ not estimated during REML calculation, V isn't positive definite or |V| less zero."))
-            return (Inf, nothing, nothing, Inf)
+            return (1e100, nothing, nothing, 1e100)
         end
         sweep!(Vp, 1:q)
         V⁻¹ = Symmetric(utriaply!(x -> -x, V))
@@ -171,7 +171,7 @@ function reml_sweep_β(lmm, data::AbstractLMMDataBlocks, θ::Vector{T}, β::Vect
     try
         logdetθ₂ = logdet(θ₂)
     catch
-        return (Inf, nothing, nothing, Inf)
+        return (1e100, nothing, nothing, 1e100)
     end
     return   θ₁ + logdetθ₂ + θ₃ + c, θ₂, θ₃ #REML, iC, θ₃
 end