/* s e t u p T Q f a c t o r i s a t i o n  */
returnValue QProblem_setupTQfactorisation(QProblem *_THIS)
/* s e t u p T Q f a c t o r i s a t i o n  */
returnValue QProblem_setupTQfactorisation(QProblem *_THIS)
{
	int i, ii;
	int nFR = QProblem_getNFR( _THIS );
	int* FR_idx;
	Indexlist_getNumberArray( Bounds_getFree( &(_THIS->bounds) ),&FR_idx );
	for( i=0; i<NVMAX*NVMAX; ++i )
		_THIS->Q[i] = 0.0;
	for( i=0; i<nFR; ++i )
	{
		ii = FR_idx[i];
		QQ(ii,i) = 1.0;
	}
	for( i=0; i<NVCMIN*NVCMIN; ++i )
		_THIS->T[i] = 0.0;
	return SUCCESSFUL_RETURN;
}