#include <FilterModalEnergy.h>

Inheritance diagram for Nektar::SolverUtils::FilterModalEnergy:

Collaboration diagram for Nektar::SolverUtils::FilterModalEnergy:

Public Member Functions
SOLVER_UTILS_EXPORT	FilterModalEnergy (const LibUtilities::SessionReaderSharedPtr &pSession, const ParamMap &pParams)

virtual SOLVER_UTILS_EXPORT	~FilterModalEnergy ()

Public Member Functions inherited from Nektar::SolverUtils::Filter
SOLVER_UTILS_EXPORT	Filter (const LibUtilities::SessionReaderSharedPtr &pSession)

virtual SOLVER_UTILS_EXPORT	~Filter ()

SOLVER_UTILS_EXPORT void	Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

SOLVER_UTILS_EXPORT void	Update (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

SOLVER_UTILS_EXPORT void	Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

SOLVER_UTILS_EXPORT bool	IsTimeDependent ()

Static Public Member Functions
static FilterSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const ParamMap &pParams)

Static Public Attributes
static std::string	className

Protected Member Functions
virtual void	v_Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

virtual void	v_Update (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

virtual void	v_Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)

virtual bool	v_IsTimeDependent ()

NekDouble	L2Error (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, unsigned int field, const NekDouble &time)

void	SetUpBaseFields (SpatialDomains::MeshGraphSharedPtr &mesh)

void	ImportFldBase (std::string pInfile, SpatialDomains::MeshGraphSharedPtr pGraph)

Private Attributes
enum MultiRegions::ProjectionType	m_projectionType

Array< OneD, MultiRegions::ExpListSharedPtr >	m_base

LibUtilities::FieldIOSharedPtr	m_fld

vector< unsigned int >	m_boundaryRegionsIdList

vector< bool >	m_boundaryRegionIsInList

unsigned int	m_index

unsigned int	m_outputFrequency

unsigned int	m_outputPlane

bool	m_isHomogeneous1D

bool	m_isHomogeneous2D

bool	m_PertEnergy

int	m_npointsZ

int	m_nproc

std::string	m_outputFile

std::string	m_EqTypeStr

std::ofstream	m_outputStream

LibUtilities::BasisSharedPtr	m_homogeneousBasis

std::string	m_BoundaryString

int	m_nplanes

int	m_NumQuadPointsError

bool	m_SingleMode

bool	m_HalfMode

bool	m_MultipleModes

bool	m_useFFT

NekDouble	m_LhomZ

bool	m_homogen_dealiasing

Friends
class	MemoryManager< FilterModalEnergy >

Additional Inherited Members
Public Types inherited from Nektar::SolverUtils::Filter
typedef std::map< std::string, std::string >	ParamMap

Protected Attributes inherited from Nektar::SolverUtils::Filter
LibUtilities::SessionReaderSharedPtr	m_session

Detailed Description

Definition at line 57 of file FilterModalEnergy.h.

Constructor & Destructor Documentation

Nektar::SolverUtils::FilterModalEnergy::FilterModalEnergy	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const ParamMap &	pParams
	)

Constructor.

Definition at line 51 of file FilterModalEnergy.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::LibUtilities::Equation::Evaluate(), m_EqTypeStr, m_fld, m_isHomogeneous1D, m_isHomogeneous2D, m_LhomZ, m_NumQuadPointsError, m_outputFile, m_outputFrequency, m_outputPlane, m_PertEnergy, and Nektar::SolverUtils::Filter::m_session.

                              :
     Filter(pSession)
 {
     ParamMap::const_iterator it;
 
     // OutputFile
     it = pParams.find("OutputFile");
     if (it == pParams.end())
     {
         m_outputFile = m_session->GetSessionName();
     }
     else
     {
         ASSERTL0(it->second.length() > 0, "Missing parameter 'OutputFile'.");
         m_outputFile = it->second;
     }
     if (!(m_outputFile.length() >= 4
           && m_outputFile.substr(m_outputFile.length() - 4) == ".mdl"))
     {
         m_outputFile += ".mdl";
     }
 
     // OutputFrequency
     it = pParams.find("OutputFrequency");
     if (it == pParams.end())
     {
         m_outputFrequency = 1;
     }
     else
     {
         LibUtilities::Equation equ(m_session, it->second);
         m_outputFrequency = floor(equ.Evaluate());
     }
 
 
     m_session->MatchSolverInfo("Homogeneous", "1D", m_isHomogeneous1D, false);
     m_session->MatchSolverInfo("Homogeneous", "2D", m_isHomogeneous2D, false);
     m_session->MatchSolverInfo("CalculatePerturbationEnergy", "True",
                                m_PertEnergy, false);
     m_session->LoadParameter  ("NumQuadPointsError", m_NumQuadPointsError, 0);
     m_EqTypeStr = m_session->GetSolverInfo("EQTYPE");
 
     // OutputPlane
     if(m_isHomogeneous1D || m_isHomogeneous2D)
     {
         m_session->LoadParameter("LZ", m_LhomZ);
 
         it = pParams.find("OutputPlane");
         if (it == pParams.end())
         {
             m_outputPlane = 0;
         }
         else
         {
             LibUtilities::Equation equ(m_session, it->second);
             m_outputPlane = floor(equ.Evaluate());
         }
     }
 
     m_fld = MemoryManager<LibUtilities::FieldIO>::
                 AllocateSharedPtr(pSession->GetComm());
 
 }

Nektar::SolverUtils::FilterModalEnergy::~FilterModalEnergy ( )

virtual

Destructor.

Definition at line 120 of file FilterModalEnergy.cpp.

 {
 
 }

Member Function Documentation

static FilterSharedPtr Nektar::SolverUtils::FilterModalEnergy::create	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const ParamMap &	pParams
	)

inlinestatic

Definition at line 63 of file FilterModalEnergy.h.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

     {
         FilterSharedPtr p = MemoryManager<FilterModalEnergy>::
                                 AllocateSharedPtr(pSession, pParams);
         //p->InitObject();
         return p;
     }

void Nektar::SolverUtils::FilterModalEnergy::ImportFldBase	(	std::string	pInfile,
		SpatialDomains::MeshGraphSharedPtr	pGraph
	)

protected

Import the base flow fld file.

Definition at line 624 of file FilterModalEnergy.cpp.

References ASSERTL0, m_base, m_fld, and Nektar::SolverUtils::Filter::m_session.

Referenced by v_Update().

 {
     std::vector<LibUtilities::FieldDefinitionsSharedPtr> FieldDef;
     std::vector<std::vector<NekDouble> > FieldData;
 
     // Get Homogeneous
     m_fld->Import(pInfile,FieldDef,FieldData);
 
     int nvar = m_session->GetVariables().size();
     if(m_session->DefinesSolverInfo("HOMOGENEOUS"))
     {
         std::string HomoStr = m_session->GetSolverInfo("HOMOGENEOUS");
     }
     // Copy FieldData into m_fields
     for (int j = 0; j < nvar; ++j)
     {
         for (int i = 0; i < FieldDef.size(); ++i)
         {
             bool flag =
                 FieldDef[i]->m_fields[j] == m_session->GetVariable(j);
 
             ASSERTL0(flag, (std::string("Order of ") + pInfile
                         + std::string(" data and that defined in "
                             "m_boundaryconditions differs")).c_str());
 
             m_base[j]->ExtractDataToCoeffs(FieldDef[i], FieldData[i],
                                            FieldDef[i]->m_fields[j],
                                            m_base[j]->UpdateCoeffs());
         }
     }
 }

NekDouble Nektar::SolverUtils::FilterModalEnergy::L2Error	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		unsigned int	field,
		const NekDouble &	time
	)

protected

Calculate the L2 norm of a given field for calculating the modal energy.

Definition at line 348 of file FilterModalEnergy.cpp.

References m_NumQuadPointsError.

Referenced by v_Update().

 {
     NekDouble L2error = -1.0;
     LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
 
     if (m_NumQuadPointsError == 0)
     {
         if (pFields[field]->GetPhysState() == false)
         {
             pFields[field]->BwdTrans(pFields[field]->GetCoeffs(),
                                      pFields[field]->UpdatePhys());
         }
     }
 
     L2error = pFields[field]->L2(pFields[field]->GetPhys());
     return L2error;
 }

void Nektar::SolverUtils::FilterModalEnergy::SetUpBaseFields ( SpatialDomains::MeshGraphSharedPtr & graphShrPtr )

protected

Setup the base fields in case of stability analyses.

Definition at line 372 of file FilterModalEnergy.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::MultiRegions::eDiscontinuous, Nektar::LibUtilities::eFourier, Nektar::LibUtilities::eFourierEvenlySpaced, Nektar::LibUtilities::eFourierHalfModeRe, Nektar::LibUtilities::eFourierSingleModeSpaced, Nektar::MultiRegions::eGalerkin, Nektar::MultiRegions::eMixed_CG_Discontinuous, m_base, m_HalfMode, m_homogen_dealiasing, m_isHomogeneous1D, m_LhomZ, m_MultipleModes, m_npointsZ, m_projectionType, Nektar::SolverUtils::Filter::m_session, m_SingleMode, and m_useFFT.

Referenced by v_Update().

 {
     int i;
     int m_expdim  = graphShrPtr->GetMeshDimension();
 
     //definition of the projection tipe:
     if(m_session->DefinesSolverInfo("PROJECTION"))
     {
         std::string ProjectStr = m_session->GetSolverInfo("PROJECTION");
 
         if ((ProjectStr == "Continuous") ||
             (ProjectStr == "Galerkin")   ||
             (ProjectStr == "CONTINUOUS") ||
             (ProjectStr == "GALERKIN"))
         {
             m_projectionType = MultiRegions::eGalerkin;
         }
         else if ((ProjectStr == "MixedCGDG") ||
                  (ProjectStr == "Mixed_CG_Discontinuous"))
         {
             m_projectionType = MultiRegions::eMixed_CG_Discontinuous;
         }
         else if(ProjectStr == "DisContinuous")
         {
             m_projectionType = MultiRegions::eDiscontinuous;
         }
         else
         {
             ASSERTL0(false, "PROJECTION value not recognised");
         }
     }
     else
     {
         cerr << "Projection type not specified in SOLVERINFO,"
                         "defaulting to continuous Galerkin" << endl;
         m_projectionType = MultiRegions::eGalerkin;
     }
 
     if (m_session->DefinesSolverInfo("ModeType"))
     {
         m_session->MatchSolverInfo("ModeType", "SingleMode",
                                    m_SingleMode, false);
         m_session->MatchSolverInfo("ModeType", "HalfMode",
                                    m_HalfMode, false);
         m_session->MatchSolverInfo("ModeType", "MultipleModes",
                                    m_MultipleModes, false);
     }
 
     m_session->MatchSolverInfo("USEFFT","FFTW", m_useFFT, false);
     m_session->MatchSolverInfo("DEALIASING", "True",
                                m_homogen_dealiasing, false);
 
     if (m_homogen_dealiasing == false)
     {
         m_session->MatchSolverInfo("DEALIASING", "On",
                                    m_homogen_dealiasing, false);
     }
 
     // Stability Analysis flags
     if (m_session->DefinesSolverInfo("ModeType"))
     {
         if (m_SingleMode)
         {
             m_npointsZ = 2;
         }
         else if (m_HalfMode)
         {
             m_npointsZ = 1;
         }
         else if (m_MultipleModes)
         {
             m_npointsZ = m_session->GetParameter("HomModesZ");
         }
         else
         {
             ASSERTL0(false, "SolverInfo ModeType not valid");
         }
     }
     else
     {
         m_npointsZ = m_session->GetParameter("HomModesZ");
     }
 
     if (m_projectionType == MultiRegions::eGalerkin ||
         m_projectionType == MultiRegions::eMixed_CG_Discontinuous)
     {
         switch (m_expdim)
         {
             case 1:
             {
                 for(i = 0; i < m_base.num_elements(); i++)
                 {
                     m_base[i] = MemoryManager<MultiRegions::ContField1D>
                         ::AllocateSharedPtr(m_session, graphShrPtr,
                                             m_session->GetVariable(0));
                 }
             }
             break;
             case 2:
             {
                 if (m_isHomogeneous1D)
                 {
                     if (m_SingleMode)
                     {
                         const LibUtilities::PointsKey PkeyZ(
                             m_npointsZ,
                             LibUtilities::eFourierSingleModeSpaced);
                         const LibUtilities::BasisKey  BkeyZ(
                             LibUtilities::eFourier,
                             m_npointsZ, PkeyZ);
 
                         for (i = 0 ; i < m_base.num_elements(); i++)
                         {
                             m_base[i] = MemoryManager<MultiRegions::
                                 ContField3DHomogeneous1D>::
                                     AllocateSharedPtr(
                                         m_session, BkeyZ, m_LhomZ,
                                         m_useFFT, m_homogen_dealiasing,
                                         graphShrPtr,
                                         m_session->GetVariable(i));
 
                             m_base[i]->SetWaveSpace(true);
                         }
                     }
                     else if (m_HalfMode)
                     {
                         //1 plane field (half mode expansion)
                         const LibUtilities::PointsKey PkeyZ(
                             m_npointsZ,
                             LibUtilities::eFourierSingleModeSpaced);
                         const LibUtilities::BasisKey  BkeyZ(
                             LibUtilities::eFourierHalfModeRe,
                             m_npointsZ,PkeyZ);
 
                         for (i = 0 ; i < m_base.num_elements(); i++)
                         {
                             m_base[i] = MemoryManager<MultiRegions::
                                 ContField3DHomogeneous1D>::
                                     AllocateSharedPtr(
                                         m_session, BkeyZ, m_LhomZ,
                                         m_useFFT, m_homogen_dealiasing,
                                         graphShrPtr,
                                         m_session->GetVariable(i));
 
                             m_base[i]->SetWaveSpace(true);
                         }
                     }
                     else
                     {
                         const LibUtilities::PointsKey PkeyZ(
                             m_npointsZ,
                             LibUtilities::eFourierEvenlySpaced);
                         const LibUtilities::BasisKey  BkeyZ(
                             LibUtilities::eFourier, m_npointsZ, PkeyZ);
 
                         for (i = 0 ; i < m_base.num_elements(); i++)
                         {
                             m_base[i] = MemoryManager<MultiRegions::
                                 ContField3DHomogeneous1D>::
                                     AllocateSharedPtr(
                                         m_session, BkeyZ, m_LhomZ,
                                         m_useFFT, m_homogen_dealiasing,
                                         graphShrPtr,
                                         m_session->GetVariable(i));
 
                             m_base[i]->SetWaveSpace(false);
                         }
                     }
                 }
                 else
                 {
                     i = 0;
                     MultiRegions::ContField2DSharedPtr firstbase =
                         MemoryManager<MultiRegions::ContField2D>::
                             AllocateSharedPtr(
                                 m_session,graphShrPtr,
                                 m_session->GetVariable(i));
 
                     m_base[0] = firstbase;
 
                     for (i = 1 ; i < m_base.num_elements(); i++)
                     {
                         m_base[i] = MemoryManager<MultiRegions::
                             ContField2D>::AllocateSharedPtr(
                                 *firstbase, graphShrPtr,
                                 m_session->GetVariable(i));
                     }
                 }
             }
             break;
             case 3:
             {
                 MultiRegions::ContField3DSharedPtr firstbase =
                     MemoryManager<MultiRegions::ContField3D>::
                         AllocateSharedPtr(m_session, graphShrPtr,
                                           m_session->GetVariable(0));
                 m_base[0] = firstbase;
                 for (i = 1 ; i < m_base.num_elements(); i++)
                 {
                     m_base[i] = MemoryManager<MultiRegions::
                     ContField3D>::AllocateSharedPtr(
                         *firstbase, graphShrPtr,
                         m_session->GetVariable(0));
                 }
             }
             break;
             default:
                 ASSERTL0(false, "Expansion dimension not recognised");
                 break;
         }
     }
     else
     {
         switch (m_expdim)
         {
             case 1:
             {
                 // need to use zero for variable as may be more base
                 // flows than variables
                 for (i = 0 ; i < m_base.num_elements(); i++)
                 {
                     m_base[i] = MemoryManager<MultiRegions::
                         DisContField1D>::AllocateSharedPtr(
                             m_session, graphShrPtr,
                             m_session->GetVariable(0));
                 }
                 break;
             }
             case 2:
             {
                 for (i = 0 ; i < m_base.num_elements(); i++)
                 {
                     m_base[i] = MemoryManager<MultiRegions::
                         DisContField2D>::AllocateSharedPtr(
                             m_session, graphShrPtr,
                             m_session->GetVariable(0));
                 }
                 break;
             }
             case 3:
                 ASSERTL0(false, "3D not set up");
             default:
                 ASSERTL0(false, "Expansion dimension not recognised");
                 break;
         }
     }
 }

void Nektar::SolverUtils::FilterModalEnergy::v_Finalise	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		const NekDouble &	time
	)

protectedvirtual

Close the output stream.

Implements Nektar::SolverUtils::Filter.

Definition at line 334 of file FilterModalEnergy.cpp.

References m_outputStream.

 {
     if (pFields[0]->GetComm()->GetRank() == 0)
     {
         m_outputStream.close();
     }
 }

void Nektar::SolverUtils::FilterModalEnergy::v_Initialise	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		const NekDouble &	time
	)

protectedvirtual

Initialize the parallel communication and the output stream.

Implements Nektar::SolverUtils::Filter.

Definition at line 128 of file FilterModalEnergy.cpp.

References m_index, m_isHomogeneous1D, m_outputFile, m_outputStream, and v_Update().

 {
     LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
 
     if (vComm->GetRank() == 0)
     {
         // Open output stream
         if(m_isHomogeneous1D)
         {
             m_outputStream.open(m_outputFile.c_str());
             m_outputStream << "# Time,  Fourier Mode, Energy ";
             m_outputStream << endl;
         }
         else
         {
             m_outputStream.open(m_outputFile.c_str());
             m_outputStream << "# Time, Energy ";
             m_outputStream << endl;
         }
 
     }
 
     m_index = 0;
     v_Update(pFields, time);
 }

bool Nektar::SolverUtils::FilterModalEnergy::v_IsTimeDependent ( )

protectedvirtual

Flag for time-dependent flows.

Implements Nektar::SolverUtils::Filter.

Definition at line 661 of file FilterModalEnergy.cpp.

 {
     return true;
 }

void Nektar::SolverUtils::FilterModalEnergy::v_Update	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		const NekDouble &	time
	)

protectedvirtual

Update the modal energy every m_outputFrequency.

Implements Nektar::SolverUtils::Filter.

Definition at line 160 of file FilterModalEnergy.cpp.

References ASSERTL0, ImportFldBase(), L2Error(), m_base, m_EqTypeStr, m_index, m_isHomogeneous1D, m_isHomogeneous2D, m_npointsZ, m_outputFrequency, m_outputStream, m_PertEnergy, Nektar::SolverUtils::Filter::m_session, Nektar::SpatialDomains::MeshGraph::Read(), SetUpBaseFields(), Vmath::Vadd(), and Vmath::Vsub().

Referenced by v_Initialise().

 {
     // Only output every m_outputFrequency
     if ((m_index++) % m_outputFrequency)
     {
         return;
     }
 
     LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
 
     // Homogeneous 1D implementation
     if (m_isHomogeneous1D)
     {
         int colrank = vComm->GetColumnComm()->GetRank();
         int nproc   = vComm->GetColumnComm()->GetSize();
         m_npointsZ  = (m_session->GetParameter("HomModesZ"));
         int locsize = m_npointsZ/nproc/2;
 
         Array<OneD, NekDouble> energy    (locsize, 0.0);
         Array<OneD, NekDouble> energy_tmp(locsize, 0.0);
         Array<OneD, NekDouble> tmp;
 
         // Calculate the energy of the perturbation for stability
         // analysis
         if (m_PertEnergy)
         {
             // Compressible Flow Solver
             if (m_EqTypeStr=="EulerCFE"   ||
                 m_EqTypeStr=="EulerADCFE" ||
                 m_EqTypeStr=="NavierStokesCFE")
             {
                 ASSERTL0(false, "Stability analysis module not "
                          "implemented for the Compressible Flow "
                          "Solver. Please remove the function BaseFlow "
                          "from your .xml file");
             }
             // Incompressible Navier-Stokes Solver
             else
             {
                 SpatialDomains::MeshGraphSharedPtr graphShrPtr =
                     SpatialDomains::MeshGraph::Read(m_session);
                 SetUpBaseFields(graphShrPtr);
                 string file = m_session->
                     GetFunctionFilename("BaseFlow", 0);
                 ImportFldBase(file, graphShrPtr);
 
                 for (int i = 0; i < pFields.num_elements()-1; ++i)
                 {
                     Vmath::Vsub(pFields[i]->GetNcoeffs(),
                                 pFields[i]->GetCoeffs(), 1,
                                 m_base [i]->GetCoeffs(), 1,
                                 pFields[i]->UpdateCoeffs(), 1);
 
                     energy_tmp = pFields[i]->HomogeneousEnergy();
                     Vmath::Vadd(locsize, energy_tmp, 1,
                                 energy, 1, energy, 1);
 
                     Vmath::Vadd(pFields[i]->GetNcoeffs(),
                                 pFields[i]->GetCoeffs(), 1,
                                 m_base[i]->GetCoeffs(), 1,
                                 pFields[i]->UpdateCoeffs(), 1);
                 }
             }
         }
         // Calculate the modal energy for general simulation
         else
         {
             // Compressible Flow Solver
             if (m_EqTypeStr=="EulerCFE"   ||
                 m_EqTypeStr=="EulerADCFE" ||
                 m_EqTypeStr=="NavierStokesCFE")
             {
                 // Extracting kinetic energy
                 for (int i = 1; i < pFields.num_elements()-1; ++i)
                 {
                     energy_tmp = pFields[i]->HomogeneousEnergy();
                     Vmath::Vadd(locsize, energy_tmp, 1,
                                 energy, 1, energy, 1);
                 }
             }
             // Incompressible Navier-Stokes Solver
             else
             {
                 // Extracting kinetic energy
                 for (int i = 0; i < pFields.num_elements()-1; ++i)
                 {
                     energy_tmp = pFields[i]->HomogeneousEnergy();
                     Vmath::Vadd(locsize, energy_tmp, 1,
                                 energy, 1, energy, 1);
                 }
             }
         }
 
         // Send to root process
         if (colrank == 0)
         {
             int j, m = 0;
 
             for (j = 0; j < energy.num_elements(); ++j, ++m)
             {
                 m_outputStream << setw(10) << time
                                << setw(5)  << m
                                << setw(18) << energy[j] << endl;
             }
 
             for (int i = 1; i < nproc; ++i)
             {
                 vComm->GetColumnComm()->Recv(i, energy);
 
                 for (j = 0; j < energy.num_elements(); ++j, ++m)
                 {
                     m_outputStream << setw(10) << time
                                    << setw(5)  << m
                                    << setw(18) << energy[j] << endl;
                 }
             }
         }
         else
         {
             vComm->GetColumnComm()->Send(0, energy);
         }
     }
     // Homogeneous 2D implementation
     else if (m_isHomogeneous2D)
     {
         ASSERTL0(false, "3D Homogeneous 2D energy "
                  "dumping not implemented yet");
     }
     // General implementation
     else
     {
         // Compressible Flow Solver
         if (m_EqTypeStr=="EulerCFE"   ||
             m_EqTypeStr=="EulerADCFE" ||
             m_EqTypeStr=="NavierStokesCFE")
         {
             // Total energy
             NekDouble energy = 0.0;
             for (int i = 1; i < pFields.num_elements()-1; ++i)
             {
                 pFields[i]->SetPhysState(true);
                 NekDouble norm = L2Error(pFields, i, time);
                 energy += norm * norm;
             }
 
             m_outputStream << setprecision(6) << time;
             m_outputStream.width(25);
             m_outputStream << setprecision(8) << 0.5*energy;
             m_outputStream << endl;
         }
         // Incompressible Navier-Stokes Solver
         else
         {
             // Kinetic energy
             NekDouble energy = 0.0;
             for (int i = 0; i < pFields.num_elements()-1; ++i)
             {
                 pFields[i]->SetPhysState(true);
                 NekDouble norm = L2Error(pFields, i, time);
                 energy += norm * norm;
             }
             m_outputStream << setprecision(6) << time;
             m_outputStream.width(25);
             m_outputStream << setprecision(8) << 0.5*energy;
             m_outputStream << endl;
         }
     }
 }