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