Home
last modified time | relevance | path

Searched refs:j (Results 1 – 25 of 1093) sorted by relevance

12345678910>>...44

/petsc/src/ksp/ksp/impls/cg/pipecg2/
H A Dpipecg2.c12 PetscInt j, n; in VecMergedDot_Private() local
22 for (j = 0; j < n; j++) { in VecMergedDot_Private()
23 sumwu += PW[j] * PetscConj(PU[j]); in VecMergedDot_Private()
24 sumru += PR[j] * PetscConj(PU[j]); in VecMergedDot_Private()
25 sumuu += PU[j] * PetscConj(PU[j]); in VecMergedDot_Private()
29 for (j = 0; j < n; j++) { in VecMergedDot_Private()
30 sumwu += PW[j] * PetscConj(PU[j]); in VecMergedDot_Private()
31 sumru += PR[j] * PetscConj(PU[j]); in VecMergedDot_Private()
32 sumuu += PR[j] * PetscConj(PR[j]); in VecMergedDot_Private()
36 for (j = 0; j < n; j++) { in VecMergedDot_Private()
[all …]
/petsc/src/benchmarks/streams/
H A DSSEVersion.c76 int BytesPerWord, j, k, size; in main() local
118 for (j = 0; j < N; j++) { in main()
119 a[j] = 1.0; in main()
120 b[j] = 2.0; in main()
121 c[j] = 0.0; in main()
124 for (j = 0; j < N; j++) a[j] = 1.0; in main()
125 for (j = 0; j < N; j++) b[j] = 2.0; in main()
126 for (j = 0; j < N; j++) c[j] = 0.0; in main()
135 for (j = 0; j < N; j++) a[j] = 2.0E0 * a[j]; in main()
161 for (j = 0; j < N; j += 8) { in main()
[all …]
H A DBasicVersion.c91 register int j, k; in main() local
96 for (j = 0; j < N; j++) { in main()
97 a[j] = 1.0; in main()
98 b[j] = 2.0; in main()
99 c[j] = 0.0; in main()
103 for (j = 0; j < N; j++) a[j] = 2.0E0 * a[j]; in main()
113 for (j = 0; j < N; j++) c[j] = a[j]; in main()
118 for (j = 0; j < N; j++) b[j] = scalar * c[j]; in main()
122 for (j = 0; j < N; j++) c[j] = a[j] + b[j]; in main()
126 for (j = 0; j < N; j++) a[j] = b[j] + scalar * c[j]; in main()
[all …]
/petsc/src/ts/tutorials/multirate/
H A Dex8.c126 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample_3WaySplit() local
145 for (j = 0; j < N + 1; j++) { in FVSample_3WaySplit()
146 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample_3WaySplit()
148 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit()
154 for (j = 0; j < N + 1; j++) { in FVSample_3WaySplit()
155 xj = xi + hm * (j - N / 2) / (PetscReal)N; in FVSample_3WaySplit()
157 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit()
163 for (j = 0; j < N + 1; j++) { in FVSample_3WaySplit()
164 xj = xi + hf * (j - N / 2) / (PetscReal)N; in FVSample_3WaySplit()
166 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit()
[all …]
H A Dex7.c172 PetscInt i, j, Mx, dof, xs, xm, sf = ctx->sf, fs = ctx->fs; in FVRHSFunction() local
194 for (j = 0; j < dof; j++) x[i * dof + j] = x[j]; in FVRHSFunction()
197 for (j = 0; j < dof; j++) x[i * dof + j] = x[(xs + xm - 1) * dof + j]; in FVRHSFunction()
208 for (j = 0; j < dof; j++) { in FVRHSFunction()
209 … r[j] = (x[i * dof + j] - x[(i - 1) * dof + j]) / (x[(i - 1) * dof + j] - x[(i - 2) * dof + j]); in FVRHSFunction()
210 min[j] = PetscMin(r[j], 2.0); in FVRHSFunction()
211 …u[j] = x[(i - 1) * dof + j] + PetscMax(0, PetscMin(min[j], alpha[0] + gamma[0] * r[j])) * (x[(i … in FVRHSFunction()
216 for (j = 0; j < dof; j++) f[(i - 1) * dof + j] -= ctx->flux[j] / hs; in FVRHSFunction()
219 for (j = 0; j < dof; j++) f[i * dof + j] += ctx->flux[j] / hs; in FVRHSFunction()
225 for (j = 0; j < dof; j++) { in FVRHSFunction()
[all …]
H A Dex6.c133 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample_2WaySplit() local
150 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
151 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample_2WaySplit()
153 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit()
159 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
160 xj = xi + hf * (j - N / 2) / (PetscReal)N; in FVSample_2WaySplit()
162 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit()
168 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
169 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample_2WaySplit()
171 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit()
[all …]
H A Dex4.c263 PetscInt i, j; in PhysicsCharacteristic_Conservative() local
267 for (j = 0; j < m; j++) Xi[i * m + j] = X[i * m + j] = (PetscScalar)(i == j); in PhysicsCharacteristic_Conservative()
477 PetscInt i, j, k, dof, xs, xm, Mx; in FVSample_2WaySplit() local
494 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
495 xj = xi + hs * (j - N / 2) / (PetscReal)N; in FVSample_2WaySplit()
497 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit()
503 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
504 xj = xi + hf * (j - N / 2) / (PetscReal)N; in FVSample_2WaySplit()
506 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit()
512 for (j = 0; j < N + 1; j++) { in FVSample_2WaySplit()
[all …]
H A Dex5.c190 PetscInt i, j, k, Mx, dof, xs, xm, len_slow; in FVRHSFunctionslow() local
212 for (j = 0; j < dof; j++) x[i * dof + j] = x[j]; in FVRHSFunctionslow()
215 for (j = 0; j < dof; j++) x[i * dof + j] = x[(xs + xm - 1) * dof + j]; in FVRHSFunctionslow()
228 for (j = 0; j < dof; j++) { in FVRHSFunctionslow()
230 jmpL = x[(i + 0) * dof + j] - x[(i - 1) * dof + j]; in FVRHSFunctionslow()
231 jmpR = x[(i + 1) * dof + j] - x[(i + 0) * dof + j]; in FVRHSFunctionslow()
233 cjmpL[k] += ctx->Rinv[k + j * dof] * jmpL; in FVRHSFunctionslow()
234 cjmpR[k] += ctx->Rinv[k + j * dof] * jmpR; in FVRHSFunctionslow()
241 for (j = 0; j < dof; j++) ctx->cslope[j] /= hx; /* rescale to a slope */ in FVRHSFunctionslow()
242 for (j = 0; j < dof; j++) { in FVRHSFunctionslow()
[all …]
/petsc/src/snes/tutorials/
H A Dex19.c187 PetscInt i, j, mx, xs, ys, xm, ym; in FormInitialGuess() local
217 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
219 x[j][i].u = 0.0; in FormInitialGuess()
220 x[j][i].v = 0.0; in FormInitialGuess()
221 x[j][i].omega = 0.0; in FormInitialGuess()
222 x[j][i].temp = (grashof > 0) * i * dx; in FormInitialGuess()
236 PetscInt xints, xinte, yints, yinte, i, j; in FormFunctionLocal() local
267 j = 0; in FormFunctionLocal()
271 f[j][i].u = x[j][i].u; in FormFunctionLocal()
272 f[j][i].v = x[j][i].v; in FormFunctionLocal()
[all …]
H A Dex30.c285 static inline PetscScalar UInterp(Field **x, PetscInt i, PetscInt j) in UInterp() argument
287 return 0.25 * (x[j][i].u + x[j + 1][i].u + x[j][i + 1].u + x[j + 1][i + 1].u); in UInterp()
290 static inline PetscScalar WInterp(Field **x, PetscInt i, PetscInt j) in WInterp() argument
292 return 0.25 * (x[j][i].w + x[j + 1][i].w + x[j][i + 1].w + x[j + 1][i + 1].w); in WInterp()
295 static inline PetscScalar PInterp(Field **x, PetscInt i, PetscInt j) in PInterp() argument
297 return 0.25 * (x[j][i].p + x[j + 1][i].p + x[j][i + 1].p + x[j + 1][i + 1].p); in PInterp()
300 static inline PetscScalar TInterp(Field **x, PetscInt i, PetscInt j) in TInterp() argument
302 return 0.25 * (x[j][i].T + x[j + 1][i].T + x[j][i + 1].T + x[j + 1][i + 1].T); in TInterp()
306 static inline PetscScalar HorizVelocity(PetscInt i, PetscInt j, AppCtx *user) in HorizVelocity() argument
314 z = (j - grid->jlid - 0.5) * grid->dz; in HorizVelocity()
[all …]
H A Dex15.c286 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormInitialGuess() local
317 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
318 temp = (PetscReal)(PetscMin(j, My - j - 1)) * hy; in FormInitialGuess()
320 if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) { in FormInitialGuess()
322 x[j][i] = 0.0; in FormInitialGuess()
326 … x[j][i] = temp1 * PetscSqrtReal(PetscMin((PetscReal)(PetscMin(i, Mx - i - 1)) * hx, temp)); in FormInitialGuess()
330 … const PetscReal xx = 2 * (PetscReal)i / (Mx - 1) - 1, yy = 2 * (PetscReal)j / (My - 1) - 1; in FormInitialGuess()
331 x[j][i] = (1 - xx * xx) * (1 - yy * yy) * xx * yy; in FormInitialGuess()
334 x[j][i] = 0.; in FormInitialGuess()
336 … const PetscReal xx = 2 * (PetscReal)i / (Mx - 1) - 1, yy = 2 * (PetscReal)j / (My - 1) - 1; in FormInitialGuess()
[all …]
H A Dex18.c111 PetscInt i, j, xs, ys, xm, ym; in FormInitialGuess() local
125 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
126 for (i = xs; i < xs + xm; i++) x[j][i] = tleft; in FormInitialGuess()
135 PetscInt i, j, mx, my, xs, ys, xm, ym; in FormFunction() local
164 for (j = ys; j < ys + ym; j++) { in FormFunction()
166 t0 = x[j][i]; in FormFunction()
168 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1) { in FormFunction()
171 tw = x[j][i - 1]; in FormFunction()
176 te = x[j][i + 1]; in FormFunction()
181 ts = x[j - 1][i]; in FormFunction()
[all …]
/petsc/src/ksp/ksp/impls/gmres/
H A Dborthog2.c38 PetscInt j; in KSPGMRESClassicalGramSchmidtOrthogonalization() local
53 for (j = 0; j <= it; j++) { in KSPGMRESClassicalGramSchmidtOrthogonalization()
54 hh[j] = 0.0; in KSPGMRESClassicalGramSchmidtOrthogonalization()
55 hes[j] = 0.0; in KSPGMRESClassicalGramSchmidtOrthogonalization()
63 for (j = 0; j <= it; j++) { in KSPGMRESClassicalGramSchmidtOrthogonalization()
64 KSPCheckDot(ksp, lhh[j]); in KSPGMRESClassicalGramSchmidtOrthogonalization()
66 lhh[j] = -lhh[j]; in KSPGMRESClassicalGramSchmidtOrthogonalization()
75 for (j = 0; j <= it; j++) { in KSPGMRESClassicalGramSchmidtOrthogonalization()
76 hh[j] -= lhh[j]; /* hh += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
77 hes[j] -= lhh[j]; /* hes += <v,vnew> */ in KSPGMRESClassicalGramSchmidtOrthogonalization()
[all …]
/petsc/src/ts/tutorials/
H A Dex26.c211 PetscInt i, j, mx, xs, ys, xm, ym; in FormInitialSolution() local
241 for (j = ys; j < ys + ym; j++) { in FormInitialSolution()
243 x[j][i].u = 0.0; in FormInitialSolution()
244 x[j][i].v = 0.0; in FormInitialSolution()
245 x[j][i].omega = 0.0; in FormInitialSolution()
246 x[j][i].temp = (grashof > 0) * i * dx; in FormInitialSolution()
260 PetscInt xints, xinte, yints, yinte, i, j; in FormIFunctionLocal() local
291 j = 0; in FormIFunctionLocal()
295 f[j][i].u = x[j][i].u; in FormIFunctionLocal()
296 f[j][i].v = x[j][i].v; in FormIFunctionLocal()
[all …]
H A Dex15.c159 PetscInt i, j, Mx, My, xs, ys, xm, ym; in FormIFunction() local
192 for (j = ys; j < ys + ym; j++) { in FormIFunction()
195 if (i == 0 || j == 0 || i == Mx - 1 || j == My - 1) { in FormIFunction()
197 f[j][i] = uarray[j][i]; /* F = U */ in FormIFunction()
199 if (i == 0 && j == 0) { /* SW corner */ in FormIFunction()
200 f[j][i] = uarray[j][i] - uarray[j + 1][i + 1]; in FormIFunction()
201 } else if (i == Mx - 1 && j == 0) { /* SE corner */ in FormIFunction()
202 f[j][i] = uarray[j][i] - uarray[j + 1][i - 1]; in FormIFunction()
203 } else if (i == 0 && j == My - 1) { /* NW corner */ in FormIFunction()
204 f[j][i] = uarray[j][i] - uarray[j - 1][i + 1]; in FormIFunction()
[all …]
/petsc/src/ksp/ksp/tests/
H A Dex10.c85 PetscInt i, j, k, i1, i2, j_1, j2, k1, k2, h1, h2, shiftx, shifty, shiftz; in GetElasticityMatrix() local
107 for (j = 0; j < m; j++) { in GetElasticityMatrix()
110 base = 2 * k * shiftz + 2 * j * shifty + 2 * i * shiftx; in GetElasticityMatrix()
212 PetscInt i, j, k, l, Ii, J; in Elastic20Stiff() local
286 for (j = 0; j < 81; j++) Ke[i][j] = 0.0; in Elastic20Stiff()
291 for (j = 0; j < 20; j++) { in Elastic20Stiff()
293 for (l = 0; l < 3; l++) Ke[3 * rmap[i] + k][3 * rmap[j] + l] = K[Ii + k][J + l]; in Elastic20Stiff()
302 for (j = 0; j < i; j++) Ke[i][j] = (Ke[i][j] + Ke[j][i]) / 2.0; in Elastic20Stiff()
312 PetscInt i, j, k, cnt; in paulsetup20() local
351 for (j = 0; j < n_int; j++) { in paulsetup20()
[all …]
/petsc/src/dm/impls/plex/
H A Dgmshlex.h14 #define SI2(p, i, j) ((i) + (SN2(p) - SN2((p) - (j)))) argument
15 #define SI3(p, i, j, k) (SI2((p) - (k), i, j) + (SN3(p) - SN3((p) - (k)))) argument
17 #define SL2(p, i, j) SL1((p) - 1, i) SL1((p) - (i), j) argument
18 #define SL3(p, i, j, k) SL1((p) - 2, i) SL1((p) - (i), j) SL1((p) - (i) - (j), k) argument
24 #define BI2(p, i, j) ((i) + (j) * BN1(p)) argument
25 #define BI3(p, i, j, k) ((i) + BI2(p, j, k) * BN1(p)) argument
27 #define BL2(p, i, j) BL1(p, i) BL1(p, j) argument
28 #define BL3(p, i, j, k) BL1(p, i) BL1(p, j) BL1(p, k) argument
70 #define loop2(i, j) SL2(p, i, j) in GmshLexOrder_TRI() argument
71 #define index(i, j) SI2(p, i, j) in GmshLexOrder_TRI() argument
[all …]
/petsc/src/snes/tests/
H A Dex20.c112 PetscInt i, j, k, xs, ys, xm, ym, zs, zm; in FormInitialGuess() local
124 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
125 for (i = xs; i < xs + xm; i++) x[k][j][i] = user->tleft; in FormInitialGuess()
135 PetscInt i, j, k, mx, my, mz, xs, ys, zs, xm, ym, zm; in FormFunction() local
167 for (j = ys; j < ys + ym; j++) { in FormFunction()
169 t0 = x[k][j][i]; in FormFunction()
171 if (i > 0 && i < mx - 1 && j > 0 && j < my - 1 && k > 0 && k < mz - 1) { in FormFunction()
174 tw = x[k][j][i - 1]; in FormFunction()
179 te = x[k][j][i + 1]; in FormFunction()
184 ts = x[k][j - 1][i]; in FormFunction()
[all …]
H A Dex69.c162 PetscInt i, j, mx, xs, ys, xm, ym; in FormInitialGuess() local
192 for (j = ys; j < ys + ym; j++) { in FormInitialGuess()
194 x[j][i].u = 0.0; in FormInitialGuess()
195 x[j][i].v = 0.0; in FormInitialGuess()
196 x[j][i].omega = 0.0; in FormInitialGuess()
197 x[j][i].temp = (grashof > 0) * i * dx; in FormInitialGuess()
211 PetscInt xints, xinte, yints, yinte, i, j; in FormFunctionLocal() local
248 j = 0; in FormFunctionLocal()
252 f[j][i].u = x[j][i].u; in FormFunctionLocal()
253 f[j][i].v = x[j][i].v; in FormFunctionLocal()
[all …]
/petsc/src/ts/tutorials/power_grid/
H A Dex6.c129 PetscInt i, j; in ini_bou() local
148 for (j = ys; j < ys + ym; j++) { in ini_bou()
149 xi = coors[j][i].x; in ini_bou()
150 yi = coors[j][i].y; in ini_bou()
151 if (i == 0 || j == 0 || i == M - 1 || j == N - 1) p[j][i] = 0.0; in ini_bou()
153 …p[j][i] = (0.5 / (PETSC_PI * sigmax * sigmay * PetscSqrtScalar(1.0 - rho * rho))) * PetscExpScalar… in ini_bou()
164 PetscErrorCode adv1(PetscScalar **p, PetscScalar y, PetscInt i, PetscInt j, PetscInt M, PetscScalar… in adv1() argument
183 *p1 = f * (p[j][i + 1] - p[j][i - 1]) / (2 * user->dx); in adv1()
188 PetscErrorCode adv2(PetscScalar **p, PetscScalar x, PetscInt i, PetscInt j, PetscInt N, PetscScalar… in adv2() argument
207 *p2 = f * (p[j + 1][i] - p[j - 1][i]) / (2 * user->dy); in adv2()
[all …]
/petsc/src/ksp/pc/impls/cp/
H A Dcp.c12 PetscInt *i, *j; /* offsets of nonzeros by column, non-zero indices by column */ member
18 PetscInt i, j, *colcnt; in PCSetUp_CP() local
32 PetscCall(PetscFree3(cp->a, cp->i, cp->j)); in PCSetUp_CP()
37 if (!cp->a) PetscCall(PetscMalloc3(aij->nz, &cp->a, cp->n + 1, &cp->i, aij->nz, &cp->j)); in PCSetUp_CP()
40 for (i = 0; i < aij->nz; i++) colcnt[aij->j[i]]++; in PCSetUp_CP()
45 for (j = aij->i[i]; j < aij->i[i + 1]; j++) { /* over columns in row */ in PCSetUp_CP()
46 cp->j[cp->i[aij->j[j]] + colcnt[aij->j[j]]] = i; in PCSetUp_CP()
47 cp->a[cp->i[aij->j[j]] + colcnt[aij->j[j]]++] = aij->a[j]; in PCSetUp_CP()
55 …for (j = cp->i[i]; j < cp->i[i + 1]; j++) cp->d[i] += cp->a[j] * cp->a[j]; /* over rows in column … in PCSetUp_CP()
65 PetscInt i, j; in PCApply_CP() local
[all …]
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresleja.c26 PetscInt i, j; in KSPAGMRESLejaCfpdMax() local
33 for (j = 0; j < nbre; j++) { in KSPAGMRESLejaCfpdMax()
34 rd = rm[i] - rm[spos[j]]; in KSPAGMRESLejaCfpdMax()
35 id = im[i] - im[spos[j]]; in KSPAGMRESLejaCfpdMax()
51 PetscInt i, pos, j; in KSPAGMRESLejaOrdering() local
57 j = 0; in KSPAGMRESLejaOrdering()
58 while (j < m) { in KSPAGMRESLejaOrdering()
59 if (im[j] != 0.0) { /* complex eigenvalue */ in KSPAGMRESLejaOrdering()
60 if (im[j] < 0.0) { /* change the order */ in KSPAGMRESLejaOrdering()
61 temp = im[j + 1]; in KSPAGMRESLejaOrdering()
[all …]
/petsc/src/vec/is/is/tests/
H A Dex13.c13 PetscInt *idx, i, j; in CreateIS() local
22 for (i = 0, j = first; i < n; i++, j += step) idx[i] = j; in CreateIS()
30 for (i = 0, j = first; i < n; i++, j += step) idx[i] = j; in CreateIS()
42 PetscInt i, j = 0, type; in main() local
57 PetscCall(CreateIS(comm, type, n, first + offset, step, &is[j])); in main()
58 j++; in main()
60 PetscCall(CreateIS(comm, type, n, first + offset, step, &is[j])); in main()
61 PetscCall(ISCopy(is[j], is[j])); in main()
62 j++; in main()
65 PetscCall(ISDuplicate(tmp, &is[j])); in main()
[all …]
/petsc/src/dm/tests/
H A Dex10k.kokkos.cxx60 for (PetscInt j = ys; j < ys + ym; j++) { /* run on host */ in main() local
62 for (PetscInt c = 0; c < dof; c++) garray[j][i][c] = 100 * j + 10 * (i + 1) + c; in main()
71 KOKKOS_LAMBDA(PetscInt j, PetscInt i, PetscInt c) /* might run on device */ in main()
72 { gview(j, i, c) = 100 * j + 10 * (i + 1) + c; }); in main()
82 for (PetscInt j = ys; j < ys + ym; j++) { in main() local
84 …cInt c = 0; c < dof; c++) garray[j][i][c] = (larray[j][i - 1][c] + larray[j][i + 1][c] + larray[j … in main()
95 …S_LAMBDA(PetscInt j, PetscInt i, PetscInt c) { gview(j, i, c) = (lview(j, i - 1, c) + lview(j, i +… in main()
133 for (PetscInt j = ys; j < ys + ym; j++) { in main() local
135 garray2[j][i].u = 100 * j + 10 * (i + 1) + 111; in main()
136 garray2[j][i].v = 100 * j + 10 * (i + 1) + 222; in main()
[all …]
/petsc/src/tao/leastsquares/impls/pounders/
H A Dgqt.c7 PetscInt i, j; in estsv() local
42 for (j = i + 1; j < n; j++) sm += PetscAbs(z[j] + wm * r[i + ldr * j]); in estsv()
59 for (j = n - 1; j >= 0; j--) { in estsv()
61 if (PetscAbs(z[j]) > PetscAbs(r[j + ldr * j])) { in estsv()
62 temp = PetscMin(0.01, PetscAbs(r[j + ldr * j] / z[j])); in estsv()
66 if (r[j + ldr * j] == 0) { in estsv()
67 z[j] = 1.0; in estsv()
69 z[j] = z[j] / r[j + ldr * j]; in estsv()
71 temp = -z[j]; in estsv()
72 PetscCall(PetscBLASIntCast(j, &blasj)); in estsv()
[all …]

12345678910>>...44