ViewVC Help
View File | Revision Log | Show Annotations | View Changeset | Root Listing
root/freemol/trunk/src/mengine/src/eangle.c
(Generate patch)
# Line 8 | Line 8
8   #include "derivs.h"
9   #include "hess.h"
10   #include "utility.h"
11 #include "gmmx.h"
11   #include "field.h"
12  
13   void eangle2a(int);
# Line 48 | Line 47
47          fprintf(pcmlogfile,"\nAngle Terms \n");
48          fprintf(pcmlogfile,"             At1      At2      At3     Angle     Thet0   Tconst      Ebend\n");
49      }
51    gmmx.hybrida = FALSE;
50  
51      for (i=0; i < angles.nang; i++)
52      {
# Line 87 | Line 85
85                         cosine = -1.0;
86                      angle = radian*acos(cosine);
87  
90                    if (angle < 0.75*angles.anat[i])
91                      if ( atom[ib].mmx_type != 80)
92                         gmmx.hybrida = TRUE;
93
88                      if (angles.angtype[i] == HARMONIC)
89                      {
90                          dt = angle - angles.anat[i];
# Line 219 | Line 213
213            deriv.dea[i][1] = 0.0;
214            deriv.dea[i][2] = 0.0;
215        }
222    gmmx.hybrida = FALSE;
216      for (i=0; i < angles.nang; i++)
217      {
218          ia = angles.i13[i][0];
# Line 260 | Line 253
253                      if (cosine < -1.0)
254                         cosine = -1.0;
255                      angle = radian*acos(cosine);
263                    if (angle < 0.75*angles.anat[i])
264                      if ( atom[ib].mmx_type != 80)
265                         gmmx.hybrida = TRUE;
256                      if (angles.angtype[i] == HARMONIC)
257                      {
258                          dt = angle - angles.anat[i];

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines