| /petsc/src/snes/tutorials/ |
| H A D | ex4.c | 119 PetscScalar f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormObjectiveLocal() local
130 xl = xr = xb = xt = xc; in FormObjectiveLocal()
145 xt = user->top[i + 1]; in FormObjectiveLocal()
146 } else xt = x[j + 1][i]; in FormObjectiveLocal()
150 d3 = (xc - xt); in FormObjectiveLocal()
217 PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionLocal() local
230 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionLocal()
248 xt = user->top[i + 1]; in FormFunctionLocal()
250 } else xt = x[j + 1][i]; in FormFunctionLocal()
257 d3 = (xc - xt); in FormFunctionLocal()
[all …]
|
| H A D | ex58.c | 149 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormGradient() local
181 xlt = xrb = xl = xr = xb = xt = xc; in FormGradient()
199 xt = ctx->top[i + 1]; in FormGradient()
201 } else xt = x[j + 1][i]; in FormGradient()
208 d3 = (xc - xt); in FormGradient()
213 d8 = (xt - xlt); in FormGradient()
277 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormJacobian() local
312 xlt = xrb = xl = xr = xb = xt = xc; in FormJacobian()
334 xt = ctx->top[i + 1]; in FormJacobian()
336 } else xt = x[j + 1][i]; in FormJacobian()
[all …]
|
| /petsc/src/ksp/pc/impls/cp/ |
| H A D | cp.c | 64 PetscScalar *b, *x, xt; in PCApply_CP() local
73 xt = 0.; in PCApply_CP()
74 … for (j = cp->i[i]; j < cp->i[i + 1]; j++) xt += cp->a[j] * b[cp->j[j]]; /* over rows in column */ in PCApply_CP()
75 xt *= cp->d[i]; in PCApply_CP()
76 x[i] = xt; in PCApply_CP()
77 …for (j = cp->i[i]; j < cp->i[i + 1]; j++) b[cp->j[j]] -= xt * cp->a[j]; /* over rows in column upd… in PCApply_CP()
80 xt = 0.; in PCApply_CP()
81 … for (j = cp->i[i]; j < cp->i[i + 1]; j++) xt += cp->a[j] * b[cp->j[j]]; /* over rows in column */ in PCApply_CP()
82 xt *= cp->d[i]; in PCApply_CP()
83 x[i] = xt; in PCApply_CP()
[all …]
|
| /petsc/src/tao/complementarity/tutorials/ |
| H A D | minsurf1.c | 142 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormConstraints() local
161 xlt = xrb = xl = xr = xb = xt = xc; in FormConstraints()
185 xt = user->top[i + 1]; in FormConstraints()
188 xt = x[row + mx]; in FormConstraints()
196 d3 = (xc - xt); in FormConstraints()
201 d8 = (xt - xlt); in FormConstraints()
264 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormJacobian() local
285 xlt = xrb = xl = xr = xb = xt = xc; in FormJacobian()
310 xt = user->top[i + 1]; in FormJacobian()
313 xt = x[row + mx]; in FormJacobian()
[all …]
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | minsurf1.c | 131 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormFunctionGradient() local
148 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient()
171 xt = user->top[i + 1]; in FormFunctionGradient()
174 xt = x[row + mx]; in FormFunctionGradient()
182 d3 = (xc - xt); in FormFunctionGradient()
187 d8 = (xt - xlt); in FormFunctionGradient()
310 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in QuadraticH() local
331 xlt = xrb = xl = xr = xb = xt = xc; in QuadraticH()
356 xt = user->top[i + 1]; in QuadraticH()
359 xt = x[row + mx]; in QuadraticH()
[all …]
|
| H A D | minsurf2.c | 170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local
211 xt = user->top[i - xs + 1]; in FormFunction()
213 xt = x[j + 1][i]; in FormFunction()
218 d3 = (xc - xt); in FormFunction()
304 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionGradient() local
327 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient()
351 xt = user->top[i - xs + 1]; in FormFunctionGradient()
354 xt = x[j + 1][i]; in FormFunctionGradient()
362 d3 = (xc - xt); in FormFunctionGradient()
367 d8 = (xt - xlt); in FormFunctionGradient()
[all …]
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 76 PetscReal xc, xl, xr, xt, xb, xlt, xrb
116 xt = xc
145 xt = top_v(1 + i - xs + 1)
148 xt = x_v(1 + row + gxm)
161 d3 = xc - xt
166 d8 = xt - xlt
308 PetscReal xc, xl, xr, xt, xb, xlt, xrb
355 xt = xc
384 xt = top_v(1 + i - xs + 1)
387 xt = x_v(1 + row + gxm)
[all …]
|
| H A D | plate2.c | 207 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionGradient() local
238 xlt = xrb = xl = xr = xb = xt = xc; in FormFunctionGradient()
262 xt = top[i - xs + 1]; in FormFunctionGradient()
265 xt = x[row + gxm]; in FormFunctionGradient()
273 d3 = (xc - xt); in FormFunctionGradient()
278 d8 = (xt - xlt); in FormFunctionGradient()
423 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormHessian() local
460 xlt = xrb = xl = xr = xb = xt = xc; in FormHessian()
485 xt = top[i - xs + 1]; in FormHessian()
488 xt = x[row + gxm]; in FormHessian()
[all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex54f.F90 | 309 PetscReal xo, xs, xt, yo, ys, yt, xsm, xsp, xtm
331 xt = -xl(1, 1) - xl(1, 2) + xl(1, 3) + xl(1, 4) + xo*s
338 xsj1 = xs*yt - xt*ys
352 xt = (xt + xt)*xsj1
362 xtm = xt*tm
364 xtp = xt*tp
|
| /petsc/src/snes/impls/gs/ |
| H A D | gssecant.c | 15 PetscReal dxt, xt, ft, ft1 = 0; in SNESComputeNGSDefaultSecant() local
74 xt = 0.; in SNESComputeNGSDefaultSecant()
92 xt += PetscRealPart(PetscSqr(x)); in SNESComputeNGSDefaultSecant()
104 if (stol * PetscSqrtReal(xt) > PetscSqrtReal(dxt)) isdone = PETSC_TRUE; in SNESComputeNGSDefaultSecant()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_mat_dense.py | 120 xt, yt = AT.createVecs()
124 y.copy(xt)
125 AT.mult(xt, yt)
131 AT.multTranspose(yt, xt)
132 self.assertTrue(xt.equal(y))
|
| H A D | test_mat_py.py | 629 xt, yt = AT.createVecs()
633 y.copy(xt)
634 AT.mult(xt, yt)
640 AT.multTranspose(yt, xt)
641 self.assertTrue(xt.equal(y))
|
| H A D | test_mat_aij.py | 364 xt, yt = AT.createVecs()
368 y.copy(xt)
369 AT.mult(xt, yt)
375 AT.multTranspose(yt, xt)
376 self.assertTrue(xt.equal(y))
|
| /petsc/src/mat/impls/shell/ |
| H A D | shell.c | 323 Vec xt; in MatZeroRows_Shell() local
331 PetscCall(MatCreateVecs(mat, &xt, NULL)); in MatZeroRows_Shell()
332 PetscCall(VecCopy(x, xt)); in MatZeroRows_Shell()
334 PetscCall(VecSetValues(xt, nc, gcols, vals, INSERT_VALUES)); /* xt = [x1, 0] */ in MatZeroRows_Shell()
336 PetscCall(VecAssemblyBegin(xt)); in MatZeroRows_Shell()
337 PetscCall(VecAssemblyEnd(xt)); in MatZeroRows_Shell()
338 PetscCall(VecAYPX(xt, -1.0, x)); /* xt = [0, x2] */ in MatZeroRows_Shell()
340 PetscCall(VecGetOwnershipRange(xt, &st, NULL)); in MatZeroRows_Shell()
341 PetscCall(VecGetLocalSize(xt, &nl)); in MatZeroRows_Shell()
342 PetscCall(VecGetArray(xt, &vals)); in MatZeroRows_Shell()
[all …]
|
| /petsc/src/snes/tutorials/network/power/ |
| H A D | power.h | 88 PetscScalar xt; /* -----------------------"-------------------------- */ member
|
| /petsc/src/ts/tutorials/ |
| H A D | ex29.c | 222 PetscErrorCode FormIFunctionLocal(DMDALocalInfo *info, PetscReal ptime, Field **x, Field **xt, Fiel… in FormIFunctionLocal() argument
249 f[j][i].u = xt[j][i].u * scale; in FormIFunctionLocal()
|
| H A D | ex27.c | 187 PetscErrorCode FormIFunctionLocal(DMDALocalInfo *info, PetscScalar ptime, Field **x, Field **xt, Fi… in FormIFunctionLocal() argument
245 f[j][i].sp[l] = xt[j][i].sp[l] * scale; in FormIFunctionLocal()
|
| /petsc/src/mat/impls/dense/seq/ |
| H A D | dense.c | 110 Vec xt; in MatZeroRowsColumns_SeqDense() local
113 PetscCall(VecDuplicate(x, &xt)); in MatZeroRowsColumns_SeqDense()
114 PetscCall(VecCopy(x, xt)); in MatZeroRowsColumns_SeqDense()
115 PetscCall(VecScale(xt, -1.0)); in MatZeroRowsColumns_SeqDense()
116 PetscCall(MatMultAdd(A, xt, b, b)); in MatZeroRowsColumns_SeqDense()
117 PetscCall(VecDestroy(&xt)); in MatZeroRowsColumns_SeqDense()
1007 PetscScalar *x, *v = mat->v, zero = 0.0, xt; in MatSOR_SeqDense() local
1029 PetscCallBLAS("BLASdotu", xt = b[i] - BLASdotu_(&bm, v + i, &bm, x, &o)); in MatSOR_SeqDense()
1030 x[i] = (1. - omega) * x[i] + (xt + v[i + i * m] * x[i]) * omega / (v[i + i * m] + shift); in MatSOR_SeqDense()
1035 PetscCallBLAS("BLASdotu", xt = b[i] - BLASdotu_(&bm, v + i, &bm, x, &o)); in MatSOR_SeqDense()
[all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | surfacesphere_bin.msh | 22 …����BQ块�?""J�ݿ.C�%��߿�����-r46�?�/���V� <r�`,鿨����c���?ˊv_ɠ�i<N�Eڿ����xte��?R]��,^п�s%�4�ٿ�…
|
| /petsc/src/ksp/pc/impls/hpddm/ |
| H A D | pchpddm.cxx | 1258 Vec x, xt = nullptr; in PCHPDDMSetUpNeumannOverlap_Private() local
1263 …PetscCallBack("PCHPDDM Neumann callback", (*data->setup)(data->aux, t, x, xt, s, data->is, data->s… in PCHPDDMSetUpNeumannOverlap_Private()
|