OpenLB 1.7
Loading...
Searching...
No Matches
euklidNorm2D.hh
Go to the documentation of this file.
1/* This file is part of the OpenLB library
2 *
3 * Copyright (C) 2019 Albert Mink, Mathias J. Krause, Lukas Baron
4 * E-mail contact: info@openlb.net
5 * The most recent release of OpenLB can be downloaded at
6 * <http://www.openlb.net/>
7 *
8 * This program is free software; you can redistribute it and/or
9 * modify it under the terms of the GNU General Public License
10 * as published by the Free Software Foundation; either version 2
11 * of the License, or (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the GNU General Public
19 * License along with this program; if not, write to the Free
20 * Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
21 * Boston, MA 02110-1301, USA.
22*/
23
24#ifndef EUKLID_NORM_2D_HH
25#define EUKLID_NORM_2D_HH
26
27#include <vector>
28#include "utilities/omath.h"
29#include <limits>
30
31#include "euklidNorm2D.h"
32#include "dynamics/lbm.h" // for computation of lattice rho and velocity
35#include "blockBaseF2D.h"
36#include "functors/genericF.h"
40
41
42namespace olb {
43
44template<typename T,typename DESCRIPTOR>
47 : SuperLatticeF2D<T,DESCRIPTOR>(f.getSuperLattice(), 1), _f(f)
48{
49 this->getName() = "l2(" + _f.getName() + ")";
50 int maxC = this->_sLattice.getLoadBalancer().size();
51 this->_blockF.reserve(maxC);
52 for (int iC = 0; iC < maxC; iC++) {
53 this->_blockF.emplace_back(new BlockEuklidNorm2D<T,DESCRIPTOR>(f.getBlockF(iC)));
54 }
55}
56
57template <typename T, typename DESCRIPTOR>
59 : BlockF2D<T>(f.getBlockStructure(),1), _f(f)
60{
61 this->getName() = "l2("+f.getName()+")";
62}
63
64template <typename T, typename DESCRIPTOR>
65bool BlockEuklidNorm2D<T,DESCRIPTOR>::operator() (T output[], const int input[])
66{
67 output[0] = T(); // flash output, since this methods adds values.
68 T data[_f.getTargetDim()];
69 _f(data,input);
70 for ( int i = 0; i < _f.getTargetDim(); ++i) {
71 output[0] += data[i]*data[i];
72 }
73 output[0] = util::sqrt(output[0]);
74 return true;
75}
76
77}
78#endif
BlockL2Norm2D returns pointwise the l2-norm, e.g. of a velocity.
BlockEuklidNorm2D(BlockF2D< T > &f)
bool operator()(T output[], const int input[]) override
has to be implemented for 'every' derived class
represents all functors that operate on a cuboid in general, mother class of BlockLatticeF,...
std::string & getName()
read and write access to name
Definition genericF.hh:51
SuperEuklidNorm2D(SuperLatticeF2D< T, DESCRIPTOR > &f)
std::vector< std::unique_ptr< BlockF2D< T > > > _blockF
Super functors may consist of several BlockF2D<W> derived functors.
BlockF2D< W > & getBlockF(int iCloc)
represents all functors that operate on a SuperLattice in general, e.g. getVelocity(),...
SuperLattice< T, DESCRIPTOR > & _sLattice
The description of a generic interface for all functor classes – header file.
This file contains indicator functions.
Wrapper functions that simplify the use of MPI.
cpu::simd::Pack< T > sqrt(cpu::simd::Pack< T > value)
Definition pack.h:100
Top level namespace for all of OpenLB.
Representation of a parallel 2D geometry – header file.