#include "FLAME.h"
int Chol_unb_var1( FLA_Obj A )
{
 FLA_Obj ATL, ATR, A00, a01, A02, 
 ABL, ABR, a10t, alpha11, a12t,
 A20, a21, A22;
 FLA_Part_2x2( A, &ATL, &ATR,
 &ABL, &ABR, 0, 0, FLA_TL );
 while ( FLA_Obj_length( ATL ) < FLA_Obj_length( A ) ){
 FLA_Repart_2x2_to_3x3( ATL, /**/ ATR, &A00, /**/ &a01, &A02,
 /* ************* */ /* ************************** */
 &a10t, /**/ &alpha11, &a12t,
 ABL, /**/ ABR, &A20, /**/ &a21, &A22,
 1, 1, FLA_BR );
 /*------------------------------------------------------------*/
 // a10t = a10t * inv( tril( A00 ) )'; % TRSV
 // a10 = inv( tril( A00 ) ) * a10; % TRSV
 FLA_Trsv( FLA_LOWER_TRIANGULAR, FLA_NO_TRANSPOSE, FLA_NONUNIT_DIAG, 
	 A00, a10t );
 // alpha11 = alpha11 - a10t * a10t';
 FLA_Dotcs( FLA_NO_CONJUGATE, FLA_MINUS_ONE, a10t, a10t, FLA_ONE, alpha11 );
 // alpha11 = sqrt( alpha11 );
 FLA_Sqrt( alpha11 );
 /*------------------------------------------------------------*/
 FLA_Cont_with_3x3_to_2x2( &ATL, /**/ &ATR, A00, a01, /**/ A02,
 a10t, alpha11, /**/ a12t,
 /* ************** */ /* ************************ */
 &ABL, /**/ &ABR, A20, a21, /**/ A22,
 FLA_TL );
 }
 return FLA_SUCCESS;
}
</div><div class="naked_ctrl">
<form action="/index.cgi/larger-text" method="get" name="gate">
<p><a href="http://altstyle.alfasado.net">AltStyle</a> によって変換されたページ <a href="https://www.cs.utexas.edu/~flame/Codes/Chol/FLAMEC/Chol_unb_var1.c">(-&gt;オリジナル)</a>
/ <label>アドレス: <input type="text" name="naked_post_url" value="https://www.cs.utexas.edu/~flame/Codes/Chol/FLAMEC/Chol_unb_var1.c" size="22" /></label> <label>モード: <select name="naked_post_mode">
<option value="default">デフォルト</option>
<option value="speech">音声ブラウザ</option>
<option value="ruby">ルビ付き</option>
<option value="contrast">配色反転</option>
<option value="larger-text" selected="selected">文字拡大</option>
<option value="mobile">モバイル</option>
</select>
<input type="submit" value="表示" />
</p>
</form>
</div>