void VU_copu_f(
vsip_cmview_f *A,
vsip_cvview_f *x,
vsip_cvview_f *y)
{
vsip_length m = vsip_cmgetcollength_f(A);
vsip_length n = vsip_cmgetrowlength_f(A);
vsip_length i,j;
vsip_cscalar_f temp;
for(i=0; i<m; i++)
for(j=0; j<n; j++){
temp = vsip_cmul_f(vsip_cvget_f(x,i),vsip_cvget_f(y,j));
vsip_cmput_f(A,i,j,vsip_cadd_f(vsip_cmget_f(A,i,j),temp));
}
return;
}
vsip_cscalar_f (vsip_csub_f)(
vsip_cscalar_f x, vsip_cscalar_f y) {/* x - y*/
return vsip_cadd_f(x, vsip_cneg_f(y)); }