Lines Matching refs:hv
241 double *hv;
276 hv = update_u ? uT : hv0;
281 scale += fabs( hv[j] = t );
290 double t = (hv[j] *= f);
295 f = hv[0];
298 hv[0] = f - g;
303 /* calc temp[0:n-i] = a[i:m,i:n]'*hv[0:m-i] */
304 icvMatrAXPY_64f( m1, n1 - 1, a + 1, lda, hv, temp + 1, 0 );
307 /* modify a: a[i:m,i:n] = a[i:m,i:n] + hv[0:m-i]*temp[0:n-i]' */
308 icvMatrAXPY_64f( m1, n1 - 1, temp + 1, 0, hv, a + 1, lda );
313 /* store -2/(hv'*hv) */
319 hv[-1] = h;
333 hv = update_v ? vT : hv0;
338 scale += fabs( hv[j] = t );
347 double t = (hv[j] *= f);
352 f = hv[0];
355 hv[0] = f - g;
357 hv[-1] = 0.;
359 /* update a[i:m:i+1:n] = a[i:m,i+1:n] + (a[i:m,i+1:n]*hv[0:m-i])*... */
360 icvMatrAXPY3_64f( m1, n1, hv, lda, a, h );
366 /* store -2/(hv'*hv) */
372 hv[-1] = h;
402 hv = u0 + (lduT + 1) * i;
403 h = i == 0 ? ku0 : hv[-1];
409 uT = hv;
410 icvMatrAXPY3_64f( lh, l-1, hv+1, lduT, uT+1, h );
412 s = hv[0] * h;
413 for( k = 0; k < l; k++ ) hv[k] *= s;
414 hv[0] += 1;
419 hv[j] = 0;
421 hv[j * lduT] = 0;
422 hv[0] = 1;
445 hv = v0 + (ldvT + 1) * i;
446 h = i == 0 ? kv0 : hv[-1];
452 vT = hv;
453 icvMatrAXPY3_64f( lh, l-1, hv+1, ldvT, vT+1, h );
455 s = hv[0] * h;
456 for( k = 0; k < l; k++ ) hv[k] *= s;
457 hv[0] += 1;
462 hv[j] = 0;
464 hv[j * ldvT] = 0;
465 hv[0] = 1;
636 float *hv;
673 hv = update_u ? uT : hv0;
678 scale += fabs( hv[j] = (float)t );
687 double t = (hv[j] = (float)(hv[j]*f));
692 f = hv[0];
695 hv[0] = (float)(f - g);
700 /* calc temp[0:n-i] = a[i:m,i:n]'*hv[0:m-i] */
701 icvMatrAXPY_32f( m1, n1 - 1, a + 1, lda, hv, temp + 1, 0 );
705 /* modify a: a[i:m,i:n] = a[i:m,i:n] + hv[0:m-i]*temp[0:n-i]' */
706 icvMatrAXPY_32f( m1, n1 - 1, temp + 1, 0, hv, a + 1, lda );
711 /* store -2/(hv'*hv) */
717 hv[-1] = (float)h;
730 hv = update_v ? vT : hv0;
735 scale += fabs( hv[j] = (float)t );
744 double t = (hv[j] = (float)(hv[j]*f));
749 f = hv[0];
752 hv[0] = (float)(f - g);
754 hv[-1] = 0.f;
756 /* update a[i:m:i+1:n] = a[i:m,i+1:n] + (a[i:m,i+1:n]*hv[0:m-i])*... */
757 icvMatrAXPY3_32f( m1, n1, hv, lda, a, h );
763 /* store -2/(hv'*hv) */
769 hv[-1] = (float)h;
799 hv = u0 + (lduT + 1) * i;
800 h = i == 0 ? ku0 : hv[-1];
806 uT = hv;
807 icvMatrAXPY3_32f( lh, l-1, hv+1, lduT, uT+1, h );
809 s = hv[0] * h;
810 for( k = 0; k < l; k++ ) hv[k] = (float)(hv[k]*s);
811 hv[0] += 1;
816 hv[j] = 0;
818 hv[j * lduT] = 0;
819 hv[0] = 1;
842 hv = v0 + (ldvT + 1) * i;
843 h = i == 0 ? kv0 : hv[-1];
849 vT = hv;
850 icvMatrAXPY3_32f( lh, l-1, hv+1, ldvT, vT+1, h );
852 s = hv[0] * h;
853 for( k = 0; k < l; k++ ) hv[k] = (float)(hv[k]*s);
854 hv[0] += 1;
859 hv[j] = 0;
861 hv[j * ldvT] = 0;
862 hv[0] = 1;
