ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-4/src/applications/staticProps/Hxy.hpp
Revision: 2751
Committed: Fri May 12 21:34:43 2006 UTC (18 years, 2 months ago) by xsun
File size: 3402 byte(s)
Log Message:
Changes to calculate undulation spectrum

File Contents

# User Rev Content
1 xsun 2751 /*
2     * Copyright (c) 2005 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. Acknowledgement of the program authors must be made in any
10     * publication of scientific results based in part on use of the
11     * program. An acceptable form of acknowledgement is citation of
12     * the article in which the program was described (Matthew
13     * A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher
14     * J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented
15     * Parallel Simulation Engine for Molecular Dynamics,"
16     * J. Comput. Chem. 26, pp. 252-271 (2005))
17     *
18     * 2. Redistributions of source code must retain the above copyright
19     * notice, this list of conditions and the following disclaimer.
20     *
21     * 3. Redistributions in binary form must reproduce the above copyright
22     * notice, this list of conditions and the following disclaimer in the
23     * documentation and/or other materials provided with the
24     * distribution.
25     *
26     * This software is provided "AS IS," without a warranty of any
27     * kind. All express or implied conditions, representations and
28     * warranties, including any implied warranty of merchantability,
29     * fitness for a particular purpose or non-infringement, are hereby
30     * excluded. The University of Notre Dame and its licensors shall not
31     * be liable for any damages suffered by licensee as a result of
32     * using, modifying or distributing the software or its
33     * derivatives. In no event will the University of Notre Dame or its
34     * licensors be liable for any lost revenue, profit or data, or for
35     * direct, indirect, special, consequential, incidental or punitive
36     * damages, however caused and regardless of the theory of liability,
37     * arising out of the use of or inability to use software, even if the
38     * University of Notre Dame has been advised of the possibility of
39     * such damages.
40     *
41     *
42     * hxy.hpp
43     * OOPSE-2.0
44     *
45     * Created by Xiuquan Sun on 05/09/06.
46     * @author Xiuquan Sun
47     * @version $Id: Hxy.hpp,v 1.1 2006-05-12 21:34:43 xsun Exp $
48     *
49     */
50     #ifndef APPLICATIONS_STATICPROPS_HXY_HPP
51     #define APPLICATIONS_STATICPROPS_HXY_HPP
52    
53     #define c_re(c) ((c)[0])
54     #define c_im(c) ((c)[1])
55    
56     #include "applications/staticProps/RadialDistrFunc.hpp"
57     namespace oopse {
58    
59     class Hxy : public StaticAnalyser {
60    
61     public:
62     Hxy(SimInfo* info, const std::string& filename, const std::string& sele, int nbins_x, int nbins_y, int nrbins);
63    
64     int getnbins() {
65     return nbins;
66     }
67    
68     int getnBinsX() {
69     return nBinsX_;
70     }
71    
72     int getnBinsY() {
73     return nBinsY_;
74     }
75    
76     virtual void process();
77    
78     private:
79    
80     virtual void printSpectrum();
81    
82    
83     Snapshot* currentSnapshot_;
84    
85     int nProcessed_;
86     std::string selectionScript_;
87     SelectionEvaluator evaluator_;
88     SelectionManager seleMan_;
89    
90     double nBinsX_;
91     double nBinsY_;
92     int nbins;
93     double dfreq;
94    
95     std::vector<double> gridZ_;
96     std::vector<int> gridsample_;
97     std::vector< std::vector<double> > bin;
98     std::vector< std::vector<int> > samples;
99     std::vector<double> sum_bin, sum_bin_sq, avg_bin, avg_bin_sq;
100     std::vector<double> errbin_sum, errbin_sum_sq, errbin, errbin_sq;
101     };
102    
103     }
104     #endif
105    
106    
107