fefactory.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2021-2024 The Ikarus Developers mueller@ibb.uni-stuttgart.de
// SPDX-License-Identifier: LGPL-3.0-or-later
 
#pragma once
#include <ikarus/finiteelements/febase.hh>
namespace Ikarus {
template <typename BH, typename SK, bool useFlat = true, bool useEigenRef = false>
struct FEFactory
{
  using Skills = SK;
 
private:
  const BH* basisHandler_;
  SK skills_;
 
public:
  template <typename SK2 = SK>
  FEFactory(const BH& basisHandler, const SK2& sk)
      : basisHandler_{&basisHandler},
        skills_{sk} {}
 
  auto operator()() {
    return std::apply(
        [&]<typename... Args>(Args&&... args) {
          // the template would not be needed,
          // when https://www.open-std.org/jtc1/sc22/wg21/docs/papers/2019/p1814r0.html
          // will be implemented in clang. It is already implemented in gcc 12.2
          typename PreFE<BH, useFlat, useEigenRef>::template FE<std::decay_t<Args>::template Skill...> fe(
              *basisHandler_, std::forward<Args>(args)...);
 
          return fe;
        },
        skills_.args);
  }
};
 
template <bool useFlat = true, bool useEigenRef = false, typename BH, typename SK>
auto makeFE(const BH& basisHandler, const SK& sk) {
  FEFactory<BH, SK, useFlat, useEigenRef> factory(basisHandler, sk);
 
  return factory();
}
 
template <bool useEigenRef = false, typename BH, typename SK>
auto makeFEWithUnTouchedBasis(const BH& basisHandler, SK&& sk) {
  FEFactory<BH, SK, false, useEigenRef> factory(basisHandler, std::forward<SK>(sk));
 
  return factory();
}
 
} // namespace Ikarus