| # | Line 9 | Line 9 | GhostBend::GhostBend( Atom &a, Atom &b, Atom &c ){ | |
|---|---|---|
| 9 | ||
| 10 | ||
| 11 | ||
| 12 | < | GhostBend::GhostBend( Atom &a, Atom &b, Atom &c ){ |
| 13 | < | |
| 14 | < | set_atoms(a, b, c); |
| 12 | > | GhostBend::GhostBend( Atom &a, Atom &b ){ |
| 13 | ||
| 14 | + | c_p_a = &a; |
| 15 | + | |
| 16 | if( !b.isDirectional() ){ | |
| 17 | ||
| 18 | // if atom b is not directional, then bad things will happen | |
| # | Line 52 | Line 52 | void GhostBend::calc_forces(){ | |
| 52 | double aF[3], bF[3], bTrq[3]; | |
| 53 | ||
| 54 | c_p_a->getPos( aR ); | |
| 55 | < | c_p_b->getPos( bR ); |
| 55 | > | atomB->getPos( bR ); |
| 56 | ||
| 57 | ||
| 58 | dx = aR[0] - bR[0]; | |
| – | Removed lines |
| + | Added lines |
| < | Changed lines |
| > | Changed lines |