(Subscript[K, e])={((B^T)(Subscript[D, s])Bif   (Subscript[ψ, h])(e)≤0(B^T)(Subscript[D, v])Bif   (Subscript[ψ, h])(e)>0)