sun radius dA fix?

palumbom · Nov 29, 2024 · b052c9c · b052c9c
1 parent 4a5b01a
commit b052c9c
Show file tree

Hide file tree

Showing 4 changed files with 4 additions and 203 deletions.
diff --git a/scripts/debug_gpu.jl b/scripts/debug_gpu.jl
diff --git a/scripts/debug_gpu_eclipse.jl b/scripts/debug_gpu_eclipse.jl
@@ -75,5 +75,6 @@ plt.imshow(wsp.dA .- Array(gpu_allocs.dA))
 plt.title("dA")
 plt.colorbar()
 plt.savefig("gpu_test.png")
+plt.show()
 plt.clf()
 #z_rot & dA
diff --git a/src/eclipse_comp.jl b/src/eclipse_comp.jl
@@ -129,7 +129,8 @@ function eclipse_compute_quantities!(disk::DiskParamsEclipse{T}, epoch::T, t::In
             # calculate area element of tile
             dϕ = step(ϕe_sub) 
             dθ = step(θe_sub) 
-            dA_sub = map(x -> calc_dA(sun_radius, getindex(x,1), dϕ, dθ), subgrid)
+            dA_sub = map(x -> calc_dA(1.0, getindex(x,1), dϕ, dθ), subgrid)
+
             # get total projected, visible area of larger tile
             dA_total_proj = dA_sub .* mu_grid
             dA_total_proj_matrix[i,j] = dA_total_proj

diff --git a/src/gpu/gpu_precomps_eclipse.jl b/src/gpu/gpu_precomps_eclipse.jl
@@ -244,7 +244,7 @@ function calc_eclipse_quantities_gpu!(wavelength, μs, z_rot, ax_codes,
                 z_rot_sub += (v_orbit_sub / c_ms)
 
                 # get projected area element
-                dA_sub = calc_dA_gpu(sun_radius, ϕc_sub, dϕ, dθ)
+                dA_sub = calc_dA_gpu(1.0, ϕc_sub, dϕ, dθ)
                 dA_sub *= μ_sub
                 dA_sum += dA_sub