settings.hpp

/******************************************************************************
 * Copyright (C) Siarhei Uzunbajakau, 2023.
 *
 * This program is free software. You can use, modify, and redistribute it under
 * the terms of the GNU Lesser General Public License as published by the Free
 * Software Foundation, either version 3 or (at your option) any later version.
 * This program is distributed without any warranty.
 *
 * Refer to COPYING.LESSER for more details.
 ******************************************************************************/
 
#ifndef SettingsSCH_H__
#define SettingsSCH_H__
 
#include "constants.hpp"
#include <deal.II/base/types.h>
 
using namespace dealii;
 
class SettingsSCH : public Constants::Physics
{
public:
  SettingsSCH(){};
 
  const unsigned int nr_threads_max = 0;
 
  const double ep_0 = permittivity_fs;
 
  double a = 0.5;
 
  double b = 1.0;
 
  double d1 = 0.2;
 
  const double rd1 = sqrt(DIMENSION__) * d1;
 
#if DIMENSION__ == 2
  const double kappa_f = ep_0 / (a * std::log(b / a));
#endif
 
#if DIMENSION__ == 3
  const double kappa_f = ep_0 * b / (a * (b - a));
#endif
 
  const types::boundary_id bid = 1;
 
  const double eps = 1e-12;
 
  const bool print_time_tables = false;
 
  const bool project_exact_solution = false;
 
  const bool log_cg_convergence = false;
};
 
using Settings = SettingsSCH;
 
#endif