= update Project.toml

Author: rvignolo

Project.toml

@@ -1,7 +1,7 @@
 name = "NeutronTransport"
 uuid = "161f4978-d3bb-4696-90f4-058341551711"
 authors = ["Ramiro Vignolo <ramirovignolo@gmail.com>"]
-version = "0.1.0"
+version = "0.2.0"
 
 [deps]
 Gridap = "56d4f2e9-7ea1-5844-9cf6-b9c51ca7ce8e"

src/NeutronTransport.jl

@@ -10,8 +10,6 @@ using Gridap.Geometry: get_face_tag, get_tag_from_name
 
 import Base: show, eltype
 
-include("cross_section.jl")
-
 abstract type TransportFormulation end
 
 struct MethodOfCharacteristics <: TransportFormulation end
@@ -23,6 +21,7 @@ const MoC = MethodOfCharacteristics
 const CP = CollisionProbability
 const SN = DiscreteOrdinates
 
+include("cross_section.jl")
 include("polar_quad.jl")
 include("quadrature.jl")
 

src/mocsolver.jl

@@ -17,7 +17,7 @@ struct MoCSolution{T<:Real,P<:MoCProblem} <: TransportSolution
     # reduced source
     q::Vector{T}
 
-    # total source (integrated in energy given per region) is used for convergence purposes
+    # total source (integrated in energy given per region) used for convergence purposes
     Q::Vector{T}
 
     # angular flux at boundary
@@ -375,9 +375,7 @@ function residual(sol::MoCSolution{T}, prob::MoCProblem) where {T}
         xs = xss[idx]
 
         old_qi = Q[i]
-
         new_qi = zero(T)
-        # old_qi = zero(T)
 
         @unpack νΣf, Σs0 = xs
         fissionable = isfissionable(xs)
@@ -393,7 +391,6 @@ function residual(sol::MoCSolution{T}, prob::MoCProblem) where {T}
         end
 
         new_qi /= keff
-        # old_qi /= keff
 
         # total scattering source
         for g in 1:NGroups, g′ in 1:NGroups
@@ -402,7 +399,6 @@ function residual(sol::MoCSolution{T}, prob::MoCProblem) where {T}
             new_qi += Σsgg′ * φ[ig′]
         end
 
-
         if old_qi > 0
             ϵ += ((new_qi - old_qi) / old_qi)^2
         end