constant.hpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: 2025 VTT Technical Research Centre of Finland Ltd
// SPDX-License-Identifier: AGPL-3.0-or-later
 
#ifndef PFC_INITIAL_CONDITIONS_CONSTANT_HPP
#define PFC_INITIAL_CONDITIONS_CONSTANT_HPP
 
#include "../field_modifier.hpp"
#include "openpfc/field/operations.hpp"
 
namespace pfc {
 
class Constant : public FieldModifier {
private:
  double m_n0;
 
public:
  Constant() = default;
 
  Constant(double n0) : m_n0(n0) {}
 
  double get_density() const { return m_n0; }
 
  void set_density(double n0) { m_n0 = n0; }
 
  void apply(Model &m, double t_unused) override {
    (void)t_unused; // Silence unused parameter warning
    // Apply constant in coordinate space over the local inbox
    // Preserves distributed behavior and keeps API consistent with new ops
    pfc::field::apply(m, get_field_name(),
                      [n0 = m_n0](const pfc::Real3 &) { return n0; });
  }
};
 
} // namespace pfc
 
#endif // PFC_INITIAL_CONDITIONS_CONSTANT_HPP