Lines Matching refs:jrow
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()
495 PetscCall(VecPlaceArray(bmat->right, xx + bs * aj[jrow])); in MatMult_BlockMat_Symmetric()
496 PetscCall(MatMultAdd(aa[jrow], bmat->right, bmat->left, bmat->left)); in MatMult_BlockMat_Symmetric()
498 jrow++; in MatMult_BlockMat_Symmetric()
502 PetscCall(VecPlaceArray(bmat->right, xx + bs * aj[jrow])); /* upper triangular part */ in MatMult_BlockMat_Symmetric()
503 PetscCall(MatMultAdd(aa[jrow], bmat->right, bmat->left, bmat->left)); in MatMult_BlockMat_Symmetric()
506 PetscCall(VecPlaceArray(bmat->right, yy + bs * aj[jrow])); /* lower triangular part */ in MatMult_BlockMat_Symmetric()
507 PetscCall(MatMultTransposeAdd(aa[jrow], bmat->middle, bmat->right, bmat->right)); in MatMult_BlockMat_Symmetric()
509 jrow++; in MatMult_BlockMat_Symmetric()