collision test cleaning
[physics.git] / src / collisionManager.cpp
1 /*
2  *  Copyright (C) 2008 Patrik Gornicz, Gornicz_P (at) hotmail (dot) com.
3  *
4  *  This program is free software: you can redistribute it and/or modify
5  *  it under the terms of the GNU General Public License as published by
6  *  the Free Software Foundation, either version 3 of the License, or
7  *  (at your option) any later version.
8  *
9  *  This program is distributed in the hope that it will be useful,
10  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
11  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  *  GNU General Public License for more details.
13  *
14  *  You should have received a copy of the GNU General Public License
15  *  along with this program.  If not, see <http://www.gnu.org/licenses/>.
16  */
17
18 #include "collisionManager.h"
19 #include "debug.h"
20
21 #include "Vector2.h"
22
23 #include "Entities/Ball.h"
24 #include "Entities/PhysicsEntity.h"
25
26 #include "CollisionInfo.h"
27
28 /// ***** Private Class Header *****
29
30 /// ***** Private Method Headers *****
31
32 void applyCollisionAt(PhysicsEntity* p1, PhysicsEntity* p2);
33 void applyCollisionAt(Ball* b1, Ball* b2);
34
35 bool getInfoAt(Ball* b1, Ball* b2, CollisionInfo* cInfo);
36
37 /// ***** Private Variables *****
38
39 /// ***** Initializers/Cleaners *****
40
41 void collision::init()
42 {
43
44 }
45 void collision::clean()
46 {
47
48 }
49
50 /// ***** Public Methods *****
51
52 void collision::update(setPhys& sp)
53 {
54     for( setPhys::iterator it1 = sp.begin();
55          it1 != sp.end();
56          it1++ )
57     {
58         for( setPhys::iterator it2 = sp.begin();
59                 it2 != sp.end();
60                 it2++ )
61         {
62             if( *it1 != *it2 )
63             {
64                 applyCollisionAt(*it1, *it2);
65             }
66         }
67     }
68
69 }
70
71 /// ***** Private Methods *****
72
73 void applyCollisionAt(PhysicsEntity* p1, PhysicsEntity* p2)
74 {
75     Ball* b1 = dynamic_cast<Ball*>(p1);
76     Ball* b2 = dynamic_cast<Ball*>(p2);
77
78     if( b1 != NULL && b2 != NULL )
79     {
80         applyCollisionAt(b1, b2);
81         return;
82     }
83
84     DPF(0, "ENTITY TYPE NOT SUPPORTED BY applyCollisionAt()!!");
85 }
86
87 void applyCollisionAt(Ball* b1, Ball* b2)
88 {
89     CollisionInfo   cInfo;
90
91     if(!getInfoAt(b1, b2, &cInfo))
92         return;
93
94     // a few values to simplify the equations
95     const Vector2& normal = cInfo.normal;
96     const Vector2& point  = cInfo.point;
97
98     float m1 = b1->mass;
99     float m2 = b2->mass;
100
101     float e = (b1->CoR + b2->CoR) / 2;
102
103     Vector2 v1 = b1->velocityRaw();
104     Vector2 v2 = b2->velocityRaw();
105
106
107     float iTop = -(e + 1) * (v1 - v2).dot(normal);
108
109     // otherwise the collision happened and we do the math the below assumes
110     // collisions have no friction
111
112     float impulse = iTop / (normal.dot(normal) * (1 / m1 + 1 / m2));
113
114     b1->applyImpulse(impulse / m1 * normal, point);
115     b2->applyImpulse(-impulse / m2 * normal, point);
116 }
117
118 bool getInfoAt(Ball* b1, Ball* b2, CollisionInfo* pcInfo)
119 {
120     // a few values to simplify the equations
121     float r1 = b1->radius;
122     float r2 = b2->radius;
123
124     Vector2 p1 = b1->positionRaw();
125     Vector2 p2 = b2->positionRaw();
126     Vector2 v1 = b1->velocityRaw();
127     Vector2 v2 = b2->velocityRaw();
128
129     // quick binding box check
130     if (p1.x - r1 > p2.x + r2
131      || p1.x + r1 < p2.x - r2
132      || p1.y - r1 > p2.y + r2
133      || p1.y + r1 < p2.y - r2)
134         return false;
135
136     // test if not touching
137     if ((p1 - p2).sqrLength() >= (r1 + r2)*(r1 + r2))
138         return false;
139
140     // test if they are moving apart in some way if they aren't it's likely
141     // that they collided last frame and are still overlapping
142
143     if ((v1 - v2).dot(p1 - p2) >= 0)
144         return false;
145
146     pcInfo->point   = p1 - (p1 - p2) * r1 / (r1 + r2);
147     pcInfo->normal  = p1 - p2;
148
149     return true;
150 }