--- trunk/OOPSE/libmdtools/GhostBend.cpp 2003/08/20 19:11:51 704 +++ trunk/OOPSE/libmdtools/GhostBend.cpp 2003/08/20 19:42:31 707 @@ -9,10 +9,10 @@ GhostBend::GhostBend( Atom &a, Atom &b, Atom &c ){ -GhostBend::GhostBend( Atom &a, Atom &b, Atom &c ){ - - set_atoms(a, b, c); +GhostBend::GhostBend( Atom &a, Atom &b ){ + c_p_a = &a; + if( !b.isDirectional() ){ // if atom b is not directional, then bad things will happen @@ -52,7 +52,7 @@ void GhostBend::calc_forces(){ double aF[3], bF[3], bTrq[3]; c_p_a->getPos( aR ); - c_p_b->getPos( bR ); + atomB->getPos( bR ); dx = aR[0] - bR[0];