Lines Matching +full:- +full:p
1 void pmat(float32_t *p,int nbrows,int nbcols) in pmat() argument
7 printf("%f ",(double)p[c+r*nbcols]); in pmat()
14 void pvec(float32_t *p,int nb) in pvec() argument
18 printf("%f ",(double)p[c]); in pvec()
23 void pvec(Q7 *p,int nb) in pvec() argument
27 printf("%f ",(double)(1.0f*p[c].v/128.0f)); in pvec()
52 float32_t *p; in _arm_mat_qr_f32() local
54 if (pSrc->numRows < pSrc->numCols) in _arm_mat_qr_f32()
59 memcpy(pOutR->pData,pSrc->pData,pSrc->numCols * pSrc->numRows*sizeof(float32_t)); in _arm_mat_qr_f32()
60 pOutR->numCols = pSrc->numCols; in _arm_mat_qr_f32()
61 pOutR->numRows = pSrc->numRows; in _arm_mat_qr_f32()
63 p = pOutR->pData; in _arm_mat_qr_f32()
66 for(col=0 ; col < pSrc->numCols; col++) in _arm_mat_qr_f32()
76 beta = arm_householder_f32(pTmpA,threshold,pSrc->numRows - col,pTmpA); in _arm_mat_qr_f32()
79 //pvec(pTmpA,pSrc->numRows-col); in _arm_mat_qr_f32()
80 //pmat(p,pSrc->numRows-col,pSrc->numCols-col); in _arm_mat_qr_f32()
84 /* v.T A(col:,col:) -> tmpb */ in _arm_mat_qr_f32()
86 pa = p; in _arm_mat_qr_f32()
89 blkCnt = (pSrc->numCols-col) >> 2; in _arm_mat_qr_f32()
98 blkCnt--; in _arm_mat_qr_f32()
100 blkCnt = (pSrc->numCols-col) & 3; in _arm_mat_qr_f32()
118 pa1 = pa0 + pSrc->numCols; in _arm_mat_qr_f32()
119 pa2 = pa1 + pSrc->numCols; in _arm_mat_qr_f32()
120 pa3 = pa2 + pSrc->numCols; in _arm_mat_qr_f32()
123 blkCnt = (pSrc->numRows-col - 1) >> 2; in _arm_mat_qr_f32()
129 blkCnt2 = (pSrc->numCols-col) >> 2; in _arm_mat_qr_f32()
153 blkCnt2--; in _arm_mat_qr_f32()
155 blkCnt2 = (pSrc->numCols-col) & 3; in _arm_mat_qr_f32()
182 pa0 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
183 pa1 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
184 pa2 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
185 pa3 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
189 blkCnt--; in _arm_mat_qr_f32()
193 for(;k<pSrc->numRows-col; k++) in _arm_mat_qr_f32()
196 blkCnt2 = (pSrc->numCols-col) >> 2; in _arm_mat_qr_f32()
206 blkCnt2--; in _arm_mat_qr_f32()
208 blkCnt2 = (pSrc->numCols-col) & 3; in _arm_mat_qr_f32()
225 //pvec(pTmpB,pSrc->numCols-col); in _arm_mat_qr_f32()
226 //printf("--\r\n"); in _arm_mat_qr_f32()
228 /* A(col:,col:) - beta v tmpb */ in _arm_mat_qr_f32()
229 pa = p; in _arm_mat_qr_f32()
230 for(j=0;j<pSrc->numRows-col; j++) in _arm_mat_qr_f32()
232 float32_t f = -beta * pTmpA[j]; in _arm_mat_qr_f32()
235 blkCnt2 = (pSrc->numCols-col) >> 2; in _arm_mat_qr_f32()
246 blkCnt2--; in _arm_mat_qr_f32()
248 blkCnt2 = (pSrc->numCols-col) & 3; in _arm_mat_qr_f32()
265 pa = p + pOutR->numCols; in _arm_mat_qr_f32()
266 for(k=0;k<pSrc->numRows-col-1; k++) in _arm_mat_qr_f32()
269 pa += pOutR->numCols; in _arm_mat_qr_f32()
272 p += 1 + pOutR->numCols; in _arm_mat_qr_f32()
280 memset(pOutQ->pData,0,sizeof(float32_t)*pOutQ->numRows*pOutQ->numRows); in _arm_mat_qr_f32()
282 pa = pOutQ->pData; in _arm_mat_qr_f32()
283 for(col=0 ; col < pOutQ->numCols; col++) in _arm_mat_qr_f32()
286 pa += pOutQ->numCols+1; in _arm_mat_qr_f32()
289 nb = pOutQ->numRows - pOutQ->numCols + 1; in _arm_mat_qr_f32()
291 pc = pOutTau + pOutQ->numCols - 1; in _arm_mat_qr_f32()
292 for(col=0 ; col < pOutQ->numCols; col++) in _arm_mat_qr_f32()
299 pos = pSrc->numRows - nb; in _arm_mat_qr_f32()
300 p = pOutQ->pData + pos + pOutQ->numCols*pos ; in _arm_mat_qr_f32()
307 /* v.T A(col:,col:) -> tmpb */ in _arm_mat_qr_f32()
310 pa = p; in _arm_mat_qr_f32()
313 blkCnt2 = (pOutQ->numRows-pos) >> 2; in _arm_mat_qr_f32()
323 blkCnt2--; in _arm_mat_qr_f32()
325 blkCnt2 = (pOutQ->numRows-pos) & 3; in _arm_mat_qr_f32()
341 pa1 = pa0 + pOutQ->numRows; in _arm_mat_qr_f32()
342 pa2 = pa1 + pOutQ->numRows; in _arm_mat_qr_f32()
343 pa3 = pa2 + pOutQ->numRows; in _arm_mat_qr_f32()
346 blkCnt = (pOutQ->numRows-pos - 1) >> 2; in _arm_mat_qr_f32()
352 blkCnt2 = (pOutQ->numRows-pos) >> 2; in _arm_mat_qr_f32()
377 blkCnt2--; in _arm_mat_qr_f32()
379 blkCnt2 = (pOutQ->numRows-pos) & 3; in _arm_mat_qr_f32()
407 pa0 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
408 pa1 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
409 pa2 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
410 pa3 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
414 blkCnt--; in _arm_mat_qr_f32()
418 for(;k<pOutQ->numRows-pos; k++) in _arm_mat_qr_f32()
421 blkCnt2 = (pOutQ->numRows-pos) >> 2; in _arm_mat_qr_f32()
432 blkCnt2--; in _arm_mat_qr_f32()
434 blkCnt2 = (pOutQ->numRows-pos) & 3; in _arm_mat_qr_f32()
451 pa = p; in _arm_mat_qr_f32()
452 beta = *pc--; in _arm_mat_qr_f32()
453 for(j=0;j<pOutQ->numRows-pos; j++) in _arm_mat_qr_f32()
455 float32_t f = -beta * pTmpA[j]; in _arm_mat_qr_f32()
458 blkCnt2 = (pOutQ->numCols-pos) >> 2; in _arm_mat_qr_f32()
469 blkCnt2--; in _arm_mat_qr_f32()
471 blkCnt2 = (pOutQ->numCols-pos) & 3; in _arm_mat_qr_f32()
523 float32_t *p; in _arm_mat_qr_f32() local
525 if (pSrc->numRows < pSrc->numCols) in _arm_mat_qr_f32()
530 memcpy(pOutR->pData,pSrc->pData,pSrc->numCols * pSrc->numRows*sizeof(float32_t)); in _arm_mat_qr_f32()
531 pOutR->numCols = pSrc->numCols; in _arm_mat_qr_f32()
532 pOutR->numRows = pSrc->numRows; in _arm_mat_qr_f32()
534 p = pOutR->pData; in _arm_mat_qr_f32()
537 for(col=0 ; col < pSrc->numCols; col++) in _arm_mat_qr_f32()
543 beta = arm_householder_f32(pTmpA,threshold,pSrc->numRows - col,pTmpA); in _arm_mat_qr_f32()
548 /* v.T A(col:,col:) -> tmpb */ in _arm_mat_qr_f32()
550 pa = p; in _arm_mat_qr_f32()
551 for(j=0;j<pSrc->numCols-col; j++) in _arm_mat_qr_f32()
560 pa1 = pa0 + pSrc->numCols; in _arm_mat_qr_f32()
561 pa2 = pa1 + pSrc->numCols; in _arm_mat_qr_f32()
562 pa3 = pa2 + pSrc->numCols; in _arm_mat_qr_f32()
565 blkCnt = (pSrc->numRows-col - 1) >> 2; in _arm_mat_qr_f32()
571 for(j=0;j<pSrc->numCols-col; j++) in _arm_mat_qr_f32()
582 pa0 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
583 pa1 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
584 pa2 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
585 pa3 += col + 3*pSrc->numCols; in _arm_mat_qr_f32()
589 blkCnt--; in _arm_mat_qr_f32()
593 for(;k<pSrc->numRows-col; k++) in _arm_mat_qr_f32()
595 for(j=0;j<pSrc->numCols-col; j++) in _arm_mat_qr_f32()
604 /* A(col:,col:) - beta v tmpb */ in _arm_mat_qr_f32()
605 pa = p; in _arm_mat_qr_f32()
606 for(j=0;j<pSrc->numRows-col; j++) in _arm_mat_qr_f32()
610 for(i=0;i<pSrc->numCols-col; i++) in _arm_mat_qr_f32()
612 *pa = *pa - f * pTmpB[i] ; in _arm_mat_qr_f32()
619 pa = p + pOutR->numCols; in _arm_mat_qr_f32()
620 for(k=0;k<pSrc->numRows-col-1; k++) in _arm_mat_qr_f32()
623 pa += pOutR->numCols; in _arm_mat_qr_f32()
626 p += 1 + pOutR->numCols; in _arm_mat_qr_f32()
634 memset(pOutQ->pData,0,sizeof(float32_t)*pOutQ->numRows*pOutQ->numRows); in _arm_mat_qr_f32()
636 pa = pOutQ->pData; in _arm_mat_qr_f32()
637 for(col=0 ; col < pOutQ->numCols; col++) in _arm_mat_qr_f32()
640 pa += pOutQ->numCols+1; in _arm_mat_qr_f32()
643 nb = pOutQ->numRows - pOutQ->numCols + 1; in _arm_mat_qr_f32()
645 pc = pOutTau + pOutQ->numCols - 1; in _arm_mat_qr_f32()
646 for(col=0 ; col < pOutQ->numCols; col++) in _arm_mat_qr_f32()
650 pos = pSrc->numRows - nb; in _arm_mat_qr_f32()
651 p = pOutQ->pData + pos + pOutQ->numCols*pos ; in _arm_mat_qr_f32()
658 /* v.T A(col:,col:) -> tmpb */ in _arm_mat_qr_f32()
661 pa = p; in _arm_mat_qr_f32()
662 for(j=0;j<pOutQ->numRows-pos; j++) in _arm_mat_qr_f32()
670 pa1 = pa0 + pOutQ->numRows; in _arm_mat_qr_f32()
671 pa2 = pa1 + pOutQ->numRows; in _arm_mat_qr_f32()
672 pa3 = pa2 + pOutQ->numRows; in _arm_mat_qr_f32()
675 blkCnt = (pOutQ->numRows-pos - 1) >> 2; in _arm_mat_qr_f32()
681 for(j=0;j<pOutQ->numRows-pos; j++) in _arm_mat_qr_f32()
692 pa0 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
693 pa1 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
694 pa2 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
695 pa3 += pos + 3*pOutQ->numRows; in _arm_mat_qr_f32()
699 blkCnt--; in _arm_mat_qr_f32()
703 for(;k<pOutQ->numRows-pos; k++) in _arm_mat_qr_f32()
705 for(j=0;j<pOutQ->numRows-pos; j++) in _arm_mat_qr_f32()
714 pa = p; in _arm_mat_qr_f32()
715 beta = *pc--; in _arm_mat_qr_f32()
716 for(j=0;j<pOutQ->numRows-pos; j++) in _arm_mat_qr_f32()
720 for(i=0;i<pOutQ->numCols-pos; i++) in _arm_mat_qr_f32()
722 *pa = *pa - f * pTmpB[i] ; in _arm_mat_qr_f32()