Lines Matching refs:bdiag
1204 const PetscInt n = A->rmap->n, *ai = a->i, *aj = a->j, *bi = b->i, *bj = b->j, *bdiag = b->diag; in MatLUFactorNumeric_SeqAIJ_Inode() local
1301 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1302 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1321 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqAIJ_Inode()
1324 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1325 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1326 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1345 pv = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1346 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1347 nz = bdiag[i] - bdiag[i + 1] - 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1360 pv = b->a + bdiag[i]; in MatLUFactorNumeric_SeqAIJ_Inode()
1376 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1377 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1406 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqAIJ_Inode()
1412 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1413 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1414 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1437 pc1 = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1438 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1439 nz = bdiag[i] - bdiag[i + 1] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1450 pc1 = b->a + bdiag[i]; /* Mark diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1458 pj = b->j + bdiag[i + 1] + 1; /* beginning of U(i,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1459 nz = bdiag[i] - bdiag[i + 1] - 1; /* num of entries in U(i,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1479 pc2 = b->a + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1480 pj = b->j + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1481 nz = bdiag[i + 1] - bdiag[i + 2] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1492 pc2 = b->a + bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1509 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1510 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1544 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqAIJ_Inode()
1552 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1553 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1554 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1578 pc1 = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1579 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1580 nz = bdiag[i] - bdiag[i + 1] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1591 pc1 = b->a + bdiag[i]; /* Mark diag[i] */ in MatLUFactorNumeric_SeqAIJ_Inode()
1602 pj = b->j + bdiag[i + 1] + 1; /* beginning of U(i,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1603 nz = bdiag[i] - bdiag[i + 1] - 1; /* num of entries in U(i,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1624 pc2 = b->a + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1625 pj = b->j + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1626 nz = bdiag[i + 1] - bdiag[i + 2] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1637 pc2 = b->a + bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1645 pj = b->j + bdiag[i + 2] + 1; /* beginning of U(i+1,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1646 nz = bdiag[i + 1] - bdiag[i + 2] - 1; /* num of entries in U(i+1,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1666 pc3 = b->a + bdiag[i + 3] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1667 pj = b->j + bdiag[i + 3] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1668 nz = bdiag[i + 2] - bdiag[i + 3] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1679 pc3 = b->a + bdiag[i + 2]; in MatLUFactorNumeric_SeqAIJ_Inode()
1696 nz = bdiag[i] - bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1697 bjtmp = bj + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1736 pv = b->a + bdiag[row]; in MatLUFactorNumeric_SeqAIJ_Inode()
1746 pj = b->j + bdiag[row + 1] + 1; /* beginning of U(row,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1747 pv = b->a + bdiag[row + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1748 nz = bdiag[row] - bdiag[row + 1] - 1; /* num of entries in U(row,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1773 pc1 = b->a + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1774 pj = b->j + bdiag[i + 1] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1775 nz = bdiag[i] - bdiag[i + 1] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1786 pc1 = b->a + bdiag[i]; /* Mark diag[i] */ in MatLUFactorNumeric_SeqAIJ_Inode()
1800 pj = b->j + bdiag[i + 1] + 1; /* beginning of U(i,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1801 nz = bdiag[i] - bdiag[i + 1] - 1; /* num of entries in U(i,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1823 pc2 = b->a + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1824 pj = b->j + bdiag[i + 2] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1825 nz = bdiag[i + 1] - bdiag[i + 2] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1836 pc2 = b->a + bdiag[i + 1]; in MatLUFactorNumeric_SeqAIJ_Inode()
1847 pj = b->j + bdiag[i + 2] + 1; /* beginning of U(i+1,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1848 nz = bdiag[i + 1] - bdiag[i + 2] - 1; /* num of entries in U(i+1,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1869 pc3 = b->a + bdiag[i + 3] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1870 pj = b->j + bdiag[i + 3] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1871 nz = bdiag[i + 2] - bdiag[i + 3] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1882 pc3 = b->a + bdiag[i + 2]; in MatLUFactorNumeric_SeqAIJ_Inode()
1890 pj = b->j + bdiag[i + 3] + 1; /* beginning of U(i+2,:) */ in MatLUFactorNumeric_SeqAIJ_Inode()
1891 nz = bdiag[i + 2] - bdiag[i + 3] - 1; /* num of entries in U(i+2,:) excluding diag */ in MatLUFactorNumeric_SeqAIJ_Inode()
1911 pc4 = b->a + bdiag[i + 4] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1912 pj = b->j + bdiag[i + 4] + 1; in MatLUFactorNumeric_SeqAIJ_Inode()
1913 nz = bdiag[i + 3] - bdiag[i + 4] - 1; /* exclude diagonal */ in MatLUFactorNumeric_SeqAIJ_Inode()
1924 pc4 = b->a + bdiag[i + 3]; in MatLUFactorNumeric_SeqAIJ_Inode()
2423 MatScalar *ibdiag, *bdiag, work[25]; in MatInvertDiagonalForSOR_SeqAIJ_Inode() local
2441 …PetscCall(PetscMalloc3(cnt, &a->inode.ibdiag, cnt, &a->inode.bdiag, A->rmap->n, &a->inode.ssor_wor… in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2446 bdiag = a->inode.bdiag; in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2451 for (k = 0; k < nodesz; k++) bdiag[cnt + k * nodesz + j] = v[diag[row + j] - j + k]; in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2453 PetscCall(PetscArraycpy(ibdiag + cnt, bdiag + cnt, nodesz * nodesz)); in MatInvertDiagonalForSOR_SeqAIJ_Inode()
2497 MatScalar *ibdiag, *bdiag, *t; in MatSOR_SeqAIJ_Inode() local
2513 bdiag = a->inode.bdiag; in MatSOR_SeqAIJ_Inode()
3560 t[row] = b[row] - bdiag[cnt++] * x[row]; in MatSOR_SeqAIJ_Inode()
3566 tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 2]; in MatSOR_SeqAIJ_Inode()
3567 tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 3]; in MatSOR_SeqAIJ_Inode()
3577 tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 3] + x3 * bdiag[cnt + 6]; in MatSOR_SeqAIJ_Inode()
3578 tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 4] + x3 * bdiag[cnt + 7]; in MatSOR_SeqAIJ_Inode()
3579 tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 8]; in MatSOR_SeqAIJ_Inode()
3591 … tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 4] + x3 * bdiag[cnt + 8] + x4 * bdiag[cnt + 12]; in MatSOR_SeqAIJ_Inode()
3592 …tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 9] + x4 * bdiag[cnt + 13… in MatSOR_SeqAIJ_Inode()
3593 …tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 6] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3594 …tmp4 = x1 * bdiag[cnt + 3] + x2 * bdiag[cnt + 7] + x3 * bdiag[cnt + 11] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3608 …tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 15] +… in MatSOR_SeqAIJ_Inode()
3609 …tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 6] + x3 * bdiag[cnt + 11] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3610 …tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 7] + x3 * bdiag[cnt + 12] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3611 …tmp4 = x1 * bdiag[cnt + 3] + x2 * bdiag[cnt + 8] + x3 * bdiag[cnt + 13] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3612 …tmp5 = x1 * bdiag[cnt + 4] + x2 * bdiag[cnt + 9] + x3 * bdiag[cnt + 14] + x4 * bdiag[cnt + 1… in MatSOR_SeqAIJ_Inode()
3812 const MatScalar *bdiag = a->inode.bdiag; in MatMultDiagonalBlock_SeqAIJ_Inode() local
3826 x[row] = b[row] * bdiag[cnt++]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3832 tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 2]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3833 tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 3]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3842 tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 3] + x3 * bdiag[cnt + 6]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3843 tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 4] + x3 * bdiag[cnt + 7]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3844 tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 8]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3855 tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 4] + x3 * bdiag[cnt + 8] + x4 * bdiag[cnt + 12]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3856 … tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 9] + x4 * bdiag[cnt + 13]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3857 …tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 6] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 14]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3858 …tmp4 = x1 * bdiag[cnt + 3] + x2 * bdiag[cnt + 7] + x3 * bdiag[cnt + 11] + x4 * bdiag[cnt + 15]; in MatMultDiagonalBlock_SeqAIJ_Inode()
3871 …tmp1 = x1 * bdiag[cnt] + x2 * bdiag[cnt + 5] + x3 * bdiag[cnt + 10] + x4 * bdiag[cnt + 15] + x… in MatMultDiagonalBlock_SeqAIJ_Inode()
3872 …tmp2 = x1 * bdiag[cnt + 1] + x2 * bdiag[cnt + 6] + x3 * bdiag[cnt + 11] + x4 * bdiag[cnt + 16]… in MatMultDiagonalBlock_SeqAIJ_Inode()
3873 …tmp3 = x1 * bdiag[cnt + 2] + x2 * bdiag[cnt + 7] + x3 * bdiag[cnt + 12] + x4 * bdiag[cnt + 17]… in MatMultDiagonalBlock_SeqAIJ_Inode()
3874 …tmp4 = x1 * bdiag[cnt + 3] + x2 * bdiag[cnt + 8] + x3 * bdiag[cnt + 13] + x4 * bdiag[cnt + 18]… in MatMultDiagonalBlock_SeqAIJ_Inode()
3875 …tmp5 = x1 * bdiag[cnt + 4] + x2 * bdiag[cnt + 9] + x3 * bdiag[cnt + 14] + x4 * bdiag[cnt + 19]… in MatMultDiagonalBlock_SeqAIJ_Inode()
3997 c->inode.bdiag = NULL; in MatDuplicate_SeqAIJ_Inode()