Main Page | Namespace List | Class Hierarchy | Alphabetical List | Class List | File List | Class Members | File Members

BondElem Class Reference

#include <ComputeBonds.h>

List of all members.

Public Types

enum  { size = 2 }
enum  { bondEnergyIndex, virialIndex, reductionDataSize }
enum  { reductionChecksumLabel = REDUCTION_BOND_CHECKSUM }

Public Member Functions

void computeForce (BigReal *, BigReal *)
int hash () const
 BondElem ()
 BondElem (AtomID atom0, const TupleSignature *sig, const BondValue *v)
 BondElem (const Bond *a, const BondValue *v)
 BondElem (AtomID atom0, AtomID atom1)
 ~BondElem ()
int operator== (const BondElem &a) const
int operator< (const BondElem &a) const

Static Public Member Functions

void loadTuplesForAtom (void *, AtomID, Molecule *)
void getMoleculePointers (Molecule *, int *, int32 ***, Bond **)
void getParameterPointers (Parameters *, const BondValue **)
void getTupleInfo (AtomSignature *sig, int *count, TupleSignature **t)
void submitReductionData (BigReal *, SubmitReduction *)

Public Attributes

AtomID atomID [size]
int localIndex [size]
TuplePatchElem * p [size]
Real scale
const BondValue * value

Static Public Attributes

int pressureProfileSlabs = 0
int pressureProfileAtomTypes = 1
BigReal pressureProfileThickness = 0
BigReal pressureProfileMin = 0

Member Enumeration Documentation

anonymous enum
 

Enumeration values:
size 

Definition at line 23 of file ComputeBonds.h.

00023 { size = 2 };

anonymous enum
 

Enumeration values:
bondEnergyIndex 
virialIndex 
reductionDataSize 

Definition at line 49 of file ComputeBonds.h.

00049 { bondEnergyIndex, TENSOR(virialIndex), reductionDataSize };

anonymous enum
 

Enumeration values:
reductionChecksumLabel 

Definition at line 50 of file ComputeBonds.h.

00050 { reductionChecksumLabel = REDUCTION_BOND_CHECKSUM };

Constructor & Destructor Documentation

BondElem::BondElem (   )  [inline]
 

Copyright (c) 1995, 1996, 1997, 1998, 1999, 2000 by The Board of Trustees of the University of Illinois. All rights reserved.

Definition at line 12 of file ComputeBonds.inl.

00012 { ; }

BondElem::BondElem (  AtomID  atom0,
const TupleSignature *  sig,
const BondValue *  v
)  [inline]
 

Definition at line 14 of file ComputeBonds.inl.

References atomID, TupleSignature::offset, TupleSignature::tupleParamType, and value.

00014                                                                                     {
00015     atomID[0] = atom0;
00016     atomID[1] = atom0 + sig->offset[0];
00017     value = &v[sig->tupleParamType];
00018 }

BondElem::BondElem (  const Bond *  a,
const BondValue *  v
)  [inline]
 

Definition at line 20 of file ComputeBonds.inl.

References bond::atom1, bond::atom2, atomID, Bond, bond::bond_type, and value.

00021   {
00022     atomID[0] = a->atom1;
00023     atomID[1] = a->atom2;
00024     value = &v[a->bond_type];
00025   }

BondElem::BondElem (  AtomID  atom0,
AtomID  atom1
)  [inline]
 

Definition at line 27 of file ComputeBonds.inl.

References atomID, and AtomID.

00028   {
00029     if (atom0 > atom1) {  // Swap end atoms so lowest is first!
00030       AtomID tmp = atom1; atom1 = atom0; atom0 = tmp; 
00031     }
00032     atomID[0] = atom0;
00033     atomID[1] = atom1;
00034   }

BondElem::~BondElem (   )  [inline]
 

Definition at line 57 of file ComputeBonds.h.

00057 {};

Member Function Documentation

void BondElem::computeForce (  BigReal *  ,
BigReal * 
) 
 

Definition at line 65 of file ComputeBonds.C.

References BigReal, DebugM, Lattice::delta(), TuplePatchElem::f, Force, BondValue::k, Patch::lattice, Vector::length(), localIndex, TuplePatchElem::p, p, CompAtom::partition, CompAtom::position, pp_clamp(), pp_reduction(), pressureProfileMin, pressureProfileSlabs, pressureProfileThickness, Real, value, TuplePatchElem::x, Vector::x, BondValue::x0, Vector::y, and Vector::z.

00067 {
00068   DebugM(1, "::computeForce() localIndex = " << localIndex[0] << " "
00069                << localIndex[1] << std::endl);
00070 
00071   BigReal r;            // Distance between atoms
00072   BigReal diff;         // difference between theta and theta0
00073   BigReal energy;       // energy from the bond
00074 
00075   //DebugM(3, "::computeForce() -- starting with bond type " << bondType << std::endl);
00076 
00077   // get the bond information
00078   Real k = value->k * scale;
00079   Real x0 = value->x0;
00080 
00081   // compute vectors between atoms and their distances
00082   const Lattice & lattice = p[0]->p->lattice;
00083   const Vector r12 = lattice.delta(p[0]->x[localIndex[0]].position,
00084                                         p[1]->x[localIndex[1]].position);
00085   r = r12.length();
00086 
00087   //  Compare it to the rest bond
00088   diff = r - x0;
00089 
00090   //  Add the energy from this bond to the total energy
00091   energy = k*diff*diff;
00092 
00093   //  Determine the magnitude of the force
00094   diff *= -2.0*k;
00095 
00096   //  Scale the force vector accordingly
00097   const Force f12 = r12 * (diff/r);
00098 
00099   //  Now add the forces to each force vector
00100   p[0]->f[localIndex[0]] += f12;
00101   p[1]->f[localIndex[1]] -= f12;
00102 
00103   DebugM(3, "::computeForce() -- ending with delta energy " << energy << std::endl);
00104   reduction[bondEnergyIndex] += energy;
00105   reduction[virialIndex_XX] += f12.x * r12.x;
00106   reduction[virialIndex_XY] += f12.x * r12.y;
00107   reduction[virialIndex_XZ] += f12.x * r12.z;
00108   reduction[virialIndex_YX] += f12.y * r12.x;
00109   reduction[virialIndex_YY] += f12.y * r12.y;
00110   reduction[virialIndex_YZ] += f12.y * r12.z;
00111   reduction[virialIndex_ZX] += f12.z * r12.x;
00112   reduction[virialIndex_ZY] += f12.z * r12.y;
00113   reduction[virialIndex_ZZ] += f12.z * r12.z;
00114 
00115   if (pressureProfileData) {
00116     BigReal z1 = p[0]->x[localIndex[0]].position.z;
00117     BigReal z2 = p[1]->x[localIndex[1]].position.z;
00118     int n1 = (int)floor((z1-pressureProfileMin)/pressureProfileThickness);
00119     int n2 = (int)floor((z2-pressureProfileMin)/pressureProfileThickness);
00120     pp_clamp(n1, pressureProfileSlabs);
00121     pp_clamp(n2, pressureProfileSlabs);
00122     int p1 = p[0]->x[localIndex[0]].partition;
00123     int p2 = p[1]->x[localIndex[1]].partition;
00124     int pn = pressureProfileAtomTypes;
00125     pp_reduction(pressureProfileSlabs,
00126                 n1, n2, 
00127                 p1, p2, pn, 
00128                 f12.x * r12.x, f12.y * r12.y, f12.z * r12.z,
00129                 pressureProfileData);
00130   } 
00131 }

void BondElem::getMoleculePointers (  Molecule *  ,
int *  ,
int32 ***  ,
Bond ** 
)  [static]
 

Definition at line 50 of file ComputeBonds.C.

References Bond, int32, and NAMD_die().

00051 {
00052 #ifdef MEM_OPT_VERSION
00053   NAMD_die("Should not be called in BondElem::getMoleculePointers in memory optimized version!");
00054 #else
00055   *count = mol->numBonds;
00056   *byatom = mol->bondsByAtom;
00057   *structarray = mol->bonds;
00058 #endif
00059 }

void BondElem::getParameterPointers (  Parameters *  ,
const BondValue ** 
)  [static]
 

Definition at line 61 of file ComputeBonds.C.

References Parameters::bond_array.

00061                                                                       {
00062   *v = p->bond_array;
00063 }

void BondElem::getTupleInfo (  AtomSignature *  sig,
int *  count,
TupleSignature **  t
)  [inline, static]
 

Definition at line 33 of file ComputeBonds.h.

References AtomSignature::bondCnt, and AtomSignature::bondSigs.

00033                                                                                  {
00034         *count = sig->bondCnt;
00035         *t = sig->bondSigs;
00036     }

int BondElem::hash (   )  const [inline]
 

Definition at line 44 of file ComputeBonds.h.

00044 { return 0x7FFFFFFF & ( (atomID[0]<<16) + (atomID[1])); }

void BondElem::loadTuplesForAtom (  void *  ,
AtomID  ,
Molecule * 
)  [static]
 

int BondElem::operator< (  const BondElem &  a  )  const [inline]
 

Definition at line 41 of file ComputeBonds.inl.

References atomID.

00042   {
00043     return (atomID[0] < a.atomID[0] ||
00044             (atomID[0] == a.atomID[0] &&
00045             (atomID[1] < a.atomID[1]) ));
00046   }

int BondElem::operator== (  const BondElem &  a  )  const [inline]
 

Definition at line 36 of file ComputeBonds.inl.

References atomID.

00037   {
00038     return (a.atomID[0] == atomID[0] && a.atomID[1] == atomID[1]);
00039   }

void BondElem::submitReductionData (  BigReal *  ,
SubmitReduction * 
)  [static]
 

Definition at line 133 of file ComputeBonds.C.

References ADD_TENSOR, SubmitReduction::item(), REDUCTION_BOND_ENERGY, REDUCTION_VIRIAL_NORMAL, and virialIndex.

00134 {
00135   reduction->item(REDUCTION_BOND_ENERGY) += data[bondEnergyIndex];
00136   ADD_TENSOR(reduction,REDUCTION_VIRIAL_NORMAL,data,virialIndex);
00137 }

Member Data Documentation

AtomID BondElem::atomID[size]
 

Definition at line 24 of file ComputeBonds.h.

Referenced by BondElem(), operator<(), and operator==().

int BondElem::localIndex[size]
 

Definition at line 25 of file ComputeBonds.h.

Referenced by computeForce().

TuplePatchElem* BondElem::p[size]
 

Definition at line 26 of file ComputeBonds.h.

Referenced by computeForce().

int BondElem::pressureProfileAtomTypes = 1 [static]
 

Definition at line 45 of file ComputeBonds.C.

BigReal BondElem::pressureProfileMin = 0 [static]
 

Definition at line 47 of file ComputeBonds.C.

Referenced by computeForce().

int BondElem::pressureProfileSlabs = 0 [static]
 

Copyright (c) 1995, 1996, 1997, 1998, 1999, 2000 by The Board of Trustees of the University of Illinois. All rights reserved.

Definition at line 44 of file ComputeBonds.C.

Referenced by computeForce().

BigReal BondElem::pressureProfileThickness = 0 [static]
 

Definition at line 46 of file ComputeBonds.C.

Referenced by computeForce().

Real BondElem::scale
 

Definition at line 27 of file ComputeBonds.h.

const BondValue* BondElem::value
 

Definition at line 47 of file ComputeBonds.h.

Referenced by BondElem(), and computeForce().

The documentation for this class was generated from the following files:
Generated on Sun Jul 6 04:07:44 2008 for NAMD by  doxygen 1.3.9.1