LatticeModel.tmpl.h 18.4 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
//======================================================================================================================
//
//  This file is part of waLBerla. waLBerla is free software: you can
//  redistribute it and/or modify it under the terms of the GNU General Public
//  License as published by the Free Software Foundation, either version 3 of
//  the License, or (at your option) any later version.
//
//  waLBerla is distributed in the hope that it will be useful, but WITHOUT
//  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
//  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
//  for more details.
//
//  You should have received a copy of the GNU General Public License along
//  with waLBerla (see COPYING.txt). If not, see <http://www.gnu.org/licenses/>.
//
//! \\author Martin Bauer <martin.bauer@fau.de>
//
//======================================================================================================================


#include "core/DataTypes.h"
#include "core/logging/Logging.h"

#include "field/GhostLayerField.h"
#include "field/SwapableCompare.h"
#include "domain_decomposition/BlockDataID.h"
#include "domain_decomposition/IBlock.h"
Martin Bauer's avatar
Martin Bauer committed
28
#include "stencil/{{stencil_name}}.h"
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52

#include "lbm/lattice_model/EquilibriumDistribution.h"
#include "lbm/field/Density.h"
#include "lbm/field/DensityAndMomentumDensity.h"
#include "lbm/field/DensityAndVelocity.h"
#include "lbm/field/PressureTensor.h"
#include "lbm/field/ShearRate.h"

#include <set>

#ifdef __GNUC__
#define RESTRICT __restrict__
#elif _MSC_VER
#define RESTRICT __restrict
#else
#define RESTRICT
#endif

#ifdef WALBERLA_CXX_COMPILER_IS_GNU
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wunused-variable"
#pragma GCC diagnostic ignored "-Wunused-parameter"
#endif

53
54
55
56
57
#ifdef WALBERLA_CXX_COMPILER_IS_CLANG
#pragma clang diagnostic push
#pragma clang diagnostic ignored "-Wunused-variable"
#pragma clang diagnostic ignored "-Wunused-parameter"
#endif
58

59
{% set lmIgnores = ('pdfs', 'pdfs_tmp') %}
60
{% set lmOffsets = ('block_offset_0', 'block_offset_1', 'block_offset_2') %}
61
62
63
64
65


// Forward declarations
namespace walberla{
namespace {{namespace}} {
Martin Bauer's avatar
Martin Bauer committed
66
   class {{class_name}};
67
68
69
}}
namespace walberla {
namespace mpi {
Martin Bauer's avatar
Martin Bauer committed
70
71
    mpi::SendBuffer & operator<< (mpi::SendBuffer & buf, const ::walberla::{{namespace}}::{{class_name}} & lm);
    mpi::RecvBuffer & operator>> (mpi::RecvBuffer & buf,       ::walberla::{{namespace}}::{{class_name}} & lm);
72
73
74
75
76
77
78
79
80
81
}}




namespace walberla {
namespace {{namespace}} {


/**
Martin Bauer's avatar
Martin Bauer committed
82
{{class_name}} was generated with lbmpy. Do not edit this file directly. Instead modify {{class_name}}.py.
83
84
85
86
87
88
89
90
For details see documentation of lbmpy.

Usage:
    - Create an instance of this lattice model class: the constructor parameters vary depending on the configure
      lattice model. A model with constant force needs a single force vector, while a model with variable forces needs
      a force field. All constructor parameters are ordered alphabetically.
    - Create a PDFField with the lattice model as template argument to store the particle distribution functions.
      Use the PDFField to get and modify macroscopic values.
Martin Bauer's avatar
Martin Bauer committed
91
    - The internal class {{class_name}}::Sweep is a functor to execute one LB time step.
92
93
94
      Stream, collide steps can be executed separately, or together in an optimized stream-pull-collide scheme

*/
Martin Bauer's avatar
Martin Bauer committed
95
class {{class_name}}
96
97
98
{

public:
Martin Bauer's avatar
Martin Bauer committed
99
100
    typedef stencil::{{stencil_name}} Stencil;
    typedef stencil::{{stencil_name}} CommunicationStencil;
101
102
103
    static const real_t w[{{Q}}];
    static const real_t wInv[{{Q}}];

104
    static const bool compressible = {% if compressible %}true{% else %}false{% endif %};
105
106
107
108
109
110
111
112
113

    class Sweep
    {
    public:
        Sweep( BlockDataID _pdfsID ) : pdfsID(_pdfsID) {};

        //void stream       ( IBlock * const block, const uint_t numberOfGhostLayersToInclude = uint_t(0) );
        void collide      ( IBlock * const block, const uint_t numberOfGhostLayersToInclude = uint_t(0) );
        void streamCollide( IBlock * const block, const uint_t numberOfGhostLayersToInclude = uint_t(0) );
114
        void stream       ( IBlock * const block, const uint_t numberOfGhostLayersToInclude = uint_t(0) );
115
116
117
118
119
120
121

        void operator() ( IBlock * const block, const uint_t numberOfGhostLayersToInclude = uint_t(0) )
        {
            streamCollide( block, numberOfGhostLayersToInclude );
        }

    private:
122
        {{stream_collide_kernel|generate_members(only_fields=True)|indent(8)}}
123
124
    };

125
126
    {{class_name}}( {{stream_collide_kernel|generate_constructor_parameters(lmIgnores+lmOffsets) }} )
        : {{ stream_collide_kernel|generate_constructor_initializer_list(lmIgnores+lmOffsets) }}, currentLevel(0)
127
128
    {};

129
    void configure( IBlock & block, StructuredBlockStorage & storage )  { configureBlock( &block, &storage ); }
130

131
132
133
    // Parameters:
    {{stream_collide_kernel|generate_members(lmIgnores)|indent(4)}}

134
private:
135
    void configureBlock(IBlock * block, StructuredBlockStorage * storage)
136
    {
137
138
139
140
141
        {{stream_collide_kernel|generate_block_data_to_field_extraction(lmIgnores+lmOffsets, no_declarations=True)|indent(8)}}
        block_offset_0 = uint32_c(storage->getBlockCellBB(*block).xMin());
        block_offset_1 = uint32_c(storage->getBlockCellBB(*block).yMin());
        block_offset_2 = uint32_c(storage->getBlockCellBB(*block).zMin());
        blockId = &block->getId();
142

Martin Bauer's avatar
Martin Bauer committed
143
        {% if refinement_scaling -%}
144
145
146
147
        const uint_t targetLevel = block->getBlockStorage().getLevel(*block);

        if( targetLevel != currentLevel )
        {
Jan Hönig's avatar
Jan Hönig committed
148
            real_t level_scale_factor(1);
149
            if( currentLevel < targetLevel )
Jan Hönig's avatar
Jan Hönig committed
150
               level_scale_factor = real_c( uint_t(1) << ( targetLevel - currentLevel ) );
151
            else // currentLevel > targetLevel
Jan Hönig's avatar
Jan Hönig committed
152
               level_scale_factor = real_t(1) / real_c( uint_t(1) << ( currentLevel - targetLevel ) );
153

Martin Bauer's avatar
Martin Bauer committed
154
            {% for scalingType, name, expression in refinement_scaling.scaling_info -%}
155
156
            {% if scalingType == 'normal' %}
            {{name}} = {{expression}};
Martin Bauer's avatar
Martin Bauer committed
157
158
            {% elif scalingType in ('field_with_f', 'field_xyz') %}
            auto it = {{name}}->{% if scalingType == 'field_with_f'%} beginWithGhostLayer(){% else %}beginWithGhostLayerXYZ(){% endif %};
159
160
161
162
163
            for( ; it != {{name}}->end(); ++it )
            {
                 auto x = it.x();
                 auto y = it.y();
                 auto z = it.z();
Martin Bauer's avatar
Martin Bauer committed
164
                 {% if scalingType == 'field_with_f' -%}
165
166
167
168
169
170
171
172
173
174
175
176
177
                 auto f = it.f();
                 {% endif -%}
                 *it = {{expression}};
            }
            {% endif -%}
            {% endfor -%}

            currentLevel = targetLevel;
        }
        {% endif -%}
    }

    // Updated by configureBlock:
Martin Bauer's avatar
Martin Bauer committed
178
    {{stream_collide_kernel|generate_block_data_to_field_extraction(lmIgnores, declarations_only=True)|indent(4)}}
179
    const IBlockID * blockId;
180
181
182
    uint_t currentLevel;

    // Backend classes can access private members:
Martin Bauer's avatar
Martin Bauer committed
183
    friend class {{class_name}}::Sweep;
184
185
186
187
188
189
190
191
    template<class LM, class Enable> friend class  EquilibriumDistribution;
    template<class LM, class Enable> friend struct Equilibrium;
    template<class LM, class Enable> friend struct internal::AdaptVelocityToForce;
    template<class LM, class Enable> friend struct Density;
    template<class LM>               friend struct DensityAndVelocity;
    template<class LM, class Enable> friend struct DensityAndMomentumDensity;
    template<class LM, class Enable> friend struct MomentumDensity;
    template<class LM, class It, class Enable> friend struct DensityAndVelocityRange;
192

Martin Bauer's avatar
Martin Bauer committed
193
194
    friend mpi::SendBuffer & ::walberla::mpi::operator<< (mpi::SendBuffer & , const {{class_name}} & );
    friend mpi::RecvBuffer & ::walberla::mpi::operator>> (mpi::RecvBuffer & ,       {{class_name}} & );
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209

};




//======================================================================================================================
//
//  Implementation of macroscopic value backend
//
//======================================================================================================================



template<>
Martin Bauer's avatar
Martin Bauer committed
210
class EquilibriumDistribution< {{class_name}}, void>
211
{
212
public:
Martin Bauer's avatar
Martin Bauer committed
213
   typedef typename {{class_name}}::Stencil Stencil;
214
215
216
217
218

   static real_t get( const stencil::Direction direction,
                      const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ),
                      real_t rho = real_t(1.0) )
   {
219
        {% if not compressible %}
220
221
        rho -= real_t(1.0);
        {% endif %}
Martin Bauer's avatar
Martin Bauer committed
222
        {{equilibrium_from_direction}}
223
224
225
226
227
228
   }

   static real_t getSymmetricPart( const stencil::Direction direction,
                                   const Vector3<real_t> & u = Vector3< real_t >(real_t(0.0)),
                                   real_t rho = real_t(1.0) )
   {
229
        {% if not compressible %}
230
231
        rho -= real_t(1.0);
        {% endif %}
232
        {{symmetric_equilibrium_from_direction}}
233
234
235
236
237
238
   }

   static real_t getAsymmetricPart( const stencil::Direction direction,
                                    const Vector3< real_t > & u = Vector3<real_t>( real_t(0.0) ),
                                    real_t rho = real_t(1.0) )
   {
239
        {% if not compressible %}
240
241
        rho -= real_t(1.0);
        {% endif %}
Martin Bauer's avatar
Martin Bauer committed
242
        {{asymmetric_equilibrium_from_direction}}
243
244
245
246
247
   }

   static std::vector< real_t > get( const Vector3< real_t > & u = Vector3<real_t>( real_t(0.0) ),
                                     real_t rho = real_t(1.0) )
   {
248
      {% if not compressible %}
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
      rho -= real_t(1.0);
      {% endif %}

      std::vector< real_t > equilibrium( Stencil::Size );
      for( auto d = Stencil::begin(); d != Stencil::end(); ++d )
      {
         equilibrium[d.toIdx()] = get(*d, u, rho);
      }
      return equilibrium;
   }
};


namespace internal {

template<>
Martin Bauer's avatar
Martin Bauer committed
265
struct AdaptVelocityToForce<{{class_name}}, void>
266
267
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
268
   static Vector3<real_t> get( FieldPtrOrIterator & it, const {{class_name}} & lm,
269
270
271
272
273
                               const Vector3< real_t > & velocity, const real_t rho )
   {
      auto x = it.x();
      auto y = it.y();
      auto z = it.z();
Martin Bauer's avatar
Martin Bauer committed
274
275
      {% if macroscopic_velocity_shift %}
      return velocity - Vector3<real_t>({{macroscopic_velocity_shift | join(",") }} {% if D == 2 %}, 0.0 {%endif %} );
276
277
278
279
280
      {% else %}
      return velocity;
      {% endif %}
   }

Martin Bauer's avatar
Martin Bauer committed
281
   static Vector3<real_t> get( const cell_idx_t x, const cell_idx_t y, const cell_idx_t z, const {{class_name}} & lm,
282
283
                               const Vector3< real_t > & velocity, const real_t rho )
   {
Martin Bauer's avatar
Martin Bauer committed
284
      {% if macroscopic_velocity_shift %}
285

Martin Bauer's avatar
Martin Bauer committed
286
      return velocity - Vector3<real_t>({{macroscopic_velocity_shift | join(",") }} {% if D == 2 %}, 0.0 {%endif %} );
287
288
289
290
291
292
293
294
295
296
      {% else %}
      return velocity;
      {% endif %}
   }
};
} // namespace internal



template<>
Martin Bauer's avatar
Martin Bauer committed
297
struct Equilibrium< {{class_name}}, void >
298
299
300
301
{

   template< typename FieldPtrOrIterator >
   static void set( FieldPtrOrIterator & it,
302
                    const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ), real_t rho = real_t(1.0) )
303
   {
304
        {%if not compressible %}
305
306
307
        rho -= real_t(1.0);
        {%endif %}

308
       {% for eqTerm in equilibrium -%}
309
       it[{{loop.index0 }}] = {{eqTerm}};
310
311
312
313
314
       {% endfor -%}
   }

   template< typename PdfField_T >
   static void set( PdfField_T & pdf, const cell_idx_t x, const cell_idx_t y, const cell_idx_t z,
315
                    const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ), real_t rho = real_t(1.0) )
316
   {
317
      {%if not compressible %}
318
319
320
      rho -= real_t(1.0);
      {%endif %}

321
322
      real_t & xyz0 = pdf(x,y,z,0);
      {% for eqTerm in equilibrium -%}
323
      pdf.getF( &xyz0, {{loop.index0 }})= {{eqTerm}};
324
325
326
327
328
329
      {% endfor -%}
   }
};


template<>
Martin Bauer's avatar
Martin Bauer committed
330
struct Density<{{class_name}}, void>
331
332
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
333
   static inline real_t get( const {{class_name}} & , const FieldPtrOrIterator & it )
334
335
336
337
   {
        {% for i in range(Q) -%}
            const real_t f_{{i}} = it[{{i}}];
        {% endfor -%}
Martin Bauer's avatar
Martin Bauer committed
338
        {{density_getters | indent(8)}}
339
340
341
342
        return rho;
   }

   template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
343
   static inline real_t get( const {{class_name}} & ,
344
345
346
347
348
349
                             const PdfField_T & pdf, const cell_idx_t x, const cell_idx_t y, const cell_idx_t z )
   {
        const real_t & xyz0 = pdf(x,y,z,0);
        {% for i in range(Q) -%}
            const real_t f_{{i}} = pdf.getF( &xyz0, {{i}});
        {% endfor -%}
Martin Bauer's avatar
Martin Bauer committed
350
        {{density_getters | indent(8)}}
351
352
353
354
355
356
        return rho;
   }
};


template<>
Martin Bauer's avatar
Martin Bauer committed
357
struct DensityAndVelocity<{{class_name}}>
358
359
{
    template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
360
    static void set( FieldPtrOrIterator & it, const {{class_name}} & lm,
Jan Hönig's avatar
Jan Hönig committed
361
                     const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ), const real_t rho_in = real_t(1.0) )
362
363
364
365
366
    {
        auto x = it.x();
        auto y = it.y();
        auto z = it.z();

Martin Bauer's avatar
Martin Bauer committed
367
        {{density_velocity_setter_macroscopic_values | indent(8)}}
368
369
370
371
        {% if D == 2 -%}
        const real_t u_2(0.0);
        {% endif %}

372
        Equilibrium<{{class_name}}>::set(it, Vector3<real_t>(u_0, u_1, u_2), rho{%if not compressible %} + real_t(1) {%endif%});
373
374
375
    }

    template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
376
    static void set( PdfField_T & pdf, const cell_idx_t x, const cell_idx_t y, const cell_idx_t z, const {{class_name}} & lm,
Jan Hönig's avatar
Jan Hönig committed
377
                     const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ), const real_t rho_in = real_t(1.0) )
378
    {
Martin Bauer's avatar
Martin Bauer committed
379
        {{density_velocity_setter_macroscopic_values | indent(8)}}
380
381
382
383
        {% if D == 2 -%}
        const real_t u_2(0.0);
        {% endif %}

384
        Equilibrium<{{class_name}}>::set(pdf, x, y, z, Vector3<real_t>(u_0, u_1, u_2), rho {%if not compressible %} + real_t(1) {%endif%});
385
386
387
388
389
    }
};


template<typename FieldIteratorXYZ >
Martin Bauer's avatar
Martin Bauer committed
390
struct DensityAndVelocityRange<{{class_name}}, FieldIteratorXYZ>
391
392
{

Martin Bauer's avatar
Martin Bauer committed
393
   static void set( FieldIteratorXYZ & begin, const FieldIteratorXYZ & end, const {{class_name}} & lm,
Jan Hönig's avatar
Jan Hönig committed
394
                    const Vector3< real_t > & u = Vector3< real_t >( real_t(0.0) ), const real_t rho_in = real_t(1.0) )
395
396
397
398
399
400
   {
        for( auto cellIt = begin; cellIt != end; ++cellIt )
        {
            const auto x = cellIt.x();
            const auto y = cellIt.y();
            const auto z = cellIt.z();
Martin Bauer's avatar
Martin Bauer committed
401
            {{density_velocity_setter_macroscopic_values | indent(12)}}
402
403
404
405
            {% if D == 2 -%}
            const real_t u_2(0.0);
            {% endif %}

406
            Equilibrium<{{class_name}}>::set(cellIt, Vector3<real_t>(u_0, u_1, u_2), rho{%if not compressible %} + real_t(1) {%endif%});
407
408
409
410
411
412
413
        }
   }
};



template<>
Martin Bauer's avatar
Martin Bauer committed
414
struct DensityAndMomentumDensity<{{class_name}}>
415
416
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
417
   static real_t get( Vector3< real_t > & momentumDensity, const {{class_name}} & lm,
418
419
420
421
422
423
424
425
426
427
                      const FieldPtrOrIterator & it )
   {
        const auto x = it.x();
        const auto y = it.y();
        const auto z = it.z();

        {% for i in range(Q) -%}
            const real_t f_{{i}} = it[{{i}}];
        {% endfor -%}

Martin Bauer's avatar
Martin Bauer committed
428
        {{momentum_density_getter | indent(8) }}
429
430
431
432
433
434
435
        {% for i in range(D) -%}
            momentumDensity[{{i}}] = md_{{i}};
        {% endfor %}
        return rho;
   }

   template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
436
   static real_t get( Vector3< real_t > & momentumDensity, const {{class_name}} & lm, const PdfField_T & pdf,
437
438
439
440
441
442
443
                      const cell_idx_t x, const cell_idx_t y, const cell_idx_t z )
   {
        const real_t & xyz0 = pdf(x,y,z,0);
        {% for i in range(Q) -%}
            const real_t f_{{i}} = pdf.getF( &xyz0, {{i}});
        {% endfor -%}

Martin Bauer's avatar
Martin Bauer committed
444
        {{momentum_density_getter | indent(8) }}
445
446
447
448
449
450
451
452
453
        {% for i in range(D) -%}
            momentumDensity[{{i}}] = md_{{i}};
        {% endfor %}
       return rho;
   }
};


template<>
Martin Bauer's avatar
Martin Bauer committed
454
struct MomentumDensity< {{class_name}}>
455
456
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
457
   static void get( Vector3< real_t > & momentumDensity, const {{class_name}} & lm, const FieldPtrOrIterator & it )
458
459
460
461
462
463
464
465
466
   {
        const auto x = it.x();
        const auto y = it.y();
        const auto z = it.z();

        {% for i in range(Q) -%}
            const real_t f_{{i}} = it[{{i}}];
        {% endfor -%}

Martin Bauer's avatar
Martin Bauer committed
467
        {{momentum_density_getter | indent(8) }}
468
469
470
471
472
473
        {% for i in range(D) -%}
            momentumDensity[{{i}}] = md_{{i}};
        {% endfor %}
   }

   template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
474
   static void get( Vector3< real_t > & momentumDensity, const {{class_name}} & lm, const PdfField_T & pdf,
475
476
477
478
479
480
481
                    const cell_idx_t x, const cell_idx_t y, const cell_idx_t z )
   {
        const real_t & xyz0 = pdf(x,y,z,0);
        {% for i in range(Q) -%}
            const real_t f_{{i}} = pdf.getF( &xyz0, {{i}});
        {% endfor -%}

Martin Bauer's avatar
Martin Bauer committed
482
        {{momentum_density_getter | indent(8) }}
483
484
485
486
487
488
489
490
        {% for i in range(D) -%}
            momentumDensity[{{i}}] = md_{{i}};
        {% endfor %}
   }
};


template<>
Martin Bauer's avatar
Martin Bauer committed
491
struct PressureTensor<{{class_name}}>
492
493
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
494
   static void get( Matrix3< real_t > & /* pressureTensor */, const {{class_name}} & /* latticeModel */, const FieldPtrOrIterator & /* it */ )
495
496
497
498
499
   {
       WALBERLA_ABORT("Not implemented");
   }

   template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
500
   static void get( Matrix3< real_t > & /* pressureTensor */, const {{class_name}} & /* latticeModel */, const PdfField_T & /* pdf */,
501
                    const cell_idx_t /* x */, const cell_idx_t /* y */, const cell_idx_t /* z */ )
502
503
504
505
506
507
508
   {
       WALBERLA_ABORT("Not implemented");
   }
};


template<>
Martin Bauer's avatar
Martin Bauer committed
509
struct ShearRate<{{class_name}}>
510
511
{
   template< typename FieldPtrOrIterator >
Martin Bauer's avatar
Martin Bauer committed
512
   static inline real_t get( const {{class_name}} & /* latticeModel */, const FieldPtrOrIterator & /* it */,
513
                             const Vector3< real_t > & /* velocity */, const real_t /* rho */)
514
515
516
517
518
519
   {
       WALBERLA_ABORT("Not implemented");
       return real_t(0.0);
   }

   template< typename PdfField_T >
Martin Bauer's avatar
Martin Bauer committed
520
   static inline real_t get( const {{class_name}} & latticeModel,
521
522
                             const PdfField_T & /* pdf */, const cell_idx_t /* x */, const cell_idx_t /* y */, const cell_idx_t /* z */,
                             const Vector3< real_t > & /* velocity */, const real_t /* rho */ )
523
524
525
526
527
   {
       WALBERLA_ABORT("Not implemented");
       return real_t(0.0);
   }

528
529
   static inline real_t get( const std::vector< real_t > & /* nonEquilibrium */, const real_t /* relaxationParam */,
                             const real_t /* rho */ = real_t(1) )
530
531
532
533
534
535
536
537
538
539
540
541
542
543
   {
       WALBERLA_ABORT("Not implemented");
       return real_t(0.0);
   }
};


} // namespace {{namespace}}
} // namespace walberla



#ifdef WALBERLA_CXX_COMPILER_IS_GNU
#pragma GCC diagnostic pop
544
545
546
547
548
#endif

#ifdef WALBERLA_CXX_COMPILER_IS_CLANG
#pragma clang diagnostic pop
#endif