Line data Source code
1 : /* +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 2 : Copyright (c) 2011-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 "Colvar.h" 23 : #include "ActionRegister.h" 24 : 25 : namespace PLMD { 26 : namespace colvar { 27 : 28 : //+PLUMEDOC COLVAR VOLUME 29 : /* 30 : Calculate the volume of the simulation box. 31 : 32 : \par Examples 33 : 34 : The following input tells plumed to print the volume of the system 35 : \plumedfile 36 : vol: VOLUME 37 : PRINT ARG=vol 38 : \endplumedfile 39 : 40 : */ 41 : //+ENDPLUMEDOC 42 : 43 : 44 : class Volume : public Colvar { 45 : 46 : public: 47 : explicit Volume(const ActionOptions&); 48 : // active methods: 49 : void calculate() override; 50 : /// Register all the keywords for this action 51 : static void registerKeywords( Keywords& keys ); 52 : }; 53 : 54 10437 : PLUMED_REGISTER_ACTION(Volume,"VOLUME") 55 : 56 9 : Volume::Volume(const ActionOptions&ao): 57 9 : PLUMED_COLVAR_INIT(ao) 58 : { 59 : std::vector<AtomNumber> atoms; 60 9 : checkRead(); 61 : 62 9 : addValueWithDerivatives(); setNotPeriodic(); 63 9 : requestAtoms(atoms); 64 9 : } 65 : 66 10 : void Volume::registerKeywords( Keywords& keys ) { 67 10 : Action::registerKeywords( keys ); 68 10 : ActionWithValue::registerKeywords( keys ); 69 10 : ActionAtomistic::registerKeywords( keys ); 70 10 : } 71 : 72 : 73 : // calculator 74 148 : void Volume::calculate() { 75 : 76 148 : double v=getBox().determinant(); 77 148 : setBoxDerivatives(-v*Tensor::identity()); 78 148 : setValue (v); 79 148 : } 80 : 81 : } 82 : } 83 : 84 : 85 :