Lines Matching refs:yi
39 Vec *yi, *yiwork, *ziwork; member
339 PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt)); in StateMatMult()
341 PetscCall(MatMult(user->JsBlock, user->yi[0], user->yiwork[0])); in StateMatMult()
345 PetscCall(MatMult(user->JsBlock, user->yi[i], user->yiwork[i])); in StateMatMult()
346 PetscCall(MatMult(user->M, user->yi[i - 1], user->ziwork[i - 1])); in StateMatMult()
360 PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt)); in StateMatMultTranspose()
364 PetscCall(MatMultTranspose(user->JsBlock, user->yi[i], user->yiwork[i])); in StateMatMultTranspose()
365 PetscCall(MatMult(user->M, user->yi[i + 1], user->ziwork[i + 1])); in StateMatMultTranspose()
371 PetscCall(MatMultTranspose(user->JsBlock, user->yi[i], user->yiwork[i])); in StateMatMultTranspose()
416 PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt)); in DesignMatMult()
419 PetscCall(VecPointwiseMult(user->uxiwork[i], user->yi[i], user->uxiwork[i])); in DesignMatMult()
420 PetscCall(VecPointwiseMult(user->uyiwork[i], user->yi[i], user->uyiwork[i])); in DesignMatMult()
436 PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt)); in DesignMatMultTranspose()
441 PetscCall(VecPointwiseMult(user->uxiwork[i], user->yi[i], user->uxiwork[i])); in DesignMatMultTranspose()
442 PetscCall(VecPointwiseMult(user->uyiwork[i], user->yi[i], user->uyiwork[i])); in DesignMatMultTranspose()
491 PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt)); in StateMatInvMult()
496 PetscCall(KSPSolve(user->solver, user->yi[0], user->yiwork[0])); in StateMatInvMult()
502 PetscCall(VecAXPY(user->yi[i], 1.0, user->ziwork[i - 1])); in StateMatInvMult()
504 PetscCall(KSPSolve(user->solver, user->yi[i], user->yiwork[i])); in StateMatInvMult()
521 PetscCall(Scatter_yi(X, user->yi, user->yi_scatter, user->nt)); in StateMatInvTransposeMult()
527 PetscCall(KSPSolveTranspose(user->solver, user->yi[i], user->yiwork[i])); in StateMatInvTransposeMult()
534 PetscCall(VecAXPY(user->yi[i], 1.0, user->ziwork[i + 1])); in StateMatInvTransposeMult()
536 PetscCall(KSPSolveTranspose(user->solver, user->yi[i], user->yiwork[i])); in StateMatInvTransposeMult()
583 PetscCall(Scatter_yi(user->y, user->yi, user->yi_scatter, user->nt)); in FormConstraints()
587 PetscCall(MatMult(user->JsBlock, user->yi[0], user->yiwork[0])); in FormConstraints()
591 PetscCall(MatMult(user->JsBlock, user->yi[i], user->yiwork[i])); in FormConstraints()
592 PetscCall(MatMult(user->M, user->yi[i - 1], user->ziwork[i - 1])); in FormConstraints()
649 PetscErrorCode Scatter_yi(Vec y, Vec *yi, VecScatter *scat, PetscInt nt) in Scatter_yi() argument
655 PetscCall(VecScatterBegin(scat[i], y, yi[i], INSERT_VALUES, SCATTER_FORWARD)); in Scatter_yi()
656 PetscCall(VecScatterEnd(scat[i], y, yi[i], INSERT_VALUES, SCATTER_FORWARD)); in Scatter_yi()
661 PetscErrorCode Gather_yi(Vec y, Vec *yi, VecScatter *scat, PetscInt nt) in Gather_yi() argument
667 PetscCall(VecScatterBegin(scat[i], yi[i], y, INSERT_VALUES, SCATTER_REVERSE)); in Gather_yi()
668 PetscCall(VecScatterEnd(scat[i], yi[i], y, INSERT_VALUES, SCATTER_REVERSE)); in Gather_yi()
676 Vec XX, YY, XXwork, YYwork, yi, uxi, ui, bc; in HyperbolicInitialize() local
876 PetscCall(VecCreate(PETSC_COMM_WORLD, &yi)); in HyperbolicInitialize()
877 PetscCall(VecSetSizes(yi, PETSC_DECIDE, user->mx * user->mx)); in HyperbolicInitialize()
878 PetscCall(VecSetFromOptions(yi)); in HyperbolicInitialize()
879 PetscCall(VecDuplicateVecs(yi, user->nt, &user->yi)); in HyperbolicInitialize()
880 PetscCall(VecDuplicateVecs(yi, user->nt, &user->yiwork)); in HyperbolicInitialize()
881 PetscCall(VecDuplicateVecs(yi, user->nt, &user->ziwork)); in HyperbolicInitialize()
883 PetscCall(VecGetOwnershipRange(user->yi[i], &lo, &hi)); in HyperbolicInitialize()
886 PetscCall(VecScatterCreate(user->y, is_from_y, user->yi[i], is_to_yi, &user->yi_scatter[i])); in HyperbolicInitialize()
1104 PetscCall(VecDestroy(&yi)); in HyperbolicInitialize()
1150 PetscCall(VecDestroyVecs(user->nt, &user->yi)); in HyperbolicDestroy()