Lines Matching refs:ilinkD
1191 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApply_FieldSplit_Schur() local
1201 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1208 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1209 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1211 PetscCall(KSPSolve(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApply_FieldSplit_Schur()
1212 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApply_FieldSplit_Schur()
1214 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1215 PetscCall(VecScale(ilinkD->y, jac->schurscale)); in PCApply_FieldSplit_Schur()
1217 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1218 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1228 PetscCall(MatMult(jac->C, ilinkA->y, ilinkD->x)); in PCApply_FieldSplit_Schur()
1229 PetscCall(VecScale(ilinkD->x, -1.)); in PCApply_FieldSplit_Schur()
1230 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1232 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1233 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1235 PetscCall(KSPSolve(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApply_FieldSplit_Schur()
1237 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApply_FieldSplit_Schur()
1238 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1240 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1241 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1245 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1246 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1247 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1249 PetscCall(KSPSolve(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApply_FieldSplit_Schur()
1251 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApply_FieldSplit_Schur()
1252 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1253 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->x)); in PCApply_FieldSplit_Schur()
1256 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1262 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1305 PetscCall(MatMult(jac->C, ilinkA->y, ilinkD->x)); in PCApply_FieldSplit_Schur()
1306 PetscCall(VecScale(ilinkD->x, -1.0)); in PCApply_FieldSplit_Schur()
1307 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1308 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_Schur()
1310 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1312 PetscCall(KSPSolve(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApply_FieldSplit_Schur()
1314 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApply_FieldSplit_Schur()
1315 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApply_FieldSplit_Schur()
1316 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1320 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->y)); in PCApply_FieldSplit_Schur()
1326 } else PetscCall(MatMultAdd(AinvB, ilinkD->y, ilinkA->y, ilinkA->y)); in PCApply_FieldSplit_Schur()
1331 PetscCall(MatMult(jac->B, ilinkD->y, ilinkA->x)); in PCApply_FieldSplit_Schur()
1338 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_Schur()
1395 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCMatApply_FieldSplit_Schur() local
1407 PetscCall(MatDenseScatter_Private(ilinkD->sctx, X, ilinkD->X, INSERT_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1412 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1414 PetscCall(KSPMatSolve(jac->kspschur, ilinkD->X, ilinkD->Y)); in PCMatApply_FieldSplit_Schur()
1416 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1417 PetscCall(MatScale(ilinkD->Y, jac->schurscale)); in PCMatApply_FieldSplit_Schur()
1418 PetscCall(MatDenseScatter_Private(ilinkD->sctx, ilinkD->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1426 PetscCall(MatMatMult(jac->C, ilinkA->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkD->X)); in PCMatApply_FieldSplit_Schur()
1427 PetscCall(MatScale(ilinkD->X, -1.0)); in PCMatApply_FieldSplit_Schur()
1428 PetscCall(MatDenseScatter_Private(ilinkD->sctx, X, ilinkD->X, ADD_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1430 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1432 PetscCall(KSPMatSolve(jac->kspschur, ilinkD->X, ilinkD->Y)); in PCMatApply_FieldSplit_Schur()
1434 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1435 PetscCall(MatDenseScatter_Private(ilinkD->sctx, ilinkD->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1439 PetscCall(MatDenseScatter_Private(ilinkD->sctx, X, ilinkD->X, INSERT_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1440 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1442 PetscCall(KSPMatSolve(jac->kspschur, ilinkD->X, ilinkD->Y)); in PCMatApply_FieldSplit_Schur()
1444 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1445 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1448 PetscCall(MatDenseScatter_Private(ilinkD->sctx, ilinkD->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1518 PetscCall(MatMatMult(jac->C, ilinkA->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkD->X)); in PCMatApply_FieldSplit_Schur()
1519 PetscCall(MatScale(ilinkD->X, -1.0)); in PCMatApply_FieldSplit_Schur()
1520 PetscCall(MatDenseScatter_Private(ilinkD->sctx, X, ilinkD->X, ADD_VALUES, SCATTER_FORWARD)); in PCMatApply_FieldSplit_Schur()
1522 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1524 PetscCall(KSPMatSolve(jac->kspschur, ilinkD->X, ilinkD->Y)); in PCMatApply_FieldSplit_Schur()
1526 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->X, ilinkD->Y, NULL)); in PCMatApply_FieldSplit_Schur()
1527 PetscCall(MatDenseScatter_Private(ilinkD->sctx, ilinkD->Y, Y, INSERT_VALUES, SCATTER_REVERSE)); in PCMatApply_FieldSplit_Schur()
1531 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->Y)); in PCMatApply_FieldSplit_Schur()
1537 PetscCall(MatMatMult(AinvB, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1543 PetscCall(MatMatMult(jac->B, ilinkD->Y, MAT_REUSE_MATRIX, PETSC_DETERMINE, &ilinkA->X)); in PCMatApply_FieldSplit_Schur()
1558 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApplyTranspose_FieldSplit_Schur() local
1566 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1573 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1574 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1576 PetscCall(KSPSolveTranspose(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1578 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1579 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1580 PetscCall(VecScale(ilinkD->y, jac->schurscale)); in PCApplyTranspose_FieldSplit_Schur()
1582 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1583 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1592 PetscCall(MatMultTranspose(jac->B, ilinkA->y, ilinkD->x)); in PCApplyTranspose_FieldSplit_Schur()
1593 PetscCall(VecScale(ilinkD->x, -1.)); in PCApplyTranspose_FieldSplit_Schur()
1594 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1596 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1597 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1599 PetscCall(KSPSolveTranspose(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1601 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1602 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1604 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1605 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1608 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1609 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1610 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1612 PetscCall(KSPSolveTranspose(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1614 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1615 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1616 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->x)); in PCApplyTranspose_FieldSplit_Schur()
1619 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1625 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1636 PetscCall(MatMultTranspose(jac->B, ilinkA->y, ilinkD->x)); in PCApplyTranspose_FieldSplit_Schur()
1637 PetscCall(VecScale(ilinkD->x, -1.0)); in PCApplyTranspose_FieldSplit_Schur()
1638 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1639 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, ADD_VALUES, SCATTER_FORWARD)); in PCApplyTranspose_FieldSplit_Schur()
1641 PetscCall(PetscLogEventBegin(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1643 PetscCall(KSPSolveTranspose(jac->kspschur, ilinkD->x, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1645 PetscCall(KSPCheckSolve(jac->kspschur, pc, ilinkD->y)); in PCApplyTranspose_FieldSplit_Schur()
1646 PetscCall(PetscLogEventEnd(KSP_Solve_FS_S, jac->kspschur, ilinkD->x, ilinkD->y, NULL)); in PCApplyTranspose_FieldSplit_Schur()
1647 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1650 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->y)); in PCApplyTranspose_FieldSplit_Schur()
1660 PetscCall(MatMultTranspose(jac->C, ilinkD->y, ilinkA->x)); in PCApplyTranspose_FieldSplit_Schur()
1667 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApplyTranspose_FieldSplit_Schur()
1862 PC_FieldSplitLink ilinkA = jac->head, ilinkD = ilinkA->next; in PCApply_FieldSplit_GKB() local
1871 PetscCall(VecScatterBegin(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_GKB()
1873 PetscCall(VecScatterEnd(ilinkD->sctx, x, ilinkD->x, INSERT_VALUES, SCATTER_FORWARD)); in PCApply_FieldSplit_GKB()
1887 PetscCall(VecScale(ilinkD->x, jac->gkbnu)); in PCApply_FieldSplit_GKB()
1888 PetscCall(MatMultAdd(jac->B, ilinkD->x, ilinkA->x, ilinkA->x)); /* q = q + nu*B*b */ in PCApply_FieldSplit_GKB()
1901 PetscCall(VecAXPBY(work1, 1.0 / nu, -1.0, ilinkD->x)); /* c = b - B'*x */ in PCApply_FieldSplit_GKB()
1926 PetscCall(VecCopy(d, ilinkD->y)); in PCApply_FieldSplit_GKB()
1927 PetscCall(VecScale(ilinkD->y, -z)); in PCApply_FieldSplit_GKB()
1960 PetscCall(VecAXPY(ilinkD->y, -z, d)); /* p = p - z*d */ in PCApply_FieldSplit_GKB()
1977 PetscCall(VecScatterBegin(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_GKB()
1978 PetscCall(VecScatterEnd(ilinkD->sctx, ilinkD->y, y, INSERT_VALUES, SCATTER_REVERSE)); in PCApply_FieldSplit_GKB()