#include<bits/stdc++.h>
using namespace std;

constexpr double eps=1e-9;

int r[3][3],s[3];
double low,mid[3],ans[3];

template<typename T>T sqr(T x)
{
	return x*x;
}

bool sing(int a)
{
	int b=(a+1)%3,c=(a+2)%3;
	double bm=r[b][c]/2.0,bm2=sqr(bm),am2=(sqr(r[a][b])+sqr(r[a][c]))/2.0-bm2,am=sqrt(am2);
	mid[a]=am*2/3;
	for(int i:{b,c})
	{
		int ab=r[a][i],ab2=sqr(ab);
		double al=(am2+ab2-bm2)/am;
		if(sqr(s[i])<al*s[a]+ab2+sqr(s[a])-eps) return false;
		ans[i]=-s[a]-al/2;
	}
	ans[a]=-s[a];
	return true;
}

double hei(double h,double x,double y)
{
	return -(sqr(h)+sqr(y)-sqr(x))/2/h;
}

bool twin(int a,int b)
{
	auto dis=[](double ad,double bd,double ab,double ac,double bc)
	{
		double cosdab=(sqr(ad)+sqr(ab)-sqr(bd))/2/ab/ad;
		double coscab=(sqr(ac)+sqr(ab)-sqr(bc))/2/ab/ac;
		double cosdac=cosdab*coscab-sqrt((1-sqr(cosdab))*(1-sqr(coscab)));
		return sqrt(sqr(ad)+sqr(ac)-2*ad*ac*cosdac);
	};
	if(s[a]+s[b]<r[a][b]||abs(s[a]-s[b])>r[a][b]) return false;
	int c=3-a-b;
	double cd=dis(s[a],s[b],r[a][b],r[a][c],r[b][c]);
	if(s[c]<cd-eps) return false;
	double dg=dis(s[a],s[b],r[a][b],mid[a],mid[b]);
	if(dg>low)
	{
		low=dg;
		ans[a]=hei(dg,mid[a],s[a]);
		ans[b]=hei(dg,mid[b],s[b]);
		ans[c]=hei(dg,mid[c],cd);
	}
	return true;
}

bool pyr()
{
	double dg=3*(sqr(s[0])+sqr(s[1])+sqr(s[2]))-sqr(r[0][1])-sqr(r[1][2])-sqr(r[0][2]);
	if(dg<=0) return false;
	dg=sqrt(dg)/3;
	for(int i=0;i<3;i++) ans[i]=hei(dg,mid[i],s[i]);
	return true;
}

int main()
{
	ios::sync_with_stdio(false),cin.tie(0);
	int T;
	cin>>T;
	cout<<fixed<<setprecision(5);
	while(T--)
	{
		cin>>s[0]>>s[1]>>s[2]>>r[1][2]>>r[0][2]>>r[0][1];
		r[1][0]=r[0][1];
		r[2][1]=r[1][2];
		r[2][0]=r[0][2];
		low=0;
		sing(0)||sing(1)||sing(2)||(twin(0,1)|twin(1,2)|twin(0,2))||pyr();
		cout<<ans[0]<<' '<<ans[1]<<' '<<ans[2]<<'\n';
	}
	return 0;
}