src/integrators/Integrator.hpp

/*
 * Copyright (C) 2000-2004  Object Oriented Parallel Simulation Engine (OOPSE) project
 * 
 * Contact: oopse@oopse.org
 * 
 * This program is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public License
 * as published by the Free Software Foundation; either version 2.1
 * of the License, or (at your option) any later version.
 * All we ask is that proper credit is given for our work, which includes
 * - but is not limited to - adding the above copyright notice to the beginning
 * of your source code files, and to any copyright notice that you may distribute
 * with programs based on this work.
 * 
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU Lesser General Public License for more details.
 * 
 * You should have received a copy of the GNU Lesser General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
 *
 */
 
 /**
  * @file Integrator.hpp
  * @author tlin
  * @date 11/08/2004
  * @time 13:25am
  * @version 1.0
  */

#ifndef INTEGRATORS_INTEGRATOR_HPP
#define INTEGRATORS_INTEGRATOR_HPP

#include "brains/ForceManager.hpp"
#include "io/DumpWriter.hpp"
#include "io/StatWriter.hpp"
#include "integrators/Velocitizer.hpp"

namespace oopse {


/**
 * @class Integrator Integrator.hpp "integrators/Integrator.hpp"
 * @brief Base class of Integrator
 * @todo document
 */
class Integrator {
    public:

        virtual ~Integrator();

        //avoid public virtual function        
        void integrate() {
            doIntegrate();
        }

        void update() {
            doUpdate();
        }

        void setForceManager(ForceManager* forceMan) {
            if (forceMan_ != forceMan && forceMan_  != NULL) {
                delete forceMan_;
            }
            forceMan_ = forceMan;
        }

        void setVelocitizer(Velocitizer* velocitizer) {
            if (velocitizer_ != velocitizer && velocitizer_  != NULL) {
                delete velocitizer_;
            }
            velocitizer_  = velocitizer;
        }
        
    protected:

        virtual void doIntegrate() = 0;
        virtual void doUpdate() {}
        
        Integrator(SimInfo* info);
        
        SimInfo* info_;
        ForceManager* forceMan_;
        bool needPotential;
        bool needStress;
        
        Velocitizer* velocitizer_;
        bool needVelocityScaling;
        double targetScalingTemp;
    
        DumpWriter*dumpWriter;
        StatWriter* statWriter;
        Thermo thermo;

        double runTime;
        double sampleTime;
        double statusTime;
        double thermalTime;
        double dt;

        Snapshot* currentSnapshot_; //During the integration, the address of currentSnapshot Will not change
        
    private:
        
        virtual double calcConservedQuantity() = 0;
        
        virtual DumpWriter* createDumpWriter() = 0;
        
        virtual StatWriter* createStatWriter() = 0;
};

    
}
#endif //INTEGRATORS_INTEGRATOR_HPP
Revision:	1820
Committed:	Thu Dec 2 00:09:35 2004 UTC (19 years, 7 months ago) by tim
File size:	3238 byte(s)
Log Message:	NVE get built
#	User	Rev	Content
1	tim	1721	/*
2			* Copyright (C) 2000-2004 Object Oriented Parallel Simulation Engine (OOPSE) project
3			*
4			* Contact: oopse@oopse.org
5			*
6			* This program is free software; you can redistribute it and/or
7			* modify it under the terms of the GNU Lesser General Public License
8			* as published by the Free Software Foundation; either version 2.1
9			* of the License, or (at your option) any later version.
10			* All we ask is that proper credit is given for our work, which includes
11			* - but is not limited to - adding the above copyright notice to the beginning
12			* of your source code files, and to any copyright notice that you may distribute
13			* with programs based on this work.
14			*
15			* This program is distributed in the hope that it will be useful,
16			* but WITHOUT ANY WARRANTY; without even the implied warranty of
17			* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18			* GNU Lesser General Public License for more details.
19			*
20			* You should have received a copy of the GNU Lesser General Public License
21			* along with this program; if not, write to the Free Software
22			* Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
23			*
24			*/
25
26			/**
27			* @file Integrator.hpp
28			* @author tlin
29			* @date 11/08/2004
30			* @time 13:25am
31			* @version 1.0
32			*/
33	gezelter	1490
34	tim	1721	#ifndef INTEGRATORS_INTEGRATOR_HPP
35			#define INTEGRATORS_INTEGRATOR_HPP
36	gezelter	1490
37	tim	1819	#include "brains/ForceManager.hpp"
38			#include "io/DumpWriter.hpp"
39			#include "io/StatWriter.hpp"
40			#include "integrators/Velocitizer.hpp"
41
42	tim	1721	namespace oopse {
43	gezelter	1490
44	tim	1819
45	tim	1721	/**
46			* @class Integrator Integrator.hpp "integrators/Integrator.hpp"
47			* @brief Base class of Integrator
48	tim	1775	* @todo document
49	tim	1721	*/
50			class Integrator {
51	gezelter	1490	public:
52
53	tim	1819	virtual ~Integrator();
54	gezelter	1490
55	tim	1819	//avoid public virtual function
56			void integrate() {
57			doIntegrate();
58			}
59	tim	1774
60	tim	1819	void update() {
61			doUpdate();
62			}
63	tim	1775
64			void setForceManager(ForceManager* forceMan) {
65	tim	1820	if (forceMan_ != forceMan && forceMan_ != NULL) {
66	tim	1775	delete forceMan_;
67			}
68			forceMan_ = forceMan;
69			}
70	tim	1819
71			void setVelocitizer(Velocitizer* velocitizer) {
72	tim	1820	if (velocitizer_ != velocitizer && velocitizer_ != NULL) {
73	tim	1819	delete velocitizer_;
74			}
75			velocitizer_ = velocitizer;
76			}
77	tim	1775
78	tim	1721	protected:
79	tim	1819
80			virtual void doIntegrate() = 0;
81			virtual void doUpdate() {}
82	tim	1721
83	tim	1819	Integrator(SimInfo* info);
84
85	tim	1775	SimInfo* info_;
86			ForceManager* forceMan_;
87	tim	1819	bool needPotential;
88			bool needStress;
89
90			Velocitizer* velocitizer_;
91			bool needVelocityScaling;
92			double targetScalingTemp;
93
94			DumpWriter*dumpWriter;
95			StatWriter* statWriter;
96			Thermo thermo;
97	tim	1775
98	tim	1819	double runTime;
99			double sampleTime;
100			double statusTime;
101			double thermalTime;
102			double dt;
103
104			Snapshot* currentSnapshot_; //During the integration, the address of currentSnapshot Will not change
105
106	tim	1775	private:
107	tim	1819
108	tim	1774	virtual double calcConservedQuantity() = 0;
109	tim	1820
110	tim	1819	virtual DumpWriter* createDumpWriter() = 0;
111	tim	1820
112	tim	1819	virtual StatWriter* createStatWriter() = 0;
113	gezelter	1490	};
114
115
116	tim	1721	}
117			#endif //INTEGRATORS_INTEGRATOR_HPP