Differences between revisions 2 and 3

Deletions are marked like this. Additions are marked like this.
Line 20: Line 20:
 l21 = Trsv_ln_unb_var2(l21, L22, l21);  l21 = Trsv_ln_unb_var1(l21, L22);

function [ L_out ] = Trinv_unb_var4( L )

  [ LTL, LTR, ...
    LBL, LBR ] = FLA_Part_2x2( L, ...
                               0, 0, 'FLA_TL' );

  while ( size( LTL, 1 ) < size( L, 1 ) )

    [ L00,  l01,      L02,  ...
      l10t, lambda11, l12t, ...
      L20,  l21,      L22 ] = FLA_Repart_2x2_to_3x3( LTL, LTR, ...
                                                     LBL, LBR, ...
                                                     1, 1, 'FLA_BR' );

    %------------------------------------------------------------%

        lambda11 = 1/lambda11;
        l21 = Trsv_ln_unb_var1(l21, L22);
        l21 = -l21*lambda11; 

    %------------------------------------------------------------%

    [ LTL, LTR, ...
      LBL, LBR ] = FLA_Cont_with_3x3_to_2x2( L00,  l01,      L02,  ...
                                             l10t, lambda11, l12t, ...
                                             L20,  l21,      L22, ...
                                             'FLA_TL' );

  end

  L_out = [ LTL, LTR
            LBL, LBR ];

return

LinearAlgebraWiki: Trinv/FLAME@lab/UnbVar4 (last edited 2008-08-06 20:29:53 by MarthaGanser)