Line data Source code
1 : /* +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 2 : Copyright (c) 2015-2023 The plumed team 3 : (see the PEOPLE file at the root of the distribution for a list of names) 4 : 5 : See http://www.plumed.org for more information. 6 : 7 : This file is part of plumed, version 2. 8 : 9 : plumed is free software: you can redistribute it and/or modify 10 : it under the terms of the GNU Lesser General Public License as published by 11 : the Free Software Foundation, either version 3 of the License, or 12 : (at your option) any later version. 13 : 14 : plumed is distributed in the hope that it will be useful, 15 : but WITHOUT ANY WARRANTY; without even the implied warranty of 16 : MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 17 : GNU Lesser General Public License for more details. 18 : 19 : You should have received a copy of the GNU Lesser General Public License 20 : along with plumed. If not, see <http://www.gnu.org/licenses/>. 21 : +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ */ 22 : #include "LandmarkSelectionBase.h" 23 : #include "core/ActionRegister.h" 24 : 25 : //+PLUMEDOC LANDMARKS LANDMARK_SELECT_STRIDE 26 : /* 27 : Select every \f$k\f$th landmark from the trajectory. 28 : 29 : \par Examples 30 : 31 : */ 32 : //+ENDPLUMEDOC 33 : 34 : namespace PLMD { 35 : namespace analysis { 36 : 37 : class SelectWithStride : public LandmarkSelectionBase { 38 : public: 39 : static void registerKeywords( Keywords& keys ); 40 : explicit SelectWithStride( const ActionOptions& ao ); 41 : void selectLandmarks() override; 42 : }; 43 : 44 10435 : PLUMED_REGISTER_ACTION(SelectWithStride,"LANDMARK_SELECT_STRIDE") 45 : 46 9 : void SelectWithStride::registerKeywords( Keywords& keys ) { 47 9 : LandmarkSelectionBase::registerKeywords( keys ); 48 9 : } 49 : 50 8 : SelectWithStride::SelectWithStride( const ActionOptions& ao ): 51 : Action(ao), 52 8 : LandmarkSelectionBase(ao) 53 : { 54 8 : } 55 : 56 10 : void SelectWithStride::selectLandmarks() { 57 10 : unsigned stride = std::floor( my_input_data->getNumberOfDataPoints() / getNumberOfDataPoints() ), max=stride*getNumberOfDataPoints(); 58 296 : for(unsigned i=0; i<max; i+=stride) selectFrame( i ); 59 10 : } 60 : 61 : } 62 : }