| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcnullspace.c | 149 Vec work1, work2; in PCBDDCNullSpaceAssembleCorrection() local
156 PetscCall(VecDuplicate(shell_ctx->fw[0], &work2)); in PCBDDCNullSpaceAssembleCorrection()
178 PetscCall(VecSetRandom(work2, NULL)); in PCBDDCNullSpaceAssembleCorrection()
179 PetscCall(MatMult(local_mat, work2, work1)); in PCBDDCNullSpaceAssembleCorrection()
182 PetscCall(VecAXPY(work1, -1., work2)); in PCBDDCNullSpaceAssembleCorrection()
198 PetscCall(VecSetRandom(work2, NULL)); in PCBDDCNullSpaceAssembleCorrection()
199 PetscCall(MatMult(local_mat, work2, work1)); in PCBDDCNullSpaceAssembleCorrection()
211 PetscCall(VecAXPY(work1, -1., work2)); in PCBDDCNullSpaceAssembleCorrection()
225 PetscCall(VecDestroy(&work2)); in PCBDDCNullSpaceAssembleCorrection()
230 Vec work1, work2, work3; in PCBDDCNullSpaceAssembleCorrection() local
[all …]
|
| /petsc/src/ksp/pc/impls/composite/ |
| H A D | composite.c | 18 Vec work2; member
39 if (next->next && !jac->work2) { /* allocate second work vector */ in PCApply_Composite_Multiplicative()
40 PetscCall(VecDuplicate(jac->work1, &jac->work2)); in PCApply_Composite_Multiplicative()
47 PetscCall(VecWAXPY(jac->work2, -1.0, jac->work1, x)); /* work2 <- x - work1 */ in PCApply_Composite_Multiplicative()
48 PetscCall(PCApply(next->pc, jac->work2, jac->work1)); /* work1 <- C work2 */ in PCApply_Composite_Multiplicative()
55 PetscCall(VecWAXPY(jac->work2, -1.0, jac->work1, x)); in PCApply_Composite_Multiplicative()
56 PetscCall(PCApply(next->pc, jac->work2, jac->work1)); in PCApply_Composite_Multiplicative()
71 if (next->next && !jac->work2) { /* allocate second work vector */ in PCApplyTranspose_Composite_Multiplicative()
72 PetscCall(VecDuplicate(jac->work1, &jac->work2)); in PCApplyTranspose_Composite_Multiplicative()
81 PetscCall(VecWAXPY(jac->work2, -1.0, jac->work1, x)); in PCApplyTranspose_Composite_Multiplicative()
[all …]
|
| /petsc/src/mat/impls/lrc/ |
| H A D | lrc.c | 9 Vec work1, work2; /* sequential vectors that hold partial products */ member
57 PetscCall(VecGetArrayWrite(Na->work2, &w2)); in MatMult_LRC_kernel()
61 PetscCall(VecRestoreArrayWrite(Na->work2, &w2)); in MatMult_LRC_kernel()
64 PetscCall(VecPointwiseMult(Na->work2, Na->c, Na->work2)); in MatMult_LRC_kernel()
75 PetscCall(MatMultAdd(Uloc, Na->work2, yl, yl)); in MatMult_LRC_kernel()
78 PetscCall(MatMult(Uloc, Na->work2, yl)); in MatMult_LRC_kernel()
110 PetscCall(VecDestroy(&Na->work2)); in MatDestroy_LRC()
285 PetscCall(VecDuplicate(Na->work1, &Na->work2)); in MatSetUp_LRC()
|
| /petsc/src/ksp/ksp/utils/schurm/ |
| H A D | schurm.c | 54 if (!Na->work2) PetscCall(MatCreateVecs(Na->A, &Na->work2, NULL)); in MatMultTranspose_SchurComplement()
56 PetscCall(KSPSolveTranspose(Na->ksp, Na->work1, Na->work2)); in MatMultTranspose_SchurComplement()
57 PetscCall(MatMultTranspose(Na->B, Na->work2, y)); in MatMultTranspose_SchurComplement()
72 if (!Na->work2) PetscCall(MatCreateVecs(Na->A, &Na->work2, NULL)); in MatMult_SchurComplement()
74 PetscCall(KSPSolve(Na->ksp, Na->work1, Na->work2)); in MatMult_SchurComplement()
75 PetscCall(MatMult(Na->C, Na->work2, y)); in MatMult_SchurComplement()
90 if (!Na->work2) PetscCall(MatCreateVecs(Na->A, &Na->work2, NULL)); in MatMultAdd_SchurComplement()
92 PetscCall(KSPSolve(Na->ksp, Na->work1, Na->work2)); in MatMultAdd_SchurComplement()
94 PetscCall(VecScale(Na->work2, -1.0)); in MatMultAdd_SchurComplement()
95 PetscCall(MatMultAdd(Na->C, Na->work2, z, z)); in MatMultAdd_SchurComplement()
[all …]
|
| H A D | schurm.h | 10 Vec work1, work2; member
|
| /petsc/src/snes/interface/ |
| H A D | snesut.c | 145 Vec snes_solution, work1, work2; in KSPMonitorSNESResidual() local
154 PetscCall(VecDuplicate(snes_solution, &work2)); in KSPMonitorSNESResidual()
157 PetscCall(SNESComputeFunction(snes, work1, work2)); in KSPMonitorSNESResidual()
158 PetscCall(VecNorm(work2, NORM_2, &snorm)); in KSPMonitorSNESResidual()
160 PetscCall(VecDestroy(&work2)); in KSPMonitorSNESResidual()
201 Vec snes_solution, work1, work2; in KSPMonitorSNESResidualDrawLG() local
211 PetscCall(VecDuplicate(snes_solution, &work2)); in KSPMonitorSNESResidualDrawLG()
214 PetscCall(SNESComputeFunction(snes, work1, work2)); in KSPMonitorSNESResidualDrawLG()
215 PetscCall(VecNorm(work2, NORM_2, &snorm)); in KSPMonitorSNESResidualDrawLG()
217 PetscCall(VecDestroy(&work2)); in KSPMonitorSNESResidualDrawLG()
|
| /petsc/src/tao/leastsquares/impls/pounders/ |
| H A D | pounders.c | 141 …P->gqt_maxits, gnorm, qmin, mfqP->Xsubproblem, &info, &its, mfqP->work, mfqP->work2, mfqP->work3)); in gqtwrap()
514 for (j = 0; j < mfqP->n; j++) mfqP->work2[j] = mfqP->Gres[j]; in modelimprove()
515 …asn, &half, mfqP->Hres, &blasn, &mfqP->Q_tmp[i * mfqP->npmax], &blas1, &one, mfqP->work2, &blas1)); in modelimprove()
516 …ot", mfqP->work[i] = BLASdot_(&blasn, &mfqP->Q_tmp[i * mfqP->npmax], &blas1, mfqP->work2, &blas1)); in modelimprove()
543 PetscCallBLAS("BLAScopy", BLAScopy_(&blasn, mfqP->work, &ione, mfqP->work2, &ione)); in affpoints()
550 …r_("R", "N", &ione, &blasn, &blask, mfqP->Q, &blasnpmax, mfqP->tau, mfqP->work2, &ione, mfqP->mwor… in affpoints()
553 PetscCallBLAS("BLASnrm2", proj = BLASnrm2_(&blasj, &mfqP->work2[mfqP->nmodelpoints], &ione)); in affpoints()
786 mfqP->work2[k] = 0.0; in TaoSolve_POUNDERS()
787 …for (l = 0; l < mfqP->n; l++) mfqP->work2[k] += mfqP->H[j + mfqP->m * (k + l * mfqP->n)] * mfqP->w… in TaoSolve_POUNDERS()
790 mfqP->C[j] += mfqP->work[i] * (mfqP->Fdiff[i + mfqP->n * j] + 0.5 * mfqP->work2[i]); in TaoSolve_POUNDERS()
[all …]
|
| H A D | pounders.h | 14 PetscReal *work2; /* (n) */ member
|
| /petsc/src/vec/is/sf/utils/ |
| H A D | sfutils.c | 600 PetscInt cum = 0, start, *work2; in PetscLayoutMapLocal() local
603 PetscCall(PetscCalloc1(N, &work2)); in PetscLayoutMapLocal()
610 if (idxs[r] >= 0) work2[r] = start + cum++; in PetscLayoutMapLocal()
611 PetscCall(PetscSFReduceBegin(sf, MPIU_INT, work2, work, MPI_REPLACE)); in PetscLayoutMapLocal()
612 PetscCall(PetscSFReduceEnd(sf, MPIU_INT, work2, work, MPI_REPLACE)); in PetscLayoutMapLocal()
613 PetscCall(PetscFree(work2)); in PetscLayoutMapLocal()
|
| /petsc/src/ksp/pc/interface/ |
| H A D | precon.c | 812 …Vec work2; /* this is expensive, but to fix requires a second work vector argument to PCApplyBAorA… in PCApplyBAorAB() local
813 PetscCall(VecDuplicate(x, &work2)); in PCApplyBAorAB()
814 PetscCall(PCDiagonalScaleRight(pc, x, work2)); in PCApplyBAorAB()
815 PetscUseTypeMethod(pc, applyBA, side, work2, y, work); in PCApplyBAorAB()
817 PetscCall(VecDestroy(&work2)); in PCApplyBAorAB()
|
| /petsc/src/ksp/pc/impls/fieldsplit/ |
| H A D | fieldsplit.c | 1864 Vec u, v, Hu, d, work1, work2; in PCApply_FieldSplit_GKB() local
1880 work2 = jac->w2; in PCApply_FieldSplit_GKB()
1908 PetscCall(MatMult(jac->B, v, work2)); /* u = H^{-1}*B*v */ in PCApply_FieldSplit_GKB()
1909 PetscCall(PetscLogEventBegin(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1910 PetscCall(KSPSolve(ksp, work2, u)); in PCApply_FieldSplit_GKB()
1912 PetscCall(PetscLogEventEnd(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1940 PetscCall(MatMult(jac->B, v, work2)); /* u <- H^{-1}*(B*v-beta*H*u) */ in PCApply_FieldSplit_GKB()
1942 PetscCall(VecAXPY(work2, -beta, Hu)); in PCApply_FieldSplit_GKB()
1943 PetscCall(PetscLogEventBegin(ilinkA->event, ksp, work2, u, NULL)); in PCApply_FieldSplit_GKB()
1944 PetscCall(KSPSolve(ksp, work2, u)); in PCApply_FieldSplit_GKB()
[all …]
|
| /petsc/src/dm/dt/dualspace/impls/lagrange/ |
| H A D | dspacelagrange.c | 661 PetscReal *workF, *workF2, *workT, *workT2, *work, *work2; in PetscLagNodeIndicesTensor() local
695 …tscCall(PetscMalloc6(MkT, &workT, MkT, &workT2, MkF, &workF, MkF, &workF2, Nk, &work, Nk, &work2)); in PetscLagNodeIndicesTensor()
734 for (d = 0; d < Nk; d++) work2[d] = work[d]; in PetscLagNodeIndicesTensor()
735 PetscCall(PetscDTAltVStar(dim, PetscAbsInt(formDegree), -1, work2, work)); in PetscLagNodeIndicesTensor()
742 PetscCall(PetscFree6(workT, workT2, workF, workF2, work, work2)); in PetscLagNodeIndicesTensor()
|