Nektar++
Public Member Functions | Static Public Member Functions | Static Public Attributes | List of all members
Nektar::FieldUtils::ProcessJacobianEnergy Class Reference

This processing module scales the input fld file. More...

#include <ProcessJacobianEnergy.h>

Inheritance diagram for Nektar::FieldUtils::ProcessJacobianEnergy:
[legend]

Public Member Functions

 ProcessJacobianEnergy (FieldSharedPtr f)
 
virtual ~ProcessJacobianEnergy ()
 
virtual void Process (po::variables_map &vm)
 Write mesh to output file. More...
 
virtual std::string GetModuleName ()
 
virtual std::string GetModuleDescription ()
 
virtual ModulePriority GetModulePriority ()
 
- Public Member Functions inherited from Nektar::FieldUtils::ProcessModule
 ProcessModule ()
 
 ProcessModule (FieldSharedPtr p_f)
 
- Public Member Functions inherited from Nektar::FieldUtils::Module
FIELD_UTILS_EXPORT Module (FieldSharedPtr p_f)
 
virtual ~Module ()=default
 
const ConfigOptionGetConfigOption (const std::string &key) const
 
FIELD_UTILS_EXPORT void RegisterConfig (std::string key, std::string value="")
 Register a configuration option with a module. More...
 
FIELD_UTILS_EXPORT void PrintConfig ()
 Print out all configuration options for a module. More...
 
FIELD_UTILS_EXPORT void SetDefaults ()
 Sets default configuration options for those which have not been set. More...
 
FIELD_UTILS_EXPORT void AddFile (std::string fileType, std::string fileName)
 
FIELD_UTILS_EXPORT void EvaluateTriFieldAtEquiSpacedPts (LocalRegions::ExpansionSharedPtr &exp, const Array< OneD, const NekDouble > &infield, Array< OneD, NekDouble > &outfield)
 

Static Public Member Functions

static std::shared_ptr< Modulecreate (FieldSharedPtr f)
 Creates an instance of this class. More...
 

Static Public Attributes

static ModuleKey className
 

Additional Inherited Members

- Public Attributes inherited from Nektar::FieldUtils::Module
FieldSharedPtr m_f
 Field object. More...
 
- Protected Member Functions inherited from Nektar::FieldUtils::Module
 Module ()
 
- Protected Attributes inherited from Nektar::FieldUtils::Module
std::map< std::string, ConfigOptionm_config
 List of configuration values. More...
 
std::set< std::string > m_allowedFiles
 List of allowed file formats. More...
 

Detailed Description

This processing module scales the input fld file.

Definition at line 46 of file ProcessJacobianEnergy.h.

Constructor & Destructor Documentation

◆ ProcessJacobianEnergy()

Nektar::FieldUtils::ProcessJacobianEnergy::ProcessJacobianEnergy ( FieldSharedPtr  f)

Definition at line 56 of file ProcessJacobianEnergy.cpp.

57  : ProcessModule(f)
58 {
59  m_config["topmodes"] =
60  ConfigOption(false, "1", "how many top modes to keep ");
61 }
std::map< std::string, ConfigOption > m_config
List of configuration values.
Definition: Module.h:233

References Nektar::FieldUtils::Module::m_config.

◆ ~ProcessJacobianEnergy()

Nektar::FieldUtils::ProcessJacobianEnergy::~ProcessJacobianEnergy ( )
virtual

Definition at line 63 of file ProcessJacobianEnergy.cpp.

64 {
65 }

Member Function Documentation

◆ create()

static std::shared_ptr<Module> Nektar::FieldUtils::ProcessJacobianEnergy::create ( FieldSharedPtr  f)
inlinestatic

Creates an instance of this class.

Definition at line 50 of file ProcessJacobianEnergy.h.

51  {
53  }
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.

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

◆ GetModuleDescription()

virtual std::string Nektar::FieldUtils::ProcessJacobianEnergy::GetModuleDescription ( )
inlinevirtual

Reimplemented from Nektar::FieldUtils::Module.

Definition at line 67 of file ProcessJacobianEnergy.h.

68  {
69  return "Processing Jacobian";
70  }

◆ GetModuleName()

virtual std::string Nektar::FieldUtils::ProcessJacobianEnergy::GetModuleName ( )
inlinevirtual

Implements Nektar::FieldUtils::Module.

Definition at line 62 of file ProcessJacobianEnergy.h.

63  {
64  return "ProcessJacobianEnergy";
65  }

◆ GetModulePriority()

virtual ModulePriority Nektar::FieldUtils::ProcessJacobianEnergy::GetModulePriority ( )
inlinevirtual

Implements Nektar::FieldUtils::Module.

Definition at line 72 of file ProcessJacobianEnergy.h.

73  {
74  return eModifyExp;
75  }

References Nektar::FieldUtils::eModifyExp.

◆ Process()

void Nektar::FieldUtils::ProcessJacobianEnergy::Process ( po::variables_map &  vm)
virtual

Write mesh to output file.

Implements Nektar::FieldUtils::Module.

Definition at line 67 of file ProcessJacobianEnergy.cpp.

68 {
69  m_f->SetUpExp(vm);
70 
71  int nfields = m_f->m_variables.size();
72  m_f->m_variables.push_back("jacenergy");
73  // Skip in case of empty partition
74  if (m_f->m_exp[0]->GetNumElmts() == 0)
75  {
76  return;
77  }
78 
79  int NumHomogeneousDir = m_f->m_numHomogeneousDir;
81 
82  if (nfields)
83  {
84  m_f->m_exp.resize(nfields + 1);
85  exp = m_f->AppendExpList(NumHomogeneousDir);
86 
87  m_f->m_exp[nfields] = exp;
88  }
89  else
90  {
91  exp = m_f->m_exp[0];
92  }
93 
94  Array<OneD, NekDouble> phys = exp->UpdatePhys();
95  Array<OneD, NekDouble> coeffs = exp->UpdateCoeffs();
96  Array<OneD, NekDouble> tmp;
97 
98  for (int i = 0; i < exp->GetExpSize(); ++i)
99  {
100  // copy Jacobian into field
101  StdRegions::StdExpansionSharedPtr Elmt = exp->GetExp(i);
102 
103  const StdRegions::StdExpansion * sep = &( *Elmt );
104  const LocalRegions::Expansion * lep = dynamic_cast<const LocalRegions::Expansion*>( sep );
105 
106  int nquad = Elmt->GetTotPoints();
107  int coeffoffset = exp->GetCoeff_Offset(i);
108  Array<OneD, const NekDouble> Jac = lep->GetMetricInfo()->GetJac( Elmt->GetPointsKeys() );
109  if ( lep->GetMetricInfo()->GetGtype() == SpatialDomains::eRegular )
110  {
111  Vmath::Fill(nquad, Jac[0], phys, 1);
112  }
113  else
114  {
115  Vmath::Vcopy(nquad, Jac, 1, phys, 1);
116  }
117 
118  if( lep->GetMetricInfo()->GetGtype() == SpatialDomains::eDeformed )
119  {
120  NekDouble jacmax = Vmath::Vmax(nquad, Jac, 1);
121  NekDouble jacmin = Vmath::Vmin(nquad, Jac, 1);
122 
123  NekDouble jacmeasure = jacmax / jacmin - 1.0;
124  Vmath::Fill(nquad, jacmeasure, phys, 1);
125  }
126  else
127  {
128  Vmath::Fill(nquad, 0.0, phys, 1);
129  }
130 
131  Elmt->FwdTrans(phys, tmp = coeffs + coeffoffset);
132  }
133  exp->BwdTrans(coeffs, phys);
134 }
FieldSharedPtr m_f
Field object.
Definition: Module.h:230
std::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
@ eRegular
Geometry is straight-sided with constant geometric factors.
@ eDeformed
Geometry is curved or has non-constant factors.
std::shared_ptr< StdExpansion > StdExpansionSharedPtr
double NekDouble
T Vmin(int n, const T *x, const int incx)
Return the minimum element in x - called vmin to avoid conflict with min.
Definition: Vmath.cpp:992
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:45
T Vmax(int n, const T *x, const int incx)
Return the maximum element in x – called vmax to avoid conflict with max.
Definition: Vmath.cpp:892
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1199

References Nektar::SpatialDomains::eDeformed, Nektar::SpatialDomains::eRegular, Vmath::Fill(), Nektar::LocalRegions::Expansion::GetMetricInfo(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::FieldUtils::Module::m_f, Vmath::Vcopy(), Vmath::Vmax(), and Vmath::Vmin().

Member Data Documentation

◆ className

ModuleKey Nektar::FieldUtils::ProcessJacobianEnergy::className
static
Initial value:
=
ModuleKey(eProcessModule, "jacobianenergy"),
"Show high frequency energy of Jacobian.")
static std::shared_ptr< Module > create(FieldSharedPtr f)
Creates an instance of this class.
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, std::string pDesc="")
Register a class with the factory.
Definition: NekFactory.hpp:200
std::pair< ModuleType, std::string > ModuleKey
Definition: Module.h:290
ModuleFactory & GetModuleFactory()
Definition: Module.cpp:49

Definition at line 54 of file ProcessJacobianEnergy.h.