-
Nils Kohl authored* HFG updated to 27578ea8116ca0920c88bc0da76f366448559a37 to support gradient. * Generated new operators (diffusion (P2), gradient, divergence) * Regenerated epsilon, full_stokes (with and without icoshell blending) after fixing trial and test swap and component passing.
Nils Kohl authored* HFG updated to 27578ea8116ca0920c88bc0da76f366448559a37 to support gradient. * Generated new operators (diffusion (P2), gradient, divergence) * Regenerated epsilon, full_stokes (with and without icoshell blending) after fixing trial and test swap and component passing.
P2ElementwiseDiffusion.hpp 13.30 KiB
/*
* Copyright (c) 2017-2024 Nils Kohl, Daniel Bauer, Fabian Böhm.
*
* This file is part of HyTeG
* (see https://i10git.cs.fau.de/hyteg/hyteg).
*
* 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 3 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, see <http://www.gnu.org/licenses/>.
*/
/*
* The entire file was generated with the HyTeG form generator.
*
* Avoid modifying this file. If buggy, consider fixing the generator itself.
*/
#pragma once
#include "core/DataTypes.h"
#include "hyteg/LikwidWrapper.hpp"
#include "hyteg/communication/Syncing.hpp"
#include "hyteg/edgedofspace/EdgeDoFMacroCell.hpp"
#include "hyteg/operators/Operator.hpp"
#include "hyteg/p2functionspace/P2Function.hpp"
#include "hyteg/primitivestorage/PrimitiveStorage.hpp"
#include "hyteg/solvers/Smoothables.hpp"
#include "hyteg/sparseassembly/SparseMatrixProxy.hpp"
#define FUNC_PREFIX
namespace hyteg {
namespace operatorgeneration {
/// Diffusion operator without coefficients.
///
/// Geometry map: IdentityMap
///
/// Weak formulation
///
/// u: trial function (space: Lagrange, degree: 2)
/// v: test function (space: Lagrange, degree: 2)
///
/// ∫ ∇u · ∇v
class P2ElementwiseDiffusion : public Operator< P2Function< walberla::float64 >, P2Function< walberla::float64 > >,
public OperatorWithInverseDiagonal< P2Function< walberla::float64 > >
{
public:
P2ElementwiseDiffusion( const std::shared_ptr< PrimitiveStorage >& storage, size_t minLevel, size_t maxLevel );
void apply( const P2Function< walberla::float64 >& src,
const P2Function< walberla::float64 >& dst,
uint_t level,
DoFType flag,
UpdateType updateType = Replace ) const;
void toMatrix( const std::shared_ptr< SparseMatrixProxy >& mat,
const P2Function< idx_t >& src, const P2Function< idx_t >& dst,
uint_t level,
DoFType flag ) const;
void computeInverseDiagonalOperatorValues();
std::shared_ptr< P2Function< walberla::float64 > > getInverseDiagonalValues() const;
protected:
private:
/// Kernel type: apply
/// - quadrature rule: Dunavant 2 | points: 3, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 183 286 0 0 0 0 0 0
void apply_macro_2D( walberla::float64* RESTRICT _data_dstEdge,
walberla::float64* RESTRICT _data_dstVertex,
walberla::float64* RESTRICT _data_srcEdge,
walberla::float64* RESTRICT _data_srcVertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
int64_t micro_edges_per_macro_edge,
walberla::float64 micro_edges_per_macro_edge_float ) const;
/// Kernel type: apply
/// - quadrature rule: Hammer-Marlowe-Stroud 1 | points: 4, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 993 1389 0 0 0 0 0 0
void apply_macro_3D( walberla::float64* RESTRICT _data_dstEdge,
walberla::float64* RESTRICT _data_dstVertex,
walberla::float64* RESTRICT _data_srcEdge,
walberla::float64* RESTRICT _data_srcVertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_0comp2,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_1comp2,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
walberla::float64 macro_vertex_coord_id_2comp2,
walberla::float64 macro_vertex_coord_id_3comp0,
walberla::float64 macro_vertex_coord_id_3comp1,
walberla::float64 macro_vertex_coord_id_3comp2,
int64_t micro_edges_per_macro_edge,
walberla::float64 micro_edges_per_macro_edge_float ) const;
/// Kernel type: toMatrix
/// - quadrature rule: Dunavant 2 | points: 3, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 147 250 0 0 0 0 0 3
void toMatrix_macro_2D( idx_t* RESTRICT _data_dstEdge,
idx_t* RESTRICT _data_dstVertex,
idx_t* RESTRICT _data_srcEdge,
idx_t* RESTRICT _data_srcVertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
std::shared_ptr< SparseMatrixProxy > mat,
int64_t micro_edges_per_macro_edge, walberla::float64 micro_edges_per_macro_edge_float ) const;
/// Kernel type: toMatrix
/// - quadrature rule: Hammer-Marlowe-Stroud 1 | points: 4, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 893 1289 0 0 0 0 0 3
void toMatrix_macro_3D( idx_t* RESTRICT _data_dstEdge,
idx_t* RESTRICT _data_dstVertex,
idx_t* RESTRICT _data_srcEdge,
idx_t* RESTRICT _data_srcVertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_0comp2,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_1comp2,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
walberla::float64 macro_vertex_coord_id_2comp2,
walberla::float64 macro_vertex_coord_id_3comp0,
walberla::float64 macro_vertex_coord_id_3comp1,
walberla::float64 macro_vertex_coord_id_3comp2,
std::shared_ptr< SparseMatrixProxy > mat,
int64_t micro_edges_per_macro_edge,
walberla::float64 micro_edges_per_macro_edge_float ) const;
/// Kernel type: computeInverseDiagonalOperatorValues
/// - quadrature rule: Dunavant 2 | points: 3, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 78 103 0 0 0 0 0 0
void computeInverseDiagonalOperatorValues_macro_2D( walberla::float64* RESTRICT _data_invDiag_Edge,
walberla::float64* RESTRICT _data_invDiag_Vertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
int64_t micro_edges_per_macro_edge,
walberla::float64 micro_edges_per_macro_edge_float ) const;
/// Kernel type: computeInverseDiagonalOperatorValues
/// - quadrature rule: Hammer-Marlowe-Stroud 1 | points: 4, degree: 2
/// - operations per element:
/// adds muls divs pows abs assignments function_calls unknown_ops
/// ------ ------ ------ ------ ----- ------------- ---------------- -------------
/// 288 425 0 0 0 0 0 0
void computeInverseDiagonalOperatorValues_macro_3D( walberla::float64* RESTRICT _data_invDiag_Edge,
walberla::float64* RESTRICT _data_invDiag_Vertex,
walberla::float64 macro_vertex_coord_id_0comp0,
walberla::float64 macro_vertex_coord_id_0comp1,
walberla::float64 macro_vertex_coord_id_0comp2,
walberla::float64 macro_vertex_coord_id_1comp0,
walberla::float64 macro_vertex_coord_id_1comp1,
walberla::float64 macro_vertex_coord_id_1comp2,
walberla::float64 macro_vertex_coord_id_2comp0,
walberla::float64 macro_vertex_coord_id_2comp1,
walberla::float64 macro_vertex_coord_id_2comp2,
walberla::float64 macro_vertex_coord_id_3comp0,
walberla::float64 macro_vertex_coord_id_3comp1,
walberla::float64 macro_vertex_coord_id_3comp2,
int64_t micro_edges_per_macro_edge,
walberla::float64 micro_edges_per_macro_edge_float ) const;
std::shared_ptr< P2Function< walberla::float64 > > invDiag_;
};
} // namespace operatorgeneration
} // namespace hyteg