Исправление сегментационной ошибки

Ищу причины

Сейчас я разбираюсь в некоторых основных причинах сегментационных сбоев в C++: неинициализированные переменные, недопустимые указатели или выход за пределы массива.

Проверка расчетов

Теперь я смотрю на ввод программы, чтобы убедиться, что переменные читаются корректно. Вычисляю стороны треугольника и проверяю неравенство треугольника.

Вычисляю радиусы

Я нахожу радиусы вписанной и описанной окружностей, высоты, медианы и синусы углов треугольника, проверяя их точность. Затем вывожу результаты.

Ищу причины

Я размышляю о причинах возникновения сегментационного сбоя, включая деление на ноль или другие неверные операции, и подозреваю, что возможные нулевые длины сторон могут быть связаны с коллинеарностью точек.

Проверяю корректность

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

Выявление проблемы

Похоже, проблема в подсчете площади треугольника: вычисления s2 и s3 оказываются нерациональными.

Использование операторов

Оператор 'std::abs' применяется к long double. Если стороны равны нулю, это может привести к делению на ноль.

Рассматриваю возможность наличия неопределенного поведения, как по мне, так и из-за множительной идентичности.

Проверяю формулу

Проверяю формулу для bisector длины, используя длины сторон. Интересно, как вычисляется результат, применяя my_sqrt.

Проверяю вычисления

Запускаю процесс выяснения правильного формулы для длины биссектора угла, используя различные параметры длины сторон и периметра.

Проверяю корректность

Анализирую функции 'my_arcsin' и 'my_sqrt', чтобы проверить корректность входных данных и обрабатывать исключительные значения корректно.

Применение квадрата суммы сторонுள்ள ஐ्ट्रங்க பயன்படுத்துவது interesting idea. This method takes the squares of the magnitudes of the double vectors and sums them.

Заменяю π на его значение и вот эту красивую формулу и несколько уточнений. Интересно, что это может изменить в расчетах.

Рассматриваю деление

Окей, я заметил, что в делении 2 / (side1 + side2) может произойти целочисленное деление. Лучше написать 2.0 для точности.

Понимание ситуации

Я разбираюсь с использованием функции std::abs для long double аргументов в C++. Стандартная функция определяет её для long double, и это не должно вызвать проблем.

Я разбираюсь с формулами для вычисления длины биссектрисы. Интересно, как конкретная формула влияет на результаты.

Реагировать на варианты

Если пользователь вводит ненумерический ввод, std::cin установит failbit, и переменные не будут инициализированы. Они могут принимать нулевые или неизмененные значения, что приводит к делению на ноль или нарушению неравенства треугольника, но не к сегментационному файлу.

Проверяю функцию возьмем аргумент функции в диапазоне от -1 до 1.

Исправляю формулы

Проверяю формулы для высот и углов в радианах. Убеждаюсь, что использование sinyc как угла, а не синуса, правильно. Это подчеркивает точность.

рекомендации Синуса угла в коде неправильно используют. Обычно у нас sin(yc), а не сам yc.

Понимание ошибки возникновения аварийной остановки

Ищу альтернативы

Я анализирую формулы для медиан и вдохновляюсь, ведь разные подходы дают новые идеи.

Проверяю возможности

Портуйте руки и поднимайте задачи. Am I passing invalid args to 'my_arcsin' or seeing 'my_sqrt' convergence issues?

Проверяю и замечаю

В функции my_sqrt я обсуждаю точность и максимальное количество итераций, но обсуждаю допустимость нуля на входе.

Создаю корректную формулу

Я нахожу, что в вычислениях для 'b1', 'b2', 'b3' формула может быть неверной, приводя к отрицательному аргументу в 'my_sqrt', что вызывает обрывок.

plaintext
Пробую понять, как улучшить формулу и вычисления углов. Изменяю подход, хотя не уверен в конечном результате. Это приводит к интересным мыслям о возможных методах.

**Проверяю формулы**

Я проверяю формулы для угла и биссектора в треугольнике. Убедился, что все выражения положительны, и 'my_sqrt' получает положительное значение.

**Программирование и ошибки**

Ищи несоответствия в коде, проверяю использование 'my_sqrt' с 'max_iter', чтобы выключить бесконечный цикл.

Предполагаю, что код работает, но возникает деление на ноль. Возможные причины: деление int на ноль, uninitialized переменные, или ошибки ввода.

Проверяю ввод и это приводит к ошибкам. Важно указывать шесть чисел, даже наоборот.

**Ищем источник проблемы**

Сегментация может происходить из-за вводов 'std::cin', переменные могут получать мусорные значения. Приводит к ошибкам в s3 и бисекторах, добавление валидации ввода поможет.



**Проверка и корректировка**

Я разбираюсь в инициализации и корректности переменных. Возможно, преобразование 'long double' в 'double' улучшит результат.

**Анализируя возможные причины**

Сейчас я разбираюсь с исключением разделения на ноль и инициализацией переменных в программе, чтобы понять, из-за чего падает программа.

**Ищем причину**

Похоже, что формулы могут быть некорректными, что ведет к ошибкам. Можем заглянуть в ячейки и формулы, которые их содержат, и провести проверку на правильность.

Сегментационная ошибка (segmentation fault) обычно возникает из-за доступа к недопустимой памяти, однако в предоставленном коде явных указателей или операций, приводящих к таким ошибкам, нет. Вместо этого, проблема, скорее всего, связана с логическими ошибками или некорректными вычислениями, которые могут привести к неопределенному поведению, особенно при работе с математическими функциями.

Ниже приведены несколько потенциальных причин и соответствующие исправления:

1. Некорректное использование арксинуса:

   • Проблема: В коде вы используете yc, ya и yb как углы в радианах, но затем пытаетесь использовать их напрямую в вычислениях, где ожидается значение синуса угла.
   • Решение: Вместо использования самих углов (ya, yb, yc) используйте их синусы (sinya, sinyb, sinyc) для дальнейших вычислений.

2. Некорректная формула для биссектрис:

   • Проблема: Формула для вычисления биссектрис в коде неверна.

   • Решение: Используйте корректные формулы для вычисления биссектрис. Один из вариантов формулы для длины биссектрисы из вершины C:

     $$b_c = \frac{2ab \cos(\frac{C}{2})}{a + b}$$

     Однако для простоты можно использовать формулу:

     $$b_c = \sqrt{ab \left(1 - \left(\frac{c}{a + b}\right)^2\right)}$$

3. Отсутствие проверки на корректность ввода:

   • Проблема: Если пользователь вводит некорректные данные, переменные могут содержать неожиданные значения.
   • Решение: Добавьте проверки после ввода данных, чтобы убедиться, что ввод прошел успешно.

4. Использование long double:

   • Проблема: В большинстве случаев double предоставляет достаточную точность, и использование long double может привести к неожиданным результатам на некоторых компиляторах или архитектурах.
   • Решение: Используйте double вместо long double для большей совместимости.

Ниже представлен исправленный код с учетом вышеуказанных рекомендаций:

cpp
#include <iostream>
#include <cmath>    // Для математических функций
#include <limits>   // Для проверки на ввод

// Функция вычисления квадрата числа
double kvadr(double a)
{
    return a * a;
}

// Функция вычисления модуля разности двух чисел
double absrazn(double a, double b)
{
    return std::abs(a - b);
}

// Функция вычисления квадратного корня (на основе метода Ньютона)
double my_sqrt(double n)
{
    if (n < 0)
    {
        std::cerr << "Отрицательный аргумент для квадратного корня.\n";
        return std::numeric_limits<double>::quiet_NaN();
    }
    if (n == 0)
    {
        return 0;
    }
    else
    {
        const double tochn = 1E-12;
        double x = n;
        int max_iter = 10000; // Максимальное количество итераций
        int iter = 0;

        while (iter < max_iter)
        {
            double nx = 0.5 * (x + n / x);
            if (absrazn(x, nx) < tochn)
                return nx;
            x = nx;
            iter++;
        }

        std::cerr << "Функция my_sqrt не сошлась\n";
        return x;
    }
}

// Функция вычисления арксинуса с проверкой
double my_arcsin(double x)
{
    if (x < -1 || x > 1)
    {
        std::cerr << "Аргумент функции my_arcsin должен быть в диапазоне [-1, 1]\n";
        return std::numeric_limits<double>::quiet_NaN();
    }

    return std::asin(x); // Используем стандартную функцию для точности и надежности
}

const double pi = 3.141592653589793;

// Основная функция
int main()
{
    double x1, y1, x2, y2, x3, y3;

    std::cout << "Введите координаты первой вершины треугольника: ";
    if (!(std::cin >> x1 >> y1))
    {
        std::cerr << "Некорректный ввод.\n";
        return 1;
    }

    std::cout << "Введите координаты второй вершины треугольника: ";
    if (!(std::cin >> x2 >> y2))
    {
        std::cerr << "Некорректный ввод.\n";
        return 1;
    }

    std::cout << "Введите координаты третьей вершины треугольника: ";
    if (!(std::cin >> x3 >> y3))
    {
        std::cerr << "Некорректный ввод.\n";
        return 1;
    }

    // Посчитаем длины сторон треугольника
    double side1 = my_sqrt(kvadr(x2 - x1) + kvadr(y2 - y1));
    double side2 = my_sqrt(kvadr(x3 - x1) + kvadr(y3 - y1));
    double side3 = my_sqrt(kvadr(x2 - x3) + kvadr(y2 - y3));

    // Проверка на существование треугольника
    if ((side1 + side2 <= side3) || (side1 + side3 <= side2) || (side3 + side2 <= side1))
    {
        std::cout << "Треугольник не существует\n";
        return 0;
    }

    // Посчитаем периметр и площадь треугольника
    double per = side1 + side2 + side3;
    double s = per / 2.0;
    double area_arg = s * (s - side1) * (s - side2) * (s - side3);
    double s1 = my_sqrt(area_arg);

    // Проверка на корректность площади
    if (std::isnan(s1))
    {
        std::cerr << "Ошибка при вычислении площади треугольника.\n";
        return 1;
    }

    // Посчитаем радиусы вписанной и описанной окружностей треугольника
    double rvpis = s1 / s;
    double ropis = (side1 * side2 * side3) / (4.0 * s1);

    // Посчитаем площадь и периметр этих окружностей
    double svpis = pi * kvadr(rvpis);
    double sopis = pi * kvadr(ropis);
    double pervpis = 2 * pi * rvpis;
    double peropis = 2 * pi * ropis;

    // Посчитаем длины высот треугольника
    double h1 = (2.0 * s1) / side1;
    double h2 = (2.0 * s1) / side2;
    double h3 = (2.0 * s1) / side3;

    // Посчитаем длины медиан треугольника
    double m1 = my_sqrt(2 * kvadr(side2) + 2 * kvadr(side3) - kvadr(side1)) / 2.0;
    double m2 = my_sqrt(2 * kvadr(side1) + 2 * kvadr(side3) - kvadr(side2)) / 2.0;
    double m3 = my_sqrt(2 * kvadr(side1) + 2 * kvadr(side2) - kvadr(side3)) / 2.0;

    // Найдем синусы углов треугольника
    double sinya = h3 / side2; // sin(A) = высота на сторону BC / сторона AB
    double sinyb = h3 / side1; // sin(B) = высота на сторону AC / сторона AB
    double sinyc = h1 / side2; // sin(C) = высота на сторону AB / сторона AC

    // Найдем арксинусы углов — сами углы
    double ya = my_arcsin(sinya);
    double yb = my_arcsin(sinyb);
    double yc = my_arcsin(sinyc);

    // Проверка на корректность углов
    if (std::isnan(ya) || std::isnan(yb) || std::isnan(yc))
    {
        std::cerr << "Ошибка при вычислении углов треугольника.\n";
        return 1;
    }

    // Посчитаем градусные меры углов
    double gradya = ya * (180.0 / pi);
    double gradyb = yb * (180.0 / pi);
    double gradyc = yc * (180.0 / pi);

    // Посчитаем длины биссектрис треугольника
    // Используем корректную формулу для биссектрисы
    auto bisector = [&](double a, double b, double c) -> double {
        return (2.0 * a * b * std::cos(my_arcsin( (2.0 * s1) / (a * b) ))) / (a + b);
    };

    double b1 = bisector(side2, side3, side1);
    double b2 = bisector(side1, side3, side2);
    double b3 = bisector(side1, side2, side3);

    // Посчитаем площадь треугольника различными способами
    double s2 = s1; // Площадь треугольника через радиус вписанной окружности
    double s3 = (side1 * side2 * std::sin(ya)) / 2.0; // Площадь через две стороны и угол между ними

    // Вывод результатов
    std::cout << "\n--- Результаты ---\n";

    std::cout << "Стороны треугольника: " << side1 << ", " << side2 << ", " << side3 << "\n";
    std::cout << "Периметр треугольника: " << per << "\n";
    std::cout << "Площадь треугольника: " << s1 << "\n";

    std::cout << "Радиус вписанной окружности: " << rvpis << "\n";
    std::cout << "Радиус описанной окружности: " << ropis << "\n";

    std::cout << "Площадь вписанной окружности: " << svpis << "\n";
    std::cout << "Площадь описанной окружности: " << sopis << "\n";

    std::cout << "Периметр вписанной окружности: " << pervpis << "\n";
    std::cout << "Периметр описанной окружности: " << peropis << "\n";

    std::cout << "Длины высот треугольника: " << h1 << ", " << h2 << ", " << h3 << "\n";
    std::cout << "Длины медиан треугольника: " << m1 << ", " << m2 << ", " << m3 << "\n";
    std::cout << "Длины биссектрис треугольника: " << b1 << ", " << b2 << ", " << b3 << "\n";

    std::cout << "Градусные меры углов треугольника: " << gradya << "°, " << gradyb << "°, " << gradyc << "°\n";
    std::cout << "Радианные меры углов треугольника: " << ya << ", " << yb << ", " << yc << " радиан\n";

    std::cout << "Площадь треугольника различными методами: " << s1 << ", " << s2 << ", " << s3 << "\n";

    return 0;
}

Изменения и пояснения:

1. Замена long double на double:

   • Использование double обеспечивает большую совместимость и обычно достаточную точность для большинства задач.

2. Добавление проверок ввода:

   • После каждого std::cin добавлены проверки, чтобы убедиться в корректности введенных данных. Если ввод некорректен, программа выводит ошибку и завершается.

3. Исправление вычисления площади s3:

   • Вместо использования самого угла yc в формуле площади, теперь используется std::sin(ya).

4. Исправление формулы для биссектрис:

   • Использована корректная формула для вычисления длины биссектрисы с учетом углов треугольника.

5. Улучшение функции my_sqrt:

   • Добавлена проверка на отрицательные аргументы, чтобы избежать неопределенного поведения при попытке вычислить квадратный корень из отрицательного числа.

6. Обработка возможных ошибок при вычислениях:

   • Добавлены проверки на NaN после вычисления площади и углов, чтобы предотвратить дальнейшие ошибки в программе.

Рекомендации:

• Тестирование: После внесения изменений рекомендуется протестировать программу с различными входными данными, включая граничные случаи (например, треугольники с нулевой площадью, равносторонние треугольники и т.д.).

• Использование стандартных функций: В некоторых случаях можно использовать стандартные функции C++ для вычисления геометрических свойств треугольника, что повысит надежность и точность программы.

• Логирование ошибок: Для более детальной отладки можно добавить дополнительное логирование, чтобы отслеживать значения переменных на различных этапах вычислений.

Если после внесения этих изменений проблема с сегментационной ошибкой сохраняется, рекомендуется использовать отладчик (например, gdb) для пошагового выполнения программы и определения точного места возникновения ошибки.