जूलिया - जटिल सरणी

Question

का असली हिस्सा प्राप्त करें मुझे वास्तविक और जटिल भागों के लिए परिवर्तनीय z::Array{Complex128,1} को दो सरणी में विभाजित करने की आवश्यकता है। एक तरीका यह है नया चर ::Array{Float64,1} बनाने के लिए और तत्व द्वारा उन्हें तत्व को भरने के लिए है:

for i = 1:size(z)[1] 
    ri[i] = z[i].re 
    ii[i] = z[i].im 
end 

वहाँ यह है कि डेटा की प्रतिलिपि, किसी भी तरह प्रगति और z की ऑफसेट जोड़ तोड़ की तरह शामिल नहीं करता है करने के लिए एक तरीका है?

Answer 1

संपादित करें: मूल संस्करण एक अनावश्यक reshape ऑपरेशन का उपयोग करता था। जैसा कि @ डीएनएफ ने टिप्पणियों में बताया, वह अनावश्यक था। जवाब के बाद से संशोधित किया गया है।

मामले ऐसे नकल कोई मुद्दा नहीं है, सिर्फ करना real(z) और imag(z) (नाम बदला real.(z) और v0.6 में imag.(z))। मैं इसे भविष्य के पाठकों की मदद करने के लिए शामिल करता हूं जिनके पास समान समस्या है, लेकिन जो प्रतिलिपि बनाने की परवाह नहीं करते हैं।

जैसा कि आप सुझाव देते हैं, आप डेटा कॉपी करने से बचने के लिए z के चरणों में हेरफेर कर सकते हैं। सीधे शब्दों में

zfl = reinterpret(Float64, z) 
zre = @view zfl[1:2:end-1] 
zim = @view zfl[2:2:end] 

संयुक्त, हम देख कोई डेटा नकल (आवंटन ढेर आबंटित सरणी दृष्टिकोण के कारण उत्पन्न कर रहे हैं, और कम से कम कर रहे हैं) नहीं है।

julia> z = Vector{Complex128}(100000); 

julia> function reimvec(z) 
      zfl = reinterpret(Float64, z) 
      zre = @view zfl[1:2:end-1] 
      zim = @view zfl[2:2:end] 
      zre, zim 
     end 
reimvec (generic function with 1 method) 

julia> @time reimvec(z); 
    0.000005 seconds (9 allocations: 400 bytes) 

हम देख सकते हैं, पर्दे के पीछे, इस तरह के एक सरणी strided है:

julia> strides(reimvec(z)[1]) 
(2,)