Skip to content

Commit

Permalink
format
Browse files Browse the repository at this point in the history
  • Loading branch information
palday committed Mar 5, 2024
1 parent b10b801 commit d46079d
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 3 deletions.
2 changes: 1 addition & 1 deletion src/linalg.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
2 changes: 1 addition & 1 deletion src/linalg/rankUpdate.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
4 changes: 3 additions & 1 deletion src/linearmixedmodel.jl
Original file line number Diff line number Diff line change
Expand Up @@ -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),
)
Expand Down

0 comments on commit d46079d

Please sign in to comment.