1

#include "PhysicsObject.hh"

     2

#include "Engine.hh"

     3

     4

#include <cmath>

     5

#include <utility>

     6

#include <queue>

     7

     8

PhysicsObject::PhysicsObject (PhysicsWorld &world, float mass, Vector position, Vector velocity, ObjectType type, 

     9

            float collision_elasticity, bool enabled) :

    22

{  

    23

    if (enabled)

    24

        enable();  

    25

}

    26

    27

PhysicsObject::~PhysicsObject (void) {

    28

//    Engine::log(DEBUG, "PhysicsObject.destructor") << this /* << ": objects.size=" << ((int) world.objects.size()) */;

    29

}

    30

    31

/**

    32

 * Player walks on floor.

    33

 */

    34

Vector PhysicsObject::walk_one_step (float partial, bool right) {

    35

    // which way we are walking

    36

    float deltaX = right ? partial : -partial;

    37

    Vector reached = this->position;

    38

    if(reached.roundToInt() == (reached+Vector(deltaX, 0)).roundToInt()) {

    39

        return reached+Vector(deltaX, 0);

    40

    }

    41

    // Is there upward ramp

    42

    if(!possibleLocation(position+Vector(deltaX, 0))) {

    43

        // Yes. Then we check n pixels up

    44

        for(int i = 1; i < 3; i++) {

    45

            if(possibleLocation(position+Vector(deltaX, -i))) {

    46

                // and when there is finally EMPTY, we can walk

    47

                reached = position+Vector(deltaX, -i);

    48

                break;

    49

            }

    50

        }

    51

    } else {

    52

        // Or downward ramp or flat

    53

        for(int i = 0; 1; i++) {

    54

    55

            // And when there is finally ground we can step on

    56

            // it. If there is no gound we still step there,

    57

            // but will fall one pixel down

    58

            if(possibleLocation(position+Vector(deltaX, i))) {

    59

                reached = position+Vector(deltaX, i);

    60

            } else {

    61

                break;

    62

            }

    63

    64

            // If the fall is big enough, set the worm in the air

    65

            if (i >= 2) {

    66

//                Vector back = walk(dt, !right);

    67

                this->inAir = true;

    68

//                this->velocity.x = right ? velocity : -velocity;

    69

                // Avoid stepping two pixels down when it starts to free fall

    70

                reached.y -= 2;

    71

//                this->velocity = (reached-back)*1000/dt;

    72

                break;

    73

            }

    74

        }

    75

    }

    76

    // And we return where we got

    77

    return reached;

    78

}

    79

void PhysicsObject::walk (TimeMS dt, bool right) {

    80

    float velocity = PLAYER_WALK_SPEED;

    81

    float walkAmount = (velocity*dt)/1000;

    82

    Vector reached = this->position;

    83

    while (walkAmount > 0 && !this->inAir) {

    84

        this->position = walk_one_step((1 < walkAmount ? 1 : walkAmount), right);

    85

        walkAmount--;

    86

    }

    87

    // TODO: Should the remaining walkAmount be handled somehow?

    88

}

    89

    90

/**

    91

 * Makes the player jump in the air.

    92

 * @param direction -1: jump left, 0: jump up, 1: jump right

    93

 */

    94

void PhysicsObject::jump (int direction) {

    95

    // Jump only if player is "on the ground"

    96

    if (!this->inAir) {

    97

 	    velocity.y = -100;

    98

        switch (direction) {

    99

            case 1:

   100

                velocity.x += 20;

   101

                break;

   102

            case -1:

   103

                velocity.x -= 20;

   104

                break;

   105

            case 0:

   106

                break;

   107

            default:

   108

                throw std::logic_error("Invalid jump direction");

   109

        }

   110

	    inAir = true;

   111

    }

   112

}

   113

   114

bool PhysicsObject::possibleLocation (Vector loc) {

   115

    for(unsigned int i = 0; i < this->shape.size(); i++) {

   116

        if(world.collides(loc+shape[i]))

   117

            return false;

   118

    }

   119

    return true;

   120

}

   121

   122

/**

   123

 * Updates object speed and position. This function organises force

   124

 * integration and collision detection.

   125

 */   

   126

void PhysicsObject::updatePosition (TimeMS dt) {

   127

    // Add gravity to the force queue

   128

    applyForce(world.gravity);

   129

   130

    // If the object (practically player) has a pivot point add

   131

    // a force towards that

   132

    if (pivot != NULL) {

   133

        Vector direction(pivot->getPosition()-position);

   134

        float magnitude = pivot->getPivotForce(this);

   135

        float length = direction.length();

   136

        if (length > 0) {

   137

            direction = direction / length * magnitude;

   138

            applyForce(direction);

   139

        }

   140

    }

   141

   142

    std::pair<Force, TimeMS> force;

   143

    std::queue<std::pair<Force, TimeMS> > newfq;

   144

    Force total;

   145

    while (!forceq.empty()) {

   146

        force = forceq.front();

   147

        if (force.second > dt) {

   148

            force.second -= dt;

   149

            newfq.push(force);

   150

        }

   151

        total += force.first;

   152

        forceq.pop();

   153

    }

   154

    forceq = newfq;

   155

   156

    // If the player has stopped and there's some ground under some of the 3 some of the 3t

   157

    // set inAir false

   158

    if (this->velocity == Vector(0,0)) {

   159

        this->inAir = !world.collides(this->position+shape[1]+Vector(0, 1))

   160

                      && !world.collides(this->position+shape[2]+Vector(0, 1))

   161

                      && !world.collides(this->position+shape[3]+Vector(0, 1));

   162

        // If, however, there's a force caused by a bomb, e.g., set it in air.

   163

        // Still, we have to be able to separate forces caused by walking attempts

   164

        // and bombs etc (+0.1 because float comparison can be dangerous)

   166

            this->inAir = true;

   167

    }

   168

   169

    if(!possibleLocation(position)) {

   170

        //if we are trapped in ground form dirtball or something

   171

        //we might want to just return and set velocity to some value

   172

        //return;

   173

    }

   174

   175

    // If the worm is not in the air make it walk,

   176

    // otherwise integrate the new position and velocity

   177

    if (!this->inAir) {

   178

        // It walks only if there's some vertical force

   179

        if (total.x != 0) {

   180

            walk(dt, total.x > 0);

   181

            this->velocity = Vector(0,0);

   182

        }   

   183

        // Now the possible walking has been done so we can return from this function.

   184

        // In walk inAir could have been set true, but that will be handled in the next tick.

   185

        return;

   186

    }

   187

   188

    if (!possibleLocation(position))

   189

        Engine::log(DEBUG, "PhysicsObject.updatePosition") << "impossible location: " << position;

   190

   191

    Vector newPosition;

   192

    Vector velAfterTick;

   193

    // Calculate new position and velocity to the given references

   194

    integrate(total, dt, newPosition, velAfterTick);

   195

    this->velocity = velAfterTick;

   196

   197

    // Collision detection

   198

    bool collided = false;

   199

    Vector collisionPoint;

   200

   201

    const Vector diffVec = newPosition - position;

   202

    const Vector unitVector = diffVec / diffVec.length();

   203

    Vector reached = position;

   204

   205

    while ((position - reached).sqrLength() < diffVec.sqrLength()) {

   206

        reached += unitVector;

   207

   208

        // Check if any of the shapes points collide

   209

        for (uint64_t i = 0; i < shape.size(); i++) {

   210

            if (world.collides(reached+shape[i])) {  // Collision

   211

   212

                collisionPoint = reached+shape[i];

   213

   214

                if (inAir)

   215

                    this->bounce(world.getNormal(reached+shape[i], reached-unitVector+shape[i]));

   216

   217

                reached = reached - unitVector; // Return to last point

   218

                collided = true;

   219

   220

                if (this->velocity.sqrLength() < PLAYER_MIN_SPEED * PLAYER_MIN_SPEED)

   221

                    this->velocity = Vector(0,0);

   222

   223

                break;

   224

            }

   225

        }

   226

   227

        if (collided)

   228

            break;

   229

    }

   230

   231

   232

    if (!possibleLocation(reached))

   233

        Engine::log(DEBUG, "PhysicsObject.updatePosition") << "impossible location: " << position << ", diffVec=" << diffVec;

   234

   235

    // In case of some float error check the final coordinate

   236

    if (!collided) {

   237

        if (!possibleLocation(newPosition)) {

   238

            newPosition = reached;

   239

        } else {

   240

            // This means everything was ok, so no need to do anything

   241

        }

   242

   243

        this->position = newPosition;

   244

   245

    } else {

   246

        newPosition = reached;

   247

        this->position = newPosition;

   248

   249

        // the following may delete this object, so it must be the last thing called

   250

        onCollision(collisionPoint, NULL);

   251

   252

        return;

   253

    }

   254

}

   255

   256

/**

   257

 * Bounces from straight wall in any direction.

   258

 * Direction given as normal of that wall

   259

 */

   260

void PhysicsObject::bounce (Vector normal) {

   261

    // normal.sqrLength can't be 0 when got from getNormal()

   262

    if (normal.sqrLength() != 0) {

   263

        Vector nvel = velocity;

   264

        // We project the velocity on normal and remove twice that much from velocity

   265

        nvel = nvel - ((2)*((nvel*normal)/(normal*normal))*normal);

   266

        velocity = nvel;

   267

        // We lose some of our speed on collision

   268

        this->velocity *= this->collision_elasticity;

   269

    }

   270

}

   271

   272

/**

   273

 * Integrates given force over time and stores new position to

   274

 * posAfterTick and new velocity to velAfterTick.

   275

 * @param force Force vector.

   276

 * @param dt The time the force is applied (<=PHYSICS_TICK_MS)

   277

 */

   278

void PhysicsObject::integrate(Force force, TimeMS dt, Vector &posAfterTick, Vector &velAfterTick) {

   279

    posAfterTick = position;

   280

    velAfterTick = velocity;

   281

    Derivative tmpd;

   282

    Derivative k1 = evaluate(force, 0, tmpd, posAfterTick, velAfterTick);

   285

    Derivative k4 = evaluate(force, dt, k3, posAfterTick, velAfterTick);

   286

   287

   288

    const Vector dxdt = (k1.dx + (k2.dx + k3.dx) * 2.0f + k4.dx) * 1.0f/6.0f;

   289

    const Vector dvdt = (k1.dv + (k2.dv + k3.dv) * 2.0f + k4.dv) * 1.0f/6.0f;

   290

   291

    //    Engine::log(DEBUG, "PhysicsObject.integrate") << "Changes: "<< dxdt << " " << dvdt << " Time: " <<dt;

   292

    posAfterTick = posAfterTick + (dxdt * dt)/1000;

   293

    velAfterTick = velAfterTick + (dvdt * dt)/1000;

   294

    //Engine::log(DEBUG, "PhysicsObject.integrate") << "velAfterTick: " << velAfterTick;

   295

}

   296

   297

Derivative PhysicsObject::evaluate(Force force, TimeMS dt, Derivative &d, const Vector &posAfterTick, const Vector &velAfterTick) {

   298

    Vector curPos = posAfterTick + (d.dx*dt)/1000;

   299

    Vector curVel = velAfterTick + (d.dv*dt)/1000;

   300

   301

    Derivative out;

   302

    out.dx = curVel;

   303

    out.dv = acceleration(force);

   304

    //Engine::log(DEBUG, "PhysicsObject.evaluate") << "Out.dx: " << out.dx;

   305

    return out;

   306

}

   307

   308

Vector PhysicsObject::acceleration(const Force &force) {

   309

    return (force/mass);

   310

}

   311

   312

void PhysicsObject::applyForce (Force force, TimeMS dt) {

   313

    // Add applied force to the queue

   314

    forceq.push(std::make_pair(force, dt));

   315

}

   316

   317

void PhysicsObject::changeAim(float da) {

   318

    this->aim += da;

   319

   320

    if (this->aim > PLAYER_AIM_MAX) this->aim = PLAYER_AIM_MAX;

   321

    if (this->aim < PLAYER_AIM_MIN) this->aim = PLAYER_AIM_MIN;

   322

    //Engine::log(DEBUG, "PhysicsObject.changeAim") << "Player aim: " << this->aim;

   323

}

   324

   325

void PhysicsObject::setFacing (FacingDirection facing) {

   326

    this->facing = facing;

   327

}

   328

   329

void PhysicsObject::updatePhysics (Vector position, Vector velocity, bool inAir, FacingDirection facing, float aim) {

   330

    this->position = position;

   331

    this->velocity = velocity;

   332

    this->inAir = inAir;

   333

    this->facing = facing;

   334

    this->aim = aim;

   335

}

   336

   337

Vector PhysicsObject::getPosition (void) const {

   338

    return position;

   339

}

   340

   341

PixelCoordinate PhysicsObject::getCoordinate (void) const {

   342

    return world.getPixelCoordinate(position);

   343

}

   344

   345

Vector PhysicsObject::getVelocity (void) const {

   346

    return velocity;

   347

}

   348

   349

FacingDirection PhysicsObject::getFacing (void) const {

   350

    return facing;

   351

}

   352

   353

float PhysicsObject::getAim (void) const {

   354

    return aim;

   355

}

   356

   357

Vector PhysicsObject::getDirection (void) const {

   358

    return facing == FACING_RIGHT ? Vector(cos(aim), -sin(aim)) : Vector(-cos(aim), -sin(aim));

   359

}

   360

   361

const std::vector<Vector>& PhysicsObject::getShape () const {

   362

    return shape;

   363

}

   364

   365

void PhysicsObject::setShape (std::vector<Vector> shape) {

   366

    this->shape = shape;

   367

}

   368

   369

void PhysicsObject::setPivot (PhysicsObject *pivot) {

   370

    this->pivot = pivot;

   371

}

   372

   373

void PhysicsObject::tick (TimeMS tick_length) {

   374

    this->updatePosition(tick_length);

   375

}

   376

   377

void PhysicsObject::enable (void) {

   378

    // only enable once until disabled

   379

    if (alive)

   380

        return;

   381

   382

    // mark as alive

   383

    alive = true;

   384

   385

    // add the world objects list

   386

    world.addPhysicsObject(this);

   387

}

   388

   389

void PhysicsObject::disable (void) {

   390

    // mark as disabled

   391

    alive = false;

   392

}

   393

   394

void PhysicsObject::destroy (void) {

   395

    // mark as disabled and for deletion

   396

    alive = false;

   397

    shouldDelete = true;

   398

}

   399

   400

bool PhysicsObject::isDestroyed (void) {

   401

    return !alive;

   402

}

   403

   404

bool PhysicsObject::removeIfDestroyed (void) {

   405

    if (!alive) {

   406

        if (shouldDelete)

   407

            delete this;

   408

   409

        return true;

   410

   411

    } else {

   412

        return false;

   413

    }

   414

}

   415

   416

float PhysicsObject::getPivotForce (PhysicsObject *bob) { 

   419

    return 0.0; 

   420

}

   421

   422

bool PhysicsObject::collides (const PhysicsObject &obj) {

   423

    const std::vector<Vector> oShape = obj.getShape();

   424

    Vector p1, p2, p3;

   425

    int8_t sign, nsign;

   426

    for (std::vector<Vector>::const_iterator i = oShape.begin(); i != oShape.end(); i++) { // For every point in other shape

   427

        p3 = *i + obj.getPosition();

   428

        sign = 0;

   429

        for (std::vector<Vector>::const_iterator j = shape.begin(); j != shape.end(); j++) {

   430

            p1 = *j + position;

   431

            if ( (j+1) == shape.end() ) p2 = *shape.begin() + position;

   432

            else p2 = *(j+1) + position;

   433

            nsign = crossProduct(p1, p2, p3);

   434

            if ( sign == 0 ) sign = nsign;

   435

            else if ( ((sign < 0) && (nsign < 0)) || ((sign > 0) && (nsign > 0)) ) continue;

   436

            else return false;

   437

        }

   438

    }

   439

    return true;

   440

}

   446

   447

int8_t crossProduct (const Vector &p1, const Vector &p2, const Vector &p3) {

   448

    float p = (p2.x - p1.x)*(p3.y - p1.y) - (p2.y - p1.y)*(p3.x - p1.x);

   449

    if (p < 0)

   450

        return -1;

   451

    else

   452

        return 1;

   453

}