text
		R.pop_back(); - в этой строчке у меня вызывается исключение, но я не понимаю с чем это может быть связано так как на каждой итерации алгоритма у меня удаляется один интервал и вместо него добавляются два новых это просто невозможно чтобы вектор опустел: static __declspec(noalias)std::vector<float> Base_LNA_1_2_Mer_AGP(const uint32_t N,
const uint32_t global_iterations, const float a, const float b, const float c, const float d, const float r,
const bool mode, const float epsilon, const float seed) noexcept

{

text
if (N == 1u) {
	const float initial_length = b - a;
	float dmax = initial_length;
	const float start_val = ShekelFunc(a, seed);
	float best_f = ShekelFunc(b, seed);
	uint32_t schetchick = 0u;

	std::pair<float, float> x_Rmax = std::make_pair(a, b);
	std::pair<float, float> y_Rmax = std::make_pair(start_val, best_f);

	std::vector<float> Extr;
	Extr.reserve((global_iterations << 1u) + 4u);
	std::vector<Interval> R;
	R.reserve(global_iterations + 1u);

	R.emplace_back(std::make_pair(a, start_val), std::make_pair(b, best_f), N);
	float Mmax = R.back().M;
	float m = r * Mmax;

	const float threshold_03 = 0.3f * initial_length;
	const float inv_threshold_03 = 1.0f / threshold_03;

	while (true) {
		const float new_point = Shag(m, x_Rmax.first, x_Rmax.second, y_Rmax.first, y_Rmax.second, N, r);
		const float new_value = ShekelFunc(new_point, seed);
		const std::pair<float, float> promejutochnaya_tochka{ new_point, new_value };

		if (new_value < best_f) {
			best_f = new_value;
			Extr.emplace_back(best_f);
			Extr.emplace_back(new_point);
		}

		std::pop_heap(R.begin(), R.end(), Compare);
		Interval promejutochny_otrezok = std::move(R.back());
		R.pop_back();
		const size_t r_size = R.size();

		const float len2 = promejutochny_otrezok.end.first - new_point;
		const float len1 = new_point - promejutochny_otrezok.start.first;
		const float dmin = len1 < len2 ? len1 : len2;

		if (++schetchick == global_iterations || dmin < epsilon) {
			Extr.emplace_back(dmin);
			Extr.emplace_back(schetchick);
			return Extr;
		}

		Interval curr(promejutochny_otrezok.start, promejutochnaya_tochka, N);
		Interval curr1(promejutochnaya_tochka, promejutochny_otrezok.end, N);
		const float currM = curr.M > curr1.M ? curr.M : curr1.M;

		if (len2 + len1 == dmax) {
			dmax = len2 > len1 ? len2 : len1;

#pragma loop(ivdep)
for (size_t i = 0u; i < r_size; ++i) {
const float len = R[i].end.first - R[i].start.first;
if (len > dmax) {
dmax = len;
}
}
}

text
		if (mode) {
			if (threshold_03 > dmax && schetchick % 3u == 0u || 10.0f * dmax < initial_length) {
				if (currM > Mmax) {
					Mmax = currM;
					m = r * Mmax;
				}
				const float progress = fmaf(-inv_threshold_03, dmax, 1.0f);
				const float alpha = fmaf(progress, progress, 1.0f);
				const float betta = 2.0f - alpha;

				const float MULTIPLIER = (1.0f / dmax) * Mmax;
				const float global_coeff = fmaf(MULTIPLIER, r, -MULTIPLIER);
				const float GLOBAL_FACTOR = betta * global_coeff;

				curr.ChangeCharacteristic(fmaf(GLOBAL_FACTOR, len1, curr.M * alpha), N);
				curr1.ChangeCharacteristic(fmaf(GLOBAL_FACTOR, len2, curr1.M * alpha), N);

#pragma loop(ivdep)
for (size_t i = 0u; i < r_size; ++i) {
const float len_item = R[i].end.first - R[i].start.first;
R[i].ChangeCharacteristic(fmaf(GLOBAL_FACTOR, len_item, R[i].M * alpha), N);
}
std::make_heap(R.begin(), R.end(), Compare);
}
else {
if (currM > Mmax) {
Mmax = currM;
m = r * Mmax;
curr.ChangeCharacteristic(m, N);
curr1.ChangeCharacteristic(m, N);
#pragma loop(ivdep)
for (size_t i = 0u; i < r_size; ++i) {
R[i].ChangeCharacteristic(m, N);
}
std::make_heap(R.begin(), R.end(), Compare);
}
else {
curr.ChangeCharacteristic(m, N);
curr1.ChangeCharacteristic(m, N);
}
}
}
else {
if (currM > Mmax) {
Mmax = currM;
m = r * Mmax;
curr.ChangeCharacteristic(m, N);
curr1.ChangeCharacteristic(m, N);
#pragma loop(ivdep)
for (size_t i = 0u; i < r_size; ++i) {
R[i].ChangeCharacteristic(m, N);
}
std::make_heap(R.begin(), R.end(), Compare);
}
else {
curr.ChangeCharacteristic(m, N);
curr1.ChangeCharacteristic(m, N);
}
}

text
		R.emplace_back(std::move(curr));
		std::push_heap(R.begin(), R.end(), Compare);
		R.emplace_back(std::move(curr1));
		std::push_heap(R.begin(), R.end(), Compare);

		const Interval& top = R.front();
		x_Rmax = std::make_pair(top.start.first, top.end.first);
		y_Rmax = std::make_pair(top.start.second, top.end.second);
	}
}