1310 |
|
} |
1311 |
|
|
1312 |
|
gBend = new GhostBend( *the_atoms[a], |
1313 |
< |
*the_atoms[b] ); |
1313 |
> |
*the_atoms[b], |
1314 |
> |
*the_atoms[c] ); |
1315 |
> |
|
1316 |
|
gBend->setConstants( currentBendType->k1, |
1317 |
|
currentBendType->k2, |
1318 |
|
currentBendType->k3, |
1328 |
|
currentBendType->k3, |
1329 |
|
currentBendType->t0 ); |
1330 |
|
bendArray[i] = qBend; |
1331 |
< |
} |
1331 |
> |
} |
1332 |
|
} |
1333 |
|
} |
1334 |
|
} |
1612 |
|
} |
1613 |
|
|
1614 |
|
info.d0 = atof( the_token ); |
1615 |
+ |
|
1616 |
+ |
info.k0=0.0; |
1617 |
|
} |
1618 |
|
else if( !strcmp( bondType, "harmonic" ) ){ |
1619 |
|
info.type = HARMONIC_BOND; |