106 using Teuchos::Workspace;
107 typedef Teuchos::ScalarTraits<Scalar> ST;
108 typedef typename Teuchos::ArrayRCP<Scalar>::iterator
iter_t;
109 typedef typename Teuchos::ArrayRCP<const Scalar>::iterator
const_iter_t;
110 Teuchos::WorkspaceStore*
wss = Teuchos::get_default_workspace_store().get();
147 if(
l_beta==ST::zero() ) {
150 for(
int j = 0;
j < subDim; ++
j )
158 else if(
l_beta==ST::one() ) {
163 for(
int j = 0;
j < subDim; ++
j )
174 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
192 if(
l_beta==ST::zero() ) {
193 if(
alpha0 == ST::one() ) {
194 if(
alpha1 == ST::one() ) {
197 for(
int j = 0;
j < subDim; ++
j )
202 (*
z0_val) = (*v0_val) + (*v1_val);
208 for(
int j = 0;
j < subDim; ++
j )
218 if(
alpha1 == ST::one() ) {
221 for(
int j = 0;
j < subDim; ++
j )
232 for(
int j = 0;
j < subDim; ++
j )
242 else if(
l_beta==ST::one() ) {
243 if(
alpha0 == ST::one() ) {
244 if(
alpha1 == ST::one() ) {
247 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
252 (*
z0_val) += (*v0_val) + (*v1_val);
258 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
268 if(
alpha1 == ST::one() ) {
271 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
282 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
293 if(
alpha0 == ST::one() ) {
294 if(
alpha1 == ST::one() ) {
297 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
308 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
318 if(
alpha1 == ST::one() ) {
321 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
332 for(
int j = 0;
j < subDim; ++
j, ++
z0_val )
348 if( beta_ == ST::zero() ) {
352 else if( beta_ != ST::one() ) {