review #1 changes

Author: cval26

rom/laghos_rom.cpp

@@ -1117,15 +1117,15 @@ void ROM_Sampler::Finalize(Array<int> &cutoff, ROM_Options& input)
 		CAROM::Matrix *tBasisE = basisE->getFirstNColumns(cutoff[2]);
 
 		// Form the energy identity and replace the last basis vector by it. 
-		Vector unitE(input.tL2size);
+		Vector unitE(tL2size);
 
 		// TODO: my understanding is that the "=" opetaror is overloaded, so
 		// that, by setting unitE = 1.0, the vector unitE contains the coeffs
 		// that result in the summed function being 1.0 on the grid points. 
 		// Is that correct?
 		unitE = 1.0;
 
-		for (int i = 0; i < input.tL2size; i++)
+		for (int i = 0; i < tL2size; i++)
 		{
 			// TODO: I believe I can access the last column in this way,
 			// without having to use pointer arithmetic, right?
@@ -4899,3 +4899,8 @@ void ROM_Operator::StepRK2Avg(Vector &S, double &t, double &dt) const
 
     t += dt;
 }
+
+HyperreductionSamplingType ROM_Operator::getSamplingType() const
+{
+	return hyperreductionSamplingType;
+}

rom/laghos_rom.hpp

@@ -1235,6 +1235,8 @@ class ROM_Operator : public TimeDependentOperator
 	void ForceIntegratorEQP_E(Vector const& v, Vector & res,
 			bool energy_conserve = false) const;
 
+	HyperreductionSamplingType getSamplingType() const;
+
     ~ROM_Operator()
     {
         delete mat_gf_coeff;

rom/laghos_solver.cpp

@@ -279,7 +279,7 @@ void LagrangianHydroOperator::SolveVelocity(const Vector &S,
         {
 			// TODO: this function doesn't have access to the hyperreduction 
 			// sampling type; how should I resolve this?
-            rom_op->ForceIntegratorEQP_FOM(rhs, hyperreduction.SamplingType == eqp_energy);
+            rom_op->ForceIntegratorEQP_FOM(rhs, rom_op->getSamplingType() == eqp_energy);
         }
         else
         {
@@ -412,7 +412,7 @@ void LagrangianHydroOperator::SolveEnergy(const Vector &S, const Vector &v,
         {
 			// TODO: this function doesn't have access to the hyperreduction 
 			// sampling type; how should I resolve this?
-            rom_op->ForceIntegratorEQP_E_FOM(v, e_rhs, hyperreduction.SamplingType == eqp_energy);
+            rom_op->ForceIntegratorEQP_E_FOM(v, e_rhs, rom_op->getSamplingType() == eqp_energy);
         }
         else
         {