--- trunk/tengDissertation/Appendix.tex	2006/06/08 15:44:14	2824
+++ trunk/tengDissertation/Appendix.tex	2006/06/09 02:20:39	2837
@@ -118,6 +118,7 @@ The Singleton pattern not only provides a mechanism to
 OOPSE}.
 
 \subsection{\label{appendixSection:singleton}Singleton}
+
 The Singleton pattern not only provides a mechanism to restrict
 instantiation of a class to one object, but also provides a global
 point of access to the object. Currently implemented as a global
@@ -127,22 +128,27 @@ static data approach in {\sc OOPSE}. {\tt IntegratorFa
 pollution.Although the singleton pattern can be implemented in
 various ways  to account for different aspects of the software
 designs, such as lifespan control \textit{etc}, we only use the
-static data approach in {\sc OOPSE}. {\tt IntegratorFactory} class
-is declared as
-\begin{lstlisting}[float,caption={[A classic Singleton design pattern implementation(I)] Declaration of {\tt IntegratorFactory} class.},label={appendixScheme:singletonDeclaration}]
+static data approach in {\sc OOPSE}. IntegratorFactory class is
+declared as
 
-  class IntegratorFactory {
-    public:
-      static IntegratorFactory* getInstance();
-    protected:
-      IntegratorFactory();
-    private:
-      static IntegratorFactory* instance_;
-  };
+\begin{lstlisting}[float,caption={[A classic Singleton design pattern implementation(I)] The declaration of of simple Singleton pattern.},label={appendixScheme:singletonDeclaration}]
+
+class IntegratorFactory {
+public:
+  static IntegratorFactory*
+  getInstance();
+protected:
+  IntegratorFactory();
+private:
+  static IntegratorFactory* instance_;
+};
+
 \end{lstlisting}
+
 The corresponding implementation is
-\begin{lstlisting}[float,caption={[A classic implementation of Singleton design pattern (II)] Implementation of {\tt IntegratorFactory} class.},label={appendixScheme:singletonImplementation}]
 
+\begin{lstlisting}[float,caption={[A classic implementation of Singleton design pattern (II)] The implementation of simple Singleton pattern.},label={appendixScheme:singletonImplementation}]
+
 IntegratorFactory::instance_ = NULL;
 
 IntegratorFactory* getInstance() {
@@ -151,176 +157,238 @@ IntegratorFactory* getInstance() {
   }
   return instance_;
 }
+
 \end{lstlisting}
-Since constructor is declared as {\tt protected}, a client can not
-instantiate {\tt IntegratorFactory} directly. Moreover, since the
-member function {\tt getInstance} serves as the only entry of access
-to {\tt IntegratorFactory}, this approach fulfills the basic
-requirement, a single instance. Another consequence of this approach
-is the automatic destruction since static data are destroyed upon
-program termination.
 
+Since constructor is declared as protected, a client can not
+instantiate IntegratorFactory directly. Moreover, since the member
+function getInstance serves as the only entry of access to
+IntegratorFactory, this approach fulfills the basic requirement, a
+single instance. Another consequence of this approach is the
+automatic destruction since static data are destroyed upon program
+termination.
+
 \subsection{\label{appendixSection:factoryMethod}Factory Method}
 
 Categoried as a creational pattern, the Factory Method pattern deals
 with the problem of creating objects without specifying the exact
 class of object that will be created. Factory Method is typically
 implemented by delegating the creation operation to the subclasses.
+Parameterized Factory pattern where factory method (
+createIntegrator member function) creates products based on the
+identifier (see List.~\ref{appendixScheme:factoryDeclaration}). If
+the identifier has been already registered, the factory method will
+invoke the corresponding creator (see List.~\ref{integratorCreator})
+which utilizes the modern C++ template technique to avoid excess
+subclassing.
 
-Registers a creator with a type identifier. Looks up the type
-identifier in the internal map. If it is found, it invokes the
-corresponding creator for the type identifier and returns its
-result.
-\begin{lstlisting}[float,caption={[The implementation of Factory pattern (I)].},label={appendixScheme:factoryDeclaration}]
-  class IntegratorCreator;
-  class IntegratorFactory {
-    public:
-      typedef std::map<string, IntegratorCreator*> CreatorMapType;
+\begin{lstlisting}[float,caption={[The implementation of Parameterized Factory pattern (I)]Source code of IntegratorFactory class.},label={appendixScheme:factoryDeclaration}]
 
-      bool registerIntegrator(IntegratorCreator* creator);
+class IntegratorFactory {
+public:
+  typedef std::map<string, IntegratorCreator*> CreatorMapType;
 
-      Integrator* createIntegrator(const string& id, SimInfo* info);
-
-    private:
-      CreatorMapType creatorMap_;
-  };
-\end{lstlisting}
-
-\begin{lstlisting}[float,caption={[The implementation of Factory pattern (II)].},label={appendixScheme:factoryDeclarationImplementation}]
-  bool IntegratorFactory::unregisterIntegrator(const string& id) {
-    return creatorMap_.erase(id) == 1;
+  bool registerIntegrator(IntegratorCreator* creator) {
+    return creatorMap_.insert(creator->getIdent(), creator).second;
   }
 
-  Integrator*
-  IntegratorFactory::createIntegrator(const string& id, SimInfo* info) {
+  Integrator* createIntegrator(const string& id, SimInfo* info) {
+    Integrator* result = NULL;
     CreatorMapType::iterator i = creatorMap_.find(id);
     if (i != creatorMap_.end()) {
-      //invoke functor to create object
-      return (i->second)->create(info);
-    } else {
-      return NULL;
+      result = (i->second)->create(info);
     }
+    return result;
   }
+
+private:
+  CreatorMapType creatorMap_;
+};
 \end{lstlisting}
 
-\begin{lstlisting}[float,caption={[The implementation of Factory pattern (III)].},label={appendixScheme:integratorCreator}]
+\begin{lstlisting}[float,caption={[The implementation of Parameterized Factory pattern (III)]Source code of creator classes.},label={appendixScheme:integratorCreator}]
 
-  class IntegratorCreator {
-  public:
+class IntegratorCreator {
+public:
     IntegratorCreator(const string& ident) : ident_(ident) {}
 
     const string& getIdent() const { return ident_; }
 
     virtual Integrator* create(SimInfo* info) const = 0;
 
-  private:
+private:
     string ident_;
-  };
+};
 
-  template<class ConcreteIntegrator>
-  class IntegratorBuilder : public IntegratorCreator {
-  public:
-    IntegratorBuilder(const string& ident) : IntegratorCreator(ident) {}
-    virtual  Integrator* create(SimInfo* info) const {
-      return new ConcreteIntegrator(info);
-    }
-  };
+template<class ConcreteIntegrator>
+class IntegratorBuilder : public IntegratorCreator {
+public:
+  IntegratorBuilder(const string& ident)
+                   : IntegratorCreator(ident) {}
+  virtual  Integrator* create(SimInfo* info) const {
+    return new ConcreteIntegrator(info);
+  }
+};
 \end{lstlisting}
 
 \subsection{\label{appendixSection:visitorPattern}Visitor}
 
-The purpose of the Visitor Pattern is to encapsulate an operation
-that you want to perform on the elements. The operation being
-performed on a structure can be switched without changing the
-interfaces  of the elements. In other words, one can add virtual
-functions into a set of classes without modifying their interfaces.
-The UML class diagram of Visitor patten is shown in
-Fig.~\ref{appendixFig:visitorUML}. {\tt Dump2XYZ} program in
-Sec.~\ref{appendixSection:Dump2XYZ} uses Visitor pattern
-extensively.
+The visitor pattern is designed to decouple the data structure and
+algorithms used upon them by collecting related operation from
+element classes into other visitor classes, which is equivalent to
+adding virtual functions into a set of classes without modifying
+their interfaces. Fig.~\ref{appendixFig:visitorUML} demonstrates the
+structure of Visitor pattern which is used extensively in {\tt
+Dump2XYZ}. In order to convert an OOPSE dump file, a series of
+distinct operations are performed on different StuntDoubles (See the
+class hierarchy in Fig.~\ref{oopseFig:hierarchy} and the declaration
+in List.~\ref{appendixScheme:element}). Since the hierarchies
+remains stable, it is easy to define a visit operation (see
+List.~\ref{appendixScheme:visitor}) for each class of StuntDouble.
+Note that using Composite pattern\cite{Gamma1994}, CompositVisitor
+manages a priority visitor list and handles the execution of every
+visitor in the priority list on different StuntDoubles.
 
 \begin{figure}
 \centering
-\includegraphics[width=\linewidth]{architecture.eps}
-\caption[The architecture of {\sc OOPSE}] {Overview of the structure
-of {\sc OOPSE}} \label{appendixFig:visitorUML}
+\includegraphics[width=\linewidth]{visitor.eps}
+\caption[The UML class diagram of Visitor patten] {The UML class
+diagram of Visitor patten.} \label{appendixFig:visitorUML}
 \end{figure}
 
-\begin{lstlisting}[float,caption={[The implementation of Visitor pattern (I)]Source code of the visitor classes.},label={appendixScheme:visitor}]
-  class BaseVisitor{
-    public:
-      virtual void visit(Atom* atom);
-      virtual void visit(DirectionalAtom* datom);
-      virtual void visit(RigidBody* rb);
-  };
-\end{lstlisting}
+%\begin{figure}
+%\centering
+%\includegraphics[width=\linewidth]{hierarchy.eps}
+%\caption[Class hierarchy for ojects in {\sc OOPSE}]{ A diagram of
+%the class hierarchy.
+%\begin{itemize}
+%\item A {\bf StuntDouble} is {\it any} object that can be manipulated by the
+%integrators and minimizers.
+%\item An {\bf Atom} is a fundamental point-particle that can be moved around during a simulation.
+%\item A {\bf DirectionalAtom} is an atom which has {\it orientational} as well as translational degrees of freedom.
+%\item A {\bf RigidBody} is a collection of {\bf Atom}s or {\bf
+%DirectionalAtom}s which behaves as a single unit.
+%\end{itemize}
+%} \label{oopseFig:hierarchy}
+%\end{figure}
 
 \begin{lstlisting}[float,caption={[The implementation of Visitor pattern (II)]Source code of the element classes.},label={appendixScheme:element}]
-  class StuntDouble {
-    public:
-      virtual void accept(BaseVisitor* v) = 0;
-  };
 
-  class Atom: public StuntDouble {
-    public:
-      virtual void accept{BaseVisitor* v*} {v->visit(this);}
-  };
+class StuntDouble { public:
+  virtual void accept(BaseVisitor* v) = 0;
+};
 
-  class DirectionalAtom: public Atom {
-    public:
-      virtual void accept{BaseVisitor* v*} {v->visit(this);}
-  };
+class Atom: public StuntDouble { public:
+  virtual void accept{BaseVisitor* v*} {
+    v->visit(this);
+  }
+};
 
-  class RigidBody: public StuntDouble {
-    public:
-      virtual void accept{BaseVisitor* v*} {v->visit(this);}
-  };
+class DirectionalAtom: public Atom { public:
+  virtual void accept{BaseVisitor* v*} {
+    v->visit(this);
+  }
+};
+
+class RigidBody: public StuntDouble { public:
+  virtual void accept{BaseVisitor* v*} {
+    v->visit(this);
+  }
+};
+
+\end{lstlisting}
+
+\begin{lstlisting}[float,caption={[The implementation of Visitor pattern (I)]Source code of the visitor classes.},label={appendixScheme:visitor}]
+
+class BaseVisitor{
+public:
+  virtual void visit(Atom* atom);
+  virtual void visit(DirectionalAtom* datom);
+  virtual void visit(RigidBody* rb);
+};
+
+class BaseAtomVisitor:public BaseVisitor{ public:
+  virtual void visit(Atom* atom);
+  virtual void visit(DirectionalAtom* datom);
+  virtual void visit(RigidBody* rb);
+};
+
+class SSDAtomVisitor:public BaseAtomVisitor{ public:
+  virtual void visit(Atom* atom);
+  virtual void visit(DirectionalAtom* datom);
+  virtual void visit(RigidBody* rb);
+};
+
+class CompositeVisitor: public BaseVisitor {
+public:
+
+  typedef list<pair<BaseVisitor*, int> > VistorListType;
+  typedef VistorListType::iterator VisitorListIterator;
+  virtual void visit(Atom* atom) {
+    VisitorListIterator i;
+    BaseVisitor* curVisitor;
+    for(i = visitorList.begin();i != visitorList.end();++i) {
+      atom->accept(*i);
+    }
+  }
+
+  virtual void visit(DirectionalAtom* datom) {
+    VisitorListIterator i;
+    BaseVisitor* curVisitor;
+    for(i = visitorList.begin();i != visitorList.end();++i) {
+      atom->accept(*i);
+    }
+  }
 
+  virtual void visit(RigidBody* rb) {
+    VisitorListIterator i;
+    std::vector<Atom*> myAtoms;
+    std::vector<Atom*>::iterator ai;
+    myAtoms = rb->getAtoms();
+    for(i = visitorList.begin();i != visitorList.end();++i) {{
+      rb->accept(*i);
+      for(ai = myAtoms.begin(); ai != myAtoms.end(); ++ai){
+        (*ai)->accept(*i);
+    }
+  }
+
+  void addVisitor(BaseVisitor* v, int priority);
+
+  protected:
+    VistorListType visitorList;
+};
+
 \end{lstlisting}
+
 \section{\label{appendixSection:concepts}Concepts}
 
 OOPSE manipulates both traditional atoms as well as some objects
 that {\it behave like atoms}.  These objects can be rigid
 collections of atoms or atoms which have orientational degrees of
-freedom.  Here is a diagram of the class heirarchy:
-
-%\begin{figure}
-%\centering
-%\includegraphics[width=3in]{heirarchy.eps}
-%\caption[Class heirarchy for StuntDoubles in {\sc oopse}-3.0]{ \\
-%The class heirarchy of StuntDoubles in {\sc oopse}-3.0. The
-%selection syntax allows the user to select any of the objects that
-%are descended from a StuntDouble.} \label{oopseFig:heirarchy}
-%\end{figure}
-
-\begin{itemize}
-\item A {\bf StuntDouble} is {\it any} object that can be manipulated by the
-integrators and minimizers.
-\item An {\bf Atom} is a fundamental point-particle that can be moved around during a simulation.
-\item A {\bf DirectionalAtom} is an atom which has {\it orientational} as well as translational degrees of freedom.
-\item A {\bf RigidBody} is a collection of {\bf Atom}s or {\bf
-DirectionalAtom}s which behaves as a single unit.
-\end{itemize}
-
-Every Molecule, Atom and DirectionalAtom in {\sc OOPSE} have their
-own names which are specified in the {\tt .md} file. In contrast,
-RigidBodies are denoted by their membership and index inside a
-particular molecule: [MoleculeName]\_RB\_[index] (the contents
-inside the brackets depend on the specifics of the simulation). The
-names of rigid bodies are generated automatically. For example, the
-name of the first rigid body in a DMPC molecule is DMPC\_RB\_0.
+freedom.  A diagram of the class hierarchy is illustrated in
+Fig.~\ref{oopseFig:hierarchy}. Every Molecule, Atom and
+DirectionalAtom in {\sc OOPSE} have their own names which are
+specified in the {\tt .md} file. In contrast, RigidBodies are
+denoted by their membership and index inside a particular molecule:
+[MoleculeName]\_RB\_[index] (the contents inside the brackets depend
+on the specifics of the simulation). The names of rigid bodies are
+generated automatically. For example, the name of the first rigid
+body in a DMPC molecule is DMPC\_RB\_0.
 
 \section{\label{appendixSection:syntax}Syntax of the Select Command}
 
-The most general form of the select command is: {\tt select {\it
-expression}}. This expression represents an arbitrary set of
-StuntDoubles (Atoms or RigidBodies) in {\sc OOPSE}. Expressions are
-composed of either name expressions, index expressions, predefined
-sets, user-defined expressions, comparison operators, within
-expressions, or logical combinations of the above expression types.
-Expressions can be combined using parentheses and the Boolean
-operators.
+{\sc OOPSE} provides a powerful selection utility to select
+StuntDoubles. The most general form of the select command is:
 
+{\tt select {\it expression}}.
+
+This expression represents an arbitrary set of StuntDoubles (Atoms
+or RigidBodies) in {\sc OOPSE}. Expressions are composed of either
+name expressions, index expressions, predefined sets, user-defined
+expressions, comparison operators, within expressions, or logical
+combinations of the above expression types. Expressions can be
+combined using parentheses and the Boolean operators.
+
 \subsection{\label{appendixSection:logical}Logical expressions}
 
 The logical operators allow complex queries to be constructed out of