如何从 Julia 中的 if 语句中打破 for 循环？

Question

我正在尝试将 python 模型化（不是很有效）转换为 Julia 代码（在这种情况下更有用），但我没有找到如何打破 for 循环以及如何打印此上下文中的值。

Python代码

for i in np.range(1, nt)
    if t[x == s] == 100:
        s_time = i*dt
        print(s_time)
        break

我试过的 Julia 代码：

for i = 1:nt
    if round.(Int, t[x.==s]) == 100
        s_time = i*dt
        println(s_time)
        break
    end
end

这里的变量并不重要，但是dt是一个值，可以让我在i的for循环中知道[y]的时间，而t[x.==s] == 100 让我知道什么时候我想打破循环。

Julia 代码不会 returns 我出现任何错误，但不会中断也不会打印任何内容。

完整代码（针对 Jupyter Notebook 用户）：

# NUMERICS
ti = 30  # surface temperature
td = 1000  # dyke temperature
tel = 30  # right lim temperature
teg = tel  # left lim temperature
x1 = 80  # central position dyke 2
x2 = 110  # central position dyke 1
d1 = 24  # dyke diameter
d2 = 10  # dyke diameter
s = 160  # probe position

yd1 = 0  # dyke 1 intrusion year
yd2 = 15  # dyke 2 intrusion year

d = 40  # rock diffusivity
size = 200  # size of the model
ttot = 30;  # total time in year

nx = 2001  # cells number
dx = size/(nx-1)  # cells size
dt = dx^2/(2.1*d)  # frame condition
nt = round(Int, ttot/dt)  # frame number
x = LinRange(0, size, nx); # size model

# INITIALISATION
indexing_1_l = round(Int, (x1 - d1/2)*nx/size)
indexing_1_r = round(Int, (x1 + d1/2)*nx/size)
indexing_2_l = round(Int, (x2 - d2/2)*nx/size)
indexing_2_r = round(Int, (x2 + d2/2)*nx/size);

t = ti * ones(nx)
t[indexing_1_l:indexing_1_r] .= td;
using Plots
gr() ;

# PLOTTING LOOP
for i = 1:nt 
    
    q = -d*(diff(t)/dx)
    dcdt = -diff(q)/dx
    t[2:end-1] = t[2:end-1] + dcdt*dt
    
    t[1] = teg
    t[end] = tel
    
    if i == round(Int, yd2/dt)
        t[indexing_2_l:indexing_2_r] .= td
    end
    
    # IF STATEMENT THAT DOESN'T WORK
    if round.(Int, t[x.==s]) == 100
        s_time = i*dt
        println("$s_time")
        break
    end

    if mod(i, 1000) == 0
        IJulia.clear_output(true)
        pt = plot(x, t, label = :false, c="red", title=" $(round.(Int, i*dt)) [y] - Dyke intrusion", 
            xlabel="Distance [m]", ylabel="Temperature [°C]", ylims=(0, 1000), size = (800, 600))
        Plots.display(pt)
    end
    
end

Answer 1

你有一个语法错误！ println(s_time)

Answer 2

这永远不会 return true 因为数组不能等于数字:

if round.(Int, t[x.==s]) == 100

看看这个

julia> [100] == 100
false

不是很清楚你想要它的意思，但也许是这样

if all(round.(Int, t[x.==s]) .== 100)

或者 any?

如何从 Julia 中的 if 语句中打破 for 循环？

How to break a for loop from inside an if statement in Julia?

python

for-loop

if-statement

break

julia