Ай, Карамба, у меня RuntimeWarning: invalid value encountered in sqrt

@InsomniacG · Регистрация: 28.05.2018

Студворк — интернет-сервис помощи студентам

Привет всем,очень нужна помощь в устранении ошибки.
Вроде все логично,программа строит нужные графики,считает все параметры абсолютно верно,но в какой-то момент выдает ошибку в вычислении функции( значения положительной половины окружности в координатах),причем зависимость я так и не смог найти...Входные данные сугубо положительные,типа float .
Буду благодарен за помощь,ибо сам не справляюсь(
Есть подозрения,что беда в np.linspace,который является аргументом Mrc ,что приводит к RuntimeWarning: invalid value encountered in sqrt
y =np.sqrt(r * r - (x - x0) * (x - x0)),
тип данных мб не тот,но float и комплекс не решают эту проблему

Python
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
     def Mrc(x, r, x0):
         y =np.sqrt(r * r - (x - x0) * (x - x0))
         return y
Dlunk = float(input())
     Dotv = float(input())
     a = float(input())
     b = float(input())
     Pkgs = float(input())
     S=a*b*N/100
     F=np.pi* Dlunk * Dotv/100
     Sit = float(Pkgs / (10*S)) 
     Pv = Pkgs / (10*F)
 
     Ki = Pv/Sit
TanFi0 =0.5*(Pv-Sit)/np.sqrt(Sit*Pv) 
     Fi0=np.arctan(TanFi0)*180/np.pi
     C0i = np.sqrt(Sit * Pv)  
     TanFi=(np.sqrt(Pv)-np.sqrt(Sit))*0.5/(np.sqrt(np.sqrt(Pv*Sit)))
     Fi = np.arctan(TanFi) * 180 / np.pi
     Ci=np.sqrt(np.sqrt(Pv*Sit))*(np.sqrt(Pv)+np.sqrt(Sit))*0.5
     Sis=Pv+np.sqrt(Sit*Pv)
     Sivs=2*Pkgs / (10*(S+F))
     Khr=Sis/Sivs
     KPuass=1.4/(Khr-2.05)
     Si3s=2*Pkgs
     Xc1 = (-1) * Sit
     R1 = Sit
     Xc2 = 0
     R2 = Sit * Pv * 2/ (Sit + Pv)
     Xc3 = (Pv - Sit) /2
     R3 = (Pv - Sit) /2
     Xc4 = (Pv + C0i) / 2
     R4 = (Pv + C0i) / 2
     Xc5=(Pv+2*C0i+Sit)/2
     R5=(Pv+2*C0i-Sit)/2
     Xc6=Pv+C0i
     R6=Pv
 
     Xc7=(Pv+C0i)*(2+np.sqrt(Pv/Sit))/2
     R7=(Pv+C0i)*np.sqrt(Pv/Sit)/2
     Xc8=(Pv+C0i)*(Pv/Sit+np.sqrt(Pv/Sit)-2)/4
     R8=3*C0i/2+((Pv-3*Sit)*Pv/Sit)/4
     Xc9 = (Pv + Sit) / 2
     R9 = Sit + Pv / 2
 
 
    X_l1 = np.linspace(Xc1-R1,Xc1+R1, 1000)#графика кругов 1-8
     Y1 = Mrc(X_l1, R1, Xc1)
     X_l2 = np.linspace(Xc2 - R2, Xc2 + R2, 100000)
     Y2 = Mrc(X_l2, R2, Xc2)
     X_l3 = np.linspace(Xc3 - R3, Xc3 + R3, 100000)
     Y3 = Mrc(X_l3, R3, Xc3)
     X_l4 = np.linspace(Xc4 - R4, Xc4 + R4, 100000)
     Y4 = Mrc(X_l4, R4, Xc4)
     X_l5 = np.linspace(Xc5 - R5, Xc5 + R5, 100000)
     Y5 = Mrc(X_l5, R5, Xc5)
     X_l8 = np.linspace(Xc8 - R8, Xc8 + R8, 100000)
     Y8 = Mrc(X_l8, R8, Xc8)
     X_l6 = np.linspace(Xc6 - R6, Xc6 + R6, 100000)
     Y6 = Mrc(X_l6, R6, Xc6)
     X_l9 = np.linspace(Xc9 - R9, Xc9 + R9, 100000)
     Y9 = Mrc(X_l9, R9, Xc9)
     X_l7 = np.linspace(Xc7 - R7, Xc7 + R7, 100000)
     Y7 = Mrc(X_l7, R7, Xc7)  # конец кругов
     X_l9 = np.linspace(Xc9 - R9, Xc9 + R9, 100000)
     Y9 = Mrc(X_l9, R9, Xc9)

@Welemir1 · 20.06.2018, 06:42

InsomniacG, ну поотслеживай, поставь перехват исключения или перед вычислением корня выводи в консоль r * r - (x - x0) * (x - x0) просто посмотреть не будет ли он отрицателен,

Сообщение от InsomniacG

Входные данные сугубо положительные,типа float .

а это вполне может быть и при положительных аргументах, если х-х0 будет больше r

@InsomniacG · 20.06.2018, 11:33 **[ТС]**

Welemir1 , я пытался через try/except (Runtime Warning), но все равно то же самое выдает. В исключении домножал на -1 подкоренное выражение.

@Welemir1 · 20.06.2018, 11:46

InsomniacG, юзай дебаггер или просто вот в лоб

Python
1
2
3
4
5
6
def Mrc(x, r, x0):
         z = r * r - (x - x0) * (x - x0)
         if z<0:
              print (z, x,r,x0)
         y =np.sqrt(x)
         return y

и в консоли потом поглядишь что там за параметры при которых падает

@InsomniacG · 20.06.2018, 12:50 **[ТС]**

К сожалению,для массивов это не работает...Traceback (most recent call last):
File "C:/project geo/MAINGT — копия.py", line 164, in <module>
Y1 = Mrc(X_l1, R1, Xc1)
File "C:/project geo/MAINGT — копия.py", line 153, in Mrc
if z < 0:
ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .

@InsomniacG 0 / 0 / 0 Регистрация: 28.05.2018 Сообщений: 70
	20.06.2018, 11:33 [ТС]
	Welemir1 , я пытался через try/except (Runtime Warning), но все равно то же самое выдает. В исключении домножал на -1 подкоренное выражение. 0

@InsomniacG 0 / 0 / 0 Регистрация: 28.05.2018 Сообщений: 70
	20.06.2018, 12:50 [ТС]
	К сожалению,для массивов это не работает...Traceback (most recent call last): File "C:/project geo/MAINGT — копия.py", line 164, in <module> Y1 = Mrc(X_l1, R1, Xc1) File "C:/project geo/MAINGT — копия.py", line 153, in Mrc if z < 0: ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all() 0