| 476 |
|
// Init the atomStruct mpi type |
| 477 |
|
|
| 478 |
|
atomStruct atomProto; // mpiPrototype |
| 479 |
< |
int atomBC[3] = {15,6,4}; // block counts |
| 479 |
> |
int atomBC[3] = {15,11,4}; // block counts |
| 480 |
|
MPI_Aint atomDspls[3]; // displacements |
| 481 |
|
MPI_Datatype atomMbrTypes[3]; // member mpi types |
| 482 |
|
|
| 1310 |
|
} |
| 1311 |
|
|
| 1312 |
|
gBend = new GhostBend( *the_atoms[a], |
| 1313 |
< |
*the_atoms[b] ); |
| 1313 |
> |
*the_atoms[b]); |
| 1314 |
> |
|
| 1315 |
|
gBend->setConstants( currentBendType->k1, |
| 1316 |
|
currentBendType->k2, |
| 1317 |
|
currentBendType->k3, |
| 1327 |
|
currentBendType->k3, |
| 1328 |
|
currentBendType->t0 ); |
| 1329 |
|
bendArray[i] = qBend; |
| 1330 |
< |
} |
| 1330 |
> |
} |
| 1331 |
|
} |
| 1332 |
|
} |
| 1333 |
|
} |
| 1611 |
|
} |
| 1612 |
|
|
| 1613 |
|
info.d0 = atof( the_token ); |
| 1614 |
+ |
|
| 1615 |
+ |
info.k0=0.0; |
| 1616 |
|
} |
| 1617 |
|
else if( !strcmp( bondType, "harmonic" ) ){ |
| 1618 |
|
info.type = HARMONIC_BOND; |