Lines Matching defs:beta
39 Scalar beta = *reinterpret_cast<Scalar*>(pbeta);
53 if(beta!=Scalar(1))
55 if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
56 else matrix(c, *m, *n, *ldc) *= beta;
236 // c = alpha*a*b + beta*c for side = 'L'or'l'
237 // c = alpha*b*a + beta*c for side = 'R'or'r
240 // std::cerr << "in symm " << *side << " " << *uplo << " " << *m << "x" << *n << " lda:" << *lda << " ldb:" << *ldb << " ldc:" << *ldc << " alpha:" << *palpha << " beta:" << *pbeta << "\n";
245 Scalar beta = *reinterpret_cast<Scalar*>(pbeta);
258 if(beta!=Scalar(1))
260 if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
261 else matrix(c, *m, *n, *ldc) *= beta;
303 // c = alpha*a*a' + beta*c for op = 'N'or'n'
304 // c = alpha*a'*a + beta*c for op = 'T'or't','C'or'c'
333 Scalar beta = *reinterpret_cast<Scalar*>(pbeta);
345 if(beta!=Scalar(1))
348 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
349 else matrix(c, *n, *n, *ldc).triangularView<Upper>() *= beta;
351 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
352 else matrix(c, *n, *n, *ldc).triangularView<Lower>() *= beta;
379 // c = alpha*a*b' + alpha*b*a' + beta*c for op = 'N'or'n'
380 // c = alpha*a'*b + alpha*b'*a + beta*c for op = 'T'or't'
387 Scalar beta = *reinterpret_cast<Scalar*>(pbeta);
400 if(beta!=Scalar(1))
403 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
404 else matrix(c, *n, *n, *ldc).triangularView<Upper>() *= beta;
406 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
407 else matrix(c, *n, *n, *ldc).triangularView<Lower>() *= beta;
444 // c = alpha*a*b + beta*c for side = 'L'or'l'
445 // c = alpha*b*a + beta*c for side = 'R'or'r
452 Scalar beta = *reinterpret_cast<Scalar*>(pbeta);
454 // std::cerr << "in hemm " << *side << " " << *uplo << " " << *m << " " << *n << " " << alpha << " " << *lda << " " << beta << " " << *ldc << "\n";
467 if(beta==Scalar(0)) matrix(c, *m, *n, *ldc).setZero();
468 else if(beta!=Scalar(1)) matrix(c, *m, *n, *ldc) *= beta;
499 // c = alpha*a*conj(a') + beta*c for op = 'N'or'n'
500 // c = alpha*conj(a')*a + beta*c for op = 'C'or'c'
524 RealScalar beta = *pbeta;
526 // std::cerr << "in herk " << *uplo << " " << *op << " " << *n << " " << *k << " " << alpha << " " << *lda << " " << beta << " " << *ldc << "\n";
540 if(beta!=RealScalar(1))
543 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
544 else matrix(c, *n, *n, *ldc).triangularView<StrictlyUpper>() *= beta;
546 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
547 else matrix(c, *n, *n, *ldc).triangularView<StrictlyLower>() *= beta;
549 if(beta!=Scalar(0))
551 matrix(c, *n, *n, *ldc).diagonal().real() *= beta;
564 // c = alpha*a*conj(b') + conj(alpha)*b*conj(a') + beta*c, for op = 'N'or'n'
565 // c = alpha*conj(a')*b + conj(alpha)*conj(b')*a + beta*c, for op = 'C'or'c'
572 RealScalar beta = *pbeta;
585 if(beta!=RealScalar(1))
588 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Upper>().setZero();
589 else matrix(c, *n, *n, *ldc).triangularView<StrictlyUpper>() *= beta;
591 if(beta==Scalar(0)) matrix(c, *n, *n, *ldc).triangularView<Lower>().setZero();
592 else matrix(c, *n, *n, *ldc).triangularView<StrictlyLower>() *= beta;
594 if(beta!=Scalar(0))
596 matrix(c, *n, *n, *ldc).diagonal().real() *= beta;