~cypheon/elfelli

3c6f8765ea525ce1d8c04b92ebd2c83f39a7bbab — Johann Rudloff 13 years ago e930745
Fix profiling switch in SConstruct. Speed up simulation by avoiding double calculation.
3 files changed, 16 insertions(+), 16 deletions(-)

M SConstruct
M src/Simulation.cpp
M src/Simulation.h
M SConstruct => SConstruct +3 -3
@@ 40,12 40,12 @@ if not conf.PkgConfig('gtkmm-2.4', '2.8'):
env.AppendUnique(CCFLAGS=['-Wall'])

if env['profiling']:
	env.AppendUnique(CPPDEFINES='PROFILING')
	env.AppendUnique(CPPDEFINES=['PROFILING'])

if env['debug']:
	env.AppendUnique(CCFLAGS=['-g', '-O0'], CPPDEFINES='DEBUG')
	env.AppendUnique(CCFLAGS=['-g', '-O0'], CPPDEFINES=['DEBUG'])
else:
	env.AppendUnique(CCFLAGS=['-O3'], CPPDEFINES='NDEBUG')
	env.AppendUnique(CCFLAGS=['-O3'], CPPDEFINES=['NDEBUG'])

paths = {"bindir": env['prefix'] + '/bin',
         "datadir": env['prefix'] + '/share/elfelli',

M src/Simulation.cpp => src/Simulation.cpp +11 -11
@@ 143,21 143,21 @@ Vec2 Simulation::force_at(const Vec2& pos, float charge)
    {
      PlateBody& plate = plates[i];

      float ax, ay, bx, by;
      float xa, ya, xb, yb;
      Vec2 v;
      ax = plate.pos_a.get_x() - pos.get_x();
      ay = plate.pos_a.get_y() - pos.get_y();
      bx = plate.pos_b.get_x() - plate.pos_a.get_x();
      by = plate.pos_b.get_y() - plate.pos_a.get_y();
      xa = plate.pos_a.get_x() - pos.get_x();
      ya = plate.pos_a.get_y() - pos.get_y();
      xb = plate.pos_b.get_x() - plate.pos_a.get_x();
      yb = plate.pos_b.get_y() - plate.pos_a.get_y();

      v.set_x((2*by*atan((by*by+ay*by+bx*bx+ax*bx)/(ax*by-ay*bx))+bx*log(by*by+2*ay*by+bx*bx+2*ax*bx+ay*ay+ax*ax)-2*by*
               atan((ay*by+ax*bx)/(ax*by-ay*bx))-log(ay*ay+ax*ax)*bx)/((2*by*by+2*bx*bx)*(atan((by*by+ay*by+bx*bx+ax*bx)/(ax*by-ay*bx))/(ax*by-ay*bx)-atan((ay*by+ax*bx)/(ax*by-ay*bx))/(ax*by-ay*bx))));

      v.set_y((-2*bx*atan((by*by+ay*by+bx*bx+ax*bx)/(ax*by-ay*bx))+by*log(by*by+2*ay*by+bx*bx+2*ax*bx+ay*ay+ax*ax)+2*bx*
               atan((ay*by+ax*bx)/(ax*by-ay*bx))-log(ay*ay+ax*ax)*by)/((2*by*by+2*bx*bx)*(atan((by*by+ay*by+bx*bx+ax*bx)/(ax*by-ay*bx))/(ax*by-ay*bx)-atan((ay*by+ax*bx)/(ax*by-ay*bx))/(ax*by-ay*bx))));
      float u = (-(2*ya*yb+2*xa*xb)*atan((yb*yb+ya*yb+xb*xb+xa*xb)/(xa*yb-xb*ya))-(xb*ya-xa*yb)*
                 log(yb*yb+2*ya*yb+ya*ya+xb*xb+2*xa*xb+xa*xa)-(-2*ya*yb-2*xa*xb)*atan((ya*yb+xa*xb)/(xa*yb-xb*ya))-xa*log(ya*ya+xa*xa)*
                 yb+xb*ya*log(ya*ya+xa*xa))/((2*yb*yb+2*xb*xb)*atan((yb*yb+ya*yb+xb*xb+xa*xb)/(xa*yb-xb*ya))+(-2*yb*yb-2*xb*xb)*atan((ya*yb+xa*xb)/(xa*yb-xb*ya)));
      v.set_x(xa+u*xb);
      v.set_y(ya+u*yb);

      float dist = v.length();
      float length = sqrt(bx*bx+by*by)/30;
      float length = sqrt(xb*xb+yb*yb)/30;
      Vec2 t = (v.normalize())/(dist*dist);

      f -= t * (charge * plate.charge/length);

M src/Simulation.h => src/Simulation.h +2 -2
@@ 105,8 105,8 @@ public:
  void add_body(const Vec2& v, float charge);
  void add_plate(const Vec2& a, const Vec2& b, float charge);

  const std::vector<Body> get_bodies() const{return bodies;};
  const std::vector<PlateBody> get_plates() const{return plates;};
  const std::vector<Body>& get_bodies() const{return bodies;};
  const std::vector<PlateBody>& get_plates() const{return plates;};

  Body& operator[](int n){return bodies[n];};
  int n_bodies(){return bodies.size();};