| /petsc/src/ksp/ksp/impls/tcqmr/ |
| H A D | tcqmrimpl.h | 14 #define r ksp->work[VEC_OFFSET + 1]
15 #define um1 ksp->work[VEC_OFFSET + 2]
16 #define u ksp->work[VEC_OFFSET + 3]
17 #define vm1 ksp->work[VEC_OFFSET + 4]
18 #define v ksp->work[VEC_OFFSET + 5]
19 #define v0 ksp->work[VEC_OFFSET + 6]
20 #define pvec1 ksp->work[VEC_OFFSET + 7]
21 #define pvec2 ksp->work[VEC_OFFSET + 8]
22 #define p ksp->work[VEC_OFFSET + 9]
23 #define y ksp->work[VEC_OFFSET + 10]
[all …]
|
| /petsc/src/ts/impls/explicit/ssp/ |
| H A D | ssp.c | 13 Vec *work; member
18 static PetscErrorCode TSSSPGetWorkVectors(TS ts, PetscInt n, Vec **work) in TSSSPGetWorkVectors() argument
25 if (ssp->nwork > 0) PetscCall(VecDestroyVecs(ssp->nwork, &ssp->work)); in TSSSPGetWorkVectors()
26 PetscCall(VecDuplicateVecs(ts->vec_sol, n, &ssp->work)); in TSSSPGetWorkVectors()
29 *work = ssp->work; in TSSSPGetWorkVectors()
34 static PetscErrorCode TSSSPRestoreWorkVectors(TS ts, PetscInt n, Vec **work) in TSSSPRestoreWorkVectors() argument
40 PetscCheck(*work == ssp->work, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Wrong work vectors checked out"); in TSSSPRestoreWorkVectors()
42 *work = NULL; in TSSSPRestoreWorkVectors()
56 Vec *work, F; in TSSSPStep_RK_2() local
61 PetscCall(TSSSPGetWorkVectors(ts, 2, &work)); in TSSSPStep_RK_2()
[all …]
|
| /petsc/src/mat/utils/ |
| H A D | getcolv.c | 167 PetscReal *work; in MatGetColumnSums() local
176 PetscCall(PetscCalloc1(n, &work)); in MatGetColumnSums()
177 PetscCall(MatGetColumnSumsRealPart(A, work)); in MatGetColumnSums()
178 for (i = 0; i < n; i++) sums[i] = work[i]; in MatGetColumnSums()
179 PetscCall(MatGetColumnSumsImaginaryPart(A, work)); in MatGetColumnSums()
180 for (i = 0; i < n; i++) sums[i] += work[i] * PETSC_i; in MatGetColumnSums()
181 PetscCall(PetscFree(work)); in MatGetColumnSums()
255 PetscReal *work; in MatGetColumnMeans() local
264 PetscCall(PetscCalloc1(n, &work)); in MatGetColumnMeans()
265 PetscCall(MatGetColumnMeansRealPart(A, work)); in MatGetColumnMeans()
[all …]
|
| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | comm.c | 62 PetscErrorCode PCTFS_giop(PetscInt *vals, PetscInt *work, PetscInt n, PetscInt *oprs) in PCTFS_giop() argument
71 …tscCheck(vals && work && oprs, PETSC_COMM_SELF, PETSC_ERR_PLIB, "PCTFS_giop() :: vals=%p, work=%p,… in PCTFS_giop()
97 …PetscCallMPI(MPI_Recv(work, n, MPIU_INT, MPI_ANY_SOURCE, MSGTAG0 + edge_not_pow_2, MPI_COMM_WORLD,… in PCTFS_giop()
98 PetscCall((*fp)(vals, work, n, oprs)); in PCTFS_giop()
109 …PetscCallMPI(MPI_Recv(work, n, MPIU_INT, MPI_ANY_SOURCE, MSGTAG2 + dest, MPI_COMM_WORLD, &status)); in PCTFS_giop()
110 PetscCall((*fp)(vals, work, n, oprs)); in PCTFS_giop()
139 PetscErrorCode PCTFS_grop(PetscScalar *vals, PetscScalar *work, PetscInt n, PetscInt *oprs) in PCTFS_grop() argument
148 …tscCheck(vals && work && oprs, PETSC_COMM_SELF, PETSC_ERR_PLIB, "PCTFS_grop() :: vals=%p, work=%p,… in PCTFS_grop()
173 …PetscCallMPI(MPI_Recv(work, n, MPIU_SCALAR, MPI_ANY_SOURCE, MSGTAG0 + edge_not_pow_2, MPI_COMM_WOR… in PCTFS_grop()
174 PetscCall((*fp)(vals, work, n, oprs)); in PCTFS_grop()
[all …]
|
| /petsc/src/ts/impls/implicit/radau5/ |
| H A D | radau5.c | 8 Vec work, workf; member
19 PetscCallAbort(PETSC_COMM_SELF, VecPlaceArray(cvode->work, Y)); in FVPOL()
27 PetscCallAbort(PETSC_COMM_SELF, TSComputeRHSFunction(ts, *X, cvode->work, cvode->workf)); in FVPOL()
31 PetscCallAbort(PETSC_COMM_SELF, VecDuplicate(cvode->work, &yydot)); in FVPOL()
33 …PetscCallAbort(PETSC_COMM_SELF, TSComputeIFunction(ts, *X, cvode->work, yydot, cvode->workf, PETSC… in FVPOL()
38 PetscCallAbort(PETSC_COMM_SELF, VecResetArray(cvode->work)); in FVPOL()
50 PetscCallAbort(PETSC_COMM_SELF, VecPlaceArray(cvode->work, Y)); in JVPOL()
51 PetscCallAbort(PETSC_COMM_SELF, VecDuplicate(cvode->work, &yydot)); in JVPOL()
55 …PetscCallAbort(PETSC_COMM_SELF, TSComputeIJacobian(ts, *X, cvode->work, yydot, 0, mat, mat, PETSC_… in JVPOL()
59 PetscCallAbort(PETSC_COMM_SELF, VecResetArray(cvode->work)); in JVPOL()
[all …]
|
| /petsc/src/mat/graphops/color/impls/minpack/ |
| H A D | color.c | 15 PetscInt *work; in MatFDColoringDegreeSequence_Minpack() local
18 PetscCall(PetscMalloc1(m, &work)); in MatFDColoringDegreeSequence_Minpack()
21 PetscCall(MINPACKdegr(&m, cja, cia, rja, ria, *seq, work)); in MatFDColoringDegreeSequence_Minpack()
23 PetscCall(PetscFree(work)); in MatFDColoringDegreeSequence_Minpack()
29 PetscInt *list, *work, clique, *seq, *coloring, n; in MatColoringApply_SL() local
62 PetscCall(PetscMalloc2(n, &list, 4 * n, &work)); in MatColoringApply_SL()
64 …all(MINPACKslo(&n, cja, cia, rja, ria, seq, list, &clique, work, work + n, work + 2 * n, work + 3 … in MatColoringApply_SL()
67 PetscCall(MINPACKseq(&n, cja, cia, rja, ria, list, coloring, &ncolors, work)); in MatColoringApply_SL()
69 PetscCall(PetscFree2(list, work)); in MatColoringApply_SL()
129 PetscInt *list, *work, *seq, *coloring, n; in MatColoringApply_LF() local
[all …]
|
| /petsc/src/sys/utils/ |
| H A D | pdisplay.c | 35 char *str, work[256]; in PetscOptionsGetenv() local
42 PetscCall(PetscStrncpy(work, "-", sizeof(work))); in PetscOptionsGetenv()
44 PetscCall(PetscStrlcat(work, name + 6, sizeof(work))); in PetscOptionsGetenv()
46 PetscCall(PetscStrlcat(work, name, sizeof(work))); in PetscOptionsGetenv()
48 PetscCall(PetscStrtolower(work)); in PetscOptionsGetenv()
50 PetscCall(PetscOptionsGetString(NULL, NULL, work, env, len, &flg)); in PetscOptionsGetenv()
67 PetscCall(PetscOptionsHasName(NULL, NULL, work, flag)); in PetscOptionsGetenv()
|
| H A D | mpitr.c | 29 double tsends, trecvs, work; in PetscMPIDump() local
44 work = petsc_irecv_ct + petsc_recv_ct; in PetscMPIDump()
45 PetscCallMPI(MPI_Reduce(&work, &trecvs, 1, MPI_DOUBLE, MPI_SUM, 0, PETSC_COMM_WORLD)); in PetscMPIDump()
46 work = petsc_isend_ct + petsc_send_ct; in PetscMPIDump()
47 PetscCallMPI(MPI_Reduce(&work, &tsends, 1, MPI_DOUBLE, MPI_SUM, 0, PETSC_COMM_WORLD)); in PetscMPIDump()
|
| /petsc/src/vec/vec/impls/mpi/mpiviennacl/ |
| H A D | mpiviennacl.cxx | 36 PetscReal sum, work = 0.0; in VecNorm_MPIViennaCL() local
40 PetscCall(VecNorm_SeqViennaCL(xin, NORM_2, &work)); in VecNorm_MPIViennaCL()
41 work *= work; in VecNorm_MPIViennaCL()
42 …PetscCallMPI(MPIU_Allreduce(&work, &sum, 1, MPIU_REAL, MPIU_SUM, PetscObjectComm((PetscObject)xin)… in VecNorm_MPIViennaCL()
46 PetscCall(VecNorm_SeqViennaCL(xin, NORM_1, &work)); in VecNorm_MPIViennaCL()
48 … PetscCallMPI(MPIU_Allreduce(&work, z, 1, MPIU_REAL, MPIU_SUM, PetscObjectComm((PetscObject)xin))); in VecNorm_MPIViennaCL()
51 PetscCall(VecNorm_SeqViennaCL(xin, NORM_INFINITY, &work)); in VecNorm_MPIViennaCL()
53 … PetscCallMPI(MPIU_Allreduce(&work, z, 1, MPIU_REAL, MPIU_MAX, PetscObjectComm((PetscObject)xin))); in VecNorm_MPIViennaCL()
67 PetscScalar sum, work; in VecDot_MPIViennaCL() local
70 PetscCall(VecDot_SeqViennaCL(xin, yin, &work)); in VecDot_MPIViennaCL()
[all …]
|
| /petsc/src/ts/impls/bdf/ |
| H A D | bdf.c | 18 Vec work[6 + 2]; member
114 PetscInt i, n = PETSC_STATIC_ARRAY_LENGTH(bdf->work); in TSBDF_Advance()
115 Vec tail = bdf->work[n - 1], tvtail = bdf->tvwork[n - 1]; in TSBDF_Advance()
120 bdf->work[i] = bdf->work[i - 1]; in TSBDF_Advance()
125 bdf->work[1] = tail; in TSBDF_Advance()
137 Vec *vecs = bdf->work; in TSBDF_VecLTE()
155 Vec *vecs = bdf->work + 1; in TSBDF_Extrapolate()
171 Vec *vecs = bdf->work; in TSBDF_Interpolate()
196 for (i = 1; i < n; i++) vecs[i] = bdf->transientvar ? bdf->tvwork[i] : bdf->work[i]; in TSBDF_PreSolve()
231 PetscCall(VecCopy(bdf->work[1], bdf->work[0])); in TSBDF_Restart()
[all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_ksp_py.py | 15 self.work = []
18 for v in self.work:
22 self.work[:] = ksp.getWorkVecs(right=2, left=None)
25 for v in self.work:
27 del self.work[:]
47 r, z = self.work
63 d = self.work[0].duplicate()
65 self.work += [d, q]
70 r, z, d, q = self.work
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | ij.c | 33 PetscInt *work, *ia, *ja, *j, i, nz, row, col; in MatToSymmetricIJ_SeqAIJ() local
39 PetscCall(PetscMalloc1(m + 1, &work)); in MatToSymmetricIJ_SeqAIJ()
62 work[i - 1] = row - shiftout; in MatToSymmetricIJ_SeqAIJ()
81 if (col != row) ja[work[col]++] = row + shiftout; in MatToSymmetricIJ_SeqAIJ()
82 ja[work[row]++] = col + shiftout; in MatToSymmetricIJ_SeqAIJ()
85 PetscCall(PetscFree(work)); in MatToSymmetricIJ_SeqAIJ()
|
| /petsc/src/ksp/pc/impls/svd/ |
| H A D | svd.c | 8 Vec diag, work; member
41 PetscScalar *a, *u, *v, *d, *work; in PCSetUp_SVD() local
59 PetscCall(MatCreateVecs(jac->A, &jac->diag, &jac->work)); in PCSetUp_SVD()
70 PetscCall(PetscMalloc1(lwork, &work)); in PCSetUp_SVD()
79 …PetscCallBLAS("LAPACKgesvd", LAPACKgesvd_("A", "A", &nb, &nb, a, &nb, d, u, &nb, v, &nb, work, &lw… in PCSetUp_SVD()
90 …PetscCallBLAS("LAPACKgesvd", LAPACKgesvd_("A", "A", &nb, &nb, a, &nb, dd, u, &nb, v, &nb, work, &l… in PCSetUp_SVD()
130 PetscCall(PetscFree(work)); in PCSetUp_SVD()
213 Vec work = jac->work, xred, yred; in PCApply_SVD() local
219 PetscCall(MatMultTranspose(jac->U, xred, work)); in PCApply_SVD()
221 PetscCall(MatMultHermitianTranspose(jac->U, xred, work)); in PCApply_SVD()
[all …]
|
| /petsc/src/ksp/ksp/impls/bcgs/qmrcgs/ |
| H A D | qmrcgs.c | 28 R = ksp->work[0]; in KSPSolve_QMRCGS()
29 P = ksp->work[1]; in KSPSolve_QMRCGS()
30 PH = ksp->work[2]; in KSPSolve_QMRCGS()
31 V = ksp->work[3]; in KSPSolve_QMRCGS()
32 D2 = ksp->work[4]; in KSPSolve_QMRCGS()
33 X2 = ksp->work[5]; in KSPSolve_QMRCGS()
34 S = ksp->work[6]; in KSPSolve_QMRCGS()
35 SH = ksp->work[7]; in KSPSolve_QMRCGS()
36 T = ksp->work[8]; in KSPSolve_QMRCGS()
37 D = ksp->work[9]; in KSPSolve_QMRCGS()
[all …]
|
| /petsc/src/ksp/ksp/impls/bcgs/pipebcgs/ |
| H A D | pipebcgs.c | 27 R = ksp->work[0]; in KSPSolve_PIPEBCGS()
28 RP = ksp->work[1]; in KSPSolve_PIPEBCGS()
29 S = ksp->work[2]; in KSPSolve_PIPEBCGS()
30 Y = ksp->work[3]; in KSPSolve_PIPEBCGS()
31 Q = ksp->work[4]; in KSPSolve_PIPEBCGS()
32 Q2 = ksp->work[5]; in KSPSolve_PIPEBCGS()
33 P2 = ksp->work[6]; in KSPSolve_PIPEBCGS()
34 R2 = ksp->work[7]; in KSPSolve_PIPEBCGS()
35 S2 = ksp->work[8]; in KSPSolve_PIPEBCGS()
36 W = ksp->work[9]; in KSPSolve_PIPEBCGS()
[all …]
|
| /petsc/src/mat/impls/normal/ |
| H A D | normm.c | 147 PetscScalar *work, *values; in MatGetDiagonal_Normal() local
153 PetscCall(PetscMalloc1(A->cmap->N, &work)); in MatGetDiagonal_Normal()
154 PetscCall(PetscArrayzero(work, A->cmap->N)); in MatGetDiagonal_Normal()
158 for (j = 0; j < nnz; j++) work[cols[j]] += mvalues[j] * mvalues[j]; in MatGetDiagonal_Normal()
161 …PetscCallMPI(MPIU_Allreduce(MPI_IN_PLACE, work, A->cmap->N, MPIU_SCALAR, MPIU_SUM, PetscObjectComm… in MatGetDiagonal_Normal()
165 PetscCall(PetscArraycpy(values, work + rstart, rend - rstart)); in MatGetDiagonal_Normal()
167 PetscCall(PetscFree(work)); in MatGetDiagonal_Normal()
267 Mat work[2]; member
290 PetscCall(MatProductNumeric(contents->work[0])); in MatProductNumeric_Normal_Dense()
292 PetscCall(MatDensePlaceArray(contents->work[1], array)); in MatProductNumeric_Normal_Dense()
[all …]
|
| /petsc/src/dm/dt/fe/tests/ |
| H A D | ex4.c | 54 Vec llM, lM, ones, work; in main() local
79 PetscCall(DMGetGlobalVector(dm, &work)); in main()
86 PetscCall(MatMult(M, ones, work)); in main()
87 PetscCall(VecViewFromOptions(work, NULL, "-mass_rowsum_view")); in main()
88 PetscCall(VecSum(work, &vals[3])); in main()
89 PetscCall(VecSet(work, 0.0)); in main()
90 PetscCall(DMLocalToGlobal(dm, llM, ADD_VALUES, work)); in main()
91 PetscCall(VecSum(work, &vals[4])); in main()
101 PetscCall(DMRestoreGlobalVector(dm, &work)); in main()
|
| /petsc/src/dm/field/impls/da/ |
| H A D | dmfieldda.c | 8 PetscScalar *work; member
18 PetscCall(PetscFree3(dafield->cornerVals, dafield->cornerCoeffs, dafield->work)); in DMFieldDestroy_DA()
96 PetscReal work[8]; in MultilinearEvaluate() local
103 for (k = 0; k < whol; k++) work[k] = etaB[k]; in MultilinearEvaluate()
105 etaB[k] = work[2 * k] * e; in MultilinearEvaluate()
106 etaB[k + half] = work[2 * k + 1]; in MultilinearEvaluate()
109 for (k = 0; k < whol; k++) work[k] = etaD[k]; in MultilinearEvaluate()
111 etaD[k + half] = work[2 * k]; in MultilinearEvaluate()
112 etaD[k] = work[2 * k + 1] * d; in MultilinearEvaluate()
141 for (k = 0; k < whol; k++) work[k] = etaB[k]; in MultilinearEvaluate()
[all …]
|
| /petsc/src/sys/webclient/ |
| H A D | client.c | 332 char work[256]; in PetscPullJSONValue() local
336 PetscCall(PetscStrncpy(work, "\"", sizeof(work))); in PetscPullJSONValue()
337 PetscCall(PetscStrlcat(work, key, sizeof(work))); in PetscPullJSONValue()
338 PetscCall(PetscStrlcat(work, "\":", sizeof(work))); in PetscPullJSONValue()
339 PetscCall(PetscStrstr(buff, work, &v)); in PetscPullJSONValue()
340 PetscCall(PetscStrlen(work, &len)); in PetscPullJSONValue()
344 work[len++ - 1] = 0; in PetscPullJSONValue()
345 PetscCall(PetscStrlcat(work, " :", sizeof(work))); in PetscPullJSONValue()
346 PetscCall(PetscStrstr(buff, work, &v)); in PetscPullJSONValue()
|
| /petsc/src/ksp/ksp/impls/tfqmr/ |
| H A D | tfqmr.c | 21 R = ksp->work[0]; in KSPSolve_TFQMR()
22 RP = ksp->work[1]; in KSPSolve_TFQMR()
23 V = ksp->work[2]; in KSPSolve_TFQMR()
24 T = ksp->work[3]; in KSPSolve_TFQMR()
25 Q = ksp->work[4]; in KSPSolve_TFQMR()
26 P = ksp->work[5]; in KSPSolve_TFQMR()
27 U = ksp->work[6]; in KSPSolve_TFQMR()
28 D = ksp->work[7]; in KSPSolve_TFQMR()
29 T1 = ksp->work[8]; in KSPSolve_TFQMR()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg/ |
| H A D | pipecg.c | 35 R = ksp->work[0]; in KSPSolve_PIPECG()
36 Z = ksp->work[1]; in KSPSolve_PIPECG()
37 P = ksp->work[2]; in KSPSolve_PIPECG()
38 N = ksp->work[3]; in KSPSolve_PIPECG()
39 W = ksp->work[4]; in KSPSolve_PIPECG()
40 Q = ksp->work[5]; in KSPSolve_PIPECG()
41 U = ksp->work[6]; in KSPSolve_PIPECG()
42 M = ksp->work[7]; in KSPSolve_PIPECG()
43 S = ksp->work[8]; in KSPSolve_PIPECG()
|
| /petsc/lib/petsc/bin/maint/abi-compliance-checker/ |
| H A D | LICENSE | 75 combined work, a derivative of the original library. The ordinary
111 "work based on the library" and a "work that uses the library". The
128 The "Library", below, refers to any such software library or work
129 which has been distributed under these terms. A "work based on the
130 Library" means either the Library or any derivative work under
131 copyright law: that is to say, a work containing the Library or a
136 "Source code" for a work means the preferred form of the work for
145 such a program is covered only if its contents constitute a work based
163 of it, thus forming a work based on the Library, and copy and
164 distribute such modifications or work under the terms of Section 1
[all …]
|
| /petsc/src/tao/linesearch/impls/armijo/ |
| H A D | armijo.c | 18 PetscCall(VecDestroy(&armP->work)); in TaoLineSearchDestroy_Armijo()
100 if (!armP->work) { in TaoLineSearchApply_Armijo()
101 PetscCall(VecDuplicate(x, &armP->work)); in TaoLineSearchApply_Armijo()
106 PetscCall(VecDestroy(&armP->work)); in TaoLineSearchApply_Armijo()
107 PetscCall(VecDuplicate(x, &armP->work)); in TaoLineSearchApply_Armijo()
195 PetscCall(VecWAXPY(armP->work, ls->step, s, x)); in TaoLineSearchApply_Armijo()
196 if (ls->bounded) PetscCall(VecMedian(ls->lower, armP->work, ls->upper, armP->work)); in TaoLineSearchApply_Armijo()
201 PetscCall(TaoLineSearchComputeObjective(ls, armP->work, f)); in TaoLineSearchApply_Armijo()
204 PetscCall(TaoLineSearchComputeObjectiveAndGTS(ls, armP->work, f, &gdx)); in TaoLineSearchApply_Armijo()
207 PetscCall(TaoLineSearchComputeObjectiveAndGradient(ls, armP->work, f, g)); in TaoLineSearchApply_Armijo()
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/pipeprcg/ |
| H A D | pipeprcg.c | 59 R = ksp->work[0]; in KSPSolve_PIPEPRCG()
60 RT = ksp->work[1]; in KSPSolve_PIPEPRCG()
61 W = ksp->work[2]; in KSPSolve_PIPEPRCG()
62 WT = ksp->work[3]; in KSPSolve_PIPEPRCG()
63 P = ksp->work[4]; in KSPSolve_PIPEPRCG()
64 S = ksp->work[5]; in KSPSolve_PIPEPRCG()
65 ST = ksp->work[6]; in KSPSolve_PIPEPRCG()
66 U = ksp->work[7]; in KSPSolve_PIPEPRCG()
67 UT = ksp->work[8]; in KSPSolve_PIPEPRCG()
|
| /petsc/src/ksp/ksp/impls/bcgs/fbcgs/ |
| H A D | fbcgs.c | 28 R = ksp->work[0]; in KSPSolve_FBCGS()
29 RP = ksp->work[1]; in KSPSolve_FBCGS()
30 V = ksp->work[2]; in KSPSolve_FBCGS()
31 T = ksp->work[3]; in KSPSolve_FBCGS()
32 S = ksp->work[4]; in KSPSolve_FBCGS()
33 P = ksp->work[5]; in KSPSolve_FBCGS()
34 S2 = ksp->work[6]; in KSPSolve_FBCGS()
35 P2 = ksp->work[7]; in KSPSolve_FBCGS()
|