Home
last modified time | relevance | path

Searched refs:nrhs (Results 1 – 25 of 61) sorted by relevance

123

/petsc/src/mat/tests/
H A Dex267.c5 PetscErrorCode TestMatrix(const char *test, Mat A, PetscInt nrhs, PetscBool inplace, PetscBool chol) in TestMatrix() argument
82 PetscCall(MatSetSizes(RHS, PETSC_DECIDE, PETSC_DECIDE, n, nrhs)); in TestMatrix()
85 for (i = 0; i < nrhs; i++) PetscCall(MatSetValue(RHS, i, i, 1.0, INSERT_VALUES)); in TestMatrix()
121 PetscInt i, n = 8, nrhs = 2; in main() local
134 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
136 PetscCheck(nrhs <= n, PETSC_COMM_WORLD, PETSC_ERR_ARG_SIZ, "Must have nrhs <= n"); in main()
156 PetscCall(TestMatrix("LU T", At, nrhs, inplace, PETSC_FALSE)); in main()
157 PetscCall(TestMatrix("LU HT", Aht, nrhs, inplace, PETSC_FALSE)); in main()
159 PetscCall(TestMatrix("Chol T", At, nrhs, inplace, PETSC_TRUE)); in main()
160 PetscCall(TestMatrix("Chol HT", Aht, nrhs, inplace, PETSC_TRUE)); in main()
[all …]
H A Dex214.c11 PetscInt m, n, nrhs, M, N, i, Istart, Iend, Ii, j, J, test; in main() local
76 nrhs = N; in main()
77 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
79 PetscCall(MatSetSizes(C, m, PETSC_DECIDE, PETSC_DECIDE, nrhs)); in main()
90 …s: nrhs %" PetscInt_FMT ", size mat %" PetscInt_FMT " x %" PetscInt_FMT "\n", solver, nrhs, M, N)); in main()
131 for (i = 0; i < nrhs; i++) { in main()
140 …intf(PETSC_COMM_SELF, " \n(2) first %" PetscInt_FMT " columns of inv(A) with dense RHS:\n", nrhs)); in main()
159 PetscCall(MatSetSizes(spRHST, nrhs, M, PETSC_DECIDE, PETSC_DECIDE)); in main()
168 for (i = 0; i < nrhs; i++) PetscCall(MatSetValues(spRHST, 1, &i, 1, &i, &v, INSERT_VALUES)); in main()
176 …ntf(PETSC_COMM_SELF, " \n(3) first %" PetscInt_FMT " columns of inv(A) with sparse RHS:\n", nrhs)); in main()
[all …]
H A Dex129.c29 PetscInt dof = 1, M = 8, m, n, nrhs; in main() local
62 nrhs = 2; in main()
63 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
64 PetscCall(ComputeRHSMatrix(m, nrhs, &RHS)); in main()
150 PetscErrorCode ComputeRHSMatrix(PetscInt m, PetscInt nrhs, Mat *C) in ComputeRHSMatrix() argument
159 PetscCall(MatSetSizes(RHS, m, PETSC_DECIDE, PETSC_DECIDE, nrhs)); in ComputeRHSMatrix()
170 if (nrhs > 1) { in ComputeRHSMatrix()
171 for (k = 1; k < nrhs; k++) { in ComputeRHSMatrix()
H A Dex1.c7 static PetscErrorCode createMatsAndVecs(PetscInt m, PetscInt n, PetscInt nrhs, PetscBool full, Mat … in createMatsAndVecs() argument
19 PetscCall(MatSetSizes(RHS, PETSC_DECIDE, PETSC_DECIDE, m, nrhs)); in createMatsAndVecs()
33 PetscCall(MatSetSizes(SOLU, PETSC_DECIDE, PETSC_DECIDE, n, nrhs)); in createMatsAndVecs()
90 PetscInt m = 15, n = 10, i, j, nrhs = 2; in main() local
105 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
111 PetscCall(createMatsAndVecs(n, n, nrhs, full, &mat, &RHS, &SOLU, &x, &y, &b)); in main()
203 PetscCall(createMatsAndVecs(m, n, nrhs, full, &mat, &RHS, &SOLU, &x, &y, &b)); in main()
H A Dex215.c10 PetscInt m, n, nsolve, nrhs; in main() local
74 nrhs = 2; in main()
75 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
77 PetscCall(MatSetSizes(C, m, PETSC_DECIDE, PETSC_DECIDE, nrhs)); in main()
H A Dex192.c13 PetscInt isolver = 0, size_schur, m, n, nfact, nsolve, nrhs; in main() local
94 nrhs = 2; in main()
95 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
97 PetscCall(MatSetSizes(C, m, PETSC_DECIDE, PETSC_DECIDE, nrhs)); in main()
142 … %d, size schur %" PetscInt_FMT ", size mat %" PetscInt_FMT "\n", solver, nrhs, symm, herm, size_s… in main()
/petsc/src/sys/classes/viewer/impls/socket/mex-scripts/
H A Dsclose.c52 PETSC_EXTERN void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) in mexFunction() argument
60 if (!nrhs) PETSC_MEX_ERROR("Needs one argument, the port"); in mexFunction()
H A Dsread.c17 PETSC_EXTERN void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) in mexFunction() argument
24 if (nrhs != 3) PETSC_MEX_ERROR("Receive requires three input arguments."); in mexFunction()
H A Dswrite.c17 PETSC_EXTERN void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) in mexFunction() argument
23 if (nrhs != 3) PETSC_MEX_ERROR("Receive requires three input arguments."); in mexFunction()
H A Dsopen.c158 PETSC_EXTERN void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) argument
166 if (!nrhs) {
/petsc/src/ksp/ksp/tests/
H A Dex26.c47 PetscInt i, its, Nx = PETSC_DECIDE, Ny = PETSC_DECIDE, nlocal, nrhs = 1; in main() local
61 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nrhs", &nrhs, NULL)); in main()
97 …PetscCall(MatCreateDense(PETSC_COMM_WORLD, nlocal, PETSC_DECIDE, fine_ctx.mx * fine_ctx.my, nrhs, … in main()
107 for (i = 0; i < nlocal * nrhs; i++) b[i] = 1.0; in main()
116 if ((flg || nrhs == 1) && !Brand) { in main()
122 for (n = 0; n < nrhs; n++) { in main()
/petsc/src/mat/impls/aij/seq/mkl_pardiso/
H A Dmkl_pardiso.c72 INT_TYPE phase, maxfct, mnum, mtype, n, nrhs, msglvl, err; member
244 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, mpardiso->schur_size, mpardiso->nrhs, B, &Bmat)); in MatMKLPardisoSolveSchur_Private()
245 PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, mpardiso->schur_size, mpardiso->nrhs, X, &Xmat)); in MatMKLPardisoSolveSchur_Private()
334 …so->phase, &mat_mkl_pardiso->n, NULL, NULL, NULL, NULL, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->i… in MatDestroy_MKL_PARDISO()
358 for (i = 0; i < mpardiso->nrhs; i++) { in MatMKLPardisoScatterSchur_Private()
366 for (i = 0; i < mpardiso->nrhs; i++) { in MatMKLPardisoScatterSchur_Private()
383 mat_mkl_pardiso->nrhs = 1; in MatSolve_MKL_PARDISO()
399 …&mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (void *)xarray, (void *)… in MatSolve_MKL_PARDISO()
404 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatSolve_MKL_PARDISO()
431 …mat_mkl_pardiso->perm, &mat_mkl_pardiso->nrhs, mat_mkl_pardiso->iparm, &mat_mkl_pardiso->msglvl, (… in MatSolve_MKL_PARDISO()
[all …]
/petsc/src/mat/impls/aij/mpi/mkl_cpardiso/
H A Dmkl_cpardiso.c91 INT_TYPE phase, maxfct, mnum, mtype, n, nrhs, msglvl, err; member
359 …so->phase, &mat_mkl_cpardiso->n, NULL, NULL, NULL, mat_mkl_cpardiso->perm, &mat_mkl_cpardiso->nrhs, in MatDestroy_MKL_CPARDISO()
383 mat_mkl_cpardiso->nrhs = 1; in MatSolve_MKL_CPARDISO()
390 …mat_mkl_cpardiso->perm, &mat_mkl_cpardiso->nrhs, mat_mkl_cpardiso->iparm, &mat_mkl_cpardiso->msglv… in MatSolve_MKL_CPARDISO()
406 mat_mkl_cpardiso->nrhs = 1; in MatForwardSolve_MKL_CPARDISO()
413 …mat_mkl_cpardiso->perm, &mat_mkl_cpardiso->nrhs, mat_mkl_cpardiso->iparm, &mat_mkl_cpardiso->msglv… in MatForwardSolve_MKL_CPARDISO()
429 mat_mkl_cpardiso->nrhs = 1; in MatBackwardSolve_MKL_CPARDISO()
436 …mat_mkl_cpardiso->perm, &mat_mkl_cpardiso->nrhs, mat_mkl_cpardiso->iparm, &mat_mkl_cpardiso->msglv… in MatBackwardSolve_MKL_CPARDISO()
463 PetscCall(MatGetSize(B, NULL, (PetscInt *)&mat_mkl_cpardiso->nrhs)); in MatMatSolve_MKL_CPARDISO()
465 if (mat_mkl_cpardiso->nrhs > 0) { in MatMatSolve_MKL_CPARDISO()
[all …]
/petsc/src/mat/tests/output/
H A Dex192_mumps_3.out1 Solving with mumps: nrhs 2, sym 1, herm 1, size schur 5, size mat 12
H A Dex192_mkl_pardiso.out1 Solving with mkl_pardiso: nrhs 2, sym 0, herm 0, size schur 5, size mat 12
H A Dex192_mumps.out1 Solving with mumps: nrhs 2, sym 0, herm 0, size schur 5, size mat 12
H A Dex192_mkl_pardiso_1.out1 Solving with mkl_pardiso: nrhs 2, sym 1, herm 0, size schur 5, size mat 12
H A Dex192_mumps_2.out1 Solving with mumps: nrhs 2, sym 1, herm 0, size schur 5, size mat 12
H A Dex192_mkl_pardiso_3.out1 Solving with mkl_pardiso: nrhs 2, sym 1, herm 1, size schur 5, size mat 12
H A Dex125_cusparse_2_cholesky-0.out1 ex125: nrhs 2
H A Dex125_superlu_dist_complex.out1 ex125: nrhs 2
H A Dex125_cusparse_2_cholesky-1.out1 ex125: nrhs 2
H A Dex125_cusparse_cholesky-0.out1 ex125: nrhs 2
H A Dex125_superlu_dist_nonsymmetric.out1 ex125: nrhs 2
H A Dex125_cusparse_cholesky-1.out1 ex125: nrhs 2

123