/* This file is part of the OpenLB library
*
* Copyright (C) 2013 Patrick Nathen, Mathias J. Krause
* E-mail contact: info@openlb.net
* The most recent release of OpenLB can be downloaded at
* <http://www.openlb.net/>
*
* This program 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 2
* of the License, or (at your option) any later version.
*
* This program 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 this program; if not, write to the Free
* Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
* Boston, MA 02110-1301, USA.
*/
#ifndef TURBULENT_F_3D_HH
#define TURBULENT_F_3D_HH
#include<vector>
#include<cmath>
#include "turbulentF3D.h"
#include "blockLatticeLocalF3D.h"
#include "dynamics/smagorinskyBGKdynamics.h"
#include "core/superLattice3D.h"
#include "core/finiteDifference.h"
#include "geometry/superGeometry3D.h"
#include "utilities/vectorHelpers.h"
#include "dynamics/lbHelpers.h" // for computation of lattice rho and velocity
namespace olb {
///////////////////////////// SuperLatticeYplus3D //////////////////////////////
template <typename T, typename DESCRIPTOR>
SuperLatticeYplus3D<T,DESCRIPTOR>::SuperLatticeYplus3D(SuperLattice3D<T,DESCRIPTOR>& sLattice,
const UnitConverter<T,DESCRIPTOR>& converter, SuperGeometry3D<T>& superGeometry,
IndicatorF3D<T>& indicator, const int material )
: SuperLatticePhysF3D<T,DESCRIPTOR>(sLattice,converter,1),
_superGeometry(superGeometry), _indicator(indicator), _material(material)
{
this->getName() = "yPlus";
}
template <typename T, typename DESCRIPTOR>
bool SuperLatticeYplus3D<T,DESCRIPTOR>::operator() (T output[], const int input[])
{
int globIC = input[0];
int locix = input[1];
int lociy = input[2];
int lociz = input[3];
output[0]=T();
if 