Fix 1D models analytic reference tests

test/ColumnTests.cpp

@@ -360,7 +360,7 @@ namespace column
 
 			// Get data from simulation
 			cadet::InternalStorageUnitOpRecorder const* const simData = drv.solution()->unitOperation(0);
-			double const* outlet = (forwardFlow ? simData->outlet() : simData->inlet());
+			double const* outlet = simData->outlet();
 
 			// Compare
 			for (unsigned int i = 0; i < simData->numDataPoints() * simData->numComponents() * simData->numInletPorts(); ++i, ++outlet)
@@ -398,7 +398,7 @@ namespace column
 
 			// Get data from simulation
 			cadet::InternalStorageUnitOpRecorder const* const simData = drv.solution()->unitOperation(0);
-			double const* outlet = (forwardFlow ? simData->outlet() : simData->inlet());
+			double const* outlet = simData->outlet();
 
 			// Compare
 			for (unsigned int i = 0; i < simData->numDataPoints() * simData->numComponents() * simData->numInletPorts(); ++i, ++outlet)
@@ -994,20 +994,20 @@ namespace column
 					const unsigned int nDataPoints = fwdData->numDataPoints() * nComp * nPorts;
 					double const* const time = drvFwd.solution()->time();
 					double const* fwdOutlet = fwdData->sensOutlet(0);
-					double const* bwdInlet = bwdData->sensInlet(0);
+					double const* bwdOutlet = bwdData->sensOutlet(0);
 
 					// Compare
 					unsigned int numPassed = 0;
-					for (unsigned int i = 0; i < nDataPoints; ++i, ++fwdOutlet, ++bwdInlet)
+					for (unsigned int i = 0; i < nDataPoints; ++i, ++fwdOutlet, ++bwdOutlet)
 					{
 						const unsigned int comp = (i % (nComp * nPorts)) % nComp;
 						const unsigned int port = (i % (nComp * nPorts)) / nComp;
 						const unsigned int timeIdx = i / (nComp * nPorts);
 
 						INFO("Time " << time[timeIdx] << " Port " << port << " Component " << comp << " time point idx " << timeIdx);
-						CHECK(*bwdInlet == makeApprox(*fwdOutlet, relTol, absTol));
+						CHECK(*bwdOutlet == makeApprox(*fwdOutlet, relTol, absTol));
 
-						const bool relativeOK = std::abs(*bwdInlet - *fwdOutlet) <= relTol * std::abs(*fwdOutlet);
+						const bool relativeOK = std::abs(*bwdOutlet - *fwdOutlet) <= relTol * std::abs(*fwdOutlet);
 						if (relativeOK)
 							++numPassed;
 					}