Moved netSigma out of the sink strength (#100).

ORNL-Fusion · Mar 10, 2023 · b679d44 · b679d44
1 parent aba660e
commit b679d44
Showing 1 changed file with 7 additions and 2 deletions.
diff --git a/xolotl/core/include/xolotl/core/network/impl/FeCrReaction.tpp b/xolotl/core/include/xolotl/core/network/impl/FeCrReaction.tpp
@@ -386,6 +386,12 @@ FeCrSinkReaction::computeRate(IndexType gridIndex, double time)
 
 	double strength = this->asDerived()->getSinkStrength() * dc;
 
+	auto clReg = cl.getRegion();
+	Composition comp = clReg.getOrigin();
+	if (comp.isOnAxis(Species::Free)) {
+		strength *= this->_clusterData->extraData.netSigma(this->_reactant);
+	}
+
 	return strength;
 }
 
@@ -440,8 +446,7 @@ FeCrSinkReaction::getSinkStrength()
 		Z = rCoal;
 
 		Z *= 2.0 * ::xolotl::core::fecrDisloAlignment *
-			((1.0 - portion) * this->getSinkBias() + portion) *
-			this->_clusterData->extraData.netSigma(this->_reactant);
+			((1.0 - portion) * this->getSinkBias() + portion);
 	}
 
 	return density * Z;