Home | History | Annotate | Download | only in blas

Lines Matching refs:Matrix

64     if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
65 else matrix(c, *m, *n, *ldc) *= beta;
175 matrix(b,*m,*n,*ldb) *= alpha;
269 Matrix<Scalar,Dynamic,Dynamic,ColMajor> tmp = matrix(b,*m,*n,*ldb);
270 matrix(b,*m,*n,*ldb).setZero();
310 if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
311 else matrix(c, *m, *n, *ldc) *= beta;
321 // FIXME add support for symmetric complex matrix
322 Matrix<Scalar,Dynamic,Dynamic,ColMajor> matA(size,size);
325 matA.triangularView<Upper>() = matrix(a,size,size,*lda);
326 matA.triangularView<Lower>() = matrix(a,size,size,*lda).transpose();
330 matA.triangularView<Lower>() = matrix(a,size,size,*lda);
331 matA.triangularView<Upper>() = matrix(a,size,size,*lda).transpose();
334 matrix(c, *m, *n, *ldc) += alpha * matA * matrix(b, *m, *n, *ldb);
336 matrix(c, *m, *n, *ldc) += alpha * matrix(b, *m, *n, *ldb) * matA;
399 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
400 else matrix(c, *n, *n, *ldc).triangularView<Upper>() *= beta;
402 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
403 else matrix(c, *n, *n, *ldc).triangularView<Lower>() *= beta;
410 // FIXME add support for symmetric complex matrix
414 matrix(c, *n, *n, *ldc).triangularView<Upper>() += alpha * matrix(a,*n,*k,*lda) * matrix(a,*n,*k,*lda).transpose();
416 matrix(c, *n, *n, *ldc).triangularView<Upper>() += alpha * matrix(a,*k,*n,*lda).transpose() * matrix(a,*k,*n,*lda);
421 matrix(c, *n, *n, *ldc).triangularView<Lower>() += alpha * matrix(a,*n,*k,*lda) * matrix(a,*n,*k,*lda).transpose();
423 matrix(c, *n, *n, *ldc).triangularView<Lower>() += alpha * matrix(a,*k,*n,*lda).transpose() * matrix(a,*k,*n,*lda);
462 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
463 else matrix(c, *n, *n, *ldc).triangularView<Upper>() *= beta;
465 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
466 else matrix(c, *n, *n, *ldc).triangularView<Lower>() *= beta;
476 matrix(c, *n, *n, *ldc).triangularView<Upper>()
477 += alpha *matrix(a, *n, *k, *lda)*matrix(b, *n, *k, *ldb).transpose()
478 + alpha*matrix(b, *n, *k, *ldb)*matrix(a, *n, *k, *lda).transpose();
481 matrix(c, *n, *n, *ldc).triangularView<Lower>()
482 += alpha*matrix(a, *n, *k, *lda)*matrix(b, *n, *k, *ldb).transpose()
483 + alpha*matrix(b, *n, *k, *ldb)*matrix(a, *n, *k, *lda).transpose();
488 matrix(c, *n, *n, *ldc).triangularView<Upper>()
489 += alpha*matrix(a, *k, *n, *lda).transpose()*matrix(b, *k, *n, *ldb)
490 + alpha*matrix(b, *k, *n, *ldb).transpose()*matrix(a, *k, *n, *lda);
492 matrix(c, *n, *n, *ldc).triangularView<Lower>()
493 += alpha*matrix(a, *k, *n, *lda).transpose()*matrix(b, *k, *n, *ldb)
494 + alpha*matrix(b, *k, *n, *ldb).transpose()*matrix(a, *k, *n, *lda);
527 if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
528 else if(beta!=Scalar(1)) matrix(c, *m, *n, *ldc) *= beta;
548 if(UPLO(*uplo)==UP) matrix(c,*m,*n,*ldc) += alpha * matrix(b,*m,*n,*ldb) * matrix(a,*n,*n,*lda).selfadjointView<Upper>();/*internal::product_selfadjoint_matrix<Scalar,DenseIndex,ColMajor,false,false, RowMajor,true,Conj, ColMajor>
607 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
608 else matrix(c, *n, *n, *ldc).triangularView<StrictlyUpper>() *= beta;
610 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
611 else matrix(c, *n, *n, *ldc).triangularView<StrictlyLower>() *= beta;
615 matrix(c, *n, *n, *ldc).diagonal().real() *= beta;
616 matrix(c, *n, *n, *ldc).diagonal().imag().setZero();
624 matrix(c, *n, *n, *ldc).diagonal().imag().setZero();
656 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
657 else matrix(c, *n, *n, *ldc).triangularView<StrictlyUpper>() *= beta;
659 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
660 else matrix(c, *n, *n, *ldc).triangularView<StrictlyLower>() *= beta;
664 matrix(c, *n, *n, *ldc).diagonal().real() *= beta;
665 matrix(c, *n, *n, *ldc).diagonal().imag().setZero();
669 matrix(c, *n, *n, *ldc).diagonal().imag().setZero();
678 matrix(c, *n, *n, *ldc).triangularView<Upper>()
679 += alpha *matrix(a, *n, *k, *lda)*matrix(b, *n, *k, *ldb).adjoint()
680 + numext::conj(alpha)*matrix(b, *n, *k, *ldb)*matrix(a, *n, *k, *lda).adjoint();
683 matrix(c, *n, *n, *ldc).triangularView<Lower>()
684 += alpha*matrix(a, *n, *k, *lda)*matrix(b, *n, *k, *ldb).adjoint()
685 + numext::conj(alpha)*matrix(b, *n, *k, *ldb)*matrix(a, *n, *k, *lda).adjoint();
690 matrix(c, *n, *n, *ldc).triangularView<Upper>()
691 += alpha*matrix(a, *k, *n, *lda).adjoint()*matrix(b, *k, *n, *ldb)
692 + numext::conj(alpha)*matrix(b, *k, *n, *ldb).adjoint()*matrix(a, *k, *n, *lda);
694 matrix(c, *n, *n, *ldc).triangularView<Lower>()
695 += alpha*matrix(a, *k, *n, *lda).adjoint()*matrix(b, *k, *n, *ldb)
696 + numext::conj(alpha)*matrix(b, *k, *n, *ldb).adjoint()*matrix(a, *k, *n, *lda);