updated pre_computed_set_argminmax with inface check

Author: WenjieXiao-2022

src/utils.jl

@@ -382,11 +382,9 @@ Base.length(storage::DeletedVertexStorage) = length(storage.storage)
 Computes the linear minimizer in the direction on the precomputed_set.
 Precomputed_set stores the vertices computed as extreme points v in each iteration.
 """
-function pre_computed_set_argminmax(pre_computed_set, direction, x)
+function pre_computed_set_argminmax(lmo, pre_computed_set, direction, x)
 	val = convert(eltype(direction), Inf)
 	valM = convert(eltype(direction), -Inf)
-	x_zeros = isa(x, AbstractMatrix) ? LinearIndices(x)[findall(iszero, x)] : findall(iszero, x)
-	# x_zeros = setdiff(collect(1:length(x)), SparseArrays.nonzeroinds(x))
 	idx = -1
 	idxM = -1
 	for i in eachindex(pre_computed_set)
@@ -395,8 +393,7 @@ function pre_computed_set_argminmax(pre_computed_set, direction, x)
 			val = temp_val
 			idx = i
 		end
-
-		if (isempty(x_zeros) || iszero(pre_computed_set[i][x_zeros])) && temp_val > valM
+		if is_inface_feasible(lmo, pre_computed_set[i], x) && temp_val > valM
 			valM = temp_val
 			idxM = i
 		end