OpenLB 1.8.1
Loading...
Searching...
No Matches
olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR > Class Template Reference

functor to compute 0.5*(f[extractDim]-f_wanted[0])^2/f_wanted^2 on a lattice More...

#include <dualFunctors3D.h>

+ Inheritance diagram for olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR >:
+ Collaboration diagram for olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR >:

Public Member Functions

 DrelativeDifferenceObjectiveComponentDf3D (SuperLattice< T, DESCRIPTOR > &sLattice, FunctorPtr< SuperF3D< T, T > > &&f, int extractDim, FunctorPtr< SuperF3D< T, T > > &&dFdF, FunctorPtr< AnalyticalF3D< T, T > > &&wantedF, FunctorPtr< SuperIndicatorF3D< T > > &&indicatorF)
 
bool operator() (T output[], const int input[])
 
- Public Member Functions inherited from olb::SuperF3D< T, T >
SuperF3D< T, W > & operator- (SuperF3D< T, W > &rhs)
 
SuperF3D< T, W > & operator+ (SuperF3D< T, W > &rhs)
 
SuperF3D< T, W > & operator* (SuperF3D< T, W > &rhs)
 
SuperF3D< T, W > & operator/ (SuperF3D< T, W > &rhs)
 
SuperStructure< T, 3 > & getSuperStructure ()
 
int getBlockFSize () const
 
BlockF3D< W > & getBlockF (int iCloc)
 
bool operator() (W output[], const int input[])
 has to be implemented for 'every' derived class
 
- Public Member Functions inherited from olb::GenericF< W, int >
virtual ~GenericF ()=default
 
int getSourceDim () const
 read only access to member variable _m
 
int getTargetDim () const
 read only access to member variable _n
 
std::string & getName ()
 read and write access to name
 
std::string const & getName () const
 read only access to name
 
bool operator() (W output[])
 wrapper that call the pure virtual operator() (T output[], const S input[]) from above
 
bool operator() (W output[], int input0)
 
bool operator() (W output[], int input0, int input1)
 
bool operator() (W output[], int input0, int input1, int input2)
 
bool operator() (W output[], int input0, int input1, int input2, int input3)
 

Additional Inherited Members

- Public Types inherited from olb::SuperF3D< T, T >
using identity_functor_type
 
- Public Types inherited from olb::GenericF< W, int >
using targetType
 
using sourceType
 
- Public Attributes inherited from olb::GenericF< W, int >
std::shared_ptr< GenericF< W, int > > _ptrCalcC
 memory management, frees resouces (calcClass)
 
- Static Public Attributes inherited from olb::SuperF3D< T, T >
static constexpr bool isSuper
 
static constexpr unsigned d
 
- Protected Member Functions inherited from olb::SuperF3D< T, T >
 SuperF3D (SuperStructure< T, 3 > &superStructure, int targetDim)
 
- Protected Member Functions inherited from olb::GenericF< W, int >
 GenericF (int targetDim, int sourceDim)
 
- Protected Attributes inherited from olb::SuperF3D< T, T >
SuperStructure< T, 3 > & _superStructure
 
std::vector< std::unique_ptr< BlockF3D< W > > > _blockF
 Super functors may consist of several BlockF3D<W> derived functors.
 

Detailed Description

template<typename T, typename DESCRIPTOR>
class olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR >

functor to compute 0.5*(f[extractDim]-f_wanted[0])^2/f_wanted^2 on a lattice

Definition at line 278 of file dualFunctors3D.h.

Constructor & Destructor Documentation

◆ DrelativeDifferenceObjectiveComponentDf3D()

template<typename T , typename DESCRIPTOR >
olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR >::DrelativeDifferenceObjectiveComponentDf3D ( SuperLattice< T, DESCRIPTOR > & sLattice,
FunctorPtr< SuperF3D< T, T > > && f,
int extractDim,
FunctorPtr< SuperF3D< T, T > > && dFdF,
FunctorPtr< AnalyticalF3D< T, T > > && wantedF,
FunctorPtr< SuperIndicatorF3D< T > > && indicatorF )

Definition at line 593 of file dualFunctors3D.hh.

599 : SuperF3D<T,T>(sLattice, DESCRIPTOR::q),
600 _f(std::forward<decltype(f)>(f), extractDim),
601 _dFdF(std::move(dFdF)),
602 _wantedF(std::forward<decltype(wantedF)>(wantedF), sLattice),
603 _indicatorF(std::move(indicatorF))
604{
605 this->getName() = "dRelativeDifferenceObjectiveDf";
606
607 SuperL2Norm3D<T> wantedFlp(_wantedF, *_indicatorF);
608 T output[1] {};
609 wantedFlp(output, nullptr);
610
611 const T globalValue = util::pow(output[0], 2);
612 const T weight = util::pow(_f.getSuperStructure().getCuboidDecomposition().getDeltaR(), 3);
613
614 for (int iC = 0; iC < sLattice.getLoadBalancer().size(); ++iC) {
615 this->_blockF.emplace_back(
616 new BlockDrelativeDifferenceObjectiveComponentDf3D<T,DESCRIPTOR>(
617 sLattice.getBlock(iC),
618 static_cast<BlockExtractComponentF3D<T>&>(_f.getBlockF(iC)),
619 _dFdF->getBlockF(iC),
620 _wantedF.getBlockF(iC),
621 _indicatorF->getBlockIndicatorF(iC),
622 globalValue,
623 weight)
624 );
625 }
626}
std::string & getName()
Definition genericF.hh:51
SuperStructure< T, 3 > & getSuperStructure()
BlockF3D< W > & getBlockF(int iCloc)
std::vector< std::unique_ptr< BlockF3D< W > > > _blockF
CuboidDecomposition< T, D > & getCuboidDecomposition()
Read and write access to cuboid geometry.
Expr pow(Expr base, Expr exp)
Definition expr.cpp:235
SuperLpNorm3D< T, W, 2 > SuperL2Norm3D
Functor that returns the L2 norm over omega of the the euklid norm of the input functor.

References olb::SuperF3D< T, T >::_blockF, olb::SuperLattice< T, DESCRIPTOR >::getBlock(), olb::SuperF3D< T, U >::getBlockF(), olb::SuperStructure< T, D >::getCuboidDecomposition(), olb::SuperStructure< T, D >::getLoadBalancer(), olb::GenericF< W, int >::getName(), olb::SuperF3D< T, U >::getSuperStructure(), and olb::util::pow().

+ Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename DESCRIPTOR >
bool olb::opti::DrelativeDifferenceObjectiveComponentDf3D< T, DESCRIPTOR >::operator() ( T output[],
const int input[] )

Definition at line 629 of file dualFunctors3D.hh.

630{
631 for (int i=0; i < DESCRIPTOR::q; ++i) {
632 dJdF[i] = T{};
633 }
634
635 auto& load = _f.getSuperStructure().getLoadBalancer();
636
637 if (load.isLocal(latticeR[0])) {
638 return this->getBlockF(load.loc(latticeR[0]))(dJdF, &latticeR[1]);
639 }
640 else {
641 return true;
642 }
643}
LoadBalancer< T > & getLoadBalancer()
Read and write access to the load balancer.

The documentation for this class was generated from the following files: