Files
Mindustry/core/src/mindustry/graphics/InverseKinematics.java
2023-03-25 11:11:06 -04:00

43 lines
1.9 KiB
Java

package mindustry.graphics;
import arc.math.*;
import arc.math.geom.*;
public class InverseKinematics{
private static final Vec2[] mat1 = {new Vec2(), new Vec2()}, mat2 = {new Vec2(), new Vec2()};
private static final Vec2 temp = new Vec2(), temp2 = new Vec2(), at1 = new Vec2();
public static boolean solve(float lengthA, float lengthB, Vec2 end, boolean side, Vec2 result){
at1.set(end).rotate(side ? 1 : -1).setLength(lengthA + lengthB).add(end.x / 2f, end.y / 2f);
return solve(lengthA, lengthB, end, at1, result);
}
/**
inputs:
@param lengthA first line segment length
@param lengthB second line segment length
@param end length of the endpoint you want to reach
@param attractor direction you want the result to be closer to (since there are usually 2 solutions)
output:
@param result a point in-between (0, 0) and (end) such that (0, 0).dst(result) == lengthA and result.dst(end) == lengthB - or in basic terms, the position of a joint between (0, 0) and end where the two lengths of segments are lengthA and lengthB
@return whether IK succeeded (this can fail if end is too far, for example)
*/
public static boolean solve(float lengthA, float lengthB, Vec2 end, Vec2 attractor, Vec2 result){
Vec2 axis = mat2[0].set(end).nor();
mat2[1].set(attractor).sub(temp2.set(axis).scl(attractor.dot(axis))).nor();
mat1[0].set(mat2[0].x, mat2[1].x);
mat1[1].set(mat2[0].y, mat2[1].y);
result.set(mat2[0].dot(end), mat2[1].dot(end));
float len = result.len();
float dist = Math.max(0, Math.min(lengthA, (len + (lengthA * lengthA - lengthB * lengthB) / len) / 2));
Vec2 src = temp.set(dist, Mathf.sqrt(lengthA * lengthA - dist * dist));
result.set(mat1[0].dot(src), mat1[1].dot(src));
return dist > 0 && dist < lengthA;
}
}