OpenLB 1.7
Loading...
Searching...
No Matches
Public Member Functions | Protected Attributes | List of all members
olb::AnalyticalFfromBlockF3D< T, W > Class Template Reference

Converts block functors to analytical functors. More...

#include <interpolationF3D.h>

+ Inheritance diagram for olb::AnalyticalFfromBlockF3D< T, W >:
+ Collaboration diagram for olb::AnalyticalFfromBlockF3D< T, W >:

Public Member Functions

 AnalyticalFfromBlockF3D (BlockF3D< W > &f, Cuboid3D< T > &cuboid)
 
bool operator() (W output[], const T physC[]) override
 trilinear interpolation on cubic lattice
 
- Public Member Functions inherited from olb::AnalyticalF< D, T, S >
AnalyticalF< D, T, S > & operator- (AnalyticalF< D, T, S > &rhs)
 
AnalyticalF< D, T, S > & operator+ (AnalyticalF< D, T, S > &rhs)
 
AnalyticalF< D, T, S > & operator* (AnalyticalF< D, T, S > &rhs)
 
AnalyticalF< D, T, S > & operator/ (AnalyticalF< D, T, S > &rhs)
 
- Public Member Functions inherited from olb::GenericF< T, S >
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
 
virtual bool operator() (T output[], const S input[])=0
 has to be implemented for 'every' derived class
 
bool operator() (T output[])
 wrapper that call the pure virtual operator() (T output[], const S input[]) from above
 
bool operator() (T output[], S input0)
 
bool operator() (T output[], S input0, S input1)
 
bool operator() (T output[], S input0, S input1, S input2)
 
bool operator() (T output[], S input0, S input1, S input2, S input3)
 

Protected Attributes

BlockF3D< W > & _f
 
Cuboid3D< T > & _cuboid
 

Additional Inherited Members

- Public Types inherited from olb::AnalyticalF< D, T, S >
using identity_functor_type = AnalyticalIdentity<D,T,S>
 
- Public Types inherited from olb::GenericF< T, S >
using targetType = T
 
using sourceType = S
 
- Public Attributes inherited from olb::GenericF< T, S >
std::shared_ptr< GenericF< T, S > > _ptrCalcC
 memory management, frees resouces (calcClass)
 
- Static Public Attributes inherited from olb::AnalyticalF< D, T, S >
static constexpr unsigned dim = D
 
- Protected Member Functions inherited from olb::AnalyticalF< D, T, S >
 AnalyticalF (int n)
 
- Protected Member Functions inherited from olb::GenericF< T, S >
 GenericF (int targetDim, int sourceDim)
 

Detailed Description

template<typename T, typename W = T>
class olb::AnalyticalFfromBlockF3D< T, W >

Converts block functors to analytical functors.

Definition at line 53 of file interpolationF3D.h.

Constructor & Destructor Documentation

◆ AnalyticalFfromBlockF3D()

template<typename T , typename W >
olb::AnalyticalFfromBlockF3D< T, W >::AnalyticalFfromBlockF3D ( BlockF3D< W > & f,
Cuboid3D< T > & cuboid )

Definition at line 176 of file interpolationF3D.hh.

178 : AnalyticalF3D<T,W>(f.getTargetDim()),
179 _f(f), _cuboid(cuboid)
180{
181 this->getName() = "fromBlockF";
182}
int getTargetDim() const
read only access to member variable _n
Definition genericF.hh:45
std::string & getName()
read and write access to name
Definition genericF.hh:51

References olb::GenericF< T, S >::getName().

+ Here is the call graph for this function:

Member Function Documentation

◆ operator()()

template<typename T , typename W >
bool olb::AnalyticalFfromBlockF3D< T, W >::operator() ( W output[],
const T physC[] )
override

trilinear interpolation on cubic lattice

Definition at line 186 of file interpolationF3D.hh.

187{
188 int latticeC[3];
189 int latticeR[3];
190 _cuboid.getFloorLatticeR(latticeR, physC);
191
192 auto& block = _f.getBlockStructure();
193 auto padding = std::min(1, block.getPadding());
194
195 if (LatticeR<3>(latticeR) >= -padding && LatticeR<3>(latticeR) < block.getExtent()+padding-1) {
196 const int& locX = latticeR[0];
197 const int& locY = latticeR[1];
198 const int& locZ = latticeR[2];
199
200 Vector<T,3> physRiC;
201 Vector<T,3> physCv(physC);
202 _cuboid.getPhysR(physRiC.data(), {locX, locY, locZ});
203
204 // compute weights
205 Vector<W,3> d = (physCv - physRiC) * (1. / _cuboid.getDeltaR());
206 Vector<W,3> e = 1. - d;
207
208 W output_tmp[_f.getTargetDim()];
209 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
210 output_tmp[iD] = W();
211 }
212
213 latticeC[0] = locX;
214 latticeC[1] = locY;
215 latticeC[2] = locZ;
216 _f(output_tmp,latticeC);
217 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
218 output[iD] += output_tmp[iD] * e[0] * e[1] * e[2];
219 output_tmp[iD] = W();
220 }
221
222 latticeC[0] = locX;
223 latticeC[1] = locY + 1;
224 latticeC[2] = locZ;
225 _f(output_tmp,latticeC);
226 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
227 output[iD] += output_tmp[iD] * e[0] * d[1] * e[2];
228 output_tmp[iD] = W();
229 }
230
231 latticeC[0] = locX + 1;
232 latticeC[1] = locY;
233 latticeC[2] = locZ;
234 _f(output_tmp,latticeC);
235 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
236 output[iD] += output_tmp[iD] * d[0] * e[1] * e[2];
237 output_tmp[iD] = W();
238 }
239
240 latticeC[0] = locX + 1;
241 latticeC[1] = locY + 1;
242 latticeC[2] = locZ;
243 _f(output_tmp,latticeC);
244 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
245 output[iD] += output_tmp[iD] * d[0] * d[1] * e[2];
246 output_tmp[iD] = W();
247 }
248
249 latticeC[0] = locX;
250 latticeC[1] = locY;
251 latticeC[2] = locZ + 1;
252 _f(output_tmp,latticeC);
253 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
254 output[iD] += output_tmp[iD] * e[0] * e[1] * d[2];
255 output_tmp[iD] = W();
256 }
257
258 latticeC[0] = locX;
259 latticeC[1] = locY + 1;
260 latticeC[2] = locZ + 1;
261 _f(output_tmp,latticeC);
262 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
263 output[iD] += output_tmp[iD] * e[0] * d[1] * d[2];
264 output_tmp[iD] = W();
265 }
266
267 latticeC[0] = locX + 1;
268 latticeC[1] = locY;
269 latticeC[2] = locZ + 1;
270 _f(output_tmp,latticeC);
271 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
272 output[iD] += output_tmp[iD] * d[0] * e[1] * d[2];
273 output_tmp[iD] = W();
274 }
275
276 latticeC[0] = locX + 1;
277 latticeC[1] = locY + 1;
278 latticeC[2] = locZ + 1;
279 _f(output_tmp,latticeC);
280 for (int iD = 0; iD < _f.getTargetDim(); ++iD) {
281 output[iD] += output_tmp[iD] * d[0] * d[1] * d[2];
282 output_tmp[iD] = W();
283 }
284
285 return true;
286 }
287 else {
288 return false;
289 }
290}
virtual BlockStructureD< 3 > & getBlockStructure() const
constexpr int d() any_platform

References olb::Vector< T, D >::data().

+ Here is the call graph for this function:

Member Data Documentation

◆ _cuboid

template<typename T , typename W = T>
Cuboid3D<T>& olb::AnalyticalFfromBlockF3D< T, W >::_cuboid
protected

Definition at line 56 of file interpolationF3D.h.

◆ _f

template<typename T , typename W = T>
BlockF3D<W>& olb::AnalyticalFfromBlockF3D< T, W >::_f
protected

Definition at line 55 of file interpolationF3D.h.


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