Home
last modified time | relevance | path

Searched refs:xx (Results 1 – 25 of 186) sorted by relevance

12345678

/petsc/src/snes/tutorials/
H A Dex1.c40 PetscScalar pfive = .5, *xx; in main() local
118 PetscCall(VecGetArray(x, &xx)); in main()
119 xx[0] = 2.0; in main()
120 xx[1] = 3.0; in main()
121 PetscCall(VecRestoreArray(x, &xx)); in main()
164 const PetscScalar *xx; in FormFunction1() local
175 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction1()
179 ff[0] = xx[0] * xx[0] + xx[0] * xx[1] - 3.0; in FormFunction1()
180 ff[1] = xx[0] * xx[1] + xx[1] * xx[1] - 6.0; in FormFunction1()
183 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction1()
[all …]
H A Dex42.c23 PetscScalar *xx; in main() local
80 PetscCall(VecGetArray(x, &xx)); in main()
81 xx[0] = -1.2; in main()
82 xx[1] = 1.0; in main()
83 PetscCall(VecRestoreArray(x, &xx)); in main()
131 const PetscScalar *xx; in FormFunction1() local
146 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction1()
150 ff[0] = -2.0 + 2.0 * xx[0] + 400.0 * xx[0] * xx[0] * xx[0] - 400.0 * xx[0] * xx[1]; in FormFunction1()
151 ff[1] = -200.0 * xx[0] * xx[0] + 200.0 * xx[1]; in FormFunction1()
154 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction1()
[all …]
H A Dex99.c134 const PetscScalar *xx; in FormFunction() local
145 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
149 ff[0] = 8. * PetscExpScalar(-4. * (xx[0] - 2.) * (xx[0] - 2.)) * (xx[0] - 2.) + 2. * xx[0]; in FormFunction()
152 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction()
159 const PetscScalar *xx; in FormJacobian() local
167 PetscCall(VecGetArrayRead(x, &xx)); in FormJacobian()
174 …0] = 8. * ((xx[0] - 2.) * (PetscExpScalar(-4. * (xx[0] - 2.) * (xx[0] - 2.)) * -8. * (xx[0] - 2.))… in FormJacobian()
181 PetscCall(VecRestoreArrayRead(x, &xx)); in FormJacobian()
H A Dex59.c38 PetscScalar *xx; in main() local
102 PetscCall(VecGetArray(x, &xx)); in main()
104 v2 = xx[i] * sperturb; in main()
107 PetscCall(VecRestoreArray(x, &xx)); in main()
129 const PetscScalar *xx; in FormFunction() local
134 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
141 if (second_order) ff[0] = d * (0.5 * d * (-xx[2] + 4. * xx[1] - 3. * xx[0]) - X0DOT); in FormFunction()
142 else ff[0] = d * (d * (xx[1] - xx[0]) - X0DOT); in FormFunction()
144 …for (i = 1; i < n - 1; i++) ff[i] = d2 * (xx[i - 1] - 2. * xx[i] + xx[i + 1]) + xx[i] * xx[i] - FF… in FormFunction()
146 ff[n - 1] = d * d * (xx[n - 1] - X1); in FormFunction()
[all …]
H A Dex2.c214 const PetscScalar *xx, *gg; in FormFunction() local
226 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
236 ff[0] = xx[0]; in FormFunction()
237 …for (i = 1; i < n - 1; i++) ff[i] = d * (xx[i - 1] - 2.0 * xx[i] + xx[i + 1]) + xx[i] * xx[i] - gg… in FormFunction()
238 ff[n - 1] = xx[n - 1] - 1.0; in FormFunction()
243 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction()
265 const PetscScalar *xx; in FormJacobian() local
273 PetscCall(VecGetArrayRead(x, &xx)); in FormJacobian()
292 A[1] = -2.0 * d + 2.0 * xx[i]; in FormJacobian()
312 PetscCall(VecRestoreArrayRead(x, &xx)); in FormJacobian()
H A Dex6.c187 const PetscScalar *xx, *gg; in FormFunction() local
199 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
209 ff[0] = xx[0]; in FormFunction()
210 …for (i = 1; i < n - 1; i++) ff[i] = d * (xx[i - 1] - 2.0 * xx[i] + xx[i + 1]) + xx[i] * xx[i] - gg… in FormFunction()
211 ff[n - 1] = xx[n - 1] - 1.0; in FormFunction()
216 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction()
238 const PetscScalar *xx; in FormJacobian() local
246 PetscCall(VecGetArrayRead(x, &xx)); in FormJacobian()
265 A[1] = -2.0 * d + 2.0 * xx[i]; in FormJacobian()
285 PetscCall(VecRestoreArrayRead(x, &xx)); in FormJacobian()
/petsc/src/snes/tests/
H A Dex17.c60 PetscScalar pfive = .5, *xx; in assembled_system() local
135 PetscCall(VecGetArray(x, &xx)); in assembled_system()
136 xx[0] = 2.0; in assembled_system()
137 xx[1] = 3.0; in assembled_system()
138 PetscCall(VecRestoreArray(x, &xx)); in assembled_system()
182 const PetscScalar *xx; in FormFunction1() local
193 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction1()
199 ff[0] = xx[0] * xx[0] + xx[0] * xx[1] - 3.0; in FormFunction1()
200 ff[1] = xx[0] * xx[1] + xx[1] * xx[1] - 6.0; in FormFunction1()
205 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction1()
[all …]
H A Dex7.c106 const PetscScalar *xx, *FF; in FormFunction() local
111 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
117 ff[0] = xx[0]; in FormFunction()
118 …for (i = 1; i < n - 1; i++) ff[i] = d * (xx[i - 1] - 2.0 * xx[i] + xx[i + 1]) + xx[i] * xx[i] - FF… in FormFunction()
119 ff[n - 1] = xx[n - 1] - 1.0; in FormFunction()
120 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction()
128 const PetscScalar *xx, *FF; in FormFunctioni() local
136 PetscCall(VecGetArrayRead(x, &xx)); in FormFunctioni()
142 *s = xx[0]; in FormFunctioni()
144 *s = xx[n - 1] - 1.0; in FormFunctioni()
[all …]
H A Dex4.c42 PetscScalar *xx; in main() local
112 PetscCall(VecGetArray(x, &xx)); in main()
113 xx[0] = 2.0; in main()
114 xx[1] = 3.0; in main()
115 PetscCall(VecRestoreArray(x, &xx)); in main()
165 const PetscScalar *xx; in FormFunction2() local
176 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction2()
182 ff[0] = PetscSinScalar(3.0 * xx[0]) + xx[0]; in FormFunction2()
183 ff[1] = xx[1]; in FormFunction2()
188 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction2()
[all …]
H A Dex5.c216 const PetscScalar *xx, *gg; in FormFunction() local
228 PetscCall(VecGetArrayRead(x, &xx)); in FormFunction()
238 ff[0] = xx[0]; in FormFunction()
239 …for (i = 1; i < n - 1; i++) ff[i] = d * (xx[i - 1] - 2.0 * xx[i] + xx[i + 1]) + xx[i] * xx[i] - gg… in FormFunction()
240 ff[n - 1] = xx[n - 1] - 1.0; in FormFunction()
245 PetscCall(VecRestoreArrayRead(x, &xx)); in FormFunction()
267 const PetscScalar *xx; in FormJacobian() local
275 PetscCall(VecGetArrayRead(x, &xx)); in FormJacobian()
294 A[1] = -2.0 * d + 2.0 * xx[i]; in FormJacobian()
315 PetscCall(VecRestoreArrayRead(x, &xx)); in FormJacobian()
/petsc/src/ksp/pc/impls/pbjacobi/
H A Dpbjacobi.c12 const PetscScalar *xx; in PCApply_PBJacobi() local
15 PetscCall(VecGetArrayRead(x, &xx)); in PCApply_PBJacobi()
19 for (i = 0; i < m; i++) yy[i] = diag[i] * xx[i]; in PCApply_PBJacobi()
23 x0 = xx[2 * i]; in PCApply_PBJacobi()
24 x1 = xx[2 * i + 1]; in PCApply_PBJacobi()
32 x0 = xx[3 * i]; in PCApply_PBJacobi()
33 x1 = xx[3 * i + 1]; in PCApply_PBJacobi()
34 x2 = xx[3 * i + 2]; in PCApply_PBJacobi()
44 x0 = xx[4 * i]; in PCApply_PBJacobi()
45 x1 = xx[4 * i + 1]; in PCApply_PBJacobi()
[all …]
/petsc/src/ksp/pc/impls/vpbjacobi/
H A Dvpbjacobi.c10 const PetscScalar *xx; in PCApply_VPBJacobi() local
17 PetscCall(VecGetArrayRead(x, &xx)); in PCApply_VPBJacobi()
23 yy[ncnt] = *diag * xx[ncnt]; in PCApply_VPBJacobi()
26 x0 = xx[ncnt]; in PCApply_VPBJacobi()
27 x1 = xx[ncnt + 1]; in PCApply_VPBJacobi()
32 x0 = xx[ncnt]; in PCApply_VPBJacobi()
33 x1 = xx[ncnt + 1]; in PCApply_VPBJacobi()
34 x2 = xx[ncnt + 2]; in PCApply_VPBJacobi()
40 x0 = xx[ncnt]; in PCApply_VPBJacobi()
41 x1 = xx[ncnt + 1]; in PCApply_VPBJacobi()
[all …]
/petsc/src/sys/classes/draw/impls/x/
H A Dxtone.c14 int lc, rc = 0, lx, rx = 0, xx, y, c; in PetscDrawInterpolatedTriangle_X() local
78 for (xx = lx; xx <= rx; xx++) { in PetscDrawInterpolatedTriangle_X()
79 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscDrawInterpolatedTriangle_X()
81 XDrawPoint(win->disp, PetscDrawXiDrawable(win), win->gc.set, xx, y); in PetscDrawInterpolatedTriangle_X()
84 for (xx = lx; xx >= rx; xx--) { in PetscDrawInterpolatedTriangle_X()
85 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscDrawInterpolatedTriangle_X()
87 XDrawPoint(win->disp, PetscDrawXiDrawable(win), win->gc.set, xx, y); in PetscDrawInterpolatedTriangle_X()
128 for (xx = lx; xx <= rx; xx++) { in PetscDrawInterpolatedTriangle_X()
129 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscDrawInterpolatedTriangle_X()
131 XDrawPoint(win->disp, PetscDrawXiDrawable(win), win->gc.set, xx, y); in PetscDrawInterpolatedTriangle_X()
[all …]
/petsc/src/mat/tests/
H A Dex53.c11 Vec xx, yy, s1, s2; in main() local
39 PetscCall(VecCreate(PETSC_COMM_WORLD, &xx)); in main()
40 PetscCall(VecSetSizes(xx, m, PETSC_DECIDE)); in main()
41 PetscCall(VecSetFromOptions(xx)); in main()
42 PetscCall(VecDuplicate(xx, &s1)); in main()
43 PetscCall(VecDuplicate(xx, &s2)); in main()
44 PetscCall(VecDuplicate(xx, &yy)); in main()
63 PetscCall(VecSetRandom(xx, rand)); in main()
64 PetscCall(MatMult(A, xx, s1)); in main()
65 PetscCall(MatMult(B, xx, s2)); in main()
[all …]
H A Dex91.c13 Vec xx, s1, s2; in main() local
85 PetscCall(VecCreateSeq(PETSC_COMM_SELF, mm, &xx)); in main()
86 PetscCall(VecDuplicate(xx, &s1)); in main()
87 PetscCall(VecDuplicate(xx, &s2)); in main()
89 PetscCall(VecSetRandom(xx, rand)); in main()
90 PetscCall(MatMult(A, xx, s1)); in main()
91 PetscCall(MatMult(sA, xx, s2)); in main()
97 PetscCall(VecDestroy(&xx)); in main()
142 PetscCall(VecCreateSeq(PETSC_COMM_SELF, mm, &xx)); in main()
143 PetscCall(VecDuplicate(xx, &s1)); in main()
[all …]
H A Dex48.c8 Vec xx, s1, s2, yy; in main() local
29 PetscCall(VecCreateSeq(PETSC_COMM_SELF, M, &xx)); in main()
30 PetscCall(VecDuplicate(xx, &s1)); in main()
31 PetscCall(VecDuplicate(xx, &s2)); in main()
32 PetscCall(VecDuplicate(xx, &yy)); in main()
108 PetscCall(VecSetRandom(xx, rdm)); in main()
110 PetscCall(MatMult(A, xx, s1)); in main()
111 PetscCall(MatMultAdd(A, xx, s2, s2)); in main()
167 PetscCall(MatForwardSolve(Fact, yy, xx)); in main()
168 PetscCall(MatBackwardSolve(Fact, xx, s1)); in main()
[all …]
H A Dex51.c12 Vec xx, s1, s2; in main() local
110 PetscCall(VecCreateSeq(PETSC_COMM_SELF, mm, &xx)); in main()
111 PetscCall(VecDuplicate(xx, &s1)); in main()
112 PetscCall(VecDuplicate(xx, &s2)); in main()
114 PetscCall(VecSetRandom(xx, rdm)); in main()
115 PetscCall(MatMult(submatA[i], xx, s1)); in main()
116 PetscCall(MatMult(submatB[i], xx, s2)); in main()
122 PetscCall(VecDestroy(&xx)); in main()
133 PetscCall(VecCreateSeq(PETSC_COMM_SELF, mm, &xx)); in main()
134 PetscCall(VecDuplicate(xx, &s1)); in main()
[all …]
/petsc/src/ts/tutorials/phasefield/
H A Dpotentials.c11 PetscReal x, hx = .1 / Mx, pause, xx[3], yy[3]; in main() local
35 xx[0] = xx[1] = xx[2] = x; in main()
39 PetscCall(PetscDrawLGAddPoint(lg, xx, yy)); in main()
55 xx[0] = xx[1] = xx[2] = x; in main()
59 PetscCall(PetscDrawLGAddPoint(lg, xx, yy)); in main()
/petsc/src/ts/tutorials/
H A Dex8.py17 xx = da.createLocalVector()
19 da.globalToLocal(x,xx)
24 if xs == 0: f[0] = xx[0]/hx; lxs+=1;
25 if lxe == mx: f[mx-1] = xx[mx-1-gxs]/hx; lxe-=1
27 f[i] = xxdot[i-gxs] + (2.0*xx[i-gxs] - xx[i-1-gxs] - xx[i+1-gxs])/hx - hx*math.exp(xx[i-gxs])
32 xx = da.createLocalVector()
33 da.globalToLocal(x,xx)
41 P.setValues([i],[i-1,i,i+1],[-1.0/hx,2.0/hx-hx*math.exp(xx[i-gxs])+shift,-1.0/hx])
H A Dex1f.F9039 PetscScalar, pointer :: xx(:), ff(:)
54 PetscCall(VecGetArrayRead(X, xx, ierr))
60 ff(row) = xx(row)
62 u = xx(row)
63 ub = xx(row - mx)
64 ul = xx(row - 1)
65 ut = xx(row + mx)
66 ur = xx(row + 1)
74 PetscCall(VecRestoreArrayRead(X, xx, ierr))
92 PetscScalar, pointer :: xx(:)
[all …]
/petsc/src/sys/classes/draw/impls/image/
H A Ddrawimage.h115 int lc, rc = 0, lx, rx = 0, xx, y, c; in PetscImageDrawTriangle() local
177 for (xx = lx; xx <= rx; xx++) { in PetscImageDrawTriangle()
178 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscImageDrawTriangle()
179 PetscImageDrawPixel(img, xx, y, c); in PetscImageDrawTriangle()
182 for (xx = lx; xx >= rx; xx--) { in PetscImageDrawTriangle()
183 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscImageDrawTriangle()
184 PetscImageDrawPixel(img, xx, y, c); in PetscImageDrawTriangle()
221 for (xx = lx; xx <= rx; xx++) { in PetscImageDrawTriangle()
222 c = (((xx - lx) * (rc_lc)) / (rx_lx) + lc) >> SHIFT_VAL; in PetscImageDrawTriangle()
223 PetscImageDrawPixel(img, xx, y, c); in PetscImageDrawTriangle()
[all …]
/petsc/src/vec/vec/impls/seq/
H A Dbvec1.c102 const PetscScalar *xx; in VecAXPBY_Seq() local
105 PetscCall(VecGetArrayRead(xin, &xx)); in VecAXPBY_Seq()
108 for (PetscInt i = 0; i < n; ++i) yy[i] = a * xx[i]; in VecAXPBY_Seq()
111 for (PetscInt i = 0; i < n; ++i) yy[i] = a * xx[i] + b * yy[i]; in VecAXPBY_Seq()
114 PetscCall(VecRestoreArrayRead(xin, &xx)); in VecAXPBY_Seq()
123 const PetscScalar *yy, *xx; in VecAXPBYPCZ_Seq() local
128 PetscCall(VecGetArrayRead(xin, &xx)); in VecAXPBYPCZ_Seq()
132 for (PetscInt i = 0; i < n; ++i) zz[i] = xx[i] + beta * yy[i] + gamma * zz[i]; in VecAXPBYPCZ_Seq()
134 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i] + zz[i]; in VecAXPBYPCZ_Seq()
136 for (PetscInt i = 0; i < n; ++i) zz[i] = alpha * xx[i] + beta * yy[i]; in VecAXPBYPCZ_Seq()
[all …]
/petsc/src/ksp/ksp/tutorials/
H A Dex14f.F9095 PetscScalar, pointer ::xx(:)
107 PetscCall(VecGetArray(X, xx, ierr))
122 xx(row) = 0.0
125 xx(row) = temp1*sqrt(min((min(i, mx - i - 1))*hx, temp))
131 PetscCall(VecRestoreArray(X, xx, ierr))
154 PetscScalar, pointer ::xx(:), ff(:)
176 PetscCall(VecGetArrayRead(localX, xx, ierr))
194 ff(rowf) = xx(row)
197 u = xx(row)
198 uxx = (two*u - xx(row - 1) - xx(row + 1))*hydhx
[all …]
/petsc/src/mat/impls/baij/seq/
H A Dbaijsolvnat1.c8 PetscErrorCode MatSolve_SeqBAIJ_1_NaturalOrdering_inplace(Mat A, Vec bb, Vec xx) in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace() argument
20 PetscCall(VecGetArray(xx, &x)); in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
56 PetscCall(VecRestoreArray(xx, &x)); in MatSolve_SeqBAIJ_1_NaturalOrdering_inplace()
61 PetscErrorCode MatForwardSolve_SeqBAIJ_1_NaturalOrdering(Mat A, Vec bb, Vec xx) in MatForwardSolve_SeqBAIJ_1_NaturalOrdering() argument
74 PetscCall(VecGetArray(xx, &x)); in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
90 PetscCall(VecRestoreArray(xx, &x)); in MatForwardSolve_SeqBAIJ_1_NaturalOrdering()
94 PetscErrorCode MatBackwardSolve_SeqBAIJ_1_NaturalOrdering(Mat A, Vec bb, Vec xx) in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering() argument
107 PetscCall(VecGetArray(xx, &x)); in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering()
121 PetscCall(VecRestoreArray(xx, &x)); in MatBackwardSolve_SeqBAIJ_1_NaturalOrdering()
125 PetscErrorCode MatSolve_SeqBAIJ_1_NaturalOrdering(Mat A, Vec bb, Vec xx) in MatSolve_SeqBAIJ_1_NaturalOrdering() argument
[all …]
/petsc/src/sys/classes/draw/utils/
H A Dlg.c76 PetscReal xx; in PetscDrawLGAddPoint() local
93 xx = (PetscReal)lg->nopts; in PetscDrawLGAddPoint()
95 xx = x[i]; in PetscDrawLGAddPoint()
97 if (xx > lg->xmax) lg->xmax = xx; in PetscDrawLGAddPoint()
98 if (xx < lg->xmin) lg->xmin = xx; in PetscDrawLGAddPoint()
102 lg->x[lg->loc] = xx; in PetscDrawLGAddPoint()
130 PetscErrorCode PetscDrawLGAddPoints(PetscDrawLG lg, PetscInt n, PetscReal *xx[], PetscReal *yy[]) in PetscDrawLGAddPoints() argument
154 x = xx[j]; in PetscDrawLGAddPoints()

12345678