Missing `*` method for `PDiagMat`

[btmit]

ERROR: MethodError: *(::PDMats.PDiagMat{Float64, Vector{Float64}}, ::Diagonal{Float64, Vector{Float64}}) is ambiguous.

Candidates:
  *(A::AbstractMatrix, D::Diagonal)
    @ LinearAlgebra /Applications/Julia-1.9.app/Contents/Resources/julia/share/julia/stdlib/v1.9/LinearAlgebra/src/diagonal.jl:247
  *(a::PDMats.PDiagMat, x::AbstractMatrix)
    @ PDMats ~/.julia/packages/PDMats/CbBv1/src/pdiagmat.jl:51

Possible fix, define
  *(::PDMats.PDiagMat, ::Diagonal)

[RomeoV]

A workaround is to change

pdmat * diagmat

to

(diagmat * pdmat)'

since they are both symmetric.