Триангуляция Делоне

[himecu]

В коде, который использует библиотеки NumPy, Matplotlib, SymPy и SciPy, реализована функция plot_inequalities, которая должна строить график области пересечения двух неравенств, используя триангуляцию Делоне. Однако, при запуске кода, триангуляция Делоне не строится, не знается в чем может быть причина?

import numpy as np
import matplotlib.pyplot as plt
from sympy import *
from scipy.spatial import Delaunay
 
def plot_inequalities(inequality1, inequality2, x_min, x_max, y_min, y_max):
    x, y = symbols('x y')
 
    try:
        inequality1_expr = sympify(inequality1)
        inequality2_expr = sympify(inequality2)
    except:
        print("Ошибка: некорректный формат неравенства.")
        return
 
    try:
        F1 = lambdify((x, y), inequality1_expr, 'numpy')
        F2 = lambdify((x, y), inequality2_expr, 'numpy')
    except:
        print("Ошибка: не удалось скомпилировать неравенства.")
        return
 
    # Создание сетки значений x и y
    x_vals = np.concatenate(([0], np.linspace(x_min, x_max, 400)))
    y_vals = np.linspace(y_min, y_max, 400)
    X, Y = np.meshgrid(x_vals, y_vals)
    # Проверка неравенств в каждой точке сетки
    inequality1_result = F1(X, Y)
    inequality2_result = F2(X, Y)
 
    # Нахождение точек пересечения неравенств
    intersection_points = []
    for i in range(len(x_vals)):
        for j in range(len(y_vals)):
            if inequality1_result[j, i] < 0 and inequality2_result[j, i] > 0:
                intersection_points.append([x_vals[i], y_vals[j]])
 
    intersection_points = np.array(intersection_points)
 
    # Триангуляция Делоне по точкам пересечения
    if len(intersection_points) >= 3:
        tri = Delaunay(intersection_points[:, :2])
 
        # Построение графика
        plt.triplot(intersection_points[:, 0], intersection_points[:, 1], tri.simplices.copy(), color='black')
 
    # Отображение графика
    plt.xlabel('x')
    plt.ylabel('y')
    plt.grid(False)
    plt.axis('equal')
    plt.xlim(x_min, x_max)
    plt.ylim(y_min, y_max)
    plt.show()
 
# Запуск программы
print("Введите неравенство снизу в формате 'F(x, y) > 0':")
inequality1 = input()
print("Введите неравенство сверху в формате 'F(x, y) < 0':")
inequality2 = input()
print("Введите границы интервала x в формате 'x_min, x_max':")
x_min, x_max = map(float, input().split(','))
print("Введите границы интервала y в формате 'y_min, y_max':")
y_min, y_max = map(float, input().split(','))
 
plot_inequalities(inequality1, inequality2, x_min, x_max, y_min, y_max)

случай неравенств y-x>0 и y-x^2<0
Снимок.PNG

[BDA]

Просто треугольники маленькие и сливаются.

[p51x]

И проверьте знаки - вводите > и <, а используете < и >.

[himecu]

BDA, хорошо, большое спасибо вам, да просто разбиений было слишком много

[himecu]

p51x, спасибо, но при изменении нарушается триангуляция. не знаете в чем может быть проблема?
Снимок.PNG

[p51x]

в том что там "разрыв" https://www.wolframalpha.com/input?i...%3C0+and+y%3Ex

[himecu]

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