function [ y_out ] = Gemv_unb_var1( y, A, x ) [ yT, ... yB ] = FLA_Part_2x1( y, ... 0, 'FLA_TOP' ); [ AT, ... AB ] = FLA_Part_2x1( A, ... 0, 'FLA_TOP' ); while ( size( yT, 1 ) < size( y, 1 ) ) [ y0, ... psi1, ... y2 ] = FLA_Repart_2x1_to_3x1( yT, ... yB, ... 1, 'FLA_BOTTOM' ); [ A0, ... a1t, ... A2 ] = FLA_Repart_2x1_to_3x1( AT, ... AB, ... 1, 'FLA_BOTTOM' ); %------------------------------------------------------------% psi1 = a1t * x + psi1; %------------------------------------------------------------% [ yT, ... yB ] = FLA_Cont_with_3x1_to_2x1( y0, ... psi1, ... y2, ... 'FLA_TOP' ); [ AT, ... AB ] = FLA_Cont_with_3x1_to_2x1( A0, ... a1t, ... A2, ... 'FLA_TOP' ); end y_out = [ yT yB ]; return
