74 #define A(_m, _n) (float *)plasma_geteltaddr(A, ((_m)-1), ((_n)-1), eltsize)
75 #define V(_m) &(V[(_m)-1])
76 #define TAU(_m) &(TAU[(_m)-1])
87 int len1, len2, t1ed, t2st;
110 for (i = ed; i >= st+1 ; i--){
115 t2st =
max(t1ed+1,J1);
118 if(len1>0)
CORE_slarfx2(
PlasmaLeft, len1 , *
V(i), (*
TAU(i)),
A(i-1, J1 ),
ELTLDD(vA, (i-1)),
A(i, J1 ),
ELTLDD(vA, i) );
119 if(len2>0)
CORE_slarfx2(
PlasmaLeft, len2 , *
V(i), (*
TAU(i)),
A(i-1, t2st),
ELTLDD(vA, (i-1)),
A(i, t2st),
ELTLDD(vA, i) );
123 for (i = ed; i >= st+1 ; i--){
127 t2st =
max(t1ed+1,J1);
130 if(len1>0)
CORE_slarfx2(
PlasmaRight, len1, (*
V(i)), (*
TAU(i)),
A(J1,i-1),
ELTLDD(vA, J1) ,
A(J1 , i),
ELTLDD(vA, J1) );
131 if(len2>0)
CORE_slarfx2(
PlasmaRight, len2, (*
V(i)), (*
TAU(i)),
A(t2st,i-1),
ELTLDD(vA, t2st),
A(t2st, i),
ELTLDD(vA, t2st) );
137 for (i = ed; i >= st+1 ; i--){
142 t2st =
max(t1ed+1,J1);
145 if(len1>0)
CORE_slarfx2(
PlasmaRight, len1, (*
V(i)), (*
TAU(i)),
A(J1,i-1),
ELTLDD(vA, J1) ,
A(J1 , i),
ELTLDD(vA, J1) );
146 if(len2>0)
CORE_slarfx2(
PlasmaRight, len2, (*
V(i)), (*
TAU(i)),
A(t2st,i-1),
ELTLDD(vA, t2st),
A(t2st, i),
ELTLDD(vA, t2st) );
150 for (i = ed; i >= st+1 ; i--){
154 t2st =
max(t1ed+1,J1);
157 if(len1>0)
CORE_slarfx2(
PlasmaLeft, len1 , *
V(i), (*
TAU(i)),
A(i-1, J1 ),
ELTLDD(vA, (i-1)),
A(i, J1 ),
ELTLDD(vA, i) );
158 if(len2>0)
CORE_slarfx2(
PlasmaLeft, len2 , *
V(i), (*
TAU(i)),
A(i-1, t2st),
ELTLDD(vA, (i-1)),
A(i, t2st),
ELTLDD(vA, i) );