fix: use Nelder-Mead optimizer with fatol=5 instead of Powell/ftol

        x0[i] = 0.5 * sum(boundsval[i])

    method = "Nelder-Mead"

    method = "Powell"

    # method = "Powell"

    # ftol is tolerance for output power error of optimized point

    result = minimize(

        args=(optimize_keys, config, copcrit, isample),

        bounds=boundsval,

        method=method,

        options={"maxiter": 50, "ftol": 100},

        options={"maxiter": 50, "fatol": 5},

    )

    return result, config