Ignore:
File:
1 edited

Legend:

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

    r84c494 rb9c847  
    1313#include "vector.hpp"
    1414#include "World.hpp"
    15 #include "Matrix.hpp"
    1615
    1716#include <iostream>
     
    3837  Vector Min;
    3938  Vector Max;
     39  int j=0;
    4040
    4141  dialog->queryEmpty(NAME, MapOfActions::getInstance().getDescription(NAME));
     
    5757    }
    5858    // set new box size
    59     Matrix domain;
     59    double * const cell_size = World::getInstance().getDomain();
     60    for (j=0;j<6;j++)
     61      cell_size[j] = 0.;
     62    j=-1;
    6063    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;
     64      j += i+1;
     65      cell_size[j] = (Max[i]-Min[i]);
    6466    }
    65     cout << "new domain is: " << domain << endl;
    66     World::getInstance().setDomain(domain);
     67    World::getInstance().setDomain(cell_size);
    6768    // translate all atoms, such that Min is aty (0,0,0)
    6869    for (vector<atom*>::iterator AtomRunner = AllAtoms.begin(); AtomRunner != AllAtoms.end(); ++AtomRunner)
Note: See TracChangeset for help on using the changeset viewer.