/dports/math/scilab/scilab-6.1.1/scilab/modules/linear_algebra/macros/ |
H A D | fullrfk.sci | 14 function [Bk,Ck]=fullrfk(A,k) 16 // Bk*Ck=A^k where Bk is full column rank and Ck full row rank. 17 // One has Range(Bk)=Range(A^k) and ker(Ck)=ker(A^k). 28 Bk=eye(n,n);Ck=Bk; 33 [Bk,Ck]=fullrf(A); 36 [Bk,Ck]=fullrf(A);B=Bk;C=Ck; 39 Bk=Bk*B;Ck=C*Ck; // Bk*Ck = A^k (Full rank factorization)
|
H A D | projspec.sci | 40 Ck=C;Bk=B;//write(%io(2),[dim,rcond(C*B)],'(7x,f3.0,6x,e9.3)'); 41 tst=rcond(Ck*Bk); 42 if size(Ck,1)==1 then tst=norm(Ck*Bk,1);end 52 Bk=Bk*B;Ck=C*Ck; // Bk*Ck = A^k (Full rank factorization) 56 M=inv((C*B)**index); // M=inv(Ck*Bk); (Alternative computation) 57 P=eye(n,n)-Bk*M*Ck; // P=eye(n,n)-Bk*inv(Ck*Bk)*Ck; 58 S=Bk*M*inv(C*B)*Ck; // S=inv(A-P-D)+P
|
/dports/math/igraph/igraph-0.9.5/vendor/glpk/misc/ |
H A D | fp2rat.c | 116 double xk, Akm1, Ak, Bkm1, Bk, ak, bk, fk, temp; in fp2rat() local 130 Bk = 1.0; in fp2rat() 145 temp = bk * Bk + ak * Bkm1; in fp2rat() 146 Bkm1 = Bk, Bk = temp; in fp2rat() 149 fk = Ak / Bk; in fp2rat() 152 DBL_DIG, Ak, DBL_DIG, Bk, DBL_DIG, fk); in fp2rat() 158 *q = Bk; in fp2rat()
|
/dports/math/glpk/glpk-5.0/src/misc/ |
H A D | fp2rat.c | 116 double xk, Akm1, Ak, Bkm1, Bk, ak, bk, fk, temp; in fp2rat() local 130 Bk = 1.0; in fp2rat() 145 temp = bk * Bk + ak * Bkm1; in fp2rat() 146 Bkm1 = Bk, Bk = temp; in fp2rat() 149 fk = Ak / Bk; in fp2rat() 152 DBL_DIG, Ak, DBL_DIG, Bk, DBL_DIG, fk); in fp2rat() 158 *q = Bk; in fp2rat()
|
/dports/math/libxsmm/libxsmm-1.16.3/samples/xgemm/ |
H A D | bf16sgemm.c | 76 int Bk = (*k)/bk; in sgemm_trup() local 83 size_t* pob = poa + (Bm*Bk); in sgemm_trup() 84 size_t* poc = pob + (Bn*Bk); in sgemm_trup() 115 for ( lk1 = 0; lk1 < Bk; ++lk1 ) { in sgemm_trup() 121 for ( lk1 = 0; lk1 < Bk; ++lk1 ) { in sgemm_trup() 149 for ( lk1 = 0; lk1 < Bk; ++lk1 ) { in sgemm_trup() 246 int Bk = (*k)/bk; in bf16sgemm_trup() local 253 size_t* pob = poa + (Bm*Bk); in bf16sgemm_trup() 254 size_t* poc = pob + (Bn*Bk); in bf16sgemm_trup() 285 for ( lk1 = 0; lk1 < Bk; ++lk1 ) { in bf16sgemm_trup() [all …]
|
/dports/devel/ppl/ppl-1.2/demos/ppl_lpsol/ |
H A D | glpk_set_d_eps.c | 218 double xk, Akm1, Ak, Bkm1, Bk, ak, bk, fk, temp; in fp2rat() local 233 Bk = 1.0; in fp2rat() 248 temp = bk * Bk + ak * Bkm1; in fp2rat() 249 Bkm1 = Bk, Bk = temp; in fp2rat() 252 fk = Ak / Bk; in fp2rat() 254 print("%.*g / %.*g = %.*g", DBL_DIG, Ak, DBL_DIG, Bk, DBL_DIG, in fp2rat() 260 *q = Bk; in fp2rat()
|
/dports/math/openblas/OpenBLAS-0.3.18/kernel/zarch/ |
H A D | gemm_vec.c | 392 C_0_0 += Ak[0] * Bk[0]; in GEBP_block_8_4() 393 C_1_0 += Ak[1] * Bk[0]; in GEBP_block_8_4() 394 C_2_0 += Ak[2] * Bk[0]; in GEBP_block_8_4() 395 C_3_0 += Ak[3] * Bk[0]; in GEBP_block_8_4() 397 C_0_1 += Ak[0] * Bk[1]; in GEBP_block_8_4() 398 C_1_1 += Ak[1] * Bk[1]; in GEBP_block_8_4() 399 C_2_1 += Ak[2] * Bk[1]; in GEBP_block_8_4() 400 C_3_1 += Ak[3] * Bk[1]; in GEBP_block_8_4() 402 C_0_2 += Ak[0] * Bk[2]; in GEBP_block_8_4() 403 C_1_2 += Ak[1] * Bk[2]; in GEBP_block_8_4() [all …]
|
/dports/textproc/mdocml/mandoc-1.14.5/regress/mdoc/Bk/ |
H A D | badarg.out_lint | 1 mandoc: badarg.in:14:5: ERROR: skipping excess arguments: Bk ... -badarg 2 mandoc: badarg.in:23:5: ERROR: skipping excess arguments: Bk ... -lines 4 mandoc: badarg.in:32:5: ERROR: skipping excess arguments: Bk ... -line 5 mandoc: badarg.in:50:12: ERROR: skipping excess arguments: Bk ... -murks 6 mandoc: badarg.in:55:2: WARNING: empty block: Bk
|
H A D | inputlines.out_ascii | 4 BBkk--iinnppuuttlliinneess - effect of input lines breaks on word keeps 7 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 10 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 13 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 16 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 19 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 22 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 25 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 28 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… 31 …BBkk--iinnppuuttlliinneess _x _x _x _x _x _x _x _x _x _x _x _x _x _x _… [all …]
|
/dports/devel/tinygo/tinygo-0.14.1/lib/CMSIS/CMSIS/DSP_Lib/Examples/arm_convolution_example/ARM/ |
H A D | arm_convolution_example_f32.c | 119 float32_t Bk[MAX_BLOCKSIZE]; /* Input B */ variable 205 arm_fill_f32(0.0, Bk, MAX_BLOCKSIZE); in main() 209 arm_copy_f32(testInputB_f32, Bk, MAX_BLOCKSIZE/2); in main() 217 arm_cfft_radix4_f32(cfft_instance_ptr, Bk); in main() 220 arm_cmplx_mult_cmplx_f32(Ak, Bk, AxB, MAX_BLOCKSIZE/2); in main()
|
/dports/devel/tinygo/tinygo-0.14.1/lib/CMSIS/CMSIS/DSP_Lib/Examples/arm_convolution_example/GCC/ |
H A D | arm_convolution_example_f32.c | 119 float32_t Bk[MAX_BLOCKSIZE]; /* Input B */ variable 205 arm_fill_f32(0.0, Bk, MAX_BLOCKSIZE); in main() 209 arm_copy_f32(testInputB_f32, Bk, MAX_BLOCKSIZE/2); in main() 217 arm_cfft_radix4_f32(cfft_instance_ptr, Bk); in main() 220 arm_cmplx_mult_cmplx_f32(Ak, Bk, AxB, MAX_BLOCKSIZE/2); in main()
|
/dports/devel/arduino-core/Arduino-b439a77/hardware/arduino/sam/system/CMSIS/CMSIS/DSP_Lib/Examples/arm_convolution_example/ |
H A D | arm_convolution_example_f32.c | 103 float32_t Bk[MAX_BLOCKSIZE]; /* Input B */ variable 189 arm_fill_f32(0.0, Bk, MAX_BLOCKSIZE); in main() 193 arm_copy_f32(testInputB_f32, Bk, MAX_BLOCKSIZE/2); in main() 201 arm_cfft_radix4_f32(cfft_instance_ptr, Bk); in main() 204 arm_cmplx_mult_cmplx_f32(Ak, Bk, AxB, MAX_BLOCKSIZE/2); in main()
|
/dports/science/nwchem-data/nwchem-7.0.2-release/src/basis/libraries/ |
H A D | sarc-dkh | 4167 Bk S 4177 Bk S 4179 Bk S 4181 Bk S 4183 Bk S 4185 Bk S 4187 Bk S 4189 Bk S 4191 Bk S 4193 Bk S [all …]
|
H A D | sarc-zora | 4176 Bk S 4186 Bk S 4188 Bk S 4190 Bk S 4192 Bk S 4194 Bk S 4196 Bk S 4198 Bk S 4200 Bk S 4202 Bk S [all …]
|
/dports/science/nwchem/nwchem-7b21660b82ebd85ef659f6fba7e1e73433b0bd0a/src/basis/libraries/ |
H A D | sarc-zora | 4176 Bk S 4186 Bk S 4188 Bk S 4190 Bk S 4192 Bk S 4194 Bk S 4196 Bk S 4198 Bk S 4200 Bk S 4202 Bk S [all …]
|
H A D | sarc-dkh | 4167 Bk S 4177 Bk S 4179 Bk S 4181 Bk S 4183 Bk S 4185 Bk S 4187 Bk S 4189 Bk S 4191 Bk S 4193 Bk S [all …]
|
/dports/audio/praat/praat-6.2.03/external/glpk/ |
H A D | glplib03.c | 529 double xk, Akm1, Ak, Bkm1, Bk, ak, bk, fk, temp; in fp2rat() local 544 Bk = 1.0; in fp2rat() 559 temp = bk * Bk + ak * Bkm1; in fp2rat() 560 Bkm1 = Bk, Bk = temp; in fp2rat() 563 fk = Ak / Bk; in fp2rat() 565 print("%.*g / %.*g = %.*g", DBL_DIG, Ak, DBL_DIG, Bk, DBL_DIG, in fp2rat() 571 *q = Bk; in fp2rat()
|
/dports/science/py-pyscf/pyscf-2.0.1/pyscf/gto/basis/ |
H A D | sarc-dkh2.dat | 4140 Bk S 4150 Bk S 4152 Bk S 4154 Bk S 4156 Bk S 4158 Bk S 4160 Bk S 4162 Bk S 4164 Bk S 4166 Bk S [all …]
|
H A D | stuttgart_dz.dat | 1054 Bk S 1059 Bk S 1061 Bk S 1063 Bk S 1065 Bk S 2075 Bk ul 2077 Bk S 2081 Bk P 2085 Bk D 2089 Bk F [all …]
|
/dports/misc/mxnet/incubator-mxnet-1.9.0/3rdparty/mkldnn/src/gpu/ocl/gemm/ |
H A D | gen9_gemm.cpp | 354 for (int64_t Bk = 0; Bk < k; Bk += block_k) { in execute_standard() local 355 int64_t size_k = k - Bk; in execute_standard() 364 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 378 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 394 float eff_beta = (Bk == 0) ? beta : 1.0f; in execute_standard() 401 bool beta0 = (beta == 0) && (Bk == 0); in execute_standard() 555 for (int64_t Bk = 0; Bk < k; Bk += block_k) { in execute_superkernel() local 556 int64_t size_k = k - Bk; in execute_superkernel() 560 auto off_a = off_a0 + (!transa ? Bk * lda : Bk); in execute_superkernel() 561 auto off_b = off_b0 + (!transb ? Bk : Bk * ldb); in execute_superkernel() [all …]
|
H A D | xe_lp_gemm.cpp | 237 for (int64_t Bk = 0; Bk < k; Bk += size_k) { in execute_standard() local 238 size_k = k - Bk; in execute_standard() 245 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 250 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 251 apply_co = !co.is_null() && !(do_scale || (Bk > 0)); in execute_standard() 255 int eff_beta = ((Bk > 0) || (!do_scale && (beta == 1.0f))) in execute_standard()
|
H A D | gen9_gemm_x8x8s32.cpp | 239 for (int64_t Bk = 0; Bk < k; Bk += size_k) { in execute_standard() local 240 size_k = k - Bk; in execute_standard() 247 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 252 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 253 apply_co = !co.is_null() && !(do_scale || (Bk > 0)); in execute_standard() 257 int eff_beta = ((Bk > 0) || (!do_scale && (beta == 1.0f))) in execute_standard()
|
/dports/math/onednn/oneDNN-2.5.1/src/gpu/ocl/gemm/ |
H A D | gen9_gemm.cpp | 354 for (int64_t Bk = 0; Bk < k; Bk += block_k) { in execute_standard() local 355 int64_t size_k = k - Bk; in execute_standard() 364 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 378 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 394 float eff_beta = (Bk == 0) ? beta : 1.0f; in execute_standard() 401 bool beta0 = (beta == 0) && (Bk == 0); in execute_standard() 555 for (int64_t Bk = 0; Bk < k; Bk += block_k) { in execute_superkernel() local 556 int64_t size_k = k - Bk; in execute_superkernel() 560 auto off_a = off_a0 + (!transa ? Bk * lda : Bk); in execute_superkernel() 561 auto off_b = off_b0 + (!transb ? Bk : Bk * ldb); in execute_superkernel() [all …]
|
H A D | gen9_gemm_x8x8s32.cpp | 239 for (int64_t Bk = 0; Bk < k; Bk += size_k) { in execute_standard() local 240 size_k = k - Bk; in execute_standard() 247 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 252 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 253 apply_co = !co.is_null() && !(do_scale || (Bk > 0)); in execute_standard() 257 int eff_beta = ((Bk > 0) || (!do_scale && (beta == 1.0f))) in execute_standard()
|
H A D | xe_lp_gemm.cpp | 234 for (int64_t Bk = 0; Bk < k; Bk += size_k) { in execute_standard() local 235 size_k = k - Bk; in execute_standard() 242 + (!transa ? (Bm + Bk * lda) : (Bk + Bm * lda)); in execute_standard() 247 + (!transb ? (Bk + Bn * ldb) : (Bn + Bk * ldb)); in execute_standard() 248 apply_co = !co.is_null() && !(do_scale || (Bk > 0)); in execute_standard() 252 int eff_beta = ((Bk > 0) || (!do_scale && (beta == 1.0f))) in execute_standard()
|