bug fixes

dominikkiese · Dec 13, 2023 · b96aa86 · b96aa86
1 parent ad5bb49
commit b96aa86
Show file tree

Hide file tree

Showing 13 changed files with 190 additions and 240 deletions.
diff --git a/src/Action/Action.jl b/src/Action/Action.jl
@@ -123,7 +123,7 @@ function get_Γ_comp_avx!(
     )           :: Nothing
 
     # init with bare value
-    @turbo temp .= a.Γ[comp].bare
+    temp .= a.Γ[comp].bare
 
     # add s channel
     if ch_s
@@ -184,7 +184,7 @@ function replace_with!(
     )  :: Nothing
 
     # replace self energy
-    @turbo a1.Σ .= a2.Σ
+    a1.Σ .= a2.Σ
 
     # replace vertices
     for i in eachindex(a1.Γ)
@@ -215,7 +215,7 @@ function mult_with!(
     )   :: Nothing
 
     # multiply self energy
-    @turbo a.Σ .*= fac
+    a.Σ .*= fac
 
     # multiply vertices
     for i in eachindex(a.Γ)
@@ -267,7 +267,7 @@ function mult_with_add_to!(
     )   :: Nothing
 
     # multiply add for the self energy
-    @turbo a1.Σ .+= fac .* a2.Σ
+    a1.Σ .+= fac .* a2.Σ
 
     # multiply add for the vertices
     for i in eachindex(a1.Γ)

diff --git a/src/Action/channel.jl b/src/Action/channel.jl
@@ -359,10 +359,10 @@ function replace_with!(
     ch2 :: Channel
     )   :: Nothing 
 
-    @turbo ch1.q1   .= ch2.q1 
-    @turbo ch1.q2_1 .= ch2.q2_1
-    @turbo ch1.q2_2 .= ch2.q2_2
-    @turbo ch1.q3   .= ch2.q3
+    ch1.q1   .= ch2.q1 
+    ch1.q2_1 .= ch2.q2_1
+    ch1.q2_2 .= ch2.q2_2
+    ch1.q3   .= ch2.q3
 
     return nothing 
 end
@@ -373,10 +373,10 @@ function mult_with!(
     fac :: Float64
     )   :: Nothing 
 
-    @turbo ch.q1   .*= fac
-    @turbo ch.q2_1 .*= fac
-    @turbo ch.q2_2 .*= fac 
-    @turbo ch.q3   .*= fac 
+    ch.q1   .*= fac
+    ch.q2_1 .*= fac
+    ch.q2_2 .*= fac 
+    ch.q3   .*= fac 
 
     return nothing 
 end 
@@ -388,10 +388,10 @@ function mult_with_add_to!(
     ch1 :: Channel
     )   :: Nothing 
 
-    @turbo ch1.q1   .+= fac .* ch2.q1 
-    @turbo ch1.q2_1 .+= fac .* ch2.q2_1
-    @turbo ch1.q2_2 .+= fac .* ch2.q2_2
-    @turbo ch1.q3   .+= fac .* ch2.q3
+    ch1.q1   .+= fac .* ch2.q1 
+    ch1.q2_1 .+= fac .* ch2.q2_1
+    ch1.q2_2 .+= fac .* ch2.q2_2
+    ch1.q3   .+= fac .* ch2.q3
 
     return nothing 
 end

diff --git a/src/Flow/BSE_s.jl b/src/Flow/BSE_s.jl
@@ -40,13 +40,13 @@ function compute_channel_s_BSE!(
     # parse result
     for i in eachindex(a2.Γ)
         if kernel == 1
-            @turbo a2.Γ[i].ch_s.q1[:, w1] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q1[:, w1] .= view(tbuff[1], i, :)
         elseif kernel == 2
-            @turbo a2.Γ[i].ch_s.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
         elseif kernel == 3
-            @turbo a2.Γ[i].ch_s.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
         else
-            @turbo a2.Γ[i].ch_s.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
         end
     end
 

diff --git a/src/Flow/BSE_t.jl b/src/Flow/BSE_t.jl
@@ -40,13 +40,13 @@ function compute_channel_t_BSE!(
     # parse result
     for i in eachindex(a2.Γ)
         if kernel == 1
-            @turbo a2.Γ[i].ch_t.q1[:, w1] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_t.q1[:, w1] .= view(tbuff[1], i, :)
         elseif kernel == 2
-            @turbo a2.Γ[i].ch_t.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_t.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
         elseif kernel == 3
-            @turbo a2.Γ[i].ch_t.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_t.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
         else
-            @turbo a2.Γ[i].ch_t.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_t.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
         end
     end
 

diff --git a/src/Flow/BSE_u.jl b/src/Flow/BSE_u.jl
@@ -40,13 +40,13 @@ function compute_channel_u_BSE!(
     # parse result
     for i in eachindex(a2.Γ)
         if kernel == 1
-            @turbo a2.Γ[i].ch_u.q1[:, w1] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_u.q1[:, w1] .= view(tbuff[1], i, :)
         elseif kernel == 2
-            @turbo a2.Γ[i].ch_u.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_u.q2_1[:, w1, w2] .= view(tbuff[1], i, :)
         elseif kernel == 3
-            @turbo a2.Γ[i].ch_u.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_u.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
         else
-            @turbo a2.Γ[i].ch_u.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_u.q3[:, w1, w2, w3] .= view(tbuff[1], i, :)
         end
     end
 

diff --git a/src/Flow/SDE.jl b/src/Flow/SDE.jl
@@ -45,9 +45,9 @@ function compute_channel_s_reduced!(
     # parse result 
     for i in eachindex(a2.Γ)
         if kernel == 1
-            @turbo a2.Γ[i].ch_s.q1[:, w1] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q1[:, w1] .= view(tbuff[1], i, :)
         else
-            @turbo a2.Γ[i].ch_s.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
+            a2.Γ[i].ch_s.q2_2[:, w1, w3] .= view(tbuff[1], i, :)
         end
     end
 

diff --git a/src/Flow/bubbles.jl b/src/Flow/bubbles.jl
@@ -106,7 +106,7 @@ function compute_corrs_kat!(
     )     :: Nothing 
 
     # reset buffer
-    @turbo corrs .= 0.0
+    corrs .= 0.0
 
     # compute boundary corrections
     @sync for i in 1 : size(corrs, 3)
@@ -168,7 +168,7 @@ function compute_corrs!(
     )     :: Nothing 
 
     # reset buffer
-    @turbo corrs .= 0.0
+    corrs .= 0.0
 
     # compute boundary corrections
     @sync for i in 1 : size(corrs, 3)