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

Comparing trunk/OOPSE-3.0/src/types/AtomType.cpp (file contents):
Revision 2097 by tim, Wed Mar 9 17:30:29 2005 UTC vs.
Revision 2204 by gezelter, Fri Apr 15 22:04:00 2005 UTC

# Line 110 | Line 110 | void AtomType::complete() {
110    }
111  
112  
113 < void AtomType::complete() {
113 >  void AtomType::complete() {
114      int isError;
115      GenericData* data;
116      
117      //notify a new LJtype atom type is created
118      if (isLennardJones()) {
119 <        data = getPropertyByName("LennardJones");
120 <        if (data != NULL) {
121 <            LJParamGenericData* ljData = dynamic_cast<LJParamGenericData*>(data);
119 >      data = getPropertyByName("LennardJones");
120 >      if (data != NULL) {
121 >        LJParamGenericData* ljData = dynamic_cast<LJParamGenericData*>(data);
122  
123 <            if (ljData != NULL) {
124 <                LJParam ljParam = ljData->getData();
123 >        if (ljData != NULL) {
124 >          LJParam ljParam = ljData->getData();
125                  
126 <                newLJtype(&atp.ident, &ljParam.sigma, &ljParam.epsilon, &ljParam.soft_pot, &isError);
126 >          newLJtype(&atp.ident, &ljParam.sigma, &ljParam.epsilon, &ljParam.soft_pot, &isError);
127  
128 <                if (isError != 0) {
129 <                    sprintf( painCave.errMsg,
130 <                           "Fortran rejected newLJtype\n");
131 <                    painCave.severity = OOPSE_ERROR;
132 <                    painCave.isFatal = 1;
133 <                    simError();          
134 <                }
128 >          if (isError != 0) {
129 >            sprintf( painCave.errMsg,
130 >                     "Fortran rejected newLJtype\n");
131 >            painCave.severity = OOPSE_ERROR;
132 >            painCave.isFatal = 1;
133 >            simError();          
134 >          }
135                  
136 <            } else {
137 <                    sprintf( painCave.errMsg,
138 <                           "Can not cast GenericData to LJParam\n");
139 <                    painCave.severity = OOPSE_ERROR;
140 <                    painCave.isFatal = 1;
141 <                    simError();          
142 <            }            
143 <        } else {
144 <            sprintf( painCave.errMsg, "Can not find Parameters for LennardJones\n");
145 <            painCave.severity = OOPSE_ERROR;
146 <            painCave.isFatal = 1;
147 <            simError();          
148 <        }
136 >        } else {
137 >          sprintf( painCave.errMsg,
138 >                   "Can not cast GenericData to LJParam\n");
139 >          painCave.severity = OOPSE_ERROR;
140 >          painCave.isFatal = 1;
141 >          simError();          
142 >        }            
143 >      } else {
144 >        sprintf( painCave.errMsg, "Can not find Parameters for LennardJones\n");
145 >        painCave.severity = OOPSE_ERROR;
146 >        painCave.isFatal = 1;
147 >        simError();          
148 >      }
149      }
150  
151      if (isElectrostatic()) {
# Line 160 | Line 160 | void AtomType::complete() {
160      }
161        
162      if (isCharge()) {
163 <        data = getPropertyByName("Charge");
164 <        if (data != NULL) {
165 <            DoubleGenericData* doubleData= dynamic_cast<DoubleGenericData*>(data);
163 >      data = getPropertyByName("Charge");
164 >      if (data != NULL) {
165 >        DoubleGenericData* doubleData= dynamic_cast<DoubleGenericData*>(data);
166  
167 <            if (doubleData != NULL) {
168 <                double charge = doubleData->getData();
169 <                setCharge(&atp.ident, &charge, &isError);
167 >        if (doubleData != NULL) {
168 >          double charge = doubleData->getData();
169 >          setCharge(&atp.ident, &charge, &isError);
170                  
171 <                if (isError != 0) {
172 <                    sprintf( painCave.errMsg,
173 <                           "Fortran rejected setCharge\n");
174 <                    painCave.severity = OOPSE_ERROR;
175 <                    painCave.isFatal = 1;
176 <                    simError();          
177 <                }
178 <            } else {
179 <                    sprintf( painCave.errMsg,
180 <                           "Can not cast GenericData to DoubleGenericData\n");
181 <                    painCave.severity = OOPSE_ERROR;
182 <                    painCave.isFatal = 1;
183 <                    simError();          
184 <            }
185 <        } else {
186 <            sprintf( painCave.errMsg, "Can not find Charge Parameters\n");
187 <            painCave.severity = OOPSE_ERROR;
188 <            painCave.isFatal = 1;
189 <            simError();          
190 <        }
171 >          if (isError != 0) {
172 >            sprintf( painCave.errMsg,
173 >                     "Fortran rejected setCharge\n");
174 >            painCave.severity = OOPSE_ERROR;
175 >            painCave.isFatal = 1;
176 >            simError();          
177 >          }
178 >        } else {
179 >          sprintf( painCave.errMsg,
180 >                   "Can not cast GenericData to DoubleGenericData\n");
181 >          painCave.severity = OOPSE_ERROR;
182 >          painCave.isFatal = 1;
183 >          simError();          
184 >        }
185 >      } else {
186 >        sprintf( painCave.errMsg, "Can not find Charge Parameters\n");
187 >        painCave.severity = OOPSE_ERROR;
188 >        painCave.isFatal = 1;
189 >        simError();          
190 >      }
191      }
192  
193      if (isEAM()) {
194 <        data = getPropertyByName("EAM");
195 <        if (data != NULL) {
196 <            EAMParamGenericData* eamData = dynamic_cast<EAMParamGenericData*>(data);
194 >      data = getPropertyByName("EAM");
195 >      if (data != NULL) {
196 >        EAMParamGenericData* eamData = dynamic_cast<EAMParamGenericData*>(data);
197  
198 <            if (eamData != NULL) {
198 >        if (eamData != NULL) {
199  
200 <                EAMParam eamParam = eamData->getData();
200 >          EAMParam eamParam = eamData->getData();
201                  
202  
203 <                newEAMtype(&eamParam.latticeConstant, &eamParam.nrho, &eamParam.drho,  &eamParam.nr, &eamParam.dr, &eamParam.rcut,
204 <                                &eamParam.rvals[0], &eamParam.rhovals[0], &eamParam.Frhovals[0], &atp.ident, &isError );
203 >          newEAMtype(&eamParam.latticeConstant, &eamParam.nrho, &eamParam.drho,  &eamParam.nr, &eamParam.dr, &eamParam.rcut,
204 >                     &eamParam.rvals[0], &eamParam.rhovals[0], &eamParam.Frhovals[0], &atp.ident, &isError );
205  
206 <                if (isError != 0) {
207 <                    sprintf( painCave.errMsg,
208 <                           "Fortran rejected newEAMtype\n");
209 <                    painCave.severity = OOPSE_ERROR;
210 <                    painCave.isFatal = 1;
211 <                    simError();          
212 <                }
213 <            } else {
214 <                    sprintf( painCave.errMsg,
215 <                           "Can not cast GenericData to EAMParam\n");
216 <                    painCave.severity = OOPSE_ERROR;
217 <                    painCave.isFatal = 1;
218 <                    simError();          
219 <            }
220 <        } else {
221 <            sprintf( painCave.errMsg, "Can not find EAM Parameters\n");
222 <            painCave.severity = OOPSE_ERROR;
223 <            painCave.isFatal = 1;
224 <            simError();          
225 <        }
206 >          if (isError != 0) {
207 >            sprintf( painCave.errMsg,
208 >                     "Fortran rejected newEAMtype\n");
209 >            painCave.severity = OOPSE_ERROR;
210 >            painCave.isFatal = 1;
211 >            simError();          
212 >          }
213 >        } else {
214 >          sprintf( painCave.errMsg,
215 >                   "Can not cast GenericData to EAMParam\n");
216 >          painCave.severity = OOPSE_ERROR;
217 >          painCave.isFatal = 1;
218 >          simError();          
219 >        }
220 >      } else {
221 >        sprintf( painCave.errMsg, "Can not find EAM Parameters\n");
222 >        painCave.severity = OOPSE_ERROR;
223 >        painCave.isFatal = 1;
224 >        simError();          
225 >      }
226      }
227  
228      
229 < }
229 >  }
230  
231 < void AtomType::addProperty(GenericData* genData) {
231 >  void AtomType::addProperty(GenericData* genData) {
232      properties_.addProperty(genData);  
233 < }
233 >  }
234  
235 < void AtomType::removeProperty(const std::string& propName) {
235 >  void AtomType::removeProperty(const std::string& propName) {
236      properties_.removeProperty(propName);  
237 < }
237 >  }
238  
239 < void AtomType::clearProperties() {
239 >  void AtomType::clearProperties() {
240      properties_.clearProperties();
241 < }
241 >  }
242  
243 < std::vector<std::string> AtomType::getPropertyNames() {
243 >  std::vector<std::string> AtomType::getPropertyNames() {
244      return properties_.getPropertyNames();  
245 < }
245 >  }
246        
247 < std::vector<GenericData*> AtomType::getProperties() {
247 >  std::vector<GenericData*> AtomType::getProperties() {
248      return properties_.getProperties();
249 < }
249 >  }
250  
251 < GenericData* AtomType::getPropertyByName(const std::string& propName) {
251 >  GenericData* AtomType::getPropertyByName(const std::string& propName) {
252      return properties_.getPropertyByName(propName);
253 < }  
253 >  }  
254   }

Diff Legend

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