Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/Actions/WorldAction/CenterOnEdgeAction.cpp

    r84c494 r623e89  
    66 */
    77
    8 #include "Helpers/MemDebug.hpp"
    9 
    108#include "Actions/WorldAction/CenterOnEdgeAction.hpp"
    119#include "atom.hpp"
     
    1311#include "vector.hpp"
    1412#include "World.hpp"
    15 #include "Matrix.hpp"
    1613
    1714#include <iostream>
     
    3835  Vector Min;
    3936  Vector Max;
     37  int j=0;
    4038
    4139  dialog->queryEmpty(NAME, MapOfActions::getInstance().getDescription(NAME));
     
    5755    }
    5856    // set new box size
    59     Matrix domain;
     57    double * const cell_size = World::getInstance().getDomain();
     58    for (j=0;j<6;j++)
     59      cell_size[j] = 0.;
     60    j=-1;
    6061    for (int i=0;i<NDIM;i++) {
    61       double tmp = Max[i]-Min[i];
    62       tmp = fabs(tmp)>=1. ? tmp : 1.0;
    63       domain.at(i,i) = tmp;
     62      j += i+1;
     63      cell_size[j] = (Max[i]-Min[i]);
    6464    }
    65     cout << "new domain is: " << domain << endl;
    66     World::getInstance().setDomain(domain);
    6765    // translate all atoms, such that Min is aty (0,0,0)
    6866    for (vector<atom*>::iterator AtomRunner = AllAtoms.begin(); AtomRunner != AllAtoms.end(); ++AtomRunner)
Note: See TracChangeset for help on using the changeset viewer.