R/postprocess.R

In rHanso: An R Implementation of Hybrid Algorithm for Non-Smooth Optimization (HANSO)

postprocess <- function(x, g, dnorm, X, G, w) {
    dist = c()
    X = as.matrix(X)
    G = as.matrix(G)
    w = as.matrix(w)
    loc = list()
    for (j in 1:ncol(X)) dist[j] = norm(as.matrix(x - X[, j]))
    evaldist = max(dist)
    mindist = min(dist)
    indx = which.min(dist)
    if (mindist == 0 & indx == 1) 
        {
        }  #nothing to do
 else if (mindist == 0 & indx > 1) {
        X[, c(1, indx)] = X[, c(indx, 1)]
        G[, c(1, indx)] = G[, c(indx, 1)]
        w[c(1, indx)] = w[c(1, indx)]
    } else {
        X = cbind(x, X)
        G = cbind(g, G)
        tmp = qpspecial(G)
        w = tmp$x
        d = tmp$d
        dnorm = norm(d)
    }
    loc$dnorm = dnorm
    loc$evaldist = evaldist
    return(list(loc = loc, X = X, G = G, w = w))
}