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; |