fix the order of the conjugation when constructing the electron row outer product (#141)

ewanwm · web-flow · commit d859825af718 · 2025-08-22T17:05:00.000+02:00
Looking more carefully at eq 5 in Barger et al, the conjugation should be swapped when building the electon row outer product
diff --git a/nuTens/propagator/const-density-solver.cpp b/nuTens/propagator/const-density-solver.cpp
@@ -48,14 +48,14 @@ void ConstDensityMatterSolver::buildElectronOuterProduct()
     if (antiNeutrino)
     {
         electronOuter =
-            Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}).conj(), mixingMatrix.getValues({0, 0, "..."})),
+            Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}), mixingMatrix.getValues({0, 0, "..."}).conj()),
                           -nuTens::constants::Groot2 * density);
     }
 
     else 
     {
         electronOuter =
-            Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}), mixingMatrix.getValues({0, 0, "..."}).conj()),
+            Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}).conj(), mixingMatrix.getValues({0, 0, "..."})),
                           nuTens::constants::Groot2 * density);
     }
 

Original file line number	Diff line number	Diff line change
`@@ -48,14 +48,14 @@ void ConstDensityMatterSolver::buildElectronOuterProduct()`
`48`	`48`	`if (antiNeutrino)`
`49`	`49`	`{`
`50`	`50`	`electronOuter =`
`51`		`- Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}).conj(), mixingMatrix.getValues({0, 0, "..."})),`
	`51`	`+ Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}), mixingMatrix.getValues({0, 0, "..."}).conj()),`
`52`	`52`	`-nuTens::constants::Groot2 * density);`
`53`	`53`	`}`
`54`	`54`
`55`	`55`	`else`
`56`	`56`	`{`
`57`	`57`	`electronOuter =`
`58`		`- Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}), mixingMatrix.getValues({0, 0, "..."}).conj()),`
	`58`	`+ Tensor::scale(Tensor::outer(mixingMatrix.getValues({0, 0, "..."}).conj(), mixingMatrix.getValues({0, 0, "..."})),`
`59`	`59`	`nuTens::constants::Groot2 * density);`
`60`	`60`	`}`
`61`	`61`