| /dports/science/hypre/hypre-2.23.0/src/seq_mv/ |
| H A D | csr_sptrans_device.c | 55 m, n, nnzA, in hypreDevice_CSRSpTransCusparse()
67 m, n, nnzA, in hypreDevice_CSRSpTransCusparse()
127 m, n, nnzA, in hypreDevice_CSRSpTransRocsparse()
167 d_jc = hypre_TAlloc(HYPRE_Int, nnzA, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpTrans()
179 d_it = d_pm + nnzA; in hypreDevice_CSRSpTrans()
180 hypreDevice_CsrRowPtrsToIndices_v2(m, nnzA, d_ia, d_it); in hypreDevice_CSRSpTrans()
184 d_jt = d_it + nnzA; in hypreDevice_CSRSpTrans()
188 HYPRE_THRUST_CALL(sequence, d_pm, d_pm + nnzA); in hypreDevice_CSRSpTrans()
190 HYPRE_THRUST_CALL(gather, d_pm, d_pm + nnzA, d_it, d_jc); in hypreDevice_CSRSpTrans()
197 d_ic = hypreDevice_CsrRowIndicesToPtrs(n, nnzA, d_jt); in hypreDevice_CSRSpTrans()
[all …]
|
| H A D | csr_spadd_device.c | 28 HYPRE_Int nnzA, in hypreDevice_CSRSpAdd() argument
47 if (nnzA == 0 && nnzB == 0) in hypreDevice_CSRSpAdd()
62 HYPRE_Int nnzT = nnzA + nnzB, nnzC; in hypreDevice_CSRSpAdd()
85 HYPRE_THRUST_CALL(gather, d_ja, d_ja + nnzA, d_ja_map, d_jt); in hypreDevice_CSRSpAdd()
89 hypre_TMemcpy(d_jt, d_ja, HYPRE_Int, nnzA, HYPRE_MEMORY_DEVICE, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpAdd()
93 HYPRE_THRUST_CALL(gather, d_jb, d_jb + nnzB, d_jb_map, d_jt + nnzA); in hypreDevice_CSRSpAdd()
107 HYPRE_THRUST_CALL( transform, d_aa, d_aa + nnzA, d_at, alpha * _1 ); in hypreDevice_CSRSpAdd()
116 HYPRE_THRUST_CALL( transform, d_ab, d_ab + nnzB, d_at + nnzA, beta * _1 ); in hypreDevice_CSRSpAdd()
120 hypreDevice_CsrRowPtrsToIndices_v2(ma, nnzA, d_ia, d_it); in hypreDevice_CSRSpAdd()
123 hypreDevice_CsrRowPtrsToIndicesWithRowNum(mb, nnzB, d_ib, d_num_b, d_it + nnzA); in hypreDevice_CSRSpAdd()
[all …]
|
| H A D | csr_spgemm_device_rocsparse.c | 19 HYPRE_Int nnzA, in hypreDevice_CSRSpGemmRocsparse() argument
39 d_a_sorted = hypre_TAlloc(HYPRE_Complex, nnzA, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmRocsparse()
41 d_ja_sorted = hypre_TAlloc(HYPRE_Int, nnzA, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmRocsparse()
55 hypre_TMemcpy(d_ja_sorted, d_ja, HYPRE_Int, nnzA, HYPRE_MEMORY_DEVICE, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmRocsparse()
56 hypre_TMemcpy(d_a_sorted, d_a, HYPRE_Complex, nnzA, HYPRE_MEMORY_DEVICE, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmRocsparse()
91 descrA, nnzA, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmRocsparse()
102 descrA, nnzA, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmRocsparse()
116 descrA, nnzA, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmRocsparse()
141 descrA, nnzA, d_a_sorted, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmRocsparse()
153 … descrA, nnzA, (float *) d_a_sorted, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmRocsparse()
|
| H A D | csr_spgemm_device_cusparse.c | 19 HYPRE_Int nnzA, in hypreDevice_CSRSpGemmCusparse() argument
36 nnzA, d_ia, d_ja, d_a, in hypreDevice_CSRSpGemmCusparse()
41 descr_A, nnzA, d_ia, d_ja, d_a, in hypreDevice_CSRSpGemmCusparse()
74 HYPRE_Int nnzA, in hypreDevice_CSRSpGemmCusparseGenericAPI() argument
90 cusparseSpMatDescr_t matA = hypre_CSRMatrixToCusparseSpMat_core(m, k, 0, nnzA, d_ia, d_ja, d_a); in hypreDevice_CSRSpGemmCusparseGenericAPI()
182 HYPRE_Int nnzA, in hypreDevice_CSRSpGemmCusparseOldAPI() argument
202 d_a_sorted = hypre_TAlloc(HYPRE_Complex, nnzA, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmCusparseOldAPI()
204 d_ja_sorted = hypre_TAlloc(HYPRE_Int, nnzA, HYPRE_MEMORY_DEVICE); in hypreDevice_CSRSpGemmCusparseOldAPI()
223 hypre_SortCSRCusparse(m, k, nnzA, descr_A, d_ia, d_ja_sorted, d_a_sorted); in hypreDevice_CSRSpGemmCusparseOldAPI()
234 descr_A, nnzA, d_ia, d_ja_sorted, in hypreDevice_CSRSpGemmCusparseOldAPI()
[all …]
|
| H A D | csr_matop_device.c | 1432 HYPRE_Int nnzA = hypre_CSRMatrixNumNonzeros(A); in hypre_CSRMatrixIntersectPattern() local
1443 HYPRE_THRUST_CALL( sequence, idx, idx + nnzA + nnzB ); in hypre_CSRMatrixIntersectPattern()
1456 nnzA + nnzB, nnzA, Cii, Cjj, idx, markA, diag_opt ); in hypre_CSRMatrixIntersectPattern()
1547 HYPRE_Int nnzA, in hypre_SortCSRCusparse() argument
1602 HYPRE_Int nnzA = hypre_CSRMatrixNumNonzeros(A); in hypre_CSRMatrixTriLowerUpperSolveCusparse() local
1626 if (nnzA <= 0) in hypre_CSRMatrixTriLowerUpperSolveCusparse()
1774 HYPRE_Int nnzA = hypre_CSRMatrixNumNonzeros(A); in hypre_CSRMatrixTriLowerUpperSolveRocsparse() local
1798 if (nnzA <= 0) in hypre_CSRMatrixTriLowerUpperSolveRocsparse()
1944 HYPRE_Int nnzA, in hypre_SortCSRRocsparse() argument
1964 d_a_tmp = hypre_TAlloc(HYPRE_Complex, nnzA, HYPRE_MEMORY_DEVICE); in hypre_SortCSRRocsparse()
[all …]
|
| H A D | seq_mv.hpp | 21 …RE_Int m, HYPRE_Int k, HYPRE_Int n, cusparseMatDescr_t descr_A, HYPRE_Int nnzA, HYPRE_Int *d_ia, H…
23 …RE_Int m, HYPRE_Int k, HYPRE_Int n, cusparseMatDescr_t descr_A, HYPRE_Int nnzA, HYPRE_Int *d_ia, H…
25 void hypre_SortCSRCusparse( HYPRE_Int n, HYPRE_Int m, HYPRE_Int nnzA, cusparseMatDescr_t descrA, co…
29 …RE_Int m, HYPRE_Int k, HYPRE_Int n, rocsparse_mat_descr descrA, HYPRE_Int nnzA, HYPRE_Int *d_ia, H…
31 void hypre_SortCSRRocsparse( HYPRE_Int n, HYPRE_Int m, HYPRE_Int nnzA, rocsparse_mat_descr descrA, …
|
| /dports/science/ergo/ergo-3.8/source/integrals/ |
| H A D | simple_sparse_mat.cc | 60 int spmat_sort_elements(i_j_val_struct* A, int nnzA) { in spmat_sort_elements() argument
61 qsort(A, nnzA, sizeof(i_j_val_struct), compare_i_j_val_structs); in spmat_sort_elements()
63 for(int k = 0; k < nnzA-1; k++) { in spmat_sort_elements()
67 for(int k = 0; k < nnzA; k++) { in spmat_sort_elements()
70 while(kk < nnzA) { in spmat_sort_elements()
81 int spmat_multiply_matrices(const i_j_val_struct* A, int nnzA, const i_j_val_struct* B, int nnzB, i… in spmat_multiply_matrices() argument
86 for(int idxA = 0; idxA < nnzA; idxA++) in spmat_multiply_matrices()
|
| /dports/math/arpack++/arpackpp-2.3.0/examples/areig/ |
| H A D | areig.h | 470 int AREig(arcomplex<ARFLOAT> EigVal[], int n, int nnzA,
715 int AREig(double EigValR[], ARFLOAT EigValI[], int n, int nnzA,
740 int AREig(float EigValR[], ARFLOAT EigValI[], int n, int nnzA,
766 int nnzA, ARFLOAT A[], int irowA[], int pcolA[],
791 int AREig(double EigValR[], ARFLOAT EigValI[], int n, int nnzA,
818 int AREig(float EigValR[], ARFLOAT EigValI[], int n, int nnzA,
846 int nnzA, ARFLOAT A[], int irowA[], int pcolA[], int nnzB,
1029 ARluSymMatrix<ARFLOAT> matrixA(n, nnzA, A, irowA, pcolA, uplo);
1054 ARluSymMatrix<ARFLOAT> matrixA(n, nnzA, A, irowA, pcolA, uplo);
1079 ARluSymMatrix<ARFLOAT> matrixA(n, nnzA, A, irowA, pcolA, uplo);
[all …]
|
| /dports/math/osqp/osqp-0.6.2/tests/update_matrices/ |
| H A D | test_update_matrices.h | 75 c_int i, nnzP, nnzA; in test_update() local
193 nnzA = data->test_solve_A->p[data->test_solve_A->n]; in test_update()
194 Ax_new_idx = (c_int*) c_malloc(nnzA * sizeof(c_int)); in test_update()
197 for (i = 0; i < nnzA; i++) { in test_update()
226 osqp_update_A(work, data->test_solve_A_new->x, OSQP_NULL, nnzA); in test_update()
252 data->test_solve_A_new->x, Ax_new_idx, nnzA); in test_update()
281 data->test_solve_A_new->x, OSQP_NULL, nnzA); in test_update()
317 c_int i, nnzP, nnzA, exitflag; in test_update_pardiso() local
406 nnzA = data->test_solve_A->p[data->test_solve_A->n]; in test_update_pardiso()
410 for (i = 0; i < nnzA; i++) { in test_update_pardiso()
[all …]
|
| /dports/math/arpack++/arpackpp-2.3.0/examples/areig/nonsym/ |
| H A D | ansymgsc.cc | 64 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
79 NonSymMatrixE(n, nnzA, valA, irowA, pcolA); in main()
85 nconv = AREig(EigValR, EigValI, EigVec, n, nnzA, valA, in main()
91 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | ansymgre.cc | 63 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
80 StiffnessMatrix(n, rho, nnzA, valA, irowA, pcolA); in main()
86 nconv = AREig(EigValR, EigValI, EigVec, n, nnzA, valA, in main()
91 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | ansymgsh.cc | 66 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
83 StiffnessMatrix(n, rho, nnzA, valA, irowA, pcolA); in main()
89 nconv = AREig(EigValR, EigValI, EigVec, n, nnzA, valA, in main()
94 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| /dports/math/arpack++/arpackpp-2.3.0/examples/areig/sym/ |
| H A D | asymgcay.cc | 66 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
83 SymmetricMatrixC(n, nnzA, valA, irowA, pcolA, uplo); in main()
89 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, pcolA, in main()
94 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | asymgbkl.cc | 66 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
83 SymmetricMatrixC(n, nnzA, valA, irowA, pcolA, uplo); in main()
89 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, pcolA, in main()
94 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | asymgreg.cc | 66 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
83 SymmetricMatrixC(n, nnzA, valA, irowA, pcolA, uplo); in main()
89 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, in main()
94 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | asymgshf.cc | 66 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
83 SymmetricMatrixC(n, nnzA, valA, irowA, pcolA, uplo); in main()
89 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, pcolA, in main()
94 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| /dports/math/arpack++/arpackpp-2.3.0/examples/areig/complex/ |
| H A D | acompgre.cc | 65 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
81 CompMatrixE(n, rho, nnzA, valA, irowA, pcolA); in main()
87 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, in main()
92 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| H A D | acompgsh.cc | 65 int nnzA, nnzB; // Number of nonzero elements in A and B. in main() local
81 CompMatrixE(n, rho, nnzA, valA, irowA, pcolA); in main()
87 nconv = AREig(EigVal, EigVec, n, nnzA, valA, irowA, pcolA, nnzB, in main()
92 Solution(nconv, n, nnzA, valA, irowA, pcolA, nnzB, in main()
|
| /dports/math/librsb/librsb-1.2.0.10/ |
| H A D | rsb_swt.c | 43 rsb_nnz_idx_t nnzA = mtxAp->nnz; in rsb__do_switch_leaf() local
51 if( nnzA >= nrA+1 && nnzA >= ncA+1 ) in rsb__do_switch_leaf()
78 if( nnzA < nrA+1 ) in rsb__do_switch_leaf()
91 rsb__do_switch_array_to_halfword_coo(JA,nnzA,0); in rsb__do_switch_leaf()
104 rsb__do_switch_array_to_halfword_coo(JA,nnzA,0); in rsb__do_switch_leaf()
142 if(nnzA) in rsb__do_switch_leaf()
145 RSB_ASSERT( mtxAp->bpntr[nrA] == nnzA ); in rsb__do_switch_leaf()
156 rsb__do_switch_array_to_halfword_coo(IA,nnzA,0); in rsb__do_switch_leaf()
168 rsb__do_switch_array_to_halfword_coo(JA,nnzA,coff); in rsb__do_switch_leaf()
189 rsb__util_coo_array_add(IA,nnzA,roff); in rsb__do_switch_leaf()
[all …]
|
| /dports/math/octave-forge-sparsersb/sparsersb-1.0.9/src/ |
| H A D | sparsersb.cc | 358 rsb_nnz_idx_t nnzA = 0; in rsboi_mtx_get_coo() local
372 rsboi_oi2ri(IA,nnzA); in rsboi_mtx_alloc_from_csc_const()
387 rsboi_oi2ri(IA,nnzA); in rsboi_mtx_alloc_from_coo_const()
388 rsboi_oi2ri(JA,nnzA); in rsboi_mtx_alloc_from_coo_const()
758 VA.resize1(nnzA + nze); in sparse_matrix_value()
759 IA.resize1(nnzA + nze); in sparse_matrix_value()
760 JA.resize1(nnzA + nze); in sparse_matrix_value()
876 VAp[2*(nnzA-1-nzi)+0]=VAp[(nnzA-1-nzi)+0], in sparse_complex_matrix_value()
1457 nnzA = 0; in print()
1470 ", nnz = "<<nnzA in print()
[all …]
|
| /dports/math/suitesparse-config/SuiteSparse-5.10.1/GraphBLAS/CUDA/templates/ |
| H A D | GB_jit_AxB_dot3_phase3_mp.cu | 110 int64_t nnzA = 0; in AxB_dot3_phase3_mp() local
131 nnzA = xend - xstart; in AxB_dot3_phase3_mp()
144 int64_t nxy = nnzA + nnzB; in AxB_dot3_phase3_mp()
152 int x_max = GB_IMIN( diag, nnzA); in AxB_dot3_phase3_mp()
166 if (( diag > 0) &&(diag < (nnzA+nnzB)) && (Ai[xcoord+xstart] == Bi[ycoord+ystart]) ) { in AxB_dot3_phase3_mp()
177 x_max = GB_IMIN( diag_end, nnzA); in AxB_dot3_phase3_mp()
192 if ( (diag_end < (nnzA +nnzB)) && (Ai[xcoord +xstart] == Bi[ycoord + ystart]) ) { in AxB_dot3_phase3_mp()
|
| /dports/math/suitesparse-klu/SuiteSparse-5.10.1/GraphBLAS/CUDA/templates/ |
| H A D | GB_jit_AxB_dot3_phase3_mp.cu | 110 int64_t nnzA = 0; in AxB_dot3_phase3_mp() local
131 nnzA = xend - xstart; in AxB_dot3_phase3_mp()
144 int64_t nxy = nnzA + nnzB; in AxB_dot3_phase3_mp()
152 int x_max = GB_IMIN( diag, nnzA); in AxB_dot3_phase3_mp()
166 if (( diag > 0) &&(diag < (nnzA+nnzB)) && (Ai[xcoord+xstart] == Bi[ycoord+ystart]) ) { in AxB_dot3_phase3_mp()
177 x_max = GB_IMIN( diag_end, nnzA); in AxB_dot3_phase3_mp()
192 if ( (diag_end < (nnzA +nnzB)) && (Ai[xcoord +xstart] == Bi[ycoord + ystart]) ) { in AxB_dot3_phase3_mp()
|
| /dports/math/suitesparse-btf/SuiteSparse-5.10.1/GraphBLAS/CUDA/templates/ |
| H A D | GB_jit_AxB_dot3_phase3_mp.cu | 110 int64_t nnzA = 0; in AxB_dot3_phase3_mp() local
131 nnzA = xend - xstart; in AxB_dot3_phase3_mp()
144 int64_t nxy = nnzA + nnzB; in AxB_dot3_phase3_mp()
152 int x_max = GB_IMIN( diag, nnzA); in AxB_dot3_phase3_mp()
166 if (( diag > 0) &&(diag < (nnzA+nnzB)) && (Ai[xcoord+xstart] == Bi[ycoord+ystart]) ) { in AxB_dot3_phase3_mp()
177 x_max = GB_IMIN( diag_end, nnzA); in AxB_dot3_phase3_mp()
192 if ( (diag_end < (nnzA +nnzB)) && (Ai[xcoord +xstart] == Bi[ycoord + ystart]) ) { in AxB_dot3_phase3_mp()
|
| /dports/math/suitesparse-amd/SuiteSparse-5.10.1/GraphBLAS/CUDA/templates/ |
| H A D | GB_jit_AxB_dot3_phase3_mp.cu | 110 int64_t nnzA = 0; in AxB_dot3_phase3_mp() local
131 nnzA = xend - xstart; in AxB_dot3_phase3_mp()
144 int64_t nxy = nnzA + nnzB; in AxB_dot3_phase3_mp()
152 int x_max = GB_IMIN( diag, nnzA); in AxB_dot3_phase3_mp()
166 if (( diag > 0) &&(diag < (nnzA+nnzB)) && (Ai[xcoord+xstart] == Bi[ycoord+ystart]) ) { in AxB_dot3_phase3_mp()
177 x_max = GB_IMIN( diag_end, nnzA); in AxB_dot3_phase3_mp()
192 if ( (diag_end < (nnzA +nnzB)) && (Ai[xcoord +xstart] == Bi[ycoord + ystart]) ) { in AxB_dot3_phase3_mp()
|
| /dports/math/suitesparse-ldl/SuiteSparse-5.10.1/GraphBLAS/CUDA/templates/ |
| H A D | GB_jit_AxB_dot3_phase3_mp.cu | 110 int64_t nnzA = 0; in AxB_dot3_phase3_mp() local
131 nnzA = xend - xstart; in AxB_dot3_phase3_mp()
144 int64_t nxy = nnzA + nnzB; in AxB_dot3_phase3_mp()
152 int x_max = GB_IMIN( diag, nnzA); in AxB_dot3_phase3_mp()
166 if (( diag > 0) &&(diag < (nnzA+nnzB)) && (Ai[xcoord+xstart] == Bi[ycoord+ystart]) ) { in AxB_dot3_phase3_mp()
177 x_max = GB_IMIN( diag_end, nnzA); in AxB_dot3_phase3_mp()
192 if ( (diag_end < (nnzA +nnzB)) && (Ai[xcoord +xstart] == Bi[ycoord + ystart]) ) { in AxB_dot3_phase3_mp()
|