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 SettingsABC_H__
#define SettingsABC_H__
 
#include "constants.hpp"
#include <deal.II/base/types.h>
 
using namespace dealii;
 
class SettingsABC : public Constants::Physics
{
public:
  SettingsABC(){};
 
  const unsigned int nr_threads_max = 0;
 
  const double ep_0 = permittivity_fs;
 
  const double a = 0.1;
 
  const double x0 = 0.1;
 
  const double R = 0.05;
 
  const double R_mid = 1.5 * (R + x0);
 
  const types::boundary_id bid_left = 3;
 
  const types::boundary_id bid_right = 1;
 
  const types::boundary_id bid_in = 1;
 
  const types::manifold_id mfid_left = 3;
 
  const types::manifold_id mfid_right = 1;
 
  const types::manifold_id mfid_infty = 2;
 
#if BC_TYPE__ == 0
  const types::boundary_id bid_infty = 0; // Neumann boundary condition on the
                                          // outer boundary.
#endif
 
#if BC_TYPE__ == 1
  const types::boundary_id bid_infty = 5; // Dirichlet boundary condition on the
                                          // outer boundary.
#endif
 
#if BC_TYPE__ == 2
  const types::boundary_id bid_infty = 2; // ABC on the outer boundary.
#endif
 
  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 = SettingsABC;
 
#endif