/dports/math/octave/octave-6.4.0/liboctave/numeric/ |
H A D | sparse-qr.cc | 340 j < b_nc; in C() 460 Matrix x (nc, b_nc); in tall_solve() 635 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in tall_solve() 728 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in wide_solve() 1105 j < b_nc; in C() 1297 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in tall_solve() 1417 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in wide_solve() 1655 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in tall_solve() 1759 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in wide_solve() 1994 octave_idx_type sz = x_nz * (b_nc - i) / b_nc; in tall_solve() [all …]
|
H A D | sparse-dmsolve.cc | 296 octave_idx_type b_nc = b.cols (); in dmsolve_permute() local 300 a.resize (dim_vector (b_nr, b_nc)); in dmsolve_permute() 304 for (octave_idx_type j = 0; j < b_nc; j++) in dmsolve_permute() 320 octave_idx_type b_nc = b.cols (); in dmsolve_permute() local 325 a = MSparse<RT> (b_nr, b_nc, b_nz); in dmsolve_permute() 333 for (octave_idx_type j = 0; j < b_nc; j++) in dmsolve_permute() 378 octave_idx_type b_nc = b.cols (); 384 if (nr == 0 || nc == 0 || b_nc == 0) 385 retval = RT (nc, b_nc, 0.0); 418 retval.resize (nc, b_nc); [all …]
|
H A D | EIG.cc | 402 F77_INT b_nc = octave::to_f77_int (b.cols ()); in init() local 404 if (n != a_nc || nb != b_nc) in init() 539 F77_INT b_nc = octave::to_f77_int (b.cols ()); in symmetric_init() local 541 if (n != a_nc || nb != b_nc) in symmetric_init() 609 F77_INT b_nc = octave::to_f77_int (b.cols ()); in init() local 611 if (n != a_nc || nb != b_nc) in init() 717 F77_INT b_nc = octave::to_f77_int (b.cols ()); in hermitian_init() local 719 if (n != a_nc || nb != b_nc) in hermitian_init()
|
H A D | fEIG.cc | 399 F77_INT b_nc = octave::to_f77_int (b.cols ()); in init() local 401 if (n != a_nc || nb != b_nc) in init() 534 F77_INT b_nc = octave::to_f77_int (b.cols ()); in symmetric_init() local 536 if (n != a_nc || nb != b_nc) in symmetric_init() 604 F77_INT b_nc = octave::to_f77_int (b.cols ()); in init() local 606 if (n != a_nc || nb != b_nc) in init() 708 F77_INT b_nc = octave::to_f77_int (b.cols ()); in hermitian_init() local 710 if (n != a_nc || nb != b_nc) in hermitian_init()
|
/dports/math/octave/octave-6.4.0/liboctave/array/ |
H A D | fCDiagMatrix.cc | 396 octave_idx_type b_nc = b.cols (); in operator *() local 399 octave::err_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc); in operator *() 401 FloatComplexDiagMatrix c (a_nr, b_nc); in operator *() 421 octave_idx_type b_nc = b.cols (); in operator *() local 426 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in operator *() 429 FloatComplexDiagMatrix c (a_nr, b_nc); in operator *() 431 octave_idx_type len = (a_nr < b_nc ? a_nr : b_nc); in operator *() 451 octave_idx_type b_nc = b.cols (); in operator *() local 456 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in operator *() 459 FloatComplexDiagMatrix c (a_nr, b_nc); in operator *() [all …]
|
H A D | MSparse.cc | 40 octave_idx_type b_nc = b.cols (); in plus_or_minus() local 42 if (a_nr != b_nr || a_nc != b_nc) in plus_or_minus() 259 octave_idx_type b_nc = b.cols (); in plus_or_minus() local 285 else if (b_nr == 1 && b_nc == 1) in plus_or_minus() 306 else if (a_nr != b_nr || a_nc != b_nc) in plus_or_minus() 391 octave_idx_type b_nc = b.cols (); in product() local 396 r = MSparse<T> (b_nr, b_nc); in product() 410 else if (b_nr == 1 && b_nc == 1) in product() 427 else if (a_nr != b_nr || a_nc != b_nc) in product() 489 octave_idx_type b_nc = b.cols (); in quotient() local [all …]
|
H A D | CDiagMatrix.cc | 412 octave_idx_type b_nc = b.cols (); in operator *() local 415 octave::err_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc); in operator *() 417 ComplexDiagMatrix c (a_nr, b_nc); in operator *() 437 octave_idx_type b_nc = b.cols (); in operator *() local 442 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in operator *() 445 ComplexDiagMatrix c (a_nr, b_nc); in operator *() 447 octave_idx_type len = (a_nr < b_nc ? a_nr : b_nc); in operator *() 467 octave_idx_type b_nc = b.cols (); in operator *() local 472 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in operator *() 475 ComplexDiagMatrix c (a_nr, b_nc); in operator *() [all …]
|
H A D | dSparse.cc | 1503 retval.resize (nc, b_nc); in utsolve() 1583 retval.resize (nc, b_nc); in utsolve() 2013 retval.resize (nc, b_nc); in utsolve() 2094 retval.resize (nc, b_nc); in utsolve() 2529 retval.resize (nc, b_nc); in ltsolve() 3082 retval.resize (nc, b_nc); in ltsolve() 4693 sz = (static_cast<double> (b_nc) - j) / b_nc in bsolve() 6165 sz = (static_cast<double> (b_nc) - j) / b_nc in fsolve() 6686 sz = (static_cast<double> (b_nc) - j) / b_nc in fsolve() 7768 r.resize (b_nr, b_nc); in min() [all …]
|
H A D | CSparse.cc | 1568 retval.resize (nc, b_nc); in utsolve() 1648 retval.resize (nc, b_nc); in utsolve() 2079 retval.resize (nc, b_nc); in utsolve() 2159 retval.resize (nc, b_nc); in utsolve() 2591 retval.resize (nc, b_nc); in ltsolve() 3141 retval.resize (nc, b_nc); in ltsolve() 4744 sz = (static_cast<double> (b_nc) - j) / b_nc in bsolve() 6180 sz = (static_cast<double> (b_nc) - j) / b_nc in fsolve() 6680 sz = (static_cast<double> (b_nc) - j) / b_nc in fsolve() 7809 r.resize (b_nr, b_nc); in min() [all …]
|
H A D | fMatrix.cc | 1206 if (nr == 0 || nc == 0 || b_nc == 0) in utsolve() 1207 retval = FloatMatrix (nc, b_nc, 0.0); in utsolve() 1309 if (nr == 0 || nc == 0 || b_nc == 0) in ltsolve() 1310 retval = FloatMatrix (nc, b_nc, 0.0); in ltsolve() 1411 if (nr == 0 || b_nc == 0) in fsolve() 1953 if (m == 0 || n == 0 || b_nc == 0) in lssolve() 1954 retval = FloatMatrix (n, b_nc, 0.0); in lssolve() 2073 retval = FloatMatrix (n, b_nc, in lssolve() 2799 retval = FloatMatrix (a_nr, b_nc); in xgemm() 2821 retval = FloatMatrix (a_nr, b_nc); in xgemm() [all …]
|
H A D | dMatrix.cc | 1196 if (nr == 0 || nc == 0 || b_nc == 0) in utsolve() 1197 retval = Matrix (nc, b_nc, 0.0); in utsolve() 1294 if (nr == 0 || nc == 0 || b_nc == 0) in ltsolve() 1295 retval = Matrix (nc, b_nc, 0.0); in ltsolve() 1927 if (m == 0 || n == 0 || b_nc == 0) in lssolve() 1928 retval = Matrix (n, b_nc, 0.0); in lssolve() 2042 retval = Matrix (n, b_nc, 0.0); in lssolve() 2788 retval = Matrix (a_nr, b_nc, 0.0); in xgemm() 2793 retval = Matrix (a_nr, b_nc); in xgemm() 2815 retval = Matrix (a_nr, b_nc); in xgemm() [all …]
|
H A D | CMatrix.cc | 1511 F77_INT b_nc = octave::to_f77_int (b.cols ()); in utsolve() local 1517 if (nr == 0 || nc == 0 || b_nc == 0) in utsolve() 1615 if (nr == 0 || nc == 0 || b_nc == 0) in ltsolve() 1713 if (nr == 0 || b_nc == 0) in fsolve() 2269 if (m == 0 || n == 0 || b_nc == 0) in lssolve() 2394 retval = ComplexMatrix (n, b_nc, in lssolve() 3341 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in xgemm() 3342 retval = ComplexMatrix (a_nr, b_nc, 0.0); in xgemm() 3352 retval = ComplexMatrix (a_nr, b_nc, 0.0); in xgemm() 3393 if (b_nc == 1 && a_nr == 1) in xgemm() [all …]
|
H A D | fCMatrix.cc | 1513 F77_INT b_nc = octave::to_f77_int (b.cols ()); in utsolve() local 1519 if (nr == 0 || nc == 0 || b_nc == 0) in utsolve() 1615 F77_INT b_nc = octave::to_f77_int (b.cols ()); in ltsolve() local 1621 if (nr == 0 || nc == 0 || b_nc == 0) in ltsolve() 1717 F77_INT b_nc = octave::to_f77_int (b.cols ()); in fsolve() local 1723 if (nr == 0 || b_nc == 0) in fsolve() 2295 if (m == 0 || n == 0 || b_nc == 0) in lssolve() 2420 retval = FloatComplexMatrix (n, b_nc, in lssolve() 3372 if (a_nr == 0 || a_nc == 0 || b_nc == 0) in xgemm() 3424 if (b_nc == 1 && a_nr == 1) in xgemm() [all …]
|
H A D | dDiagMatrix.cc | 310 octave_idx_type b_nc = b.cols (); in operator *() local 313 octave::err_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc); in operator *() 315 DiagMatrix c (a_nr, b_nc); in operator *()
|
H A D | fDiagMatrix.cc | 288 octave_idx_type b_nc = b.cols (); in operator *() local 291 octave::err_nonconformant ("operator *", a_nr, a_nc, b_nr, b_nc); in operator *() 293 FloatDiagMatrix c (a_nr, b_nc); in operator *()
|
/dports/science/quantum-espresso/q-e-qe-6.7.0/GWW/simple_ip/ |
H A D | diagonalization.f90 | 29 COMPLEX(kind=DP), DIMENSION(:,:,:), ALLOCATABLE :: b_nc , csca_nc local 75 allocate(b_nc(sh%nkb,sh%npol,sh%ntot_e)) 77 call trilinear_parallel_nc(kptns,ik,sh,b_nc) 79 b_nc(1:sh%nkb,1:sh%npol,1:sh%ntot_e) = sh%beck_nc(1:sh%nkb, 1:sh%npol, 1:sh%ntot_e,ik) 107 sh%deeq_nc(ih,jh,na,1)*b_nc(jkb,1,j)+ & 108 sh%deeq_nc(ih,jh,na,2)*b_nc(jkb,2,j) 110 sh%deeq_nc(ih,jh,na,3)*b_nc(jkb,1,j)+& 111 sh%deeq_nc(ih,jh,na,4)*b_nc(jkb,2,j) 120 call ZGEMM ('C', 'N', sh%ntot_e, sh%ntot_e, sh%npol*sh%nkb, ( 1.d0, 0.d0 ) , b_nc, & 155 deallocate(b_nc)
|
/dports/math/octave-forge-communications/communications-1.2.3/src/ |
H A D | galois.cc | 492 int b_nc = b.cols (); in elem_pow() local 512 for (int j = 0; j < b_nc; j++) in elem_pow() 519 else if (b_nr == 1 && b_nc == 1) in elem_pow() 557 int b_nc = b.cols (); in elem_pow() local 559 if (b_nr == 1 && b_nc == 1) in elem_pow() 562 if (a_nr != b_nr || a_nc != b_nc) in elem_pow() 767 int b_nc = b.cols (); in operator *() local 1246 int b_nc = b.cols (); in solve() local 1427 int b_nc = b.cols (); in xdiv() local 1430 if (a_nc != b_nc) in xdiv() [all …]
|
/dports/math/octave/octave-6.4.0/libinterp/corefcn/ |
H A D | sparse-xpow.cc | 415 octave_idx_type b_nc = b.cols (); in elem_xpow() local 420 if (nr != b_nr || nc != b_nc) in elem_xpow() 421 octave::err_nonconformant ("operator .^", nr, nc, b_nr, b_nc); in elem_xpow() 517 octave_idx_type b_nc = b.cols (); in elem_xpow() local 522 if (nr != b_nr || nc != b_nc) in elem_xpow() 523 octave::err_nonconformant ("operator .^", nr, nc, b_nr, b_nc); in elem_xpow() 657 octave_idx_type b_nc = b.cols (); in elem_xpow() local 662 if (nr != b_nr || nc != b_nc) in elem_xpow() 663 octave::err_nonconformant ("operator .^", nr, nc, b_nr, b_nc); in elem_xpow() 724 octave_idx_type b_nc = b.cols (); in elem_xpow() local [all …]
|
H A D | xpow.cc | 803 octave_idx_type b_nc = b.cols (); in elem_xpow() local 805 if (nr != b_nr || nc != b_nc) in elem_xpow() 883 octave_idx_type b_nc = b.cols (); in elem_xpow() local 885 if (nr != b_nr || nc != b_nc) in elem_xpow() 1021 octave_idx_type b_nc = b.cols (); in elem_xpow() local 1023 if (nr != b_nr || nc != b_nc) in elem_xpow() 1071 if (nr != b_nr || nc != b_nc) in elem_xpow() 2151 if (nr != b_nr || nc != b_nc) in elem_xpow() 2231 if (nr != b_nr || nc != b_nc) in elem_xpow() 2330 if (nr != b_nr || nc != b_nc) in elem_xpow() [all …]
|
H A D | sylvester.cc | 78 octave_idx_type b_nc = arg_b.columns (); variable 99 if (b_nr != b_nc)
|
H A D | sparse-xdiv.cc | 62 octave_idx_type b_nc = b.cols (); in mx_leftdiv_conform() local 64 octave::err_nonconformant (R"(operator \)", a_nr, a_nc, b_nr, b_nc); in mx_leftdiv_conform() 91 octave_idx_type b_nc = b.cols (); in mx_div_conform() local 93 if (a_nc != b_nc) in mx_div_conform() 98 octave::err_nonconformant ("operator /", a_nr, a_nc, b_nr, b_nc); in mx_div_conform()
|
H A D | xdiv.cc | 69 octave_idx_type b_nc = b.cols (); in mx_leftdiv_conform() local 71 octave::err_nonconformant (R"(operator \)", a_nr, a_nc, b_nr, b_nc); in mx_leftdiv_conform() 90 octave_idx_type b_nc = b.cols (); in mx_div_conform() local 92 if (a_nc != b_nc) in mx_div_conform() 97 octave::err_nonconformant ("operator /", a_nr, a_nc, b_nr, b_nc); in mx_div_conform()
|
H A D | qz.cc | 275 F77_INT b_nc = octave::to_f77_int (args(1).columns ()); variable 277 if (nn != b_nc || nn != b_nr)
|
/dports/math/octave-forge-sparsersb/sparsersb-1.0.9/src/ |
H A D | sparsersb.cc | 1623 const octave_idx_type nrhs = b_nc; in scale_rows() 1637 const octave_idx_type nrhs = b_nc; in scale_rows() 1659 const octave_idx_type nrhs = b_nc; in rsboi_spmm() 1674 const octave_idx_type nrhs = b_nc; in rsboi_spmm() 1697 const octave_idx_type nrhs = b_nc; in rsboi_spmtm() 1712 const octave_idx_type nrhs = b_nc; in rsboi_spmtm() 1735 const octave_idx_type nrhs = b_nc; in rsboi_spmhm() 1775 const octave_idx_type nrhs = b_nc; in rsboi_spmm() 1815 const octave_idx_type nrhs = b_nc; in rsboi_spmtm() 1852 const octave_idx_type nrhs = b_nc; in rsboi_spmhm() [all …]
|
/dports/math/m4ri/m4ri-release-20200125/m4ri/ |
H A D | brilliantrussian.c | 1064 rci_t const b_nc = B->ncols; in _mzd_mul_m4rm() local 1066 if (b_nc < m4ri_radix-10 || a_nr < 16) { in _mzd_mul_m4rm() 1086 rci_t const klog = round(0.75 * log2_floor(MIN(MIN(a_nr,a_nc),b_nc))); in _mzd_mul_m4rm() 1103 Talign[z] = mzd_init(__M4RI_TWOPOW(k), b_nc+m4ri_radix); in _mzd_mul_m4rm() 1104 …T[z] = mzd_init_window(Talign[z], 0, c_align*m4ri_radix, Talign[z]->nrows, b_nc + c_align*m4ri_rad… in _mzd_mul_m4rm() 1106 T[z] = mzd_init(__M4RI_TWOPOW(k), b_nc); in _mzd_mul_m4rm()
|