Merge pull request #347 from Neuroblox/sDCM_cleanup

sDCM code cleanup

Project.toml

@@ -27,6 +27,7 @@ JLD2 = "033835bb-8acc-5ee8-8aae-3f567f8a3819"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LogExpFunctions = "2ab3a3ac-af41-5b50-aa03-7779005ae688"
 MAT = "23992714-dd62-5051-b70f-ba57cb901cac"
+MKL = "33e6dc65-8f57-5167-99aa-e5a354878fb2"
 MLJ = "add582a8-e3ab-11e8-2d5e-e98b27df1bc7"
 MetaGraphs = "626554b9-1ddb-594c-aa3c-2596fe9399a5"
 ModelingToolkit = "961ee093-0014-501f-94e3-6117800e7a78"

src/Neuroblox.jl

@@ -14,6 +14,7 @@ ForwardDiff.can_dual(::Type{Complex{Float64}}) = true
 using ChainRules: _eigen_norm_phase_fwd!
 
 using LinearAlgebra
+using MKL
 using ToeplitzMatrices: Toeplitz
 using ExponentialUtilities: exponential!
 
@@ -85,32 +86,6 @@ abstract type BloxConnectMultiComplex <: BloxConnection end
 # dictionary type for Blox parameters
 Para_dict = Dict{Symbol, Union{<: Real, Num}}
 
-# struct types for Variational Laplace
-mutable struct VLState
-    iter::Int
-    v::Float64             # log ascent rate
-    F::Vector{Float64}
-    dF::Vector{Float64}
-    λ::Vector{Float64}
-    ϵ_θ::Vector{Float64}
-    reset_state::Vector{Any}
-    μθ_po::Vector{Float64}
-    Σθ_po::Matrix{Float64}
-    dFdθ::Vector{Float64}
-    dFdθθ::Matrix{Float64}
-end
-
-struct VLSetup
-    model_at_x0
-    y_csd::Array{Complex}
-    tolerance::Float64
-    systemnums::Vector{Int}
-    systemvecs::Vector{Vector{Float64}}
-    systemmatrices::Vector{Matrix{Float64}}
-    Q::Matrix{Complex}
-end
-
-
 include("utilities/spectral_tools.jl")
 include("utilities/learning_tools.jl")
 include("utilities/bold_methods.jl")