Differences between revisions 3 and 4
| Deletions are marked like this. | Additions are marked like this. |
| Line 1: | Line 1: |
| {{{#!frame align=float:left,width=10%,background=yellow | {{{#!frame align=float:left,width=10%,position=relative,background=yellow |
| Line 6: | Line 6: |
|
~-[[Include(Gemv/FLAME@lab/UnbVar1)]]-~ |
[[Include(Gemv/FLAME@lab/UnbVar1)]] |
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