Reenable ABC Model

src/MetagraphOptimization.jl

@@ -110,8 +110,6 @@ import Base.delete!
 import Base.insert!
 import Base.collect
 
-include("QEDprocesses_patch.jl")
-
 include("devices/interface.jl")
 include("task/type.jl")
 include("node/type.jl")
@@ -171,7 +169,6 @@ include("models/print.jl")
 
 include("models/physics_models/interface.jl")
 
-#=
 include("models/physics_models/abc/types.jl")
 include("models/physics_models/abc/particle.jl")
 include("models/physics_models/abc/compute.jl")
@@ -179,8 +176,8 @@ include("models/physics_models/abc/create.jl")
 include("models/physics_models/abc/properties.jl")
 include("models/physics_models/abc/parse.jl")
 include("models/physics_models/abc/print.jl")
-=#
 
+include("models/physics_models/qed/utility.jl")
 include("models/physics_models/qed/types.jl")
 include("models/physics_models/qed/particle.jl")
 include("models/physics_models/qed/diagrams.jl")

src/QEDprocesses_patch.jl

@@ -1,12 +0,0 @@
-# patch QEDprocesses
-# see issue https://github.com/QEDjl-project/QEDprocesses.jl/issues/77
-@inline function QEDprocesses.number_particles(
-    proc_def::QEDbase.AbstractProcessDefinition, ::Type{PS}
-) where {
-    DIR<:QEDbase.ParticleDirection,
-    PT<:QEDbase.AbstractParticleType,
-    EL<:AbstractFourMomentum,
-    PS<:ParticleStateful{DIR,PT,EL},
-}
-    return QEDprocesses.number_particles(proc_def, DIR(), PT())
-end

src/models/physics_models/abc/compute.jl

@@ -85,12 +85,6 @@ Linearly many FLOP with growing data.
 """
 function compute(::ComputeTaskABC_Sum, data...)::Float64
     return sum_kbn([data...])
-
-    #=s = 0.0im
-    for d in data
-        s += d
-    end
-    return s=#
 end
 
 function compute(::ComputeTaskABC_Sum, data::AbstractArray)::Float64

src/models/physics_models/qed/utility.jl

@@ -0,0 +1,14 @@
+using QEDprocesses
+
+# add type overload for number_particles function
+@inline function QEDprocesses.number_particles(
+    proc_def::QEDbase.AbstractProcessDefinition,
+    ::Type{PS},
+) where {
+    DIR <: QEDbase.ParticleDirection,
+    PT <: QEDbase.AbstractParticleType,
+    EL <: AbstractFourMomentum,
+    PS <: ParticleStateful{DIR, PT, EL},
+}
+    return QEDprocesses.number_particles(proc_def, DIR(), PT())
+end

src/node/validate.jl

@@ -29,12 +29,6 @@ function is_valid_node(graph::DAG, node::Node)
         @assert is_valid(graph, node.nodeSplit)
     end=#
 
-    # every child must be in some input of the task
-    for child in node.children
-        str = Symbol(to_var_name(child.id))
-        @assert (str in task(node).t1_inputs) || (str in task(node).t2_inputs) "$str was not in any of the tasks' inputs\nt1_inputs: $(task(node).t1_inputs)\nt2_inputs: $(task(node).t2_inputs)"
-    end
-
     return true
 end