Блог

Я новичок в tcl, и я не мог понять, как перевести этот код с python на TCL

 import numpy as np

g0 = 7.88e12
Eox = np.array([155473, 15573, 1553, 1557473, 5473, 473, 1573, 19553])
E1 = 0.55e6
m = 0.7
fot= 1

D = float(input("rad dose"))

Fy = np.array([(abs(Eox)/(abs(Eox) E1))**m])
Not = np.array([g0*D*Fy*fot])
Not[Not>6.8e18] = 6.8e18
Nit = (1.7e4*D) 1e10
if Nit > 5e12:
    Nit = 5e12

print(Not)
print(Nit)
 

Tcl не обрабатывает списки значений, как numpy (это ближе к стандартному Python), без расширения. Я помню, что была такая штука, но я ею не пользовался и сейчас не могу вспомнить название. Поэтому я буду использовать стандартный Tcl. Ближайшей аналогией для массива numpy в стандартном Tcl является список значений (хотя это больше похоже на список или кортеж Python, чем на что-либо другое).

Я переведу вам эти четыре строки:

 Eox = np.array([155473, 15573, 1553, 1557473, 5473, 473, 1573, 19553])

Fy = np.array([(abs(Eox)/(abs(Eox) E1))**m])
Not = np.array([g0*D*Fy*fot])
Not[Not>6.8e18] = 6.8e18
 

Они становятся:

 set Eox {155473 15573 1553 1557473 5473 473 1573 19553}

set Fy [lmap EoxVal $Eox {expr {
    (abs($EoxVal) / abs($EoxVal   $E1)) ** $m
}}]
set Not [lmap FyVal $Fy {expr {
    $g0 * $D * $FyVal * $fot
}}]
set Not [lmap NotVal $Not {expr {
    ($NotVal > 6.8e18) ? 6.8e18 : $NotVal
}}]
 

Использование lmap делает то же самое, что numpy делает за кулисами для вас в этих массивах. Однако последние две команды можно объединить:

 set Not [lmap FyVal $Fy {expr {
    min($g0 * $D * $FyVal * $fot, 6.8e18)
}}]
 

(Кроме того, я разделил каждую lmap/expr команду на несколько строк для большей ясности. Вы можете написать более компактный код, но его труднее читать. Очистить код-очень хороший план для любого кода, который вы собираетесь хранить дольше 5 минут.)

Остальные строки кода являются чисто базовыми. set , puts , gets и expr будет делать то, что вам нужно. И подозрительно очевидное if .