Differences between revisions 2 and 3
| Deletions are marked like this. | Additions are marked like this. |
| Line 3: | Line 3: |
| function [ b_out ] = Trsv_ln_unb_var1( b, L, x ) | function [ b_out ] = Trsv_ln_unb_var1( b, L ) |
| Line 12: | Line 12: |
|
[ xT, ... xB ] = FLA_Part_2x1( x, ... 0, 'FLA_TOP' ); |
|
| Line 30: | Line 26: |
|
[ x0, ... chi1, ... x2 ] = FLA_Repart_2x1_to_3x1( xT, ... xB, ... 1, 'FLA_BOTTOM' ); |
|
| Line 56: | Line 46: |
|
[ xT, ... xB ] = FLA_Cont_with_3x1_to_2x1( x0, ... chi1, ... x2, ... 'FLA_TOP' ); |
function [ b_out ] = Trsv_ln_unb_var1( b, L )
[ bT, ...
bB ] = FLA_Part_2x1( b, ...
0, 'FLA_TOP' );
[ LTL, LTR, ...
LBL, LBR ] = FLA_Part_2x2( L, ...
0, 0, 'FLA_TL' );
while ( size( bT, 1 ) < size( b, 1 ) )
[ b0, ...
beta1, ...
b2 ] = FLA_Repart_2x1_to_3x1( bT, ...
bB, ...
1, 'FLA_BOTTOM' );
[ L00, l01, L02, ...
l10t, lambda11, l12t, ...
L20, l21, L22 ] = FLA_Repart_2x2_to_3x3( LTL, LTR, ...
LBL, LBR, ...
1, 1, 'FLA_BR' );
%------------------------------------------------------------%
beta1 = beta1 - l10t * b0;
beta1 = beta1 / lambda11;
%------------------------------------------------------------%
[ bT, ...
bB ] = FLA_Cont_with_3x1_to_2x1( b0, ...
beta1, ...
b2, ...
'FLA_TOP' );
[ LTL, LTR, ...
LBL, LBR ] = FLA_Cont_with_3x3_to_2x2( L00, l01, L02, ...
l10t, lambda11, l12t, ...
L20, l21, L22, ...
'FLA_TL' );
end
b_out = [ bT
bB ];
return