ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/group/trunk/OOPSE/libmdtools/EAM_FF.cpp
(Generate patch)

Comparing trunk/OOPSE/libmdtools/EAM_FF.cpp (file contents):
Revision 627 by chuckv, Wed Jul 16 21:49:59 2003 UTC vs.
Revision 1261 by gezelter, Fri Jun 11 14:14:10 2004 UTC

# Line 1 | Line 1
1 < #include <cstdlib>
2 < #include <cstdio>
3 < #include <cstring>
1 > #include <stdlib.h>
2 > #include <stdio.h>
3 > #include <string.h>
4  
5   #include <iostream>
6   using namespace std;
# Line 30 | Line 30 | namespace EAM_NS{
30      double mass;
31      double lattice_constant;
32      double eam_drho;  // The distance between each of the points indexed by rho.
33 +    double eam_rcut;  // The cutoff radius for eam.
34      double eam_dr;    // The distance between each of the rho points.    
35      int eam_nrho;  // Number of points indexed by rho
36      int eam_nr;    // The number of points based on r (Both Phi(r) and Rho(r)).
36 –    int eam_rcut;  // The cutoff radius for eam.
37      int eam_ident; // Atomic number
38      int ident;
39      int last;      //  0  -> default
# Line 41 | Line 41 | namespace EAM_NS{
41    } atomStruct;
42  
43    int parseAtom( char *lineBuffer, int lineNum, atomStruct &info, char *eamPotFile );
44 <  int parseEAM( atomStruct &info, char *eamPotFile );
44 >  int parseEAM( atomStruct &info, char *eamPotFile, double **eam_rvals,
45 >                double **eam_rhovals, double **eam_Frhovals);
46   #ifdef IS_MPI
47    
48    MPI_Datatype mpiAtomStructType;
# Line 75 | Line 76 | namespace EAM_NS{
76      void add( atomStruct &info, double *the_eam_rvals,
77                double *the_eam_rhovals,double *the_eam_Frhovals ){
78  
78 –      int i;
79 –
79        // check for duplicates
80        
81        if( !strcmp( info.name, name ) ){
82          sprintf( painCave.errMsg,
83 <                 "Duplicate LJ atom type \"%s\" found in "
84 <                 "the LJ_FF param file./n",
83 >                 "Duplicate EAM atom type \"%s\" found in "
84 >                 "the EAM_FF param file./n",
85                   name );
86          painCave.isFatal = 1;
87          simError();
# Line 135 | Line 134 | namespace EAM_NS{
134      double eam_drho; // The distance between each of the points indexed by rho.
135      int eam_nr;   // The number of points based on r (Both Phi(r) and Rho(r)).
136      double eam_dr;   // The distance between each of the rho points.
137 <    int eam_rcut; // The cutoff radius for eam.
137 >    double eam_rcut; // The cutoff radius for eam.
138  
139      double *eam_rvals;    // Z of r values
140      double *eam_rhovals;  // rho of r values
# Line 150 | Line 149 | namespace EAM_NS{
149  
150   }
151  
152 < using namespace LJ_NS;
152 > using namespace EAM_NS;
153  
154   //****************************************************************
155   // begins the actual forcefield stuff.  
156   //****************************************************************
157  
158 + EAM_FF::EAM_FF() {
159 +  EAM_FF("");
160 + }
161  
162 < EAM_FF::EAM_FF(){
162 > EAM_FF::EAM_FF(char* the_variant){
163  
164    char fileName[200];
165    char* ffPath_env = "FORCE_PARAM_PATH";
166    char* ffPath;
167    char temp[200];
166 –  char errMsg[1000];
168  
169    headAtomType = NULL;
170    currentAtomType = NULL;
171  
172 +  // Set eamRcut to 0.0
173 +  eamRcut = 0.0;
174 +
175    // do the funtion wrapping
176    wrapMeFF( this );
177  
# Line 178 | Line 182 | EAM_FF::EAM_FF(){
182    // Init the atomStruct mpi type
183  
184    atomStruct atomProto; // mpiPrototype
185 <  int atomBC[3] = {15,4,6};  // block counts
185 >  int atomBC[3] = {15,5,5};  // block counts
186    MPI_Aint atomDspls[3];           // displacements
187    MPI_Datatype atomMbrTypes[3];    // member mpi types
188  
# Line 200 | Line 204 | EAM_FF::EAM_FF(){
204    if( worldRank == 0 ){
205   #endif
206      
207 <    // generate the force file name
207 >    // generate the force file name  
208 >
209 >    strcpy( fileName, "EAM" );
210 >
211 >    if (strlen(the_variant) > 0) {
212 >      has_variant = 1;
213 >      strcpy( variant, the_variant);
214 >      strcat( fileName, ".");
215 >      strcat( fileName, variant );
216 >
217 >      sprintf( painCave.errMsg,
218 >               "Using %s variant of EAM force field.\n",
219 >               variant );
220 >      painCave.severity = OOPSE_INFO;
221 >      painCave.isFatal = 0;
222 >      simError();
223 >    }
224 >    strcat( fileName, ".frc");
225 >
226 >    //fprintf( stderr,"Trying to open %s\n", fileName );
227      
205 –    strcpy( fileName, "EAM_FF.frc" );
206 –    //    fprintf( stderr,"Trying to open %s\n", fileName );
207 –    
228      // attempt to open the file in the current directory first.
229      
230      frcFile = fopen( fileName, "r" );
# Line 227 | Line 247 | EAM_FF::EAM_FF(){
247        frcFile = fopen( fileName, "r" );
248        
249        if( frcFile == NULL ){
250 <        
250 >        
251          sprintf( painCave.errMsg,
252 <                 "Error opening the force field parameter file: %s\n"
253 <                 "Have you tried setting the FORCE_PARAM_PATH environment "
254 <                 "vairable?\n",
252 >                 "Error opening the force field parameter file:\n"
253 >                 "\t%s\n"
254 >                 "\tHave you tried setting the FORCE_PARAM_PATH environment "
255 >                 "variable?\n",
256                   fileName );
257 +        painCave.severity = OOPSE_ERROR;
258          painCave.isFatal = 1;
259          simError();
260        }
261      }
262 +
263      
264   #ifdef IS_MPI
265    }
# Line 263 | Line 286 | EAM_FF::~EAM_FF(){
286   #endif // is_mpi
287   }
288  
289 < void EAM_FF::initForceField( int ljMixRule ){
289 >
290 > void EAM_FF::calcRcut( void ){
291    
292 + #ifdef IS_MPI
293 +  double tempEamRcut = eamRcut;
294 +  MPI_Allreduce( &tempEamRcut, &eamRcut, 1, MPI_DOUBLE, MPI_MAX,
295 +                 MPI_COMM_WORLD);
296 + #endif  //is_mpi
297 +  entry_plug->setDefaultRcut(eamRcut);
298 + }
299 +
300 +
301 + void EAM_FF::initForceField( int ljMixRule ){
302    initFortran( ljMixRule, 0 );
303   }
304  
# Line 283 | Line 317 | void EAM_FF::cleanMe( void ){
317   #endif // is_mpi
318   }
319  
320 +
321   void EAM_FF::readParams( void ){
322  
323    atomStruct info;
324    info.last = 1; // initialize last to have the last set.
325                   // if things go well, last will be set to 0
326  
292 –  int i;
327    int identNum;
328    double *eam_rvals;    // Z of r values
329    double *eam_rhovals;  // rho of r values
# Line 305 | Line 339 | void EAM_FF::readParams( void ){
339      // read in the atom types.
340  
341      headAtomType = new LinkedAtomType;
342 <    
342 >
343      fastForward( "AtomTypes", "eam atom readParams" );
344  
345      // we are now at the AtomTypes section.
# Line 335 | Line 369 | void EAM_FF::readParams( void ){
369          // the parser returns 0 if the line was blank
370          if( parseAtom( readLine, lineNum, info, eamPotFile ) ){
371            parseEAM(info,eamPotFile, &eam_rvals,
372 <                   &eam_rhovals, &eam_Frhovals)){
372 >                   &eam_rhovals, &eam_Frhovals);
373            info.ident = identNum;
374            headAtomType->add( info, eam_rvals,
375                               eam_rhovals,eam_Frhovals );
# Line 395 | Line 429 | void EAM_FF::readParams( void ){
429      MPIcheckPoint();
430  
431      headAtomType = new LinkedAtomType;
432 <    recieveFrcStruct( &info, mpiAtomStructType );
432 >    receiveFrcStruct( &info, mpiAtomStructType );
433  
434      while( !info.last ){
435        
# Line 421 | Line 455 | void EAM_FF::readParams( void ){
455        
456        MPIcheckPoint();
457  
458 <      recieveFrcStruct( &info, mpiAtomStructType );
458 >      receiveFrcStruct( &info, mpiAtomStructType );
459  
460  
461      }
# Line 437 | Line 471 | void EAM_FF::readParams( void ){
471    int isDipole = 0;
472    int isSSD = 0;
473    int isGB = 0;
474 <  int isEam = 1;
474 >  int isEAM = 1;
475 >  int isCharge = 0;
476    double dipole = 0.0;
477 +  double charge = 0.0;
478 +  double eamSigma = 0.0;
479 +  double eamEpslon = 0.0;
480    
481    currentAtomType = headAtomType->next;
482    while( currentAtomType != NULL ){
# Line 451 | Line 489 | void EAM_FF::readParams( void ){
489                   &isDipole,
490                   &isGB,
491                   &isEAM,
492 <                 &(currentAtomType->epslon),
493 <                 &(currentAtomType->sigma),
492 >                 &isCharge,
493 >                 &eamEpslon,
494 >                 &eamSigma,
495 >                 &charge,
496                   &dipole,
497                   &isError );
498        if( isError ){
# Line 467 | Line 507 | void EAM_FF::readParams( void ){
507    }
508        
509    entry_plug->useLJ = 0;
510 <
510 >  entry_plug->useEAM = 1;
511    // Walk down again and send out EAM type
512    currentAtomType = headAtomType->next;
513    while( currentAtomType != NULL ){
514      
515      if( currentAtomType->name[0] != '\0' ){
516        isError = 0;
517 +
518        newEAMtype( &(currentAtomType->lattice_constant),
519                    &(currentAtomType->eam_nrho),
520                    &(currentAtomType->eam_drho),
521                    &(currentAtomType->eam_nr),
522                    &(currentAtomType->eam_dr),
523 +                  &(currentAtomType->eam_rcut),
524                    currentAtomType->eam_rvals,
525                    currentAtomType->eam_rhovals,
526                    currentAtomType->eam_Frhovals,
527                    &(currentAtomType->eam_ident),
528                    &isError);
529 +
530        if( isError ){
531          sprintf( painCave.errMsg,
532                   "Error initializing the \"%s\" atom type in fortran EAM\n",
# Line 503 | Line 546 | void EAM_FF::readParams( void ){
546    MPIcheckPoint();
547   #endif // is_mpi
548  
549 +
550 +
551   }
552  
553  
554   void EAM_FF::initializeAtoms( int nAtoms, Atom** the_atoms ){
555    
556    int i;
557 <
557 >  
558    // initialize the atoms
559    
515 –
516 –  Atom* thisAtom;
517 –
560    for( i=0; i<nAtoms; i++ ){
561      
562      currentAtomType = headAtomType->find( the_atoms[i]->getType() );
# Line 528 | Line 570 | void EAM_FF::initializeAtoms( int nAtoms, Atom** the_a
570      
571      the_atoms[i]->setMass( currentAtomType->mass );
572      the_atoms[i]->setIdent( currentAtomType->ident );
531 –    the_atoms[i]->setEAM();
573  
574 +    if (eamRcut < currentAtomType->eam_rcut) eamRcut = currentAtomType->eam_rcut;
575 +    
576    }
577   }
578  
579 < void LJ_FF::initializeBonds( int nBonds, Bond** BondArray,
579 > void EAM_FF::initializeBonds( int nBonds, Bond** BondArray,
580                               bond_pair* the_bonds ){
581    
582      if( nBonds ){
# Line 657 | Line 700 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
700                       double **eam_rvals,
701                       double **eam_rhovals,
702                       double **eam_Frhovals){
703 <  
703 >  double* myEam_rvals;
704 >  double* myEam_rhovals;
705 >  double* myEam_Frhovals;
706 >
707    char* ffPath_env = "FORCE_PARAM_PATH";
708    char* ffPath;
709    char* the_token;
710 +  char* eam_eof_test;
711    FILE *eamFile;
712 +  const int BUFFERSIZE = 3000;
713 +
714    char temp[200];
715    int linenumber;
716    int nReadLines;
717 +  char eam_read_buffer[BUFFERSIZE];
718  
719 +
720    int i,j;
721  
722    linenumber = 0;
# Line 674 | Line 725 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
725    eamFile = fopen( eamPotFile, "r" );
726    
727    
728 <  if( frcFile == NULL ){
728 >  if( eamFile == NULL ){
729      
730        // next see if the force path enviorment variable is set
731      
# Line 693 | Line 744 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
744  
745      
746      
747 <    if( frcFile == NULL ){
747 >    if( eamFile == NULL ){
748        
749        sprintf( painCave.errMsg,
750                 "Error opening the EAM force parameter file: %s\n"
751                 "Have you tried setting the FORCE_PARAM_PATH environment "
752 <               "vairable?\n",
752 >               "variable?\n",
753                 eamPotFile );
754        painCave.isFatal = 1;
755        simError();
# Line 706 | Line 757 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
757    }
758  
759    // First line is a comment line, read and toss it....
760 <  eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile);
760 >  eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
761    linenumber++;
762 <  if(eof_test == NULL){
762 >  if(eam_eof_test == NULL){
763      sprintf( painCave.errMsg,
764               "error in reading commment in %s\n", eamPotFile);
765      painCave.isFatal = 1;
# Line 718 | Line 769 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
769  
770  
771    // The Second line contains atomic number, atomic mass and a lattice constant
772 <  eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile);
772 >  eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
773    linenumber++;
774 <  if(eof_test == NULL){
774 >  if(eam_eof_test == NULL){
775      sprintf( painCave.errMsg,
776               "error in reading Identifier line in %s\n", eamPotFile);
777      painCave.isFatal = 1;
# Line 730 | Line 781 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
781  
782  
783      
784 <  if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
784 >  if ( (the_token = strtok( eam_read_buffer, " \n\t,;")) == NULL){
785      sprintf( painCave.errMsg,
786               "Error parseing EAM ident  line in %s\n", eamPotFile );
787      painCave.isFatal = 1;
# Line 756 | Line 807 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
807    info.lattice_constant = atof( the_token);
808  
809    // Next line is nrho, drho, nr, dr and rcut
810 <  eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile)
811 <  if(eof_test == NULL){
810 >  eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
811 >  if(eam_eof_test == NULL){
812      sprintf( painCave.errMsg,
813               "error in reading number of points line in %s\n", eamPotFile);
814      painCave.isFatal = 1;
815      simError();
816    }
817  
818 <  if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
818 >  if ( (the_token = strtok( eam_read_buffer, " \n\t,;")) == NULL){
819      sprintf( painCave.errMsg,
820               "Error parseing EAM nrho: line in %s\n", eamPotFile );
821      painCave.isFatal = 1;
# Line 806 | Line 857 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
857    info.eam_rcut = atof( the_token);
858  
859  
860 +
861 +
862 +
863    // Ok now we have to allocate point arrays and read in number of points
864    // Index the arrays for fortran, starting at 1
865 <  *eam_Frhovals = new double[info.nrho];
866 <  *eam_rvals    = new double[info.nr];
867 <  *eam_rhovals  = new double[info.nr];
865 >  myEam_Frhovals = new double[info.eam_nrho];
866 >  myEam_rvals    = new double[info.eam_nr];
867 >  myEam_rhovals  = new double[info.eam_nr];
868  
869    // Parse F of rho vals.
870  
871    // Assume for now that we have a complete number of lines
872 <  nReadLines = int(info.nrho/5);
872 >  nReadLines = int(info.eam_nrho/5);
873 >  
874  
875 +
876    for (i=0;i<nReadLines;i++){
877      j = i*5;
878 <    
878 >
879      // Read next line
880 <    eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile);
880 >    eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
881      linenumber++;
882 <    if(eof_test == NULL){
882 >    if(eam_eof_test == NULL){
883        sprintf( painCave.errMsg,
884                 "error in reading EAM file %s at line %d\n",
885                 eamPotFile,linenumber);
# Line 833 | Line 889 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
889      
890      // Parse 5 values on each line into array
891      // Value 1
892 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
892 >    if ( (the_token = strtok( eam_read_buffer, " \n\t,;")) == NULL){
893        sprintf( painCave.errMsg,
894                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
895        painCave.isFatal = 1;
896        simError();
897      }
898 +
899 +    myEam_Frhovals[j+0] = atof( the_token );
900      
843 –    *eam_Frhovals[j+0] = atof( the_token );
844 –    
901      // Value 2
902      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
903        sprintf( painCave.errMsg,
# Line 849 | Line 905 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
905        painCave.isFatal = 1;
906        simError();
907      }
908 <  
909 <    *eam_Frhovals[j+1] = atof( the_token );
908 >
909 >    myEam_Frhovals[j+1] = atof( the_token );
910      
911      // Value 3
912      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 859 | Line 915 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
915        painCave.isFatal = 1;
916        simError();
917      }
918 +
919 +    myEam_Frhovals[j+2] = atof( the_token );
920      
863 –    *eam_Frhovals[j+2] = atof( the_token );
864 –    
921      // Value 4
922      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
923        sprintf( painCave.errMsg,
# Line 869 | Line 925 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
925        painCave.isFatal = 1;
926        simError();
927      }
928 <    
929 <    *eam_Frhovals[j+3] = atof( the_token );
928 >
929 >    myEam_Frhovals[j+3] = atof( the_token );
930  
931      // Value 5
932      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 879 | Line 935 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
935        painCave.isFatal = 1;
936        simError();
937      }
938 +
939 +    myEam_Frhovals[j+4] = atof( the_token );
940      
883 –    *eam_Frhovals[j+4] = atof( the_token );
884 –    
941    }
942    // Parse Z of r vals
943    
944    // Assume for now that we have a complete number of lines
945 <  nReadLines = int(info.nr/5);
945 >  nReadLines = int(info.eam_nr/5);
946  
947    for (i=0;i<nReadLines;i++){
948      j = i*5;
949  
950      // Read next line
951 <    eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile);
951 >    eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
952      linenumber++;
953 <    if(eof_test == NULL){
953 >    if(eam_eof_test == NULL){
954        sprintf( painCave.errMsg,
955                 "error in reading EAM file %s at line %d\n",
956                 eamPotFile,linenumber);
# Line 904 | Line 960 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
960      
961      // Parse 5 values on each line into array
962      // Value 1
963 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
963 >    if ( (the_token = strtok( eam_read_buffer, " \n\t,;")) == NULL){
964        sprintf( painCave.errMsg,
965                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
966        painCave.isFatal = 1;
967        simError();
968      }
969      
970 <    *eam_rvals[j+0] = atof( the_token );
970 >    myEam_rvals[j+0] = atof( the_token );
971  
972      // Value 2
973      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 921 | Line 977 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
977        simError();
978      }
979    
980 <    *eam_rvals[j+1] = atof( the_token );
980 >    myEam_rvals[j+1] = atof( the_token );
981  
982      // Value 3
983      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 931 | Line 987 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
987        simError();
988      }
989    
990 <    *eam_rvals[j+2] = atof( the_token );
990 >    myEam_rvals[j+2] = atof( the_token );
991  
992      // Value 4
993      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 941 | Line 997 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
997        simError();
998      }
999    
1000 <    *eam_rvals[j+3] = atof( the_token );
1000 >    myEam_rvals[j+3] = atof( the_token );
1001  
1002      // Value 5
1003      if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
# Line 951 | Line 1007 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
1007        simError();
1008      }
1009    
1010 <    *eam_rvals[j+4] = atof( the_token );
1010 >    myEam_rvals[j+4] = atof( the_token );
1011  
1012    }
1013    // Parse rho of r vals
# Line 962 | Line 1018 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
1018      j = i*5;
1019  
1020      // Read next line
1021 <    eof_test = fgets(read_buffer, sizeof(read_buffer),eamFile);
1021 >    eam_eof_test = fgets(eam_read_buffer, sizeof(eam_read_buffer),eamFile);
1022      linenumber++;
1023 <    if(eof_test == NULL){
1023 >    if(eam_eof_test == NULL){
1024        sprintf( painCave.errMsg,
1025                 "error in reading EAM file %s at line %d\n",
1026                 eamPotFile,linenumber);
# Line 974 | Line 1030 | int EAM_NS::parseEAM(atomStruct &info, char *eamPotFil
1030    
1031      // Parse 5 values on each line into array
1032      // Value 1
1033 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
1033 >    if ( (the_token = strtok( eam_read_buffer, " \n\t,;")) == NULL){
1034        sprintf( painCave.errMsg,
1035                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
1036        painCave.isFatal = 1;
1037        simError();
1038      }
1039    
1040 <    *eam_rhovals[j+0] = atof( the_token );
1040 >    myEam_rhovals[j+0] = atof( the_token );
1041  
1042      // Value 2
1043 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
1043 >    if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
1044        sprintf( painCave.errMsg,
1045                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
1046        painCave.isFatal = 1;
1047        simError();
1048      }
1049    
1050 <    *eam_rhovals[j+1] = atof( the_token );
1050 >    myEam_rhovals[j+1] = atof( the_token );
1051  
1052      // Value 3
1053 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
1053 >    if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
1054        sprintf( painCave.errMsg,
1055                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
1056        painCave.isFatal = 1;
1057        simError();
1058      }
1059    
1060 <    *eam_rhovals[j+2] = atof( the_token );
1060 >    myEam_rhovals[j+2] = atof( the_token );
1061  
1062      // Value 4
1063 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
1063 >    if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
1064        sprintf( painCave.errMsg,
1065                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
1066        painCave.isFatal = 1;
1067        simError();
1068      }
1069    
1070 <    *eam_rhovals[j+3] = atof( the_token );
1070 >    myEam_rhovals[j+3] = atof( the_token );
1071  
1072      // Value 5
1073 <    if ( (the_token = strtok( read_buffer, " \n\t,;")) == NULL){
1073 >    if ( (the_token = strtok( NULL, " \n\t,;")) == NULL){
1074        sprintf( painCave.errMsg,
1075                 "Error parseing EAM nrho: line in %s\n", eamPotFile );
1076        painCave.isFatal = 1;
1077        simError();
1078      }
1079    
1080 <    *eam_rhovals[j+4] = atof( the_token );
1081 <
1080 >    myEam_rhovals[j+4] = atof( the_token );
1081 >
1082    }
1083 +  *eam_rvals = myEam_rvals;
1084 +  *eam_rhovals = myEam_rhovals;
1085 +  *eam_Frhovals = myEam_Frhovals;
1086  
1087    fclose(eamFile);
1088    return 0;

Diff Legend

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