Solucionado?

parent da417964
...@@ -2,7 +2,7 @@ using LatticeGPU ...@@ -2,7 +2,7 @@ using LatticeGPU
using CUDA using CUDA
using TimerOutputs using TimerOutputs
@timeit "fA_fP test" let @timeit "fA_fP test" begin
function bndpropagator(U, dpar::DiracParam{T}, dws::DiracWorkspace,lp::SpaceParm{4,6,1,D}, maxiter::Int64, tol::Float64, c::Int64, s::Int64) where {T,D} function bndpropagator(U, dpar::DiracParam{T}, dws::DiracWorkspace,lp::SpaceParm{4,6,1,D}, maxiter::Int64, tol::Float64, c::Int64, s::Int64) where {T,D}
...@@ -45,7 +45,7 @@ end ...@@ -45,7 +45,7 @@ end
function fP_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1.0e-16) function fP_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1.0e-16)
@timeit "fP inversion (x12)" let @timeit "fP inversion (x12)" begin
lp = SpaceParm{4}(size,(4,4,4,4),1,(0,0,0,0,0,0)); lp = SpaceParm{4}(size,(4,4,4,4),1,(0,0,0,0,0,0));
exptheta = exp.(im.*theta./lp.iL); exptheta = exp.(im.*theta./lp.iL);
...@@ -74,7 +74,7 @@ end end ...@@ -74,7 +74,7 @@ end end
end end
@timeit "fP analitical solution" let @timeit "fP analitical solution" begin
#THEORETICAL SOLUTION: hep-lat/9606016 eq (2.33) #THEORETICAL SOLUTION: hep-lat/9606016 eq (2.33)
...@@ -99,7 +99,7 @@ end ...@@ -99,7 +99,7 @@ end
function fA_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1.0e-16) function fA_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1.0e-16)
@timeit "fA inversion (x12)" let @timeit "fA inversion (x12)" begin
lp = SpaceParm{4}(size,(4,4,4,4),1,(0,0,0,0,0,0)); lp = SpaceParm{4}(size,(4,4,4,4),1,(0,0,0,0,0,0));
exptheta = exp.(im.*theta./lp.iL); exptheta = exp.(im.*theta./lp.iL);
...@@ -129,7 +129,7 @@ function fA_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1 ...@@ -129,7 +129,7 @@ function fA_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1
end end
#THEORETICAL SOLUTION: hep-lat/9606016 eq (2.32) #THEORETICAL SOLUTION: hep-lat/9606016 eq (2.32)
@timeit "fA analitical solution" let @timeit "fA analitical solution" begin
res_th = zeros(lp.iL[4]) res_th = zeros(lp.iL[4])
pp3 = ntuple(i -> theta[i]/lp.iL[i],3) pp3 = ntuple(i -> theta[i]/lp.iL[i],3)
......
...@@ -65,8 +65,6 @@ end ...@@ -65,8 +65,6 @@ end
end end
end
#compute Sum{x} D^-1(x|y)P(y) #compute Sum{x} D^-1(x|y)P(y)
@timeit "Solving propagator" begin @timeit "Solving propagator" begin
...@@ -91,11 +89,13 @@ if dif > 1.0e-15 ...@@ -91,11 +89,13 @@ if dif > 1.0e-15
error("Dwpl test for s=",s,", c=",c," failed with difference: ",dif,"\n") error("Dwpl test for s=",s,", c=",c," failed with difference: ",dif,"\n")
end end
end
return dif return dif
end end
begin begin
dif = 0.0 dif = 0.0
for i in 1:3 for j in 1:4 for i in 1:3 for j in 1:4
......
...@@ -2,7 +2,7 @@ using CUDA,LatticeGPU, TimerOutputs ...@@ -2,7 +2,7 @@ using CUDA,LatticeGPU, TimerOutputs
#Check that Dw ( (DwdagDw)^{-1} g5 Dw g5 ) psi = psi for random fields #Check that Dw ( (DwdagDw)^{-1} g5 Dw g5 ) psi = psi for random fields
@timeit "Rand solver test" let @timeit "Rand solver test" begin
function pfrandomize!(f,lp::SpaceParm,dpar::DiracParam) function pfrandomize!(f,lp::SpaceParm,dpar::DiracParam)
...@@ -32,7 +32,7 @@ function krnl_assign_pf!(f::AbstractArray{T}, p ,lp::SpaceParm) where {T} ...@@ -32,7 +32,7 @@ function krnl_assign_pf!(f::AbstractArray{T}, p ,lp::SpaceParm) where {T}
end end
@timeit "Generate random fields" let @timeit "Generate random fields" begin
lp = SpaceParm{4}((16,16,16,16), (4,4,4,4), 0, (0,0,0,0,0,0)) lp = SpaceParm{4}((16,16,16,16), (4,4,4,4), 0, (0,0,0,0,0,0))
gp = GaugeParm{Float64}(SU3{Float64}, 6.0, 1.0) gp = GaugeParm{Float64}(SU3{Float64}, 6.0, 1.0)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment