ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE-2.0/src/visitors/LipidTransVisitor.cpp
(Generate patch)

Comparing trunk/OOPSE-2.0/src/visitors/LipidTransVisitor.cpp (file contents):
Revision 2053 by tim, Fri Feb 18 23:07:32 2005 UTC vs.
Revision 2204 by gezelter, Fri Apr 15 22:04:00 2005 UTC

# Line 1 | Line 1
1 < /*
1 > /*
2   * Copyright (c) 2005 The University of Notre Dame. All Rights Reserved.
3   *
4   * The University of Notre Dame grants you ("Licensee") a
# Line 40 | Line 40
40   */
41  
42   #include "visitors/LipidTransVisitor.hpp"
43 + #include "utils/simError.h"
44  
45   namespace oopse {
46 < void LipidTransVisitor::update() {
46 >  LipidTransVisitor::LipidTransVisitor(SimInfo* info, const std::string& originSeleScript, const std::string& refSeleScript)
47 >    : BaseVisitor(), info_(info), originEvaluator_(info), originSeleMan_(info), refEvaluator_(info), refSeleMan_(info), refSd_(NULL) {
48  
49 <        Vector3d ref = refSd_->getPos();
50 <        origin_ = originDatom_->getPos();
51 <        Vector3d v1 =  ref - origin_;
52 <        info_->getSnapshotManager()->getCurrentSnapshot()->wrapVector(v1);
53 <        Vector3d zaxis = originDatom_->getElectroFrame().getColumn(2);
54 <        Vector3d xaxis = cross(v1, zaxis);
55 <        Vector3d yaxis = cross(zaxis, xaxis);
49 >
50 >      visitorName = "LipidTransVisitor";
51 >    
52 >      originEvaluator_.loadScriptString(originSeleScript);            
53 >      if (!originEvaluator_.isDynamic()) {  
54 >        originSeleMan_.setSelectionSet(originEvaluator_.evaluate());
55 >        if (originSeleMan_.getSelectionCount() == 1) {
56 >          int i;
57 >          originDatom_ = dynamic_cast<DirectionalAtom*>(originSeleMan_.beginSelected(i));
58 >          if (originDatom_ ==  NULL) {
59 >            sprintf(painCave.errMsg, "LipidTransVisitor: origin selection must select an directional atom");
60 >            painCave.isFatal = 1;
61 >            simError();                  
62 >          }
63 >        } else {
64 >          sprintf(painCave.errMsg, "LipidTransVisitor: origin selection must select an directional atom");
65 >          painCave.isFatal = 1;
66 >          simError();                  
67 >            
68 >        }
69 >      }
70  
71 <        xaxis.normalize();
72 <        yaxis.normalize();
73 <        zaxis.normalize();
71 >      refEvaluator_.loadScriptString(refSeleScript);
72 >      if (!refEvaluator_.isDynamic()) {  
73 >        refSeleMan_.setSelectionSet(refEvaluator_.evaluate());
74 >        if (refSeleMan_.getSelectionCount() == 1) {
75 >          int i;
76 >          refSd_ = refSeleMan_.beginSelected(i);
77 >            
78 >        } else {
79 >          //error
80 >            
81 >        }
82 >      }
83 >    
84 >    }
85 >
86 >  void LipidTransVisitor::update() {
87 >
88 >    Vector3d ref = refSd_->getPos();
89 >    origin_ = originDatom_->getPos();
90 >    Vector3d v1 =  ref - origin_;
91 >    info_->getSnapshotManager()->getCurrentSnapshot()->wrapVector(v1);
92 >    Vector3d zaxis = originDatom_->getElectroFrame().getColumn(2);
93 >    Vector3d xaxis = cross(v1, zaxis);
94 >    Vector3d yaxis = cross(zaxis, xaxis);
95 >
96 >    xaxis.normalize();
97 >    yaxis.normalize();
98 >    zaxis.normalize();
99          
100 <        rotMat_.setRow(0, xaxis);
101 <        rotMat_.setRow(1, yaxis);
102 <        rotMat_.setRow(2, zaxis);
100 >    rotMat_.setRow(0, xaxis);
101 >    rotMat_.setRow(1, yaxis);
102 >    rotMat_.setRow(2, zaxis);
103  
104  
105 < }
105 >  }
106  
107 < void LipidTransVisitor::internalVisit(StuntDouble *sd) {
107 >  void LipidTransVisitor::internalVisit(StuntDouble *sd) {
108      GenericData *                     data;
109      AtomData *                        atomData;
110      AtomInfo *                        atomInfo;
# Line 72 | Line 113 | void LipidTransVisitor::internalVisit(StuntDouble *sd)
113      data = sd->getPropertyByName("ATOMDATA");
114  
115      if (data != NULL) {
116 <        atomData = dynamic_cast<AtomData *>(data);
116 >      atomData = dynamic_cast<AtomData *>(data);
117  
118 <        if (atomData == NULL)
119 <            return;
118 >      if (atomData == NULL)
119 >        return;
120      } else
121 <        return;
121 >      return;
122  
123      Snapshot* currSnapshot = info_->getSnapshotManager()->getCurrentSnapshot();
124      
125      for( atomInfo = atomData->beginAtomInfo(i); atomInfo; atomInfo = atomData->nextAtomInfo(i) ) {
126  
127 <        Vector3d tmp= atomInfo->pos - origin_;
128 <        //we don't need to do the wrapping anymore, since LipidTransVisitor will be execute after
129 <        //WrappingVistor and ReplicateVisitor
130 <        atomInfo->pos = rotMat_ * tmp;
90 <        atomInfo->dipole = rotMat_ * atomInfo->dipole;
127 >      Vector3d tmp= atomInfo->pos - origin_;
128 >      currSnapshot->wrapVector(tmp);
129 >      atomInfo->pos = rotMat_ * tmp;;
130 >      atomInfo->dipole = rotMat_ * atomInfo->dipole;
131      }
132 < }
132 >  }
133  
134 < const std::string LipidTransVisitor::toString() {
134 >  const std::string LipidTransVisitor::toString() {
135      char        buffer[65535];
136      std::string result;
137  
# Line 111 | Line 151 | const std::string LipidTransVisitor::toString() {
151      result += buffer;
152  
153      return result;
154 < }
154 >  }
155  
156   }

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines