From dbf2cb3e1a87152d01b07a88f58e570e315d3720 Mon Sep 17 00:00:00 2001
From: Shikhar Kumar <shikhark@mit.edu>
Date: Sat, 10 Nov 2018 13:58:23 -0500
Subject: [PATCH] Minor fixes to 2group solver

---
 src/cmfd_solver.cpp | 6 +-----
 1 file changed, 1 insertion(+), 5 deletions(-)

diff --git a/src/cmfd_solver.cpp b/src/cmfd_solver.cpp
index 07fdd3c3a..1a2d87091 100644
--- a/src/cmfd_solver.cpp
+++ b/src/cmfd_solver.cpp
@@ -1,5 +1,3 @@
-//TODO remove
-#include <iostream>
 #include <vector>
 
 #include "openmc/cmfd_solver.h"
@@ -12,8 +10,6 @@ namespace openmc {
 // Global variables
 //==============================================================================
 
-// TODO check which variables actually necessary
-
 std::vector<int> indptr;
 
 std::vector<int> indices;
@@ -179,7 +175,7 @@ int cmfd_linsolver_2g(double* A_data, double* b, double* x, double tol) {
         for (int icol = d1idx + 2; icol < indptr[irow + 1]; icol++)
           tmp1 += A_data[icol] * x[indices[icol]];
         for (int icol = d2idx + 1; icol < indptr[irow + 2]; icol++)
-          tmp1 += A_data[icol] * x[indices[icol]];
+          tmp2 += A_data[icol] * x[indices[icol]];
 
         // Adjust with RHS vector
         tmp1 = b[irow] - tmp1;