Searched refs:jrow (Results 1 – 3 of 3) sorted by relevance
| /petsc/src/mat/impls/blockmat/seq/ |
| H A D | blockmat.c | 439 PetscInt *aj, i, *ii, jrow, m = A->rmap->n / A->rmap->bs, bs = A->rmap->bs, n, j; in MatMult_BlockMat() local
454 jrow = ii[i]; in MatMult_BlockMat()
456 n = ii[i + 1] - jrow; in MatMult_BlockMat()
458 PetscCall(VecPlaceArray(bmat->right, xx + bs * aj[jrow])); in MatMult_BlockMat()
459 PetscCall(MatMultAdd(aa[jrow], bmat->right, bmat->left, bmat->left)); in MatMult_BlockMat()
461 jrow++; in MatMult_BlockMat()
474 PetscInt *aj, i, *ii, jrow, m = A->rmap->n / A->rmap->bs, bs = A->rmap->bs, n, j; in MatMult_BlockMat_Symmetric() local
489 jrow = ii[i]; in MatMult_BlockMat_Symmetric()
490 n = ii[i + 1] - jrow; in MatMult_BlockMat_Symmetric()
494 if (aj[jrow] == i) { in MatMult_BlockMat_Symmetric()
[all …]
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 220 PetscInt jrow = ii[i]; in MatMult_MatMultAdd_SeqMAIJ_Template() local
221 const PetscInt n = ii[i + 1] - jrow; in MatMult_MatMultAdd_SeqMAIJ_Template()
226 for (PetscInt j = 0; j < n; ++j, ++jrow) { in MatMult_MatMultAdd_SeqMAIJ_Template()
227 const PetscScalar v_jrow = v[jrow]; in MatMult_MatMultAdd_SeqMAIJ_Template()
228 const PetscInt N_idx_jrow = N * idx[jrow]; in MatMult_MatMultAdd_SeqMAIJ_Template()
351 PetscInt n, i, jrow, j, dof = b->dof, k; in MatMult_SeqMAIJ_N() local
362 jrow = ii[i]; in MatMult_SeqMAIJ_N()
363 n = ii[i + 1] - jrow; in MatMult_SeqMAIJ_N()
366 for (k = 0; k < dof; k++) sums[k] += v[jrow] * x[dof * idx[jrow] + k]; in MatMult_SeqMAIJ_N()
367 jrow++; in MatMult_SeqMAIJ_N()
[all …]
|
| /petsc/src/mat/impls/kaij/ |
| H A D | kaij.c | 660 PetscInt n, i, jrow, j, l, p = b->p, q = b->q, k; in MatMultAdd_SeqKAIJ() local
678 jrow = ii[i]; in MatMultAdd_SeqKAIJ()
679 n = ii[i + 1] - jrow; in MatMultAdd_SeqKAIJ()
682 for (k = 0; k < p; k++) sums[k] += v[jrow + j] * x[q * idx[jrow + j] + k]; in MatMultAdd_SeqKAIJ()
688 jrow = ii[i]; in MatMultAdd_SeqKAIJ()
689 n = ii[i + 1] - jrow; in MatMultAdd_SeqKAIJ()
693 for (l = 0; l < q; l++) sums[k] += v[jrow + j] * t[k + l * p] * x[q * idx[jrow + j] + l]; in MatMultAdd_SeqKAIJ()
|