L2 errors calculation

Author: cval26

rom/laghos.cpp

@@ -2129,28 +2129,28 @@ int main(int argc, char *argv[])
 						cout << "Window " << romOptions.window << ": initial romS norm " << romS.Norml2() << endl;
 					}
 
-					// Recompute and print the energy information after the
-					// change of basis, before any more timestepping takes
-					// place.
-					if (romOptions.hyperreductionSamplingType == eqp_energy)
-					{
-                        timeLoopTimer.Stop();
-						basis[romOptions.window]->LiftROMtoFOM(romS, *S);
-
-						double energy_total, energy_diff;
-						energy_total = oper->InternalEnergy(*e_gf) +
-							oper->KineticEnergy(*v_gf);
-						energy_diff	= energy_total - energy_init;
-
-						if (mpi.Root())
-						{
-							cout << "\tE_tot = " << scientific << setprecision(5)
-								<< energy_total
-								<< ",\tE_diff = " << scientific << setprecision(5)
-								<< energy_diff << endl; 
-						}
-                        timeLoopTimer.Start();
-					}
+					//if (romOptions.hyperreductionSamplingType == eqp_energy)
+					//{
+                    //    // Recompute and print the energy information after the
+                    //    // change of basis, before any more timestepping takes
+                    //    // place.
+                    //    timeLoopTimer.Stop();
+					//	basis[romOptions.window]->LiftROMtoFOM(romS, *S);
+
+					//	double energy_total, energy_diff;
+					//	energy_total = oper->InternalEnergy(*e_gf) +
+					//		oper->KineticEnergy(*v_gf);
+					//	energy_diff	= energy_total - energy_init;
+
+					//	if (mpi.Root())
+					//	{
+					//		cout << "\tE_tot = " << scientific << setprecision(5)
+					//			<< energy_total
+					//			<< ",\tE_diff = " << scientific << setprecision(5)
+					//			<< energy_diff << endl; 
+					//	}
+                    //    timeLoopTimer.Start();
+					//}
 
                     delete romOper[romOptions.window-1];
 
@@ -2297,6 +2297,35 @@ int main(int argc, char *argv[])
         outfile_tw_steps.close();
     }
 
+    if (rom_online && romOptions.hyperreductionSamplingType == eqp_energy)
+    {
+        if (myid == 0) cout << endl;
+        if (myid == 0) cout << "Final energy calculation." << endl;
+       
+        double energy_total, energy_diff;
+        energy_total = oper->InternalEnergy(*e_gf) +
+            oper->KineticEnergy(*v_gf);
+        energy_diff	= energy_total - energy_init;
+        
+        if (mpi.Root())
+        {
+            cout << "\tE_tot = " << scientific << setprecision(5)
+                << energy_total
+                << ",\tE_diff = " << scientific << setprecision(5)
+                << energy_diff << endl; 
+        }
+            
+        if (myid == 0) cout << endl;
+        if (myid == 0) cout << "Solution errors calculation." << endl;
+       
+        PrintDiffParGridFunction(normtype, myid, testing_parameter_outputPath
+                + "/Sol_Position" + romOptions.basisIdentifier, x_gf);
+        PrintDiffParGridFunction(normtype, myid, testing_parameter_outputPath
+                + "/Sol_Velocity" + romOptions.basisIdentifier, v_gf);
+        PrintDiffParGridFunction(normtype, myid, testing_parameter_outputPath
+                + "/Sol_Energy" + romOptions.basisIdentifier, e_gf);
+    }
+
     if (romOptions.use_sample_mesh)
     {
         if (romOptions.GramSchmidt && !spaceTime)