Method circumcenterOf

src/main/java/field/linalg/BaseMath.java:244–267

(Vec3 A, Vec3 B, Vec3 C)

	}

	static public Vec3 circumcenterOf(Vec3 A, Vec3 B, Vec3 C) {
		Vec3 r = new Vec3();
		if (A.z == 0.0 && B.z == 0.0 && C.z == 0.0) {
			double u = ((A.x - B.x) * (A.x + B.x) + (A.y - B.y) * (A.y + B.y)) / 2.0;
			double v = ((B.x - C.x) * (B.x + C.x) + (B.y - C.y) * (B.y + C.y)) / 2.0;
			double den = (A.x - B.x) * (B.y - C.y) - (B.x - C.x) * (A.y - B.y);
			r.set(0, (float) ((u * (B.y - C.y) - v * (A.y - B.y)) / den));
			r.set(1, (float) ((v * (A.x - B.x) - u * (B.x - C.x)) / den));
			r.set(2, (float) 0.0);
		} else {
			Vec3 BmA = Vec3.sub(B, A, new Vec3());
			Vec3 CmA = Vec3.sub(C, A, new Vec3());

			double BC = BmA.dot(CmA);
			double B2 = BmA.lengthSquared(), C2 = CmA.lengthSquared();
			double den = 2.0 * (B2 * C2 - BC * BC);
			double s = C2 * (B2 - BC) / den;
			double t = B2 * (C2 - BC) / den;

			Vec3.fma(BmA, (float) s, A, r);
			Vec3.fma(CmA, (float) t, r, r);
		}
		return r;
	}

	static public Vec2 circumcenterOf(Vec2 A, Vec2 B, Vec2 C) {
		Vec2 r = new Vec2();