Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/molecule.cpp

    r2ba827 r46d958  
    88#include <boost/bind.hpp>
    99
     10#include "World.hpp"
    1011#include "atom.hpp"
    1112#include "bond.hpp"
     
    3031 * Initialises molecule list with correctly referenced start and end, and sets molecule::last_atom to zero.
    3132 */
    32 molecule::molecule(const periodentafel * const teil) : elemente(teil), start(new atom), end(new atom),
     33molecule::molecule(const periodentafel * const teil) : elemente(teil), start(World::get()->createAtom()), end(World::get()->createAtom()),
    3334  first(new bond(start, end, 1, -1)), last(new bond(start, end, 1, -1)), MDSteps(0), AtomCount(0),
    3435  BondCount(0), ElementCount(0), NoNonHydrogen(0), NoNonBonds(0), NoCyclicBonds(0), BondDistance(0.),
     
    6061  delete(first);
    6162  delete(last);
    62   delete(end);
    63   delete(start);
     63  end->getWorld()->destroyAtom(end);
     64  start->getWorld()->destroyAtom(start);
    6465};
    6566
     
    135136  OBSERVE;
    136137  if (pointer != NULL) {
    137     atom *walker = new atom(pointer);
     138    atom *walker = pointer->clone();
    138139    walker->Name = Malloc<char>(strlen(pointer->Name) + 1, "atom::atom: *Name");
    139140    strcpy (walker->Name, pointer->Name);
     
    242243  switch(TopBond->BondDegree) {
    243244    case 1:
    244       FirstOtherAtom = new atom();    // new atom
     245      FirstOtherAtom = World::get()->createAtom();    // new atom
    245246      FirstOtherAtom->type = elemente->FindElement(1);  // element is Hydrogen
    246247      FirstOtherAtom->v.CopyVector(&TopReplacement->v); // copy velocity
     
    299300
    300301      // create the two Hydrogens ...
    301       FirstOtherAtom = new atom();
    302       SecondOtherAtom = new atom();
     302      FirstOtherAtom = World::get()->createAtom();
     303      SecondOtherAtom = World::get()->createAtom();
    303304      FirstOtherAtom->type = elemente->FindElement(1);
    304305      SecondOtherAtom->type = elemente->FindElement(1);
     
    354355    case 3:
    355356      // take the "usual" tetraoidal angle and add the three Hydrogen in direction of the bond (height of the tetraoid)
    356       FirstOtherAtom = new atom();
    357       SecondOtherAtom = new atom();
    358       ThirdOtherAtom = new atom();
     357      FirstOtherAtom = World::get()->createAtom();
     358      SecondOtherAtom = World::get()->createAtom();
     359      ThirdOtherAtom = World::get()->createAtom();
    359360      FirstOtherAtom->type = elemente->FindElement(1);
    360361      SecondOtherAtom->type = elemente->FindElement(1);
     
    475476    MDSteps++;
    476477  for(i=0;i<NumberOfAtoms;i++){
    477     Walker = new atom;
     478    Walker = World::get()->createAtom();
    478479    getline(xyzfile,line,'\n');
    479480    istringstream *item = new istringstream(line);
Note: See TracChangeset for help on using the changeset viewer.