OOPSE/libmdtools/AbstractClasses.hpp

#ifndef _ABSTRACT_CLASSES_H_
#define _ABSTRACT_CLASSES_H_

#include "Atom.hpp"
#ifdef IS_MPI

#include "mpi.h"
#endif
class Constraint{

public: 
  Constraint(){}
  ~Constraint(){}
  
  int    get_a()             {return a;}
  void   set_a(int index)    {a = index;}
  int    get_b()             {return b;}
  void   set_b(int index)    {b = index;}
  double get_dsqr()          {return dsqr;}
  void   set_dsqr(double ds) {dsqr = ds;}

private:
  int a; /* index of constrained atom a */
  int b; /* index of constrained atom b */
  double dsqr; /* the square of the constraint distance */
};

class SRI{

public:
  SRI(){ c_potential_E = 0.0; }
  virtual ~SRI() {}

  virtual void calc_forces() = 0;

  double get_potential(){ return c_potential_E; }
  virtual int is_constrained() = 0;
  virtual Constraint *get_constraint() = 0;
  virtual void constrain(double bond_distance) = 0;

  virtual void printMe( void ) = 0;

protected:
  double c_potential_E;
};


class BaseIntegrator{
  
public:
  BaseIntegrator(){}
  virtual ~BaseIntegrator(){}
  
  virtual void integrate( void ) = 0;
  virtual double  getConservedQuantity(void) = 0;  
};
/*
template<typename T> class AbstractTemplateIntegrator : public T
{
//  public:
//    AbstractTemplateIntegrator() {}
//    virtual ~AbstractTemplateIntegrator() {}
    
//    virtual void integrate(void) {}
};

typedef AbstractTemplateIntegrator<BaseIntegrator> TemplateIntegrator; 
*/
#endif
Revision:	763
Committed:	Mon Sep 15 16:52:02 2003 UTC (20 years, 9 months ago) by tim
File size:	1439 byte(s)
Log Message:	add conserved quantity to statWriter fix bug of vector wrapping at NPTi
#	User	Rev	Content
1	mmeineke	377	#ifndef _ABSTRACT_CLASSES_H_
2			#define _ABSTRACT_CLASSES_H_
3
4			#include "Atom.hpp"
5			#ifdef IS_MPI
6
7			#include "mpi.h"
8			#endif
9			class Constraint{
10
11			public:
12			Constraint(){}
13			~Constraint(){}
14
15			int get_a() {return a;}
16			void set_a(int index) {a = index;}
17			int get_b() {return b;}
18			void set_b(int index) {b = index;}
19	mmeineke	424	double get_dsqr() {return dsqr;}
20	mmeineke	377	void set_dsqr(double ds) {dsqr = ds;}
21
22			private:
23			int a; /* index of constrained atom a */
24			int b; /* index of constrained atom b */
25			double dsqr; /* the square of the constraint distance */
26			};
27
28			class SRI{
29
30			public:
31			SRI(){ c_potential_E = 0.0; }
32			virtual ~SRI() {}
33
34			virtual void calc_forces() = 0;
35
36			double get_potential(){ return c_potential_E; }
37			virtual int is_constrained() = 0;
38			virtual Constraint *get_constraint() = 0;
39			virtual void constrain(double bond_distance) = 0;
40
41			virtual void printMe( void ) = 0;
42
43			protected:
44			double c_potential_E;
45			};
46
47
48	mmeineke	540	class BaseIntegrator{
49	mmeineke	377
50			public:
51	mmeineke	540	BaseIntegrator(){}
52			virtual ~BaseIntegrator(){}
53	mmeineke	377
54			virtual void integrate( void ) = 0;
55	tim	763	virtual double getConservedQuantity(void) = 0;
56	mmeineke	377	};
57	tim	645	/*
58			template<typename T> class AbstractTemplateIntegrator : public T
59			{
60			// public:
61			// AbstractTemplateIntegrator() {}
62			// virtual ~AbstractTemplateIntegrator() {}
63
64			// virtual void integrate(void) {}
65			};
66	mmeineke	377
67	tim	645	typedef AbstractTemplateIntegrator<BaseIntegrator> TemplateIntegrator;
68			*/
69	mmeineke	377	#endif