12#include <Kokkos_Core.hpp>
20template<
typename ValueType>
23 const std::shared_ptr<GeometryScheme> geometryScheme,
28 const auto& exec = surfaceField.
exec();
30 const auto [owner, neighbour, surfFaceCells] =
34 const auto [phif, phi, phiBCValue, nonOrthDeltaCoeffs] =
spans(
38 geometryScheme->nonOrthDeltaCoeffs().internalField()
46 KOKKOS_LAMBDA(
const size_t facei) {
47 phif[facei] = nonOrthDeltaCoeffs[facei] * (phi[neighbour[facei]] - phi[owner[facei]]);
53 {nInternalFaces, phif.size()},
54 KOKKOS_LAMBDA(
const size_t facei) {
55 auto faceBCI = facei - nInternalFaces;
56 auto own =
static_cast<size_t>(surfFaceCells[faceBCI]);
58 phif[facei] = nonOrthDeltaCoeffs[facei] * (phiBCValue[faceBCI] - phi[own]);
63template<
typename ValueType>
73 : Base(exec, mesh), geometryScheme_(
GeometryScheme::readOrCreate(mesh)) {};
76 : Base(exec, mesh), geometryScheme_(
GeometryScheme::readOrCreate(mesh)) {};
78 static std::string
name() {
return "uncorrected"; }
80 static std::string
doc() {
return "Uncorrected interpolation"; }
82 static std::string
schema() {
return "none"; }
93 return geometryScheme_->nonOrthDeltaCoeffs();
96 std::unique_ptr<FaceNormalGradientFactory<ValueType>>
clone()
const override
98 return std::make_unique<Uncorrected>(*
this);
103 const std::shared_ptr<GeometryScheme> geometryScheme_;
const NeoFOAM::Field< T > & value() const
Get the view storing the computed values from the boundary condition.
const labelField & faceCells() const
Get the field of face cells.
Represents an unstructured mesh in NeoFOAM.
const BoundaryMesh & boundaryMesh() const
Get the boundary mesh.
size_t nInternalFaces() const
Get the number of internal faces in the mesh.
const labelField & faceOwner() const
Get the field of face owner cells.
const labelField & faceNeighbour() const
Get the field of face neighbour cells.
const UnstructuredMesh & mesh() const
Returns a const reference to the unstructured mesh object.
const Executor & exec() const
Returns a const reference to the executor object.
const Field< ValueType > & internalField() const
Returns a const reference to the internal field.
const BoundaryFields< ValueType > & boundaryField() const
Returns a const reference to the boundary field.
Represents a surface field in a finite volume method.
Uncorrected(const Executor &exec, const UnstructuredMesh &mesh, Input input)
Uncorrected(const Executor &exec, const UnstructuredMesh &mesh)
static std::string schema()
virtual void faceNormalGrad(const VolumeField< ValueType > &volField, SurfaceField< ValueType > &surfaceField) const override
std::unique_ptr< FaceNormalGradientFactory< ValueType > > clone() const override
static std::string name()
virtual const SurfaceField< scalar > & deltaCoeffs() const override
Represents a volume field in a finite volume method.
A template class for registering derived classes with a base class.
void computeFaceNormalGrad(const VolumeField< ValueType > &volField, const std::shared_ptr< GeometryScheme > geometryScheme, SurfaceField< ValueType > &surfaceField)
std::variant< Dictionary, TokenList > Input
void parallelFor(const Executor &exec, std::pair< size_t, size_t > range, Kernel kernel)
std::variant< SerialExecutor, CPUExecutor, GPUExecutor > Executor
auto spans(Args &... fields)