如何在Julia中初始化结构数组

我想知道我应该怎么放的,而不是xxxxx有1000个项目阵列结构初始化,其中a从去20003000(即数组手段的指标1a2000,数组手段的指标2a2001,等等),并b始终零。

struct MyStruct
    a
    b
end

myArray = Vector{MyStruct}( xxxxx , 1000)

我知道我可以进行循环并单独分配值,我只是想知道Julia中是否有更快的方法。

雅各布·尼森

只需在单个帖子中收集其他用户的答案和评论即可:

不,没有比在Julia中循环更快的了

与其他脚本语言(如Python和R)不同,Julia中的循环速度很快。实际上,其他“矢量化”操作(例如广播)是根据Julia循环本身实现的。因此,一种快速的解决方案可能是:

function initialize_vector(range::AbstractRange)
    v = Vector{MyStruct}(undef, length(range))
    @inbounds for i in eachindex(range)
        v[i] = MyStruct(range[i], 0)
    end
    return v
end

广播既方便又方便

广播几乎或有时与循环一样快,并且通常更简洁和方便。在这种情况下,initialize_vector上面的函数可以写成:

initialize_vector(range::AbstractRange) = MyStruct.(range, 0)

基准测试表明,这两个功能的速度几乎相同。

记住在结构中键入字段以获得更快的代码

Julia依靠对类型的准确推断来创建快速,专业的代码。如果类型的MyStruct.aMyStruct.b可以是任何东西,它通常是无法推断到底是什么样的操作应在执行MyStruct即使在这种情况下,编译器也可以推断类型为Int,每个类型都MyStruct必须包含对堆分配的Ints的引用,而不是堆栈分配的。因此,只需更改即可获得10倍的加速

struct MyStruct
    a
    b
end

struct MyStruct
    a::Int
    b::Int
end

如果您想要类型Mystruct.aMyStruct.b能够变化,可以创建一个参量MyStruct,如下所示:

struct MyStruct{T}
    a::T
    b::T
end

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章