cylindrical_transformation_parameters.hpp

Go to the documentation of this file.

/*
 * Copyright (C) 2010-2022 The ESPResSo project
 *
 * This file is part of ESPResSo.
 *
 * ESPResSo is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 *
 * ESPResSo is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */
 
#pragma once
 
#include <cmath>
#include <limits>
#include <stdexcept>
#include <string>
 
#include <utils/math/orthonormal_vec.hpp>
 
namespace Utils {
 
/**
 * @brief A class to hold and validate parameters for a cylindrical coordinate
 * transformations.
 */
class CylindricalTransformationParameters {
public:
  CylindricalTransformationParameters() = default;
  /**
   * @brief Main constructor.
   *
   * @param center The origin of the cylindrical coordinates.
   * @param axis The "z"-axis. Must be normalized.
   * @param orientation The axis along which phi = 0. Must be normalized and
   * orthogonal to axis.
   */
  CylindricalTransformationParameters(Utils::Vector3d const &center,
                                      Utils::Vector3d const &axis,
                                      Utils::Vector3d const &orientation)
      : m_center(center), m_axis(axis), m_orientation(orientation) {
    validate();
  }
  /**
   * @brief if you only provide center and axis, an orientation will be
   * generated automatically such that it is orthogonal to axis
   */
  CylindricalTransformationParameters(Utils::Vector3d const &center,
                                      Utils::Vector3d const &axis)
      : m_center(center), m_axis(axis),
        m_orientation(Utils::calc_orthonormal_vector(axis)) {}
 
  Utils::Vector3d center() const { return m_center; }
  Utils::Vector3d axis() const { return m_axis; }
  Utils::Vector3d orientation() const { return m_orientation; }
 
private:
  void validate() const {
    auto constexpr eps = 10. * std::numeric_limits<double>::epsilon();
    if (std::fabs(m_orientation * m_axis) > eps) {
      throw std::runtime_error(
          "CylindricalTransformationParameters: Axis and orientation must be "
          "orthogonal. Scalar product is " +
          std::to_string(m_orientation * m_axis));
    }
    if (std::fabs(m_axis.norm() - 1.) > eps) {
      throw std::runtime_error("CylindricalTransformationParameters: Axis must "
                               "be normalized. Norm is " +
                               std::to_string(m_axis.norm()));
    }
    if (std::fabs(m_orientation.norm() - 1.) > eps) {
      throw std::runtime_error("CylindricalTransformationParameters: "
                               "orientation must be normalized. Norm is " +
                               std::to_string(m_orientation.norm()));
    }
  }
 
  Utils::Vector3d const m_center{};
  Utils::Vector3d const m_axis{0., 0., 1.};
  Utils::Vector3d const m_orientation{1., 0., 0.};
};
 
} // namespace Utils