| /petsc/src/ts/tutorials/output/ |
| H A D | ex24_2.out | 1 0 t= 0.0e+00 dt= 1.0e-05 f= 2.0e+06 df= 2.6e+05 it=( 0, 0) 2 1 t= 1.0e-05 dt= 1.1e-05 f= 1.5e+06 df= 2.1e+05 it=( 3, 3) 3 2 t= 2.1e-05 dt= 1.5e-05 f= 1.2e+06 df= 1.7e+05 it=( 3, 3) 4 3 t= 3.6e-05 dt= 2.0e-05 f= 8.5e+05 df= 1.4e+05 it=( 3, 3) 5 4 t= 5.6e-05 dt= 2.9e-05 f= 6.0e+05 df= 1.0e+05 it=( 3, 3) 6 5 t= 8.6e-05 dt= 4.4e-05 f= 4.0e+05 df= 7.7e+04 it=( 3, 3) 7 6 t= 1.3e-04 dt= 6.9e-05 f= 2.4e+05 df= 5.4e+04 it=( 3, 3) 8 7 t= 2.0e-04 dt= 1.2e-04 f= 1.4e+05 df= 3.5e+04 it=( 3, 3) 9 8 t= 3.2e-04 dt= 2.1e-04 f= 6.9e+04 df= 2.1e+04 it=( 3, 3) 10 9 t= 5.3e-04 dt= 4.3e-04 f= 3.0e+04 df= 1.1e+04 it=( 4, 4) [all …]
|
| H A D | ex24_1.out | 1 0 t=0.0000e+00 dt=1.0000e-03 f=1.3005e+05 df=7.5724e+04 it=( 0, 0) 2 1 t=1.0000e-03 dt=1.1000e-03 f=2.5359e+04 df=2.2799e+04 it=( 1, 1) 3 2 t=2.1000e-03 dt=3.9248e-03 f=5.5935e+03 df=7.0287e+03 it=( 1, 1) 4 3 t=6.0248e-03 dt=1.4459e-02 f=1.5992e+03 df=2.0987e+03 it=( 1, 1) 5 4 t=2.0484e-02 dt=5.0883e-02 f=6.4634e+02 df=6.5600e+02 it=( 1, 1) 6 5 t=7.1367e-02 dt=6.6206e-02 f=1.9371e+02 df=5.5459e+02 it=( 1, 1) 7 6 t=1.3757e-01 dt=2.0251e-01 f=5.3877e+01 df=1.9944e+02 it=( 1, 1) 8 7 t=3.4008e-01 dt=6.1534e-01 f=1.0661e+01 df=7.2201e+01 it=( 1, 1) 9 8 t=9.5542e-01 dt=3.6083e+00 f=2.1794e+00 df=1.3544e+01 it=( 1, 1) 10 9 t=4.5637e+00 dt=1.8537e+01 f=1.5598e+00 df=2.9000e+00 it=( 1, 1) [all …]
|
| /petsc/src/tao/unconstrained/impls/bmrm/ |
| H A D | bmrm.c | 9 static PetscErrorCode solve(TAO_DF *df) in solve() argument 15 PetscReal *x = df->x, *a = df->a, b = df->b, *l = df->l, *u = df->u, tol = df->tol; in solve() 16 PetscReal *tempv = df->tempv, *y = df->y, *g = df->g, *d = df->d, *Qd = df->Qd; in solve() 17 PetscReal *xplus = df->xplus, *tplus = df->tplus, *sk = df->sk, *yk = df->yk; in solve() 18 PetscReal **Q = df->Q, *f = df->f, *t = df->t; in solve() 19 PetscInt dim = df->dim, *ipt = df->ipt, *ipt2 = df->ipt2, *uv = df->uv; in solve() 37 project(dim, a, b, tempv, l, u, x, &lam_ext, df); in solve() 58 project(dim, a, b, y, l, u, tempv, &lam_ext, df); in solve() 84 for (innerIter = 1; innerIter <= df->maxPGMIter; innerIter++) { in solve() 89 project(dim, a, b, tempv, l, u, y, &lam_ext, df); in solve() [all …]
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_bucket.c | 40 DMSwarmDataField df; in DMSwarmDataFieldCreate() local 43 PetscCall(PetscNew(&df)); in DMSwarmDataFieldCreate() 44 PetscCall(PetscStrallocpy(registration_function, &df->registration_function)); in DMSwarmDataFieldCreate() 45 PetscCall(PetscStrallocpy(name, &df->name)); in DMSwarmDataFieldCreate() 46 df->atomic_size = size; in DMSwarmDataFieldCreate() 47 df->L = L; in DMSwarmDataFieldCreate() 48 df->bs = 1; in DMSwarmDataFieldCreate() 50 PetscCall(PetscMalloc(size * L, &df->data)); in DMSwarmDataFieldCreate() 51 PetscCall(PetscMemzero(df->data, size * L)); in DMSwarmDataFieldCreate() 52 *DF = df; in DMSwarmDataFieldCreate() [all …]
|
| /petsc/src/ts/tests/output/ |
| H A D | ex13_2.out | 29 df(1.) = 3. (reconstructed 3.) 47 df(9.9) = 294.03 (reconstructed 294.03) 66 df(3.) = 27. (reconstructed 27.) 76 df(1.1) = 3.63 (reconstructed 3.63) 81 df(1.1) = 3.63 (reconstructed 3.63) 97 df(5.6) = 94.08 (reconstructed 94.08) 148 df(5.6) = 94.08 (reconstructed 94.08) 162 df(1.1) = 3.63 (reconstructed 3.63) 165 df(1.1) = 3.63 (reconstructed 3.63) 173 df(3.) = 27. (reconstructed 27.) [all …]
|
| H A D | ex13_1.out | 29 df(1.) = 3. (reconstructed 3.) 47 df(9.9) = 294.03 (reconstructed 294.03) 66 df(3.) = 27. (reconstructed 27.) 76 df(1.1) = 3.63 (reconstructed 3.63) 81 df(1.1) = 3.63 (reconstructed 3.63) 97 df(5.6) = 94.08 (reconstructed 94.08) 148 df(5.6) = 94.08 (reconstructed 94.08) 162 df(1.1) = 3.63 (reconstructed 3.63) 165 df(1.1) = 3.63 (reconstructed 3.63) 173 df(3.) = 27. (reconstructed 27.) [all …]
|
| H A D | ex13_3.out | 35 df(1.) = 3. (reconstructed 3.) 55 df(9.9) = 294.03 (reconstructed 294.03) 76 df(3.) = 27. (reconstructed 27.) 88 df(1.1) = 3.63 (reconstructed 3.63) 93 df(1.1) = 3.63 (reconstructed 3.63) 109 df(5.6) = 94.08 (reconstructed 94.08) 162 df(5.6) = 94.08 (reconstructed 94.08) 178 df(1.1) = 3.63 (reconstructed 3.63) 181 df(1.1) = 3.63 (reconstructed 3.63) 191 df(3.) = 27. (reconstructed 27.) [all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/tests/ |
| H A D | solve_performance.c | 11 Vec x, g, dx, df, p; in main() local 27 PetscCall(VecDuplicate(x, &df)); in main() 48 PetscCall(VecSetRandom(df, rand)); in main() 49 PetscCall(VecDot(dx, df, &dot)); in main() 60 PetscCall(VecAXPY(g, fscale, df)); in main() 67 PetscCall(VecAXPY(g, fscale, df)); in main() 80 PetscCall(VecDestroy(&df)); in main()
|
| H A D | ex3.c | 51 Vec dx, df; in TestUsability() local 55 PetscCall(VecDuplicate(pattern, &df)); in TestUsability() 60 PetscCall(MatMult(B, dx, df)); in TestUsability() 61 PetscCall(MatSolve(B, dx, df)); in TestUsability() 65 PetscCall(VecCopy(dx, df)); in TestUsability() 66 PetscCall(MatLMVMUpdate(B, dx, df)); in TestUsability() 67 PetscCall(MatMult(B, dx, df)); in TestUsability() 68 PetscCall(MatSolve(B, dx, df)); in TestUsability() 71 PetscCall(VecDestroy(&df)); in TestUsability()
|
| H A D | ex1.c | 113 static PetscErrorCode SecantTest(Mat B, Vec dx, Vec df, PetscBool is_hermitian, PetscBool test_inve… in SecantTest() argument 120 PetscCall(VecDuplicate(df, &B_x)); in SecantTest() 122 PetscCall(VecAXPY(B_x, -1.0, df)); in SecantTest() 124 PetscCall(VecNorm(df, NORM_2, &scale)); in SecantTest() 130 PetscCall(VecAXPY(B_x, -1.0, df)); in SecantTest() 136 PetscCall(VecDuplicate(df, &B_x)); in SecantTest() 137 PetscCall(MatSolve(B, df, B_x)); in SecantTest() 146 PetscCall(MatSolveHermitianTranspose(B, df, B_x)); in SecantTest() 320 Vec x, dx, f, x_prev, f_prev, df; in TestUpdate() local 332 df = dfs[iter]; in TestUpdate() [all …]
|
| /petsc/lib/petsc/bin/ |
| H A D | petsc_tas_analysis.py | 327 df = pd.read_csv(fileName) 328 Nf = getNfCSV(df) 329 nProcs = int(df.columns.tolist()[25]) 354 SNESSolveFilter = (df['Event Name'] == 'SNESSolve') 356 TSStepFilter = (df['Event Name'] == 'TSStep') 357 MatLUFactorFilter = ((df['Event Name'] == 'MatLUFactorNum') 358 | (df['Event Name'] == 'MatLUFactorSym')) 359 ConvEstErrorFilter = (df['Event Name'] == 'ConvEst Error') 360 rankFilter = (df['Rank'] == 0) 365 …if (df.loc[SNESSolveFilter & (df['Stage Name'] == 'ConvEst Refinement Level 0') & rankFilter, 'Tim… [all …]
|
| /petsc/src/dm/dt/fe/interface/ |
| H A D | fe.c | 2599 for (PetscInt df = 0; df < _dE; ++df) s += G[df] * tBDJ[df]; \ 2617 for (PetscInt df = 0; df < _dE; ++df) s += tBDI[df] * G[df]; \ 2635 for (PetscInt df = 0; df < (_dE); ++df) { \ 2636 … for (PetscInt dg = 0; dg < (_dE); ++dg) s += tBDI[df] * G[df * (_dE) + dg] * tBDJ[dg]; \ 2665 … for (PetscInt df = 0; df < cdim; ++df) tmpBasisDerI[fidx * dE + df] = basisDerI[fidx * cdim + df]; in PetscFEUpdateElementMat_Internal() local 2724 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasi… in PetscFEUpdateElementMat_Internal() 2754 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ … in PetscFEUpdateElementMat_Internal() 2784 for (PetscInt df = 0; df < dE; ++df) { in PetscFEUpdateElementMat_Internal() 2785 …; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * dE + df) * … in PetscFEUpdateElementMat_Internal() 2815 PetscInt f, fc, g, gc, df, dg; in PetscFEUpdateElementMat_Hybrid_Internal() local [all …]
|
| /petsc/src/ts/utils/ |
| H A D | dmplexts.c | 407 Vec du, uhat, uhat_t, r, rhat, df; in DMTSCheckJacobian() local 418 PetscCall(VecDuplicate(u, &df)); in DMTSCheckJacobian() 419 PetscCall(MatMult(J, du, df)); in DMTSCheckJacobian() 434 PetscCall(VecAXPBYPCZ(rhat, -1.0, -h, 1.0, r, df)); in DMTSCheckJacobian() 443 PetscCall(VecDestroy(&df)); in DMTSCheckJacobian()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex48.c | 192 static PetscScalar poissonBracket(PetscInt dim, const PetscScalar df[], const PetscScalar dg[]) in poissonBracket() argument 194 PetscScalar ret = df[0] * dg[1] - df[1] * dg[0]; in poissonBracket() 196 ret += df[1] * dg[2] - df[2] * dg[1]; in poissonBracket() 197 ret += df[2] * dg[0] - df[0] * dg[2]; in poissonBracket()
|
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 943 Vec du, uhat, r, rhat, df; in DMSNESCheckJacobian() local 954 PetscCall(VecDuplicate(u, &df)); in DMSNESCheckJacobian() 955 PetscCall(MatMult(J, du, df)); in DMSNESCheckJacobian() 968 PetscCall(VecAXPBYPCZ(rhat, -1.0, -h, 1.0, r, df)); in DMSNESCheckJacobian() 976 PetscCall(VecDestroy(&df)); in DMSNESCheckJacobian()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | wedge_cylinder.cas | 255 (13 (8 df 108 3 4)
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 2561 PetscReal phi, sinp, cosp, dc, df, x, y, xc, yc; in snapToCylinder() local 2570 df = PetscAbsReal(dis / sinp); in snapToCylinder() 2575 df = PetscAbsReal(dis / cosp); in snapToCylinder() 2579 f0[0] = xc + (u[0] - xc) * (1.0 - dc) / (df - dc); in snapToCylinder() 2580 f0[1] = yc + (u[1] - yc) * (1.0 - dc) / (df - dc); in snapToCylinder()
|
| /petsc/doc/manual/ |
| H A D | ts.md | 335 If the Jacobian $\frac{dh}{du} \frac{df}{dz}$ is non-singular, then we have precisely a semi-explic…
|