Петля внутри другого цикла в R

Question

У меня проблема с результатами цикла в функции цикла. Он подсчитывает внутренний цикл только один раз и выбирает наилучшее решение для первого raw, а затем останавливается.

Я хотел бы запомнить лучшее решение для каждой строки матрицы zmienne. Что я делаю не так?

schaffer <- function(xx) 
{x1 <- xx[1] 
x2 <- xx[2] 

fact1 <- (sin(x1^2-x2^2))^2 - 0.5 
fact2 <- (1 + 0.001*(x1^2+x2^2))^2 

y <- 0.5 + fact1/fact2 
return(y) 
} 

gradient_descent <- function(func, step, niter) { 

    N <- 3 #N- number of random points 
    zmienne <- matrix(runif(N*2, min = -100, max = 100), N, 2) 
    print(zmienne) 
    h = 0.001; 
    iter_count = 0; 

    for (i in 1:N) { 
    x_0 <- zmienne[i,] 
    x_n = x_0; 

    for (j in 1:niter) { 
     func_grad = (func(x_n+h) - func(x_n))/h; 
     if (abs(func_grad) < 0.0001) { break; } 
     x_n = x_n - step * func_grad; 
     print(x_n) 
     iter_count = iter_count + 1 
    } 
    } 
return(list(iterations = niter, best_value = func_grad, best_state = x_n, x0=x_0)) 
    } 

solution_m1 <- gradient_descent(schaffer, 0.1, 20) 
solution_m1 

Answer 1

Возврат не должен находиться внутри внутреннего контура, но в конце функции.

Answer 2

Я думаю, что это то, что вы хотите:

gradient_descent <- function(func, step, niter) { 

    N <- 3 #N- number of random points 
    zmienne <- matrix(runif(N*2, min = -100, max = 100), N, 2) 
    print(zmienne) 
    h = 0.001; 
    iter_count = 0; 
    best.vals <- NULL 
    for (i in 1:N) { 
    x_0 <- zmienne[i,] 
    x_n = x_0; 

    for (j in 1:niter) { 
     func_grad = (func(x_n+h) - func(x_n))/h; 
     if (abs(func_grad) < 0.0001) { break; } 
     x_n = x_n - step * func_grad; 
     print(x_n) 
     iter_count = iter_count + 1 
    } 
    best.vals <- c(best.vals, func_grad) 
    } 

    return(list(iterations = iter_count, best_value = best.vals, best_state = x_n, x0=x_0)) 
} 

solution_m1 <- gradient_descent(schaffer, 0.1, 20) 
solution_m1