ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/freemol/trunk/src/mengine/src/pcm7.c
(Generate patch)
# Line 28 | Line 28
28      } solvent;
29  
30   int Missing_constants;
31 FILE *errfile;
31      
33 FILE * fopen_path ( char * , char * , char * ) ;
32   void message_alert(char *,char *);
33   void pcm7_min(void);
34   int isbond(int, int);
# Line 88 | Line 86
86   void hdel(int);
87   void set_atomtypes(int);
88   void type(void);
91 void generate_bonds(void);
89   void zero_data(void);
90   void xlogp(float *);
91  
# Line 98 | Line 95
95    int i,j,nRet;
96      char string[30];
97      double etot;
98 <    if (minim_control.field == MM3)
102 <    {
103 <        // copy mm3 types into type file and retype
104 <        pot.use_hbond = FALSE;
105 <        hbondreset();
106 <        hdel(1);
107 <        strcpy(string,"mm3.prm");
108 <        zero_data();
109 <        read_datafiles(string);
110 <        set_field();
111 <        default_intype = MM3;
112 <        default_outtype = MM3;
113 <        set_atomtypes(MM3);
114 <        generate_bonds();
115 <        if (minim_control.added_const)
116 <           get_added_const();
117 <    } else if (minim_control.field == MMFF94)
98 >    if (minim_control.field == MMFF94)
99      {
100          // copy mm3 types into type file and retype
101          pot.use_hbond = FALSE;
102          pot.use_picalc = FALSE;
122        hbondreset();
123        pireset();
103          hdel(1);
104          set_field();
126        default_intype = MMFF94;
127        default_outtype = MMFF94;
105          set_atomtypes(MMFF94);
129        generate_bonds();
106          if (minim_control.added_const)
107             get_added_const();
108      } else
# Line 140 | Line 116
116          }
117          type();
118          set_field();
143        generate_bonds();
119      }
120  
121              
# Line 180 | Line 155
155      }
156   /*  assign local geometry potential function parameters  */
157      Missing_constants = FALSE;
183    //     errfile = fopen_path(pcwindir,"pcmod.err","w");
158       if (pot.use_bond || pot.use_strbnd)  nRet = kbond();
159       if (nRet == FALSE) // use end_calc to free memory
160       {
187       //         message_alert("Parameters missing. See pcmod.err for information","Error");
188         //         fclose(errfile);
161           energies.total = 9999.;
162           return FALSE;
163       }
# Line 199 | Line 171
171  
172         if (Missing_constants == TRUE)
173         {
202         //           message_alert("Parameters missing. See pcmod.err for information","Error");
203         //           fclose(errfile);
174             energies.total = 9999.0;
175             return (FALSE);
206       } else
207       {
208         //           fclose(errfile);
209         //           remove("pcmod.err");
176         }    
177         if (minim_values.iprint == TRUE)
178            etot = print_energy();

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines