Home
last modified time | relevance | path

Searched refs:tt (Results 1 – 25 of 41) sorted by relevance

12

/petsc/src/ksp/pc/impls/gamg/
H A Dgeo.c129 PetscInt jj, tid, tt, idx, nselected_2; in triangulateAndFormProl() local
308 for (tt = 0; tt < 3; tt++) { in triangulateAndFormProl()
309 PetscInt cid2 = mid.trianglelist[3 * tid + tt]; in triangulateAndFormProl()
311 AA[tt][0] = coords[lid2]; in triangulateAndFormProl()
312 AA[tt][1] = coords[data_stride + lid2]; in triangulateAndFormProl()
313 AA[tt][2] = 1.0; in triangulateAndFormProl()
314 clids[tt] = cid2; /* store for interp */ in triangulateAndFormProl()
317 for (tt = 0; tt < 3; tt++) alpha[tt] = (PetscScalar)fcoord[tt]; in triangulateAndFormProl()
324 for (tt = 0, idx = 0; tt < 3; tt++) { in triangulateAndFormProl()
325 … if (PetscRealPart(alpha[tt]) > (1.0 + EPS) || PetscRealPart(alpha[tt]) < -EPS) have = PETSC_FALSE; in triangulateAndFormProl()
[all …]
/petsc/src/ksp/ksp/impls/gmres/fgmres/
H A Dfgmres.c60 PetscReal hapbnd, tt; in KSPFGMRESCycle() local
134 PetscCall(VecNorm(VEC_VV(loc_it + 1), NORM_2, &tt)); in KSPFGMRESCycle()
135 KSPCheckNorm(ksp, tt); in KSPFGMRESCycle()
137 *HH(loc_it + 1, loc_it) = tt; in KSPFGMRESCycle()
138 *HES(loc_it + 1, loc_it) = tt; in KSPFGMRESCycle()
141 hapbnd = PetscAbsScalar((tt) / *RS(loc_it)); in KSPFGMRESCycle()
144 if (tt > hapbnd) { in KSPFGMRESCycle()
146 PetscCall(VecScale(VEC_VV(loc_it + 1), 1.0 / tt)); in KSPFGMRESCycle()
259 PetscScalar tt; in KSPFGMRESBuildSoln() local
283 tt = *RS(k); in KSPFGMRESBuildSoln()
[all …]
/petsc/src/ksp/ksp/impls/gmres/pipefgmres/
H A Dpipefgmres.c51 PetscReal hapbnd, tt; in KSPPIPEFGMRESCycle() local
165 tt = PetscRealPart(lhh[loc_it + 1]); in KSPPIPEFGMRESCycle()
184 for (k = 0; k <= loc_it; k++) tt -= PetscAbsScalar(lhh[k]) * PetscAbsScalar(lhh[k]); in KSPPIPEFGMRESCycle()
186 if ((tt < 0.0) && tt > -PETSC_SMALL) tt = 0.0; in KSPPIPEFGMRESCycle()
187 if (tt < 0.0) { in KSPPIPEFGMRESCycle()
194 … "Restart due to square root breakdown at it = %" PetscInt_FMT ", tt=%g\n", ksp->its, (double)tt)); in KSPPIPEFGMRESCycle()
197 tt = PetscSqrtReal(tt); in KSPPIPEFGMRESCycle()
201 hh[loc_it + 1] = tt; in KSPPIPEFGMRESCycle()
202 hes[loc_it + 1] = tt; in KSPPIPEFGMRESCycle()
215 if (tt) PetscCall(VecScale(PREVEC(loc_it + 1), 1.0 / tt)); in KSPPIPEFGMRESCycle()
[all …]
/petsc/src/ksp/ksp/impls/gmres/lgmres/
H A Dlgmres.c92 PetscReal hapbnd, tt; in KSPLGMRESCycle() local
189 PetscCall(VecNorm(VEC_VV(loc_it + 1), NORM_2, &tt)); in KSPLGMRESCycle()
191 *HH(loc_it + 1, loc_it) = tt; in KSPLGMRESCycle()
192 *HES(loc_it + 1, loc_it) = tt; in KSPLGMRESCycle()
195 …hapbnd = PetscAbsScalar(tt / *GRS(loc_it)); /* GRS(loc_it) contains the res_norm from the last ite… in KSPLGMRESCycle()
197 if (tt > hapbnd) { in KSPLGMRESCycle()
198 tmp = 1.0 / tt; in KSPLGMRESCycle()
201 …cInfo(ksp, "Detected happy breakdown, current hapbnd = %g tt = %g\n", (double)hapbnd, (double)tt)); in KSPLGMRESCycle()
355 PetscScalar tt; in KSPLGMRESBuildSoln() local
401 tt = *GRS(k); in KSPLGMRESBuildSoln()
[all …]
/petsc/src/vec/is/is/utils/
H A Disblock.c30 PetscInt tt, gid1, *nidx; in ISCompressIndicesGeneral() local
70 PetscCall(PetscHMapIGetWithDefault(gid1_lid1, ival + 1, 0, &tt)); in ISCompressIndicesGeneral()
71 if (!tt) { in ISCompressIndicesGeneral()
88 PetscHashIterGetVal(gid1_lid1, tpos, tt); in ISCompressIndicesGeneral()
89 …PetscCheck(tt-- <= isz, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "index greater than array-dim"); in ISCompressIndicesGeneral()
90 nidx[tt] = gid1 - 1; in ISCompressIndicesGeneral()
/petsc/src/ksp/ksp/tutorials/
H A Dex66.c31 PetscScalar uu, tt; member
50 user.tt = 1.0; in main()
67 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
76 tt = user->tt; in ComputeRHS()
84 …i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(tt * pi * ((PetscReal)j… in ComputeRHS()
H A Dex50.c35 PetscScalar uu, tt; member
54 user.tt = 1.0; in main()
71 PetscScalar Hx, Hy, pi, uu, tt; in ComputeRHS() local
80 tt = user->tt; in ComputeRHS()
88 …i] = -PetscCosScalar(uu * pi * ((PetscReal)i + 0.5) * Hx) * PetscCosScalar(tt * pi * ((PetscReal)j… in ComputeRHS()
/petsc/src/ksp/ksp/impls/gmres/
H A Dgmres.c91 PetscReal res, hapbnd, tt; in KSPGMRESCycle() local
140 PetscCall(VecNormalize(VEC_VV(it + 1), &tt)); in KSPGMRESCycle()
141 KSPCheckNorm(ksp, tt); in KSPGMRESCycle()
144 *HH(it + 1, it) = tt; in KSPGMRESCycle()
145 *HES(it + 1, it) = tt; in KSPGMRESCycle()
148 hapbnd = PetscAbsScalar(tt / *GRS(it)); in KSPGMRESCycle()
150 if (tt < hapbnd) { in KSPGMRESCycle()
151 …sp, "Detected happy ending, current hapbnd = %14.12e tt = %14.12e\n", (double)hapbnd, (double)tt)); in KSPGMRESCycle()
301 PetscScalar tt; in KSPGMRESBuildSoln() local
324 tt = *GRS(k); in KSPGMRESBuildSoln()
[all …]
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h73 … x[0], z = x[1], i1func, i2func, i3func, ks, es, pi4pow, sqrt_1s, r2, rp2, r2prp2, zmzp, zmzp2, tt; in LandauTensor2D() local
87 tt = 1. / (1 + s); in LandauTensor2D()
91 …ellipticE(2 * s * tt, &es); /* 44 flops * 2 + 75 = 163 flops including 2 logs, 1 sqrt, 1 pow, 21 m… in LandauTensor2D()
93 ellipticK(2 * s * tt, &ks); /* 44 flops + 75 in rest, 21 mult */ in LandauTensor2D()
132 PetscReal fact, u1u2, diff[3], udiff2, u12, u22, wsq, rsq, tt; in LandauTensor3DRelativistic() local
147 …tt = 2. * u1u2 * (1. - g1 * g2) + (u12 * u22 + u1u2 * u1u2) / c02; // these two terms are about t… in LandauTensor3DRelativistic()
148 wsq = udiff2 + tt; in LandauTensor3DRelativistic()
/petsc/src/ts/utils/dmplexlandau/tutorials/
H A Dex2.c153 …PetscScalar user[2] = {0., ctx->charges[0]}, q[LANDAU_MAX_SPECIES], tt[LANDAU_MAX_SPECIES], … in testSpitzer() local
172 …PetscCall(DMPlexComputeIntegralFEM(plexe, XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, 0)], tt, NU… in testSpitzer()
173 J = -ctx->n_0 * ctx->v_0 * PetscRealPart(tt[0]); in testSpitzer()
178 …PetscCall(DMPlexComputeIntegralFEM(plexi, XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, 1)], tt, NU… in testSpitzer()
179 J += -ctx->n_0 * ctx->v_0 * PetscRealPart(tt[0]); in testSpitzer()
185 …PetscCall(DMPlexComputeIntegralFEM(plexe, XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, 0)], tt, NU… in testSpitzer()
186 n_e = PetscRealPart(tt[0]) * ctx->n_0; in testSpitzer()
195 …PetscCall(DMPlexComputeIntegralFEM(plexe, XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, 0)], tt, NU… in testSpitzer()
196 vz = ctx->n_0 * PetscRealPart(tt[0]) / n_e; /* non-dimensional */ in testSpitzer()
202 …PetscCall(DMPlexComputeIntegralFEM(plexe, XsubArray[LAND_PACK_IDX(ctx->batch_view_idx, 0)], tt, NU… in testSpitzer()
[all …]
H A Dex1.c259 PetscScalar user[2] = {0, 0}, tt[1]; in Monitor() local
265 PetscCall(DMPlexComputeIntegralFEM(dm, Xloc, tt, NULL)); in Monitor()
266 n = PetscRealPart(tt[0]); in Monitor()
269 PetscCall(DMPlexComputeIntegralFEM(dm, Xloc, tt, NULL)); in Monitor()
270 user[0] = vz_0 = PetscRealPart(tt[0]) / n; in Monitor()
274 PetscCall(DMPlexComputeIntegralFEM(dm, Xloc, tt, ctx)); in Monitor()
275 energy = PetscRealPart(tt[0]) * ctx->v_0 * ctx->v_0 * m_s / n / 3; // scale? in Monitor()
280 PetscCall(DMPlexComputeIntegralFEM(dm, Xloc, tt, ctx)); in Monitor()
281 e_par = PetscRealPart(tt[0]) * ctx->v_0 * ctx->v_0 * m_s / n; in Monitor()
286 PetscCall(DMPlexComputeIntegralFEM(dm, Xloc, tt, ctx)); in Monitor()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/ode/
H A Drober.py84 tt = np.asarray([v[1] for v in history]) variable
92 np.diff(tt),
99 pylab.semilogx(tt, xx[:, i], 'rgb'[i])
H A Dorego.py114 tt = np.asarray([v[1] for v in history]) variable
123 np.diff(tt),
130 pylab.semilogy(tt, xx[:, i], 'rgb'[i])
/petsc/src/ksp/ksp/impls/gmres/dgmres/
H A Ddgmres.c107 PetscReal res_norm, res, hapbnd, tt; in KSPDGMRESCycle() local
163 PetscCall(VecNormalize(VEC_VV(it + 1), &tt)); in KSPDGMRESCycle()
165 *HH(it + 1, it) = tt; in KSPDGMRESCycle()
166 *HES(it + 1, it) = tt; in KSPDGMRESCycle()
169 hapbnd = PetscAbsScalar(tt / *GRS(it)); in KSPDGMRESCycle()
171 if (tt < hapbnd) { in KSPDGMRESCycle()
172 …cInfo(ksp, "Detected happy breakdown, current hapbnd = %g tt = %g\n", (double)hapbnd, (double)tt)); in KSPDGMRESCycle()
321 PetscScalar tt; in KSPDGMRESBuildSoln() local
339 tt = *GRS(k); in KSPDGMRESBuildSoln()
340 for (j = k + 1; j <= it; j++) tt = tt - *HH(k, j) * nrs[j]; in KSPDGMRESBuildSoln()
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Dadjointreads.m60 % tt=senmask.*obj;
61 % tt(abs(tt)==0)=NaN;
62 % plot(xg,tt,'ks-','Markersize',10); drawnow;
/petsc/src/ksp/ksp/impls/gmres/agmres/
H A Dagmresorthog.c126 PetscReal c, s, rho, Ajj, val, tt, old; in KSPAGMRESRoddec() local
162 …PetscCallBLAS("BLASdot", tt = tloc[j] * BLASdot_(&len, &Qloc[j * nloc + j], &pas, &Qloc[k * nloc +… in KSPAGMRESRoddec()
163 …PetscCallBLAS("BLASaxpy", BLASaxpy_(&len, &tt, &Qloc[j * nloc + j], &pas, &Qloc[k * nloc + j], &pa… in KSPAGMRESRoddec()
247 PetscReal c, s, rho, zp, zq, yd = 0.0, tt; in KSPAGMRESRodvec() local
302 tt = tloc[j] * (BLASdot_(&blen, &Qloc[dpt], &pas, &zloc[j], &pas)); in KSPAGMRESRodvec()
303 PetscCallBLAS("BLASaxpy", BLASaxpy_(&blen, &tt, &Qloc[dpt], &pas, &zloc[j], &pas)); in KSPAGMRESRodvec()
/petsc/src/ts/tutorials/multirate/
H A Dex3.c86 PetscReal error, tt; in main() local
165 PetscCall(TSGetTime(ts, &tt)); in main()
166 PetscCall(sol_true(tt, Utrue)); in main()
H A Dex2.c89 PetscScalar error, tt; in main() local
176 PetscCall(TSGetTime(ts, &tt)); in main()
177 PetscCall(sol_true(tt, Utrue, &ctx)); in main()
H A Dex1.c89 PetscReal error, tt; in main() local
176 PetscCall(TSGetTime(ts, &tt)); in main()
177 PetscCall(sol_true(tt, Utrue, &ctx)); in main()
H A Dex3fastslowsplit.c117 PetscScalar error, tt; in main() local
204 PetscCall(TSGetTime(ts, &tt)); in main()
205 PetscCall(sol_true(tt, Utrue)); in main()
/petsc/src/sys/classes/random/tutorials/
H A Dex2.c155 double vv, tt; in readData() local
156 …PetscCheck(fscanf(fd, "%s%lf%lf", temp, &vv, &tt) == 3, PETSC_COMM_SELF, PETSC_ERR_FILE_UNEXPECTED… in readData()
158 t[i] = tt; in readData()
/petsc/src/mat/graphops/coarsen/impls/hem/
H A Dhem.c553 PetscReal max_e = 0., tt; in MatCoarsenApply_HEM_private() local
564 if ((tt = PetscRealPart(ap[jj])) > threshold && lidj != lid) { in MatCoarsenApply_HEM_private()
565 if (tt > max_e) max_e = tt; in MatCoarsenApply_HEM_private()
575 if ((tt = PetscRealPart(ap[jj])) > threshold) { in MatCoarsenApply_HEM_private()
576 if (tt > max_e) max_e = tt; in MatCoarsenApply_HEM_private()
627 PetscReal tt; in MatCoarsenApply_HEM_private() local
636 if ((tt = PetscRealPart(ap[jj])) > threshold && lidj != lid) { in MatCoarsenApply_HEM_private()
641 Edges[nEdges].weight = tt; in MatCoarsenApply_HEM_private()
652 if ((tt = PetscRealPart(ap[jj])) > threshold) { in MatCoarsenApply_HEM_private()
656 Edges[nEdges].weight = tt; in MatCoarsenApply_HEM_private()
[all …]
/petsc/src/ts/impls/explicit/rk/
H A Dmrk.c35 PetscReal tt, t; in TSInterpolate_RK_MultirateNonsplit() local
44 for (j = 0, tt = t; j < p; j++, tt *= t) { in TSInterpolate_RK_MultirateNonsplit()
45 for (i = 0; i < s; i++) b[i] += h * B[i * p + j] * tt; in TSInterpolate_RK_MultirateNonsplit()
238 PetscReal tt, t; in TSInterpolate_RK_MultirateSplit() local
260 for (j = 0, tt = t; j < p; j++, tt *= t) { in TSInterpolate_RK_MultirateSplit()
261 for (i = 0; i < s; i++) b[i] += h * B[i * p + j] * tt; in TSInterpolate_RK_MultirateSplit()
/petsc/src/ksp/ksp/impls/gmres/pgmres/
H A Dpgmres.c201 PetscScalar tt; in KSPPGMRESBuildSoln() local
219 tt = *RS(k); in KSPPGMRESBuildSoln()
220 for (j = k + 1; j <= it; j++) tt -= *HH(k, j) * nrs[j]; in KSPPGMRESBuildSoln()
221 nrs[k] = tt / *HH(k, k); in KSPPGMRESBuildSoln()
/petsc/src/tao/bound/tutorials/
H A Djbearing2.c218 PetscReal tt, f1, f2; in FormFunctionGradient() local
297 tt = 0; in FormFunctionGradient()
298 for (kk = 0; kk < k; kk++) tt += v[kk] * x[col[kk]]; in FormFunctionGradient()
300 g[row] = tt; in FormFunctionGradient()

12