revise para_type

src/AsymptoticBound/CramerRao.jl

@@ -4,18 +4,23 @@ abstract type SLD <: AbstractDtype end
 abstract type RLD <: AbstractDtype end
 abstract type LLD <: AbstractDtype end
 
-struct QFIM_Obj{D} <: AbstractObj
+struct QFIM_Obj{P,D} <: AbstractObj
     W::AbstractMatrix
     eps::Number
 end
 
-struct CFIM_obj <: AbstractObj
+struct CFIM_obj{P} <: AbstractObj
     M::AbstractVecOrMat
     W::AbstractMatrix
     eps::Number
 end
 
-struct HCRB <: AbstractObj
+struct HCRB{P} <: AbstractObj
     W::AbstractMatrix
     eps::Number
 end
+
+QFIM_Obj(W, eps, syms::Symbol...) = QFIM_Obj{eval.(syms)...}(W,eps)
+CFIM_Obj(M, W, eps, syms::Symbol...) = CFIM_Obj{eval.(syms)...}(M,W,eps)
+HCRB_Obj(W, eps, syms::Symbol...) = HCRB_Obj{eval.(syms)...}(W,eps)
+

src/AsymptoticBound/CramerRaoBase.jl

@@ -59,6 +59,7 @@ function SLD_liouville(
     ∂ρ_∂x::Vector{Matrix{T}};
     eps = eps_default,
 ) where {T<:Complex}
+
     (x -> SLD_liouville(ρ, x; eps = eps)).(∂ρ_∂x)
 end
 
@@ -187,6 +188,7 @@ function QFIM(
     dtype = :SLD,
     eps = eps_default,
 ) where {T<:Complex}
+
     Symbol("QFIM_" * string(dtype))(ρ, dρ; eps = eps)
 end
 
@@ -197,6 +199,7 @@ function QFIM(
     dtype = :SLD,
     eps = eps_default,
 ) where {T<:Complex}
+
     Symbol("QFIM_" * string(dtype))(ρ, dρ; eps = eps)
 end
 
@@ -211,6 +214,7 @@ function QFIM(
     dtype = :SLD,
     eps = eps_default,
 ) where {T<:Complex}
+
     F = QFIM(ρ, dρ; dtype = dtype, eps = eps)
     LD = eval(dtype)(ρ, dρ; eps = eps)
     return F, LD
@@ -224,6 +228,7 @@ function QFIM(
     dtype = :SLD,
     eps = eps_default,
 ) where {T<:Complex}
+
     F = QFIM(ρ, dρ; dtype = dtype, eps = eps)
     LD = eval(dtype)(ρ, dρ; eps = eps)
     return F, LD
@@ -346,11 +351,13 @@ function QFIM_Gauss(R̄::V, dR̄::VV, D::M, dD::VM) where {V,VV,M,VM<:AbstractVe
             i = 1:quad_num, j = 1:quad_num
         ] for k = 1:para_num
     ]
+
     S, cs = Williamson_form(C)
     Gs = G_Gauss(S, dC, cs)
     F = [
         tr(Gs[i] * dC[j]) + transpose(dR̄[i]) * inv(C) * dR̄[j] for i = 1:para_num,
         j = 1:para_num
     ]
+
     F |> real
 end

src/dynamics.jl

@@ -1,10 +1,6 @@
 abstract type AbstractDynamics end
 abstract type AbstractDynamicsData end
 
-abstract type AbstracParaType end
-abstract type single_para <: AbstracParaType end
-abstract type multi_para <: AbstracParaType end
-
 abstract type AbstractStateType end
 abstract type ket <: AbstractStateType end
 abstract type dm <: AbstractStateType end

src/objective.jl

@@ -1,5 +1,8 @@
 abstract type AbstractObj end
 
+abstract type AbstracParaType end
+abstract type single_para <: AbstracParaType end
+abstract type multi_para <: AbstracParaType end
 
 include("AsymptoticBound/CramerRao.jl")
 include("AsymptoticBound/Holevo")