From d46079df0de741c0c6a7b86166794e73849a7807 Mon Sep 17 00:00:00 2001 From: Phillip Alday Date: Tue, 5 Mar 2024 11:29:31 -0600 Subject: [PATCH] format --- src/linalg.jl | 2 +- src/linalg/rankUpdate.jl | 2 +- src/linearmixedmodel.jl | 4 +++- 3 files changed, 5 insertions(+), 3 deletions(-) diff --git a/src/linalg.jl b/src/linalg.jl index 630ef0299..91e628153 100644 --- a/src/linalg.jl +++ b/src/linalg.jl @@ -18,7 +18,7 @@ function LinearAlgebra.mul!( αbnz = α * bnz[ib] jj = brv[ib] for ia in nzrange(A, j) - C[arv[ia], jj] = muladd(anz[ia], αbnz, C[arv[ia], jj]) + C[arv[ia], jj] = muladd(anz[ia], αbnz, C[arv[ia], jj]) end end end diff --git a/src/linalg/rankUpdate.jl b/src/linalg/rankUpdate.jl index 1b81dd35e..e6f1aacd6 100644 --- a/src/linalg/rankUpdate.jl +++ b/src/linalg/rankUpdate.jl @@ -22,7 +22,7 @@ function MixedModels.rankUpdate!( Cdiag = C.data.diag Adiag = A.diag @inbounds for idx in eachindex(Cdiag, Adiag) - Cdiag[idx] = muladd(β,Cdiag[idx], α * abs2(Adiag[idx])) + Cdiag[idx] = muladd(β, Cdiag[idx], α * abs2(Adiag[idx])) end return C end diff --git a/src/linearmixedmodel.jl b/src/linearmixedmodel.jl index 92cf3b786..3827f41ec 100644 --- a/src/linearmixedmodel.jl +++ b/src/linearmixedmodel.jl @@ -767,7 +767,9 @@ function StatsAPI.leverage(m::LinearMixedModel{T}) where {T} z = trm.z stride = size(z, 1) mul!( - view(rhs2, muladd((trm.refs[i] - 1), stride, rhsoffset) .+ Base.OneTo(stride)), + view( + rhs2, muladd((trm.refs[i] - 1), stride, rhsoffset) .+ Base.OneTo(stride) + ), adjoint(trm.λ), view(z, :, i), )