#include<bits/stdc++.h>
using ll = long long;
const int N = 1005;
using std::cin, std::cout;
using db = long double;
struct vec2 {
    db x, y;
    db norm() const {
        return x * x + y * y;
    }
};
vec2 operator + (vec2 x, vec2 y) { return {x.x + y.x, x.y + y.y}; }
vec2 operator - (vec2 x, vec2 y) { return {x.x - y.x, x.y - y.y}; }
vec2 operator * (vec2 x, db y) { return {x.x * y, x.y * y}; }
vec2 r90(vec2 x) { return {-x.y, x.x}; }
db operator * (vec2 x, vec2 y) { return x.x * y.y - x.y * y.x; }
db operator % (vec2 x, vec2 y) { return x.x * y.x + x.y * y.y; }
vec2 get() {
    int x, y; cin >> x >> y;
    return {(db)x,db(y)};
}
const db eps = 1e-12;
struct func {
    vec2 o;
    db c;
    db prod;
    bool init(vec2 a, vec2 b) {
        o = r90(a - b);
        if(o.x < 0) a.x *= -1, b.x *= -1, std::swap(a, b);
        if(o.y < 0) a.y *= -1, b.y *= -1, std::swap(a, b);
        o = r90(a - b);
        if(fabs(o.x) < eps) return 0;
        if(fabs(o.y) < eps) return 0;
        o = o * (1 / sqrt(o.norm()));
        c = o % a;
        prod = o.norm() / o.x / o.y;
        return 1;
    }
    db calc(db x, db y) const {
        db v = std::max<db>(0, x * o.x + y * o. y - c);
        return v * v * prod;
    }
} pb, fsy;
db f(db x, db y) {
    return pb.calc(x, y) + fsy.calc(x, y);
}
db S;
db g(db x, db y) {
    db r = f(x, y);
    return (4 * x * y - r) / (S + r);
}
db ly, ry;
db findmax(db a, db c, auto f) {
    auto g = [&](db l, db r) {
        return l + (r - l) * (std::numbers::phi_v<db> - 1);
    };
    db b = g(a, c), bv = f(b);
    for(int i = 0;i < 45;++i) {
        db x = g(a, b), xv = f(x);
        if(xv > bv) {
            c = b, b = x, bv = xv;
        } else {
            a = c, c = x;
        }
    }
    return bv;
}
int main() {
    std::ios::sync_with_stdio(false), cin.tie(0);
    int T; cin >> T;
    for(int i = 0;i < T;++i) {
        vec2 a[4];
        vec2 s = {};
        for(auto & x : a) x = get(), s = s + x;
        if((a[1] - a[0]) * (a[2] - a[0]) < 0) {
            std::reverse(a, a + 4);
        }
        s = s * 0.25;
        for(auto & x : a) x = x - s, x = x * (1 / sqrt(x.norm()));
        S = (a[1] - a[0]) * (a[2] - a[0]);
        if(!pb.init(a[0], a[1])) {
            puts("1");
            continue;
        }
        if(!fsy.init(a[1], a[2])) {
            puts("1");
            continue;
        }
        db l = 0, r = a[0].x;
        ry = a[0].y;
        for(auto & x : a) {
            r = std::max(r, x.x);
            ry = std::max(ry, x.y);
        }
        printf("%.12Lf\n", findmax(l, r, [&](db x) { return findmax(0, ry, [&](db y) { return g(x, y); });}));
    }

}