Lines Matching refs:ilinkA
1111 PC_FieldSplitLink ilinkA = jac->head; in PCSetUpOnBlocks_FieldSplit_Schur() local
1112 KSP kspA = ilinkA->ksp, kspUpper = jac->kspupper; in PCSetUpOnBlocks_FieldSplit_Schur()
1134 PetscCall(VecGetArrayAndMemType(ilinkA->x, &array, &mtype)); in PCSetUpOnBlocks_FieldSplit_Schur()
1135 PetscCall(VecRestoreArrayAndMemType(ilinkA->x, &array)); in PCSetUpOnBlocks_FieldSplit_Schur()
1179 PC_FieldSplitLink ilinkA = jac->head; in PCSetUpOnBlocks_FieldSplit_GKB() local
1180 KSP ksp = ilinkA->ksp; in PCSetUpOnBlocks_FieldSplit_GKB()
1191 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApply_FieldSplit_Schur() local
1192 KSP kspA = ilinkA->ksp, kspLower = kspA, kspUpper = jac->kspupper; in PCApply_FieldSplit_Schur()
1200 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1202 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1203 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1204 PetscCall(KSPSolve(kspA, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1205 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1206 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1207 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1216 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1222 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1223 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1224 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1225 PetscCall(KSPSolve(kspA, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1226 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1227 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1228 PetscCall(MatMult(jac->C, ilinkA->y, ilinkD->x)); in PCApply_FieldSplit_Schur()
1231 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1239 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1253 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->x)); in PCApply_FieldSplit_Schur()
1254 PetscCall(VecScale(ilinkA->x, -1.)); in PCApply_FieldSplit_Schur()
1255 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1257 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1258 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1259 PetscCall(KSPSolve(kspA, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1260 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1261 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1263 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1264 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1270 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1271 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1272 PetscCall(PetscLogEventBegin(KSP_Solve_FS_L, kspLower, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1294 PetscCall(VecCopy(ilinkA->x, c)); in PCApply_FieldSplit_Schur()
1297 … PetscCall(VecCopy(c, ilinkA->y)); // retrieve the solution as the last column of the composed Mat in PCApply_FieldSplit_Schur()
1302 if (N == P) PetscCall(KSPSolve(kspLower, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1303 PetscCall(KSPCheckSolve(kspLower, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1304 PetscCall(PetscLogEventEnd(KSP_Solve_FS_L, kspLower, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1305 PetscCall(MatMult(jac->C, ilinkA->y, ilinkD->x)); in PCApply_FieldSplit_Schur()
1320 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->y)); in PCApply_FieldSplit_Schur()
1321 PetscCall(VecAXPY(ilinkA->x, -1.0, ilinkA->y)); in PCApply_FieldSplit_Schur()
1322 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1323 PetscCall(KSPSolve(kspA, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1324 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1325 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1326 } else PetscCall(MatMultAdd(AinvB, ilinkD->y, ilinkA->y, ilinkA->y)); in PCApply_FieldSplit_Schur()
1328 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_Schur()
1329 PetscCall(KSPSolve(kspA, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_Schur()
1330 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApply_FieldSplit_Schur()
1331 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->x)); in PCApply_FieldSplit_Schur()
1332 PetscCall(PetscLogEventBegin(KSP_Solve_FS_U, kspUpper, ilinkA->x, ilinkA->z, NULL)); in PCApply_FieldSplit_Schur()
1333 PetscCall(KSPSolve(kspUpper, ilinkA->x, ilinkA->z)); in PCApply_FieldSplit_Schur()
1334 PetscCall(KSPCheckSolve(kspUpper, pc, ilinkA->z)); in PCApply_FieldSplit_Schur()
1335 PetscCall(PetscLogEventEnd(KSP_Solve_FS_U, kspUpper, ilinkA->x, ilinkA->z, NULL)); in PCApply_FieldSplit_Schur()
1336 PetscCall(VecAXPY(ilinkA->y, -1.0, ilinkA->z)); in PCApply_FieldSplit_Schur()
1339 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1340 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1395 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCMatApply_FieldSplit_Schur() local
1396 KSP kspA = ilinkA->ksp, kspLower = kspA, kspUpper = jac->kspupper; in PCMatApply_FieldSplit_Schur()
1406 PetscCall(MatDenseScatter_Private(ilinkA->sctx, X, ilinkA->X, INSERT_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1408 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1409 PetscCall(KSPMatSolve(kspA, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1410 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1411 PetscCall(MatDenseScatter_Private(ilinkA->sctx, ilinkA->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1422 PetscCall(MatDenseScatter_Private(ilinkA->sctx, X, ilinkA->X, INSERT_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1423 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1424 PetscCall(KSPMatSolve(kspA, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1425 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1426 PetscCall(MatMatMult(jac->C, ilinkA->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkD->X)); in PCMatApply_FieldSplit_Schur()
1429 PetscCall(MatDenseScatter_Private(ilinkA->sctx, ilinkA->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1445 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1446 PetscCall(MatScale(ilinkA->X, -1.0)); in PCMatApply_FieldSplit_Schur()
1447 PetscCall(MatDenseScatter_Private(ilinkA->sctx, X, ilinkA->X, ADD_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1449 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1450 PetscCall(KSPMatSolve(kspA, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1451 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1452 PetscCall(MatDenseScatter_Private(ilinkA->sctx, ilinkA->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1458 PetscCall(MatDenseScatter_Private(ilinkA->sctx, X, ilinkA->X, INSERT_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1459 PetscCall(PetscLogEventBegin(KSP_Solve_FS_L, kspLower, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1508 PetscCall(MatCopy(ilinkA->X, C, SAME_NONZERO_PATTERN)); in PCMatApply_FieldSplit_Schur()
1511 …PetscCall(MatCopy(C, ilinkA->Y, SAME_NONZERO_PATTERN)); // retrieve solutions as last columns of t… in PCMatApply_FieldSplit_Schur()
1516 if (N == P) PetscCall(KSPMatSolve(kspLower, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1517 PetscCall(PetscLogEventEnd(KSP_Solve_FS_L, kspLower, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1518 PetscCall(MatMatMult(jac->C, ilinkA->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkD->X)); in PCMatApply_FieldSplit_Schur()
1531 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1532 PetscCall(MatAXPY(ilinkA->X, -1.0, ilinkA->Y, SAME_NONZERO_PATTERN)); in PCMatApply_FieldSplit_Schur()
1533 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1534 PetscCall(KSPMatSolve(kspA, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1535 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1537 PetscCall(MatMatMult(AinvB, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1538 PetscCall(MatAXPY(ilinkA->Y, 1.0, ilinkA->X, SAME_NONZERO_PATTERN)); in PCMatApply_FieldSplit_Schur()
1541 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->X, ilinkA->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1542 PetscCall(KSPMatSolve(kspA, ilinkA->X, ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1543 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1544 if (!ilinkA->Z) PetscCall(MatDuplicate(ilinkA->X, MAT_DO_NOT_COPY_VALUES, &ilinkA->Z)); in PCMatApply_FieldSplit_Schur()
1545 PetscCall(PetscLogEventBegin(KSP_Solve_FS_U, kspUpper, ilinkA->X, ilinkA->Z, NULL)); in PCMatApply_FieldSplit_Schur()
1546 PetscCall(KSPMatSolve(kspUpper, ilinkA->X, ilinkA->Z)); in PCMatApply_FieldSplit_Schur()
1547 PetscCall(PetscLogEventEnd(KSP_Solve_FS_U, kspUpper, ilinkA->X, ilinkA->Z, NULL)); in PCMatApply_FieldSplit_Schur()
1548 PetscCall(MatAXPY(ilinkA->Y, -1.0, ilinkA->Z, SAME_NONZERO_PATTERN)); in PCMatApply_FieldSplit_Schur()
1550 PetscCall(MatDenseScatter_Private(ilinkA->sctx, ilinkA->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1558 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApplyTranspose_FieldSplit_Schur() local
1559 KSP kspA = ilinkA->ksp, kspLower = kspA, kspUpper = jac->kspupper; in PCApplyTranspose_FieldSplit_Schur()
1565 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1567 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1568 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1569 PetscCall(KSPSolveTranspose(kspA, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1570 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1571 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1572 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1581 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1586 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1587 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1588 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1589 PetscCall(KSPSolveTranspose(kspA, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1590 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1591 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1592 PetscCall(MatMultTranspose(jac->B, ilinkA->y, ilinkD->x)); in PCApplyTranspose_FieldSplit_Schur()
1595 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1603 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1616 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->x)); in PCApplyTranspose_FieldSplit_Schur()
1617 PetscCall(VecScale(ilinkA->x, -1.)); in PCApplyTranspose_FieldSplit_Schur()
1618 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1620 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1621 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1622 PetscCall(KSPSolveTranspose(kspA, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1623 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1624 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1626 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1627 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1630 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1631 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1632 PetscCall(PetscLogEventBegin(KSP_Solve_FS_U, kspUpper, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1633 PetscCall(KSPSolveTranspose(kspUpper, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1634 PetscCall(KSPCheckSolve(kspUpper, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1635 PetscCall(PetscLogEventEnd(KSP_Solve_FS_U, kspUpper, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1636 PetscCall(MatMultTranspose(jac->B, ilinkA->y, ilinkD->x)); in PCApplyTranspose_FieldSplit_Schur()
1650 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1651 PetscCall(VecAXPY(ilinkA->x, -1.0, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1652 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1653 PetscCall(KSPSolveTranspose(kspA, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1654 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1655 PetscCall(PetscLogEventEnd(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1657 PetscCall(PetscLogEventBegin(ilinkA->event, kspA, ilinkA->x, ilinkA->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1658 PetscCall(KSPSolveTranspose(kspA, ilinkA->x, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1659 PetscCall(KSPCheckSolve(kspA, pc, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1660 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->x)); in PCApplyTranspose_FieldSplit_Schur()
1661 PetscCall(PetscLogEventBegin(KSP_Solve_FS_L, kspLower, ilinkA->x, ilinkA->z, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1662 PetscCall(KSPSolveTranspose(kspLower, ilinkA->x, ilinkA->z)); in PCApplyTranspose_FieldSplit_Schur()
1663 PetscCall(KSPCheckSolve(kspLower, pc, ilinkA->z)); in PCApplyTranspose_FieldSplit_Schur()
1664 PetscCall(PetscLogEventEnd(KSP_Solve_FS_L, kspLower, ilinkA->x, ilinkA->z, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1665 PetscCall(VecAXPY(ilinkA->y, -1.0, ilinkA->z)); in PCApplyTranspose_FieldSplit_Schur()
1668 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1669 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1862 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApply_FieldSplit_GKB() local
1863 KSP ksp = ilinkA->ksp; in PCApply_FieldSplit_GKB()
1870 PetscCall(VecScatterBegin(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_GKB()
1872 PetscCall(VecScatterEnd(ilinkA->sctx, x, ilinkA->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_GKB()
1888 PetscCall(MatMultAdd(jac->B, ilinkD->x, ilinkA->x, ilinkA->x)); /* q = q + nu*B*b */ in PCApply_FieldSplit_GKB()
1896 PetscCall(PetscLogEventBegin(ilinkA->event, ksp, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_GKB()
1897 PetscCall(KSPSolve(ksp, ilinkA->x, ilinkA->y)); in PCApply_FieldSplit_GKB()
1898 PetscCall(KSPCheckSolve(ksp, pc, ilinkA->y)); in PCApply_FieldSplit_GKB()
1899 PetscCall(PetscLogEventEnd(ilinkA->event, ksp, ilinkA->x, ilinkA->y, NULL)); in PCApply_FieldSplit_GKB()
1900 PetscCall(MatMultHermitianTranspose(jac->B, ilinkA->y, work1)); in PCApply_FieldSplit_GKB()
1909 PetscCall(PetscLogEventBegin(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1912 PetscCall(PetscLogEventEnd(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1925 PetscCall(VecAXPY(ilinkA->y, z, u)); in PCApply_FieldSplit_GKB()
1943 PetscCall(PetscLogEventBegin(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1946 PetscCall(PetscLogEventEnd(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1959 PetscCall(VecAXPY(ilinkA->y, z, u)); /* r = r + z*u */ in PCApply_FieldSplit_GKB()
1961 PetscCall(MatMult(jac->H, ilinkA->y, Hu)); /* ||u||_H = u'*H*u */ in PCApply_FieldSplit_GKB()
1962 PetscCall(VecDot(Hu, ilinkA->y, &nrmz2)); in PCApply_FieldSplit_GKB()
1975 PetscCall(VecScatterBegin(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_GKB()
1976 PetscCall(VecScatterEnd(ilinkA->sctx, ilinkA->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_GKB()