summaryrefslogtreecommitdiff
path: root/servers/physics_2d/joints_2d_sw.cpp
diff options
context:
space:
mode:
authorPouleyKetchoupp <pouleyketchoup@gmail.com>2021-10-18 12:24:30 -0700
committerPouleyKetchoupp <pouleyketchoup@gmail.com>2021-10-18 17:01:10 -0700
commitcc39dca9f7d960d1bb137f1dcbbf1da5cec8a505 (patch)
tree555c5721844576a85183de9f1c33cd19c671f73c /servers/physics_2d/joints_2d_sw.cpp
parent5bb3dbbedd4f32974eef36ffc83bbe29abb65ab1 (diff)
Rename Godot Physics classes from *SW to Godot*
Also moved MT physics server wrappers to the main servers folder, since they don't have to be implementation specific.
Diffstat (limited to 'servers/physics_2d/joints_2d_sw.cpp')
-rw-r--r--servers/physics_2d/joints_2d_sw.cpp486
1 files changed, 0 insertions, 486 deletions
diff --git a/servers/physics_2d/joints_2d_sw.cpp b/servers/physics_2d/joints_2d_sw.cpp
deleted file mode 100644
index b46397b8e6..0000000000
--- a/servers/physics_2d/joints_2d_sw.cpp
+++ /dev/null
@@ -1,486 +0,0 @@
-/*************************************************************************/
-/* joints_2d_sw.cpp */
-/*************************************************************************/
-/* This file is part of: */
-/* GODOT ENGINE */
-/* https://godotengine.org */
-/*************************************************************************/
-/* Copyright (c) 2007-2021 Juan Linietsky, Ariel Manzur. */
-/* Copyright (c) 2014-2021 Godot Engine contributors (cf. AUTHORS.md). */
-/* */
-/* Permission is hereby granted, free of charge, to any person obtaining */
-/* a copy of this software and associated documentation files (the */
-/* "Software"), to deal in the Software without restriction, including */
-/* without limitation the rights to use, copy, modify, merge, publish, */
-/* distribute, sublicense, and/or sell copies of the Software, and to */
-/* permit persons to whom the Software is furnished to do so, subject to */
-/* the following conditions: */
-/* */
-/* The above copyright notice and this permission notice shall be */
-/* included in all copies or substantial portions of the Software. */
-/* */
-/* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, */
-/* EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF */
-/* MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.*/
-/* IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY */
-/* CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, */
-/* TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE */
-/* SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. */
-/*************************************************************************/
-
-#include "joints_2d_sw.h"
-
-#include "space_2d_sw.h"
-
-//based on chipmunk joint constraints
-
-/* Copyright (c) 2007 Scott Lembcke
- *
- * Permission is hereby granted, free of charge, to any person obtaining a copy
- * of this software and associated documentation files (the "Software"), to deal
- * in the Software without restriction, including without limitation the rights
- * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
- * copies of the Software, and to permit persons to whom the Software is
- * furnished to do so, subject to the following conditions:
- *
- * The above copyright notice and this permission notice shall be included in
- * all copies or substantial portions of the Software.
- *
- * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
- * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
- * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
- * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
- * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
- * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
- * SOFTWARE.
- */
-
-void Joint2DSW::copy_settings_from(Joint2DSW *p_joint) {
- set_self(p_joint->get_self());
- set_max_force(p_joint->get_max_force());
- set_bias(p_joint->get_bias());
- set_max_bias(p_joint->get_max_bias());
- disable_collisions_between_bodies(p_joint->is_disabled_collisions_between_bodies());
-}
-
-static inline real_t k_scalar(Body2DSW *a, Body2DSW *b, const Vector2 &rA, const Vector2 &rB, const Vector2 &n) {
- real_t value = 0.0;
-
- {
- value += a->get_inv_mass();
- real_t rcn = (rA - a->get_center_of_mass()).cross(n);
- value += a->get_inv_inertia() * rcn * rcn;
- }
-
- if (b) {
- value += b->get_inv_mass();
- real_t rcn = (rB - b->get_center_of_mass()).cross(n);
- value += b->get_inv_inertia() * rcn * rcn;
- }
-
- return value;
-}
-
-static inline Vector2
-relative_velocity(Body2DSW *a, Body2DSW *b, Vector2 rA, Vector2 rB) {
- Vector2 sum = a->get_linear_velocity() - (rA - a->get_center_of_mass()).orthogonal() * a->get_angular_velocity();
- if (b) {
- return (b->get_linear_velocity() - (rB - b->get_center_of_mass()).orthogonal() * b->get_angular_velocity()) - sum;
- } else {
- return -sum;
- }
-}
-
-static inline real_t
-normal_relative_velocity(Body2DSW *a, Body2DSW *b, Vector2 rA, Vector2 rB, Vector2 n) {
- return relative_velocity(a, b, rA, rB).dot(n);
-}
-
-bool PinJoint2DSW::setup(real_t p_step) {
- dynamic_A = (A->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
- dynamic_B = (B->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
-
- if (!dynamic_A && !dynamic_B) {
- return false;
- }
-
- Space2DSW *space = A->get_space();
- ERR_FAIL_COND_V(!space, false);
-
- rA = A->get_transform().basis_xform(anchor_A);
- rB = B ? B->get_transform().basis_xform(anchor_B) : anchor_B;
-
- real_t B_inv_mass = B ? B->get_inv_mass() : 0.0;
-
- Transform2D K1;
- K1[0].x = A->get_inv_mass() + B_inv_mass;
- K1[1].x = 0.0f;
- K1[0].y = 0.0f;
- K1[1].y = A->get_inv_mass() + B_inv_mass;
-
- Transform2D K2;
- K2[0].x = A->get_inv_inertia() * rA.y * rA.y;
- K2[1].x = -A->get_inv_inertia() * rA.x * rA.y;
- K2[0].y = -A->get_inv_inertia() * rA.x * rA.y;
- K2[1].y = A->get_inv_inertia() * rA.x * rA.x;
-
- Transform2D K;
- K[0] = K1[0] + K2[0];
- K[1] = K1[1] + K2[1];
-
- if (B) {
- Transform2D K3;
- K3[0].x = B->get_inv_inertia() * rB.y * rB.y;
- K3[1].x = -B->get_inv_inertia() * rB.x * rB.y;
- K3[0].y = -B->get_inv_inertia() * rB.x * rB.y;
- K3[1].y = B->get_inv_inertia() * rB.x * rB.x;
-
- K[0] += K3[0];
- K[1] += K3[1];
- }
-
- K[0].x += softness;
- K[1].y += softness;
-
- M = K.affine_inverse();
-
- Vector2 gA = rA + A->get_transform().get_origin();
- Vector2 gB = B ? rB + B->get_transform().get_origin() : rB;
-
- Vector2 delta = gB - gA;
-
- bias = delta * -(get_bias() == 0 ? space->get_constraint_bias() : get_bias()) * (1.0 / p_step);
-
- return true;
-}
-
-inline Vector2 custom_cross(const Vector2 &p_vec, real_t p_other) {
- return Vector2(p_other * p_vec.y, -p_other * p_vec.x);
-}
-
-bool PinJoint2DSW::pre_solve(real_t p_step) {
- // Apply accumulated impulse.
- if (dynamic_A) {
- A->apply_impulse(-P, rA);
- }
- if (B && dynamic_B) {
- B->apply_impulse(P, rB);
- }
-
- return true;
-}
-
-void PinJoint2DSW::solve(real_t p_step) {
- // compute relative velocity
- Vector2 vA = A->get_linear_velocity() - custom_cross(rA - A->get_center_of_mass(), A->get_angular_velocity());
-
- Vector2 rel_vel;
- if (B) {
- rel_vel = B->get_linear_velocity() - custom_cross(rB - B->get_center_of_mass(), B->get_angular_velocity()) - vA;
- } else {
- rel_vel = -vA;
- }
-
- Vector2 impulse = M.basis_xform(bias - rel_vel - Vector2(softness, softness) * P);
-
- if (dynamic_A) {
- A->apply_impulse(-impulse, rA);
- }
- if (B && dynamic_B) {
- B->apply_impulse(impulse, rB);
- }
-
- P += impulse;
-}
-
-void PinJoint2DSW::set_param(PhysicsServer2D::PinJointParam p_param, real_t p_value) {
- if (p_param == PhysicsServer2D::PIN_JOINT_SOFTNESS) {
- softness = p_value;
- }
-}
-
-real_t PinJoint2DSW::get_param(PhysicsServer2D::PinJointParam p_param) const {
- if (p_param == PhysicsServer2D::PIN_JOINT_SOFTNESS) {
- return softness;
- }
- ERR_FAIL_V(0);
-}
-
-PinJoint2DSW::PinJoint2DSW(const Vector2 &p_pos, Body2DSW *p_body_a, Body2DSW *p_body_b) :
- Joint2DSW(_arr, p_body_b ? 2 : 1) {
- A = p_body_a;
- B = p_body_b;
- anchor_A = p_body_a->get_inv_transform().xform(p_pos);
- anchor_B = p_body_b ? p_body_b->get_inv_transform().xform(p_pos) : p_pos;
-
- p_body_a->add_constraint(this, 0);
- if (p_body_b) {
- p_body_b->add_constraint(this, 1);
- }
-}
-
-//////////////////////////////////////////////
-//////////////////////////////////////////////
-//////////////////////////////////////////////
-
-static inline void
-k_tensor(Body2DSW *a, Body2DSW *b, Vector2 r1, Vector2 r2, Vector2 *k1, Vector2 *k2) {
- // calculate mass matrix
- // If I wasn't lazy and wrote a proper matrix class, this wouldn't be so gross...
- real_t k11, k12, k21, k22;
- real_t m_sum = a->get_inv_mass() + b->get_inv_mass();
-
- // start with I*m_sum
- k11 = m_sum;
- k12 = 0.0f;
- k21 = 0.0f;
- k22 = m_sum;
-
- r1 -= a->get_center_of_mass();
- r2 -= b->get_center_of_mass();
-
- // add the influence from r1
- real_t a_i_inv = a->get_inv_inertia();
- real_t r1xsq = r1.x * r1.x * a_i_inv;
- real_t r1ysq = r1.y * r1.y * a_i_inv;
- real_t r1nxy = -r1.x * r1.y * a_i_inv;
- k11 += r1ysq;
- k12 += r1nxy;
- k21 += r1nxy;
- k22 += r1xsq;
-
- // add the influnce from r2
- real_t b_i_inv = b->get_inv_inertia();
- real_t r2xsq = r2.x * r2.x * b_i_inv;
- real_t r2ysq = r2.y * r2.y * b_i_inv;
- real_t r2nxy = -r2.x * r2.y * b_i_inv;
- k11 += r2ysq;
- k12 += r2nxy;
- k21 += r2nxy;
- k22 += r2xsq;
-
- // invert
- real_t determinant = k11 * k22 - k12 * k21;
- ERR_FAIL_COND(determinant == 0.0);
-
- real_t det_inv = 1.0f / determinant;
- *k1 = Vector2(k22 * det_inv, -k12 * det_inv);
- *k2 = Vector2(-k21 * det_inv, k11 * det_inv);
-}
-
-static _FORCE_INLINE_ Vector2
-mult_k(const Vector2 &vr, const Vector2 &k1, const Vector2 &k2) {
- return Vector2(vr.dot(k1), vr.dot(k2));
-}
-
-bool GrooveJoint2DSW::setup(real_t p_step) {
- dynamic_A = (A->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
- dynamic_B = (B->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
-
- if (!dynamic_A && !dynamic_B) {
- return false;
- }
-
- Space2DSW *space = A->get_space();
- ERR_FAIL_COND_V(!space, false);
-
- // calculate endpoints in worldspace
- Vector2 ta = A->get_transform().xform(A_groove_1);
- Vector2 tb = A->get_transform().xform(A_groove_2);
-
- // calculate axis
- Vector2 n = -(tb - ta).orthogonal().normalized();
- real_t d = ta.dot(n);
-
- xf_normal = n;
- rB = B->get_transform().basis_xform(B_anchor);
-
- // calculate tangential distance along the axis of rB
- real_t td = (B->get_transform().get_origin() + rB).cross(n);
- // calculate clamping factor and rB
- if (td <= ta.cross(n)) {
- clamp = 1.0f;
- rA = ta - A->get_transform().get_origin();
- } else if (td >= tb.cross(n)) {
- clamp = -1.0f;
- rA = tb - A->get_transform().get_origin();
- } else {
- clamp = 0.0f;
- //joint->r1 = cpvsub(cpvadd(cpvmult(cpvperp(n), -td), cpvmult(n, d)), a->p);
- rA = ((-n.orthogonal() * -td) + n * d) - A->get_transform().get_origin();
- }
-
- // Calculate mass tensor
- k_tensor(A, B, rA, rB, &k1, &k2);
-
- // compute max impulse
- jn_max = get_max_force() * p_step;
-
- // calculate bias velocity
- //cpVect delta = cpvsub(cpvadd(b->p, joint->r2), cpvadd(a->p, joint->r1));
- //joint->bias = cpvclamp(cpvmult(delta, -joint->constraint.biasCoef*dt_inv), joint->constraint.maxBias);
-
- Vector2 delta = (B->get_transform().get_origin() + rB) - (A->get_transform().get_origin() + rA);
-
- real_t _b = get_bias();
- gbias = (delta * -(_b == 0 ? space->get_constraint_bias() : _b) * (1.0 / p_step)).limit_length(get_max_bias());
-
- correct = true;
- return true;
-}
-
-bool GrooveJoint2DSW::pre_solve(real_t p_step) {
- // Apply accumulated impulse.
- if (dynamic_A) {
- A->apply_impulse(-jn_acc, rA);
- }
- if (dynamic_B) {
- B->apply_impulse(jn_acc, rB);
- }
-
- return true;
-}
-
-void GrooveJoint2DSW::solve(real_t p_step) {
- // compute impulse
- Vector2 vr = relative_velocity(A, B, rA, rB);
-
- Vector2 j = mult_k(gbias - vr, k1, k2);
- Vector2 jOld = jn_acc;
- j += jOld;
-
- jn_acc = (((clamp * j.cross(xf_normal)) > 0) ? j : j.project(xf_normal)).limit_length(jn_max);
-
- j = jn_acc - jOld;
-
- if (dynamic_A) {
- A->apply_impulse(-j, rA);
- }
- if (dynamic_B) {
- B->apply_impulse(j, rB);
- }
-}
-
-GrooveJoint2DSW::GrooveJoint2DSW(const Vector2 &p_a_groove1, const Vector2 &p_a_groove2, const Vector2 &p_b_anchor, Body2DSW *p_body_a, Body2DSW *p_body_b) :
- Joint2DSW(_arr, 2) {
- A = p_body_a;
- B = p_body_b;
-
- A_groove_1 = A->get_inv_transform().xform(p_a_groove1);
- A_groove_2 = A->get_inv_transform().xform(p_a_groove2);
- B_anchor = B->get_inv_transform().xform(p_b_anchor);
- A_groove_normal = -(A_groove_2 - A_groove_1).normalized().orthogonal();
-
- A->add_constraint(this, 0);
- B->add_constraint(this, 1);
-}
-
-//////////////////////////////////////////////
-//////////////////////////////////////////////
-//////////////////////////////////////////////
-
-bool DampedSpringJoint2DSW::setup(real_t p_step) {
- dynamic_A = (A->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
- dynamic_B = (B->get_mode() > PhysicsServer2D::BODY_MODE_KINEMATIC);
-
- if (!dynamic_A && !dynamic_B) {
- return false;
- }
-
- rA = A->get_transform().basis_xform(anchor_A);
- rB = B->get_transform().basis_xform(anchor_B);
-
- Vector2 delta = (B->get_transform().get_origin() + rB) - (A->get_transform().get_origin() + rA);
- real_t dist = delta.length();
-
- if (dist) {
- n = delta / dist;
- } else {
- n = Vector2();
- }
-
- real_t k = k_scalar(A, B, rA, rB, n);
- n_mass = 1.0f / k;
-
- target_vrn = 0.0f;
- v_coef = 1.0f - Math::exp(-damping * (p_step)*k);
-
- // Calculate spring force.
- real_t f_spring = (rest_length - dist) * stiffness;
- j = n * f_spring * (p_step);
-
- return true;
-}
-
-bool DampedSpringJoint2DSW::pre_solve(real_t p_step) {
- // Apply spring force.
- if (dynamic_A) {
- A->apply_impulse(-j, rA);
- }
- if (dynamic_B) {
- B->apply_impulse(j, rB);
- }
-
- return true;
-}
-
-void DampedSpringJoint2DSW::solve(real_t p_step) {
- // compute relative velocity
- real_t vrn = normal_relative_velocity(A, B, rA, rB, n) - target_vrn;
-
- // compute velocity loss from drag
- // not 100% certain this is derived correctly, though it makes sense
- real_t v_damp = -vrn * v_coef;
- target_vrn = vrn + v_damp;
- Vector2 j = n * v_damp * n_mass;
-
- if (dynamic_A) {
- A->apply_impulse(-j, rA);
- }
- if (dynamic_B) {
- B->apply_impulse(j, rB);
- }
-}
-
-void DampedSpringJoint2DSW::set_param(PhysicsServer2D::DampedSpringParam p_param, real_t p_value) {
- switch (p_param) {
- case PhysicsServer2D::DAMPED_SPRING_REST_LENGTH: {
- rest_length = p_value;
- } break;
- case PhysicsServer2D::DAMPED_SPRING_DAMPING: {
- damping = p_value;
- } break;
- case PhysicsServer2D::DAMPED_SPRING_STIFFNESS: {
- stiffness = p_value;
- } break;
- }
-}
-
-real_t DampedSpringJoint2DSW::get_param(PhysicsServer2D::DampedSpringParam p_param) const {
- switch (p_param) {
- case PhysicsServer2D::DAMPED_SPRING_REST_LENGTH: {
- return rest_length;
- } break;
- case PhysicsServer2D::DAMPED_SPRING_DAMPING: {
- return damping;
- } break;
- case PhysicsServer2D::DAMPED_SPRING_STIFFNESS: {
- return stiffness;
- } break;
- }
-
- ERR_FAIL_V(0);
-}
-
-DampedSpringJoint2DSW::DampedSpringJoint2DSW(const Vector2 &p_anchor_a, const Vector2 &p_anchor_b, Body2DSW *p_body_a, Body2DSW *p_body_b) :
- Joint2DSW(_arr, 2) {
- A = p_body_a;
- B = p_body_b;
- anchor_A = A->get_inv_transform().xform(p_anchor_a);
- anchor_B = B->get_inv_transform().xform(p_anchor_b);
-
- rest_length = p_anchor_a.distance_to(p_anchor_b);
-
- A->add_constraint(this, 0);
- B->add_constraint(this, 1);
-}