ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
script_interface/cluster_analysis/Cluster.hpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2010-2022 The ESPResSo project
3 * Copyright (C) 2002,2003,2004,2005,2006,2007,2008,2009,2010
4 * Max-Planck-Institute for Polymer Research, Theory Group
5 *
6 * This file is part of ESPResSo.
7 *
8 * ESPResSo is free software: you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation, either version 3 of the License, or
11 * (at your option) any later version.
12 *
13 * ESPResSo is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the GNU General Public License
19 * along with this program. If not, see <http://www.gnu.org/licenses/>.
20 */
21
22#pragma once
23
24#include "core/BoxGeometry.hpp"
26
29
30#include <memory>
31#include <string>
32#include <vector>
33
34namespace ScriptInterface {
35namespace ClusterAnalysis {
36
37class Cluster : public AutoParameters<Cluster> {
38public:
39 Cluster() = default;
40 Variant do_call_method(std::string const &method,
41 VariantMap const &parameters) override {
42 if (method == "particle_ids") {
43 return m_cluster->particles;
44 }
45 if (method == "particles") {
46 return m_particle_list.lock()->call_method(
47 "by_ids", {{"id_selection", m_cluster->particles}});
48 }
49 if (method == "size") {
50 return (int)m_cluster->particles.size();
51 }
52 if (method == "longest_distance") {
53 return m_cluster->longest_distance();
54 }
55 if (method == "radius_of_gyration") {
56 return m_cluster->radius_of_gyration();
57 }
58 if (method == "fractal_dimension") {
59 auto const [df, mean_sq_residual] =
60 m_cluster->fractal_dimension(get_value<double>(parameters.at("dr")));
61 return std::vector<double>({df, mean_sq_residual});
62 }
63 if (method == "center_of_mass") {
64 return m_cluster->center_of_mass();
65 }
66 return {};
67 }
68 void set_cluster(std::shared_ptr<::ClusterAnalysis::Cluster> const &c) {
69 m_cluster = c;
70 }
71 void set_particle_list(std::weak_ptr<Particles::ParticleList> const &handle) {
72 m_particle_list = handle;
73 }
74
75private:
76 std::shared_ptr<::ClusterAnalysis::Cluster> m_cluster;
77 std::weak_ptr<Particles::ParticleList> m_particle_list;
78};
79
80} /* namespace ClusterAnalysis */
81} /* namespace ScriptInterface */
Bind parameters in the script interface.
void set_particle_list(std::weak_ptr< Particles::ParticleList > const &handle)
void set_cluster(std::shared_ptr<::ClusterAnalysis::Cluster > const &c)
Variant do_call_method(std::string const &method, VariantMap const &parameters) override
T get_value(Variant const &v)
Extract value of specific type T from a Variant.
std::unordered_map< std::string, Variant > VariantMap
Definition Variant.hpp:69
boost::make_recursive_variant< None, bool, int, std::size_t, double, std::string, ObjectRef, Utils::Vector3b, Utils::Vector3i, Utils::Vector2d, Utils::Vector3d, Utils::Vector4d, std::vector< int >, std::vector< double >, std::vector< boost::recursive_variant_ >, std::unordered_map< int, boost::recursive_variant_ >, std::unordered_map< std::string, boost::recursive_variant_ > >::type Variant
Possible types for parameters.
Definition Variant.hpp:67