Line data    Source code

       1            0 : // Distributed under the MIT License.
       2              : // See LICENSE.txt for details.
       3              : 
       4              : #pragma once
       5              : 
       6              : #include <cstddef>
       7              : #include <pup.h>
       8              : #include <string>
       9              : 
      10              : #include "NumericalAlgorithms/Strahlkorper/InitialShape.hpp"
      11              : #include "Options/Context.hpp"
      12              : #include "Options/String.hpp"
      13              : #include "Utilities/Serialization/CharmPupable.hpp"
      14              : #include "Utilities/TMPL.hpp"
      15              : 
      16              : namespace ylm {
      17              : template <typename Frame>
      18              : class Strahlkorper;
      19              : 
      20              : namespace InitialShapes {
      21              : /*!
      22              :  * \ingroup SurfacesGroup
      23              :  * \brief An initial Strahlkorper shape read from an H5 file.
      24              :  */
      25              : template <typename Frame>
      26            1 : class FromFile : public InitialShape<Frame> {
      27              :  public:
      28            0 :   struct H5Filename {
      29            0 :     using type = std::string;
      30            0 :     static constexpr Options::String help = {
      31              :         "H5 file containing the Strahlkorper coefficients"};
      32              :   };
      33            0 :   struct SubfileName {
      34            0 :     using type = std::string;
      35            0 :     static constexpr Options::String help = {
      36              :         "Subfile (without leading slash or .dat extension) within "
      37              :         "the H5 file that contains the Strahlkorper coefficients"};
      38              :   };
      39            0 :   struct Time {
      40            0 :     using type = double;
      41            0 :     static constexpr Options::String help = {
      42              :         "Time at which to read the Strahlkorper coefficients"};
      43              :   };
      44            0 :   struct TimeEpsilon {
      45            0 :     using type = double;
      46            0 :     static constexpr Options::String help = {
      47              :         "Tolerance for matching the requested time to read the Strahlkorper "
      48              :         "coefficients"};
      49              :   };
      50            0 :   struct CheckFrame {
      51            0 :     using type = bool;
      52            0 :     static constexpr Options::String help = {
      53              :         "Whether to check that the frame in the file matches the requested "
      54              :         "frame"};
      55              :   };
      56              : 
      57            0 :   using options =
      58              :       tmpl::list<H5Filename, SubfileName, Time, TimeEpsilon, CheckFrame>;
      59            0 :   static constexpr Options::String help = {
      60              :       "Construct a Strahlkorper from coefficients read from an H5 file."};
      61            0 :   static std::string name() { return "FromFile"; }
      62              : 
      63            0 :   FromFile() = default;
      64            0 :   FromFile(std::string h5_filename, std::string subfile_name, double time,
      65              :            double time_epsilon, bool check_frame);
      66              : 
      67              :   /// \cond
      68              :   explicit FromFile(CkMigrateMessage* msg);
      69              :   using PUP::able::register_constructor;
      70              :   WRAPPED_PUPable_decl_template(FromFile);  // NOLINT
      71              :   /// \endcond
      72              : 
      73            1 :   Strahlkorper<Frame> strahlkorper(
      74              :       size_t l_max, const Options::Context& context) const override;
      75              : 
      76            0 :   void pup(PUP::er& p) override;
      77              : 
      78              :  private:
      79            0 :   std::string h5_filename_{};
      80            0 :   std::string subfile_name_{};
      81            0 :   double time_{};
      82            0 :   double time_epsilon_{};
      83            0 :   bool check_frame_{};
      84              : };
      85              : }  // namespace InitialShapes
      86              : }  // namespace ylm