#ifndef POLYCOEFF_H
#define POLYCOEFF_H

#include "mfem.hpp"
#include <cmath>

namespace polycoeff
{
  /**
   * @brief Computes the xi coefficient function.
   * 
   * @param x Input vector.
   * @return double The computed xi coefficient.
   */
    double nonlinearSourceCoeff(const mfem::Vector &x);

  /**
   * @brief Computes the vector xi coefficient function.
   * 
   * @param x Input vector.
   * @param v Output vector to store the computed xi coefficient.
   */
    void diffusionCoeff(const mfem::Vector &x, mfem::Vector &v);
} // namespace polyCoeff

#endif // POLYCOEFF_H