我正在研究生物系统中的进化模拟.我必须求解多项式方程,找到根(u*X ^ 3 - N*p*r*X ^ 2 - N*p*X ^ 2 + K ^ 2*u*X - N*K ^ 2*p ),其中u和K是常数,N是常数阵列,p,r是演化参数.基本上,对于每一代人口中的每个人,我需要进行以下计算(长度(N)>>长度(p)):
for i = 1:length(p)
for j = 1:length(N)
X[j,i] = mean(fzeros(S -> u*S^3 - p[i]*N[j]*r[i]*S^2 - p[i]*N[j]*S^2 + K^2*u*S - p[i]*K^2*N[j], 0, Inf) )
end
end
Run Code Online (Sandbox Code Playgroud)
我知道我可以通过使用不同的核心为不同的个体求解方程来并行代码,甚至在每个个体中我也可以并行求解每个X [j,i].我想知道处理这种情况的最佳做法/快速方法是什么?是否有可能以更快的方式解决单个方程?
我最近在Julia转移了我的代码.我想知道如何在命令行中执行Julia代码?
我知道Julia代码可以通过运行一次来编译.
但问题是我需要为集群上的模拟模型进行参数扫描,在那里我只能使用命令行 - 而不是REPL.
在群集上运行模拟复制的最佳实践是什么?