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