Commit b61fb1d2 authored by Douglas Serson's avatar Douglas Serson

Add new field convert module for evaluating the average value of variables on the domain

parent 30986d02
......@@ -18,6 +18,7 @@ SET(FieldUtilsHeaders
OutputModules/OutputXml.h
ProcessModules/ProcessAddCompositeID.h
ProcessModules/ProcessAddFld.h
ProcessModules/ProcessAvgValue.h
ProcessModules/ProcessBoundaryExtract.h
ProcessModules/ProcessCombineAvg.h
ProcessModules/ProcessConcatenateFld.h
......@@ -72,6 +73,7 @@ SET(FieldUtilsSources
OutputModules/OutputXml.cpp
ProcessModules/ProcessAddCompositeID.cpp
ProcessModules/ProcessAddFld.cpp
ProcessModules/ProcessAvgValue.cpp
ProcessModules/ProcessBoundaryExtract.cpp
ProcessModules/ProcessCombineAvg.cpp
ProcessModules/ProcessConcatenateFld.cpp
......
////////////////////////////////////////////////////////////////////////////////
//
// File: ProcessAvgValue.cpp
//
// For more information, please see: http://www.nektar.info/
//
// The MIT License
//
// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
// Department of Aeronautics, Imperial College London (UK), and Scientific
// Computing and Imaging Institute, University of Utah (USA).
//
// License for the specific language governing rights and limitations under
// Permission is hereby granted, free of charge, to any person obtaining a
// copy of this software and associated documentation files (the "Software"),
// to deal in the Software without restriction, including without limitation
// the rights to use, copy, modify, merge, publish, distribute, sublicense,
// and/or sell copies of the Software, and to permit persons to whom the
// Software is furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included
// in all copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
// DEALINGS IN THE SOFTWARE.
//
// Description: Compute the average value of each field.
//
////////////////////////////////////////////////////////////////////////////////
#include <iostream>
#include <string>
using namespace std;
#include "ProcessAvgValue.h"
#include <LibUtilities/BasicUtils/SharedArray.hpp>
namespace Nektar
{
namespace FieldUtils
{
ModuleKey ProcessAvgValue::className =
GetModuleFactory().RegisterCreatorFunction(
ModuleKey(eProcessModule, "avgvalue"),
ProcessAvgValue::create,
"compute the average value of each field.");
ProcessAvgValue::ProcessAvgValue(FieldSharedPtr f) : ProcessModule(f)
{
}
ProcessAvgValue::~ProcessAvgValue()
{
}
void ProcessAvgValue::Process(po::variables_map &vm)
{
int nfields = m_f->m_variables.size();
int spacedim = m_f->m_graph->GetMeshDimension() + m_f->m_numHomogeneousDir;
int npoints = m_f->m_exp[0]->GetNpoints();
// Calculate volume (or area)
Array<OneD, NekDouble> ones(npoints, 1.0);
NekDouble scale = m_f->m_exp[0]->Integral(ones);
// Output volume
string name[3] = {"length", "area", "volume"};
cout << "Domain " << name[spacedim-1] << " : " << scale << endl;
// Calculate integral and average of each field
for(int i = 0; i < nfields; ++i)
{
NekDouble integral = m_f->m_exp[0]->Integral(m_f->m_exp[i]->GetPhys());
if (m_f->m_comm->GetRank() == 0)
{
cout << "Integral (variable " << m_f->m_variables[i]
<< ") : " << integral << endl;
cout << "Average value (variable " << m_f->m_variables[i]
<< ") : " << integral / scale << endl;
}
}
}
}
}
////////////////////////////////////////////////////////////////////////////////
//
// File: ProcessAvgValue.h
//
// For more information, please see: http://www.nektar.info/
//
// The MIT License
//
// Copyright (c) 2006 Division of Applied Mathematics, Brown University (USA),
// Department of Aeronautics, Imperial College London (UK), and Scientific
// Computing and Imaging Institute, University of Utah (USA).
//
// License for the specific language governing rights and limitations under
// Permission is hereby granted, free of charge, to any person obtaining a
// copy of this software and associated documentation files (the "Software"),
// to deal in the Software without restriction, including without limitation
// the rights to use, copy, modify, merge, publish, distribute, sublicense,
// and/or sell copies of the Software, and to permit persons to whom the
// Software is furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included
// in all copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
// OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
// THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
// FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
// DEALINGS IN THE SOFTWARE.
//
// Description: Compute the average value of each field.
//
////////////////////////////////////////////////////////////////////////////////
#ifndef FIELDUTILS_PROCESSAVGVALUE
#define FIELDUTILS_PROCESSAVGVALUE
#include "../Module.h"
namespace Nektar
{
namespace FieldUtils
{
/**
* @brief This processing module computes the average value of each field.
*
*/
class ProcessAvgValue : public ProcessModule
{
public:
/// Creates an instance of this class
static std::shared_ptr<Module> create(FieldSharedPtr f)
{
return MemoryManager<ProcessAvgValue>::AllocateSharedPtr(f);
}
static ModuleKey className;
ProcessAvgValue(FieldSharedPtr f);
virtual ~ProcessAvgValue();
/// Write mesh to output file.
virtual void Process(po::variables_map &vm);
virtual std::string GetModuleName()
{
return "ProcessAvgValue";
}
virtual std::string GetModuleDescription()
{
return "Evaluating average value";
}
virtual ModulePriority GetModulePriority()
{
return eModifyExp;
}
};
}
}
#endif
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment