| /petsc/src/tao/unconstrained/tutorials/output/ |
| H A D | rosenbrock1_20.out | 5 sty: 0., fy: 1., dgy: -2. 8 sty: 1., fy: 99., dgy: 396. 13 sty: 0., fy: 0.773446, dgy: -0.154432 18 sty: 0., fy: 0.648599, dgy: -0.340628 21 sty: 1., fy: 0.865442, dgy: 1.37753 26 sty: 0., fy: 0.560405, dgy: -0.260357 31 sty: 0., fy: 0.379185, dgy: -0.125226 36 sty: 0., fy: 0.280758, dgy: -0.17334 41 sty: 0., fy: 0.187322, dgy: -0.878193 44 sty: 1., fy: 51.7671, dgy: 195.182 [all …]
|
| H A D | rosenbrock1_20_alt.out | 5 sty: 0., fy: 1., dgy: -2. 8 sty: 1., fy: 99., dgy: 396. 13 sty: 0., fy: 0.773446, dgy: -0.154432 18 sty: 0., fy: 0.648599, dgy: -0.340628 21 sty: 1., fy: 0.865443, dgy: 1.37753 26 sty: 0., fy: 0.560405, dgy: -0.260357 31 sty: 0., fy: 0.379185, dgy: -0.125226 36 sty: 0., fy: 0.280758, dgy: -0.173338 41 sty: 0., fy: 0.187321, dgy: -0.878193 44 sty: 1., fy: 51.7692, dgy: 195.19 [all …]
|
| /petsc/src/tao/linesearch/impls/morethuente/ |
| H A D | morethuente.c | 10 …, PetscReal *stx, PetscReal *fx, PetscReal *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, Pets… 29 …rintf(ls->viewer, "sty: %g, fy: %g, dgy: %g\n", (double)mt->sty, (double)mt->fy, (double)mt->dgy)); in TaoLineSearchMonitor_MT() 38 PetscReal dgx, dgy, dg, dg2, fx, fy, stx, sty, dgtest; in TaoLineSearchApply_MT() local 109 fy = finit; in TaoLineSearchApply_MT() 164 mt->fy = fy; in TaoLineSearchApply_MT() 230 fym = fy - sty * dgtest; in TaoLineSearchApply_MT() 240 fy = fym + sty * dgtest; /* gradient values */ in TaoLineSearchApply_MT() 245 PetscCall(Tao_mcstep(ls, &stx, &fx, &dgx, &sty, &fy, &dgy, &ls->step, f, &dg)); in TaoLineSearchApply_MT() 363 …, PetscReal *stx, PetscReal *fx, PetscReal *dx, PetscReal *sty, PetscReal *fy, PetscReal *dy, Pets… in Tao_mcstep() argument 468 theta = 3 * (*fp - *fy) / (*sty - *stp) + *dy + *dp; in Tao_mcstep() [all …]
|
| H A D | morethuente.h | 11 PetscReal sty, fy, dgy; member
|
| /petsc/src/ts/impls/implicit/sundials/ |
| H A D | sundials.c | 14 static PetscErrorCode TSPrecond_Sundials_Petsc(realtype tn, N_Vector y, N_Vector fy, booleantype jo… in TSPrecond_Sundials_Petsc() argument 41 static int TSPrecond_Sundials_Private(realtype tn, N_Vector y, N_Vector fy, booleantype jok, boolea… in TSPrecond_Sundials_Private() argument 43 …return (int)TSPrecond_Sundials_Petsc(tn, y, fy, jok, jcurPtr, _gamma, P_data, vtemp1, vtemp2, vtem… in TSPrecond_Sundials_Private() 49 static PetscErrorCode TSPSolve_Sundials_Petsc(realtype tn, N_Vector y, N_Vector fy, N_Vector r, N_V… in TSPSolve_Sundials_Petsc() argument 73 static int TSPSolve_Sundials_Private(realtype tn, N_Vector y, N_Vector fy, N_Vector r, N_Vector z, … in TSPSolve_Sundials_Private() argument 75 return (int)TSPSolve_Sundials_Petsc(tn, y, fy, r, z, _gamma, delta, lr, P_data, vtemp); in TSPSolve_Sundials_Private()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex43.c | 74 PetscScalar fy[GAUSS_POINTS]; member 458 …ouble)PetscRealPart(_coefficients[j][i].fx[p]), (double)PetscRealPart(_coefficients[j][i].fy[p]))); in DMDAViewCoefficientsGnuplot2d() 661 … void FormMomentumRhsQ1(PetscScalar Fe[], PetscScalar coords[], PetscScalar fx[], PetscScalar fy[]) in FormMomentumRhsQ1() argument 683 Fe[NSD * i + 1] += fac * Ni_p[i] * fy[p]; in FormMomentumRhsQ1() 898 prop_fy = props[ej][ei].fy; in AssembleF_Stokes() 1254 …element_props[j][i].fy[p] = PetscSinReal(opts_nz * PETSC_PI * coord_y) * PetscCosReal(1.0 * PETSC_… in solve_stokes_2d_coupled() 1279 element_props[j][i].fy[p] = 0.0; in solve_stokes_2d_coupled() 1285 element_props[j][i].fy[p] = -1.0; in solve_stokes_2d_coupled() 1310 element_props[j][i].fy[p] = 0.0; in solve_stokes_2d_coupled() 1316 element_props[j][i].fy[p] = -1.0; in solve_stokes_2d_coupled() [all …]
|
| H A D | ex49.c | 84 PetscScalar E, nu, fx, fy; member 93 PetscScalar fy[GAUSS_POINTS]; member 406 (double)PetscRealPart(_coefficients[j][i].fy[p]))); in DMDAViewCoefficientsGnuplot2d() 481 … void FormMomentumRhsQ1(PetscScalar Fe[], PetscScalar coords[], PetscScalar fx[], PetscScalar fy[]) in FormMomentumRhsQ1() argument 503 Fe[NSD * i + 1] += fac * Ni_p[i] * fy[p]; in FormMomentumRhsQ1() 678 prop_fy = props[ej][ei].fy; in AssembleF_Elasticity() 846 element_props[j][i].fy[p] = 0.0; in solve_elasticity_2d() 877 element_props[j][i].fy[p] = 0.0; in solve_elasticity_2d()
|
| H A D | ex42.c | 29 PetscScalar fy[GAUSS_POINTS]; member 823 …sQ13D(PetscScalar Fe[], PetscScalar coords[], PetscScalar fx[], PetscScalar fy[], PetscScalar fz[]) in FormMomentumRhsQ13D() argument 845 Fe[NSD * i + 1] -= fac * Ni_p[i] * fy[p]; in FormMomentumRhsQ13D() 1107 prop_fy = props->fy; in AssembleF_Stokes() 1876 … cell->fy[p] = -PetscSinReal(2.2 * PETSC_PI * coord_y) * PetscCosReal(1.0 * PETSC_PI * coord_x); in solve_stokes_3d_coupled() 1906 cell->fy[p] = Fm[1]; in solve_stokes_3d_coupled() 1930 … cell->fy[p] = -PetscSinReal(3.0 * PETSC_PI * coord_y) * PetscCosReal(1.0 * PETSC_PI * coord_x); in solve_stokes_3d_coupled() 1953 cell->fy[p] = 0.0; in solve_stokes_3d_coupled() 1995 cell->fy[p] = -PetscSinReal(2.2 * PETSC_PI * yp) * PetscCosReal(1.0 * PETSC_PI * xp); in solve_stokes_3d_coupled()
|
| H A D | ex70.c | 301 … void LForm_MomentumRHS(PetscScalar Fe[], PetscScalar coords[], PetscScalar fx[], PetscScalar fy[]) in LForm_MomentumRHS() argument 320 Fe[NSD * i + 1] -= fac * Ni_p[i] * fy[p]; in LForm_MomentumRHS()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcreate.c | 1570 PetscInt c, f, fx, fy, fz, e, ex, ey, ez; in DMPlexCreateCubeMesh_Internal() local 1579 for (fy = 0; fy < numYEdges; ++fy) { in DMPlexCreateCubeMesh_Internal() 1581 PetscInt cell = (fz * numYEdges + fy) * numXEdges + fx; in DMPlexCreateCubeMesh_Internal() 1582 PetscInt faceB = firstZFace + (fy * numXEdges + fx) * numZVertices + fz; in DMPlexCreateCubeMesh_Internal() 1583 … PetscInt faceT = firstZFace + (fy * numXEdges + fx) * numZVertices + ((fz + 1) % numZVertices); in DMPlexCreateCubeMesh_Internal() 1584 PetscInt faceF = firstYFace + (fz * numXEdges + fx) * numYVertices + fy; in DMPlexCreateCubeMesh_Internal() 1585 … PetscInt faceK = firstYFace + (fz * numXEdges + fx) * numYVertices + ((fy + 1) % numYVertices); in DMPlexCreateCubeMesh_Internal() 1586 PetscInt faceL = firstXFace + (fz * numYEdges + fy) * numXVertices + fx; in DMPlexCreateCubeMesh_Internal() 1587 … PetscInt faceR = firstXFace + (fz * numYEdges + fy) * numXVertices + ((fx + 1) % numXVertices); in DMPlexCreateCubeMesh_Internal() 1602 …if (bdY != DM_BOUNDARY_NONE && fy == numYEdges - 1 && cutLabel) PetscCall(DMLabelSetValue(cutLabel… in DMPlexCreateCubeMesh_Internal() [all …]
|
| /petsc/src/dm/impls/stag/tutorials/ |
| H A D | ex2.c | 67 static PetscScalar fy(PetscScalar x, PetscScalar y) in fy() function 337 valRhs = fy(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem()
|
| H A D | ex3.c | 92 static PetscScalar fy(PetscScalar x, PetscScalar y, PetscScalar z) in fy() function 1203 …valRhs = fy(arrCoord[ez][ey][ex][icuy[0]], arrCoord[ez][ey][ex][icuy[1]], arrCoord[ez][ey][ex][icu… in CreateSystem()
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex18.c | 69 static PetscScalar fy(PetscScalar x, PetscScalar y) in fy() function 465 valRhs = dv * fy(cArrX[ex][icenter], cArrY[ey][iprev]); in CreateSystem()
|