ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/OpenMD/trunk/src/flucq/FluctuatingChargeParameters.cpp
Revision: 1908
Committed: Fri Jul 19 21:25:45 2013 UTC (11 years, 9 months ago) by gezelter
File size: 3642 byte(s)
Log Message:
Added infrastructure for region-constrained fluctuating charges.

File Contents

# User Rev Content
1 gezelter 1731 /*
2     * Copyright (c) 2012 The University of Notre Dame. All Rights Reserved.
3     *
4     * The University of Notre Dame grants you ("Licensee") a
5     * non-exclusive, royalty free, license to use, modify and
6     * redistribute this software in source and binary code form, provided
7     * that the following conditions are met:
8     *
9     * 1. Redistributions of source code must retain the above copyright
10     * notice, this list of conditions and the following disclaimer.
11     *
12     * 2. Redistributions in binary form must reproduce the above copyright
13     * notice, this list of conditions and the following disclaimer in the
14     * documentation and/or other materials provided with the
15     * distribution.
16     *
17     * This software is provided "AS IS," without a warranty of any
18     * kind. All express or implied conditions, representations and
19     * warranties, including any implied warranty of merchantability,
20     * fitness for a particular purpose or non-infringement, are hereby
21     * excluded. The University of Notre Dame and its licensors shall not
22     * be liable for any damages suffered by licensee as a result of
23     * using, modifying or distributing the software or its
24     * derivatives. In no event will the University of Notre Dame or its
25     * licensors be liable for any lost revenue, profit or data, or for
26     * direct, indirect, special, consequential, incidental or punitive
27     * damages, however caused and regardless of the theory of liability,
28     * arising out of the use of or inability to use software, even if the
29     * University of Notre Dame has been advised of the possibility of
30     * such damages.
31     *
32     * SUPPORT OPEN SCIENCE! If you use OpenMD or its source code in your
33     * research, please cite the appropriate papers when you publish your
34     * work. Good starting points are:
35     *
36     * [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005).
37     * [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006).
38 gezelter 1879 * [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008).
39 gezelter 1731 * [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010).
40     * [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011).
41     */
42    
43     #include <iostream>
44     #include <stdlib.h>
45     #include <string.h>
46    
47     #include "flucq/FluctuatingChargeParameters.hpp"
48    
49     namespace OpenMD {
50     FluctuatingChargeParameters::FluctuatingChargeParameters() {
51     DefineOptionalParameterWithDefaultValue(Propagator, "propagator", "NVT");
52     DefineOptionalParameterWithDefaultValue(Friction, "friction", 1600.0);
53     DefineOptionalParameterWithDefaultValue(Tolerance, "tolerance", 1.0e-6);
54     DefineOptionalParameterWithDefaultValue(MaxIterations, "maxIterations", 100);
55     DefineOptionalParameterWithDefaultValue(TargetTemp, "targetTemp", 1.0e-6);
56     DefineOptionalParameterWithDefaultValue(TauThermostat, "tauThermostat", 10.0);
57 gezelter 1766 DefineOptionalParameterWithDefaultValue(DragCoefficient, "dragCoefficient", 0.01);
58 gezelter 1908 DefineOptionalParameterWithDefaultValue(ConstrainRegions, "constrainRegions", false);
59 gezelter 1731 }
60    
61     FluctuatingChargeParameters::~FluctuatingChargeParameters() {
62     }
63    
64     void FluctuatingChargeParameters::validate() {
65     CheckParameter(Propagator, isEqualIgnoreCase("NVT") || isEqualIgnoreCase("Langevin") || isEqualIgnoreCase("Minimizer") || isEqualIgnoreCase("Exact") );
66     CheckParameter(Friction, isNonNegative());
67     CheckParameter(Tolerance, isPositive());
68     CheckParameter(MaxIterations, isPositive());
69     CheckParameter(TargetTemp, isNonNegative());
70     CheckParameter(TauThermostat, isPositive());
71 gezelter 1766 CheckParameter(DragCoefficient, isPositive());
72 gezelter 1731 }
73    
74     }