Home
last modified time | relevance | path

Searched refs:xt (Results 1 – 20 of 20) sorted by relevance

/petsc/src/snes/tutorials/
H A Dex4.c119 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 Dex58.c149 …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 Dcp.c64 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 Dminsurf1.c142 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 Dminsurf1.c131 …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 Dminsurf2.c170 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 Dplate2f.F9076 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 Dplate2.c207 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 Dex54f.F90309 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 Dgssecant.c15 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 Dtest_mat_dense.py120 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 Dtest_mat_py.py629 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 Dtest_mat_aij.py364 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 Dshell.c323 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 Dpower.h88 PetscScalar xt; /* -----------------------"-------------------------- */ member
/petsc/src/ts/tutorials/
H A Dex29.c222 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 Dex27.c187 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 Ddense.c110 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 Dsurfacesphere_bin.msh22 …࿦����BQ块�?""J�ݿ.C�%��߿�����-r46�?�/���V� <r�`,鿨����c���?ˊv_ɠ�i<N�Eڿ����xte��?R]��,^п�s%�4�ٿ�…
/petsc/src/ksp/pc/impls/hpddm/
H A Dpchpddm.cxx1258 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()