| /petsc/src/tao/leastsquares/impls/brgn/ ! |
| H A D | brgn.c | 7 TAO_BRGN *gn; in GNHessianProd() local
10 PetscCall(MatShellGetContext(H, &gn)); in GNHessianProd()
11 PetscCall(MatMult(gn->subsolver->ls_jac, in, gn->r_work)); in GNHessianProd()
12 PetscCall(MatMultTranspose(gn->subsolver->ls_jac, gn->r_work, out)); in GNHessianProd()
13 switch (gn->reg_type) { in GNHessianProd()
15 PetscCall(MatMult(gn->Hreg, in, gn->x_work)); in GNHessianProd()
16 PetscCall(VecAXPY(out, gn->lambda, gn->x_work)); in GNHessianProd()
19 PetscCall(VecAXPY(out, gn->lambda, in)); in GNHessianProd()
22 PetscCall(VecAXPY(out, gn->lambda, in)); in GNHessianProd()
26 if (gn->D) { in GNHessianProd()
[all …]
|
| /petsc/src/snes/tutorials/ ! |
| H A D | ex18.c | 297 PetscScalar dn, ds, de, dw, an, as, ae, aw, bn, bs, be, bw, gn, gs, ge, gw; in FormJacobian() local
358 gn = coef * bn * (tn - t0); in FormJacobian()
366 v[2] = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
372 v[4] = -hxdhy * (dn + gn); in FormJacobian()
400 gn = coef * bn * (tn - t0); in FormJacobian()
402 v[0] = hxdhy * (dn - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
408 v[2] = -hxdhy * (dn + gn); in FormJacobian()
427 gn = coef * bn * (tn - t0); in FormJacobian()
432 v[1] = hxdhy * (ds + dn + gs - gn) + hydhx * (dw + de + gw - ge); in FormJacobian()
438 v[3] = -hxdhy * (dn + gn); in FormJacobian()
[all …]
|
| /petsc/src/snes/tests/ ! |
| H A D | ex20.c | 452 PetscScalar ***x, bn, bs, be, bw, bu, bd, gn, gs, ge, gw, gu, gd; in FormJacobian() local
516 gn = coef * bn * (tn - t0); in FormJacobian()
547 …v[3] = hzhxdhy * (ds + dn + gs - gn) + hyhzdhx * (dw + de + gw - ge) + hxhydhz * (dd + du + gd -… in FormJacobian()
555 v[5] = -hzhxdhy * (dn + gn); in FormJacobian()
585 gn = coef * bn * (tn - t0); in FormJacobian()
599 v[0] = hzhxdhy * (dn - gn) + hyhzdhx * (dw + de + gw - ge) + hxhydhz * (du - gu); in FormJacobian()
607 v[2] = -hzhxdhy * (dn + gn); in FormJacobian()
637 … v[1] = hzhxdhy * (dn - gn) + hyhzdhx * (dw + de + gw - ge) + hxhydhz * (dd + du + gd - gu); in FormJacobian()
645 v[3] = -hzhxdhy * (dn + gn); in FormJacobian()
668 v[1] = hzhxdhy * (dn - gn) + hyhzdhx * (dw + de + gw - ge) + hxhydhz * (dd + gd); in FormJacobian()
[all …]
|
| /petsc/src/dm/impls/da/ ! |
| H A D | grglvis.c | 97 PetscInt sx, sy, sz, gsx, gsy, gsz, ist, jst, kst, gm, gn, gp; in DMDASampleGLVisFields_Private() local
107 PetscCall(DMDAGetGhostCorners(da, &gsx, &gsy, &gsz, &gm, &gn, &gp)); in DMDASampleGLVisFields_Private()
127 i = ke * gm * gn + je * gm + ie; in DMDASampleGLVisFields_Private()
344 PetscInt gsx, gsy, gsz, gm, gn, gp, kst, jst, ist; in DMDAView_GLVis_ASCII() local
471 gn = jen; in DMDAView_GLVis_ASCII()
480 PetscCall(DMDAGetGhostCorners(da, &gsx, &gsy, &gsz, &gm, &gn, &gp)); in DMDAView_GLVis_ASCII()
494 i = ke * gm * gn + je * gm + ie; in DMDAView_GLVis_ASCII()
|
| /petsc/src/dm/dt/interface/ ! |
| H A D | dtds.c | 3740 PetscInt Nf, Nfn, fn, gn; in PetscDSSelectEquations() local
3765 for (gn = 0; gn < numFields; ++gn) { in PetscDSSelectEquations()
3766 const PetscInt g = fields ? fields[gn] : gn; in PetscDSSelectEquations()
3775 PetscCall(PetscDSSetJacobian(newprob, fn, gn, g0, g1, g2, g3)); in PetscDSSelectEquations()
3776 PetscCall(PetscDSSetJacobianPreconditioner(newprob, fn, gn, g0p, g1p, g2p, g3p)); in PetscDSSelectEquations()
3777 PetscCall(PetscDSSetBdJacobian(newprob, fn, gn, g0Bd, g1Bd, g2Bd, g3Bd)); in PetscDSSelectEquations()
|
| /petsc/src/dm/impls/plex/cgns/ ! |
| H A D | plexcgns2.c | 2126 PetscInt gn = node_l2g[n]; in DMView_PlexCGNS() local
2127 if (gn < nStart || nEnd <= gn) continue; in DMView_PlexCGNS()
2128 x[gn - nStart] = X[n * coord_dim + d]; in DMView_PlexCGNS()
2472 PetscInt gn = cgv->node_l2g[n]; in VecView_Plex_Local_CGNS() local
2473 if (gn < cgv->nStart || cgv->nEnd <= gn) continue; in VecView_Plex_Local_CGNS()
2474 cgv->nodal_field[gn - cgv->nStart] = v[off + c]; in VecView_Plex_Local_CGNS()
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ ! |
| H A D | petscsnes.pxi | 441 cdef object gn = toReal(gnorm)
447 reason = converged(Snes, it, (xn, gn, fn), *args, **kargs)
|
| /petsc/src/dm/impls/plex/hdf5/ ! |
| H A D | plexhdf5.c | 1507 PetscInt *gspoints, n = 0, gn, p; in DMPlexLabelsView_HDF5_Internal() local
1517 for (gn = 0, p = 0; p < n; ++p) in DMPlexLabelsView_HDF5_Internal()
1518 if (gpoint[spoints[p]] >= 0) ++gn; in DMPlexLabelsView_HDF5_Internal()
1519 PetscCall(PetscMalloc1(gn, &gspoints)); in DMPlexLabelsView_HDF5_Internal()
1520 for (gn = 0, p = 0; p < n; ++p) in DMPlexLabelsView_HDF5_Internal()
1521 if (gpoint[spoints[p]] >= 0) gspoints[gn++] = gpoint[spoints[p]]; in DMPlexLabelsView_HDF5_Internal()
1523 …PetscCall(ISCreateGeneral(PetscObjectComm((PetscObject)dm), gn, gspoints, PETSC_OWN_POINTER, &glob… in DMPlexLabelsView_HDF5_Internal()
|
| /petsc/src/dm/impls/plex/ ! |
| H A D | plexcreate.c | 383 PetscInt gn; in DMPlexCreateCoordinateSpace() local
386 PetscCall(VecGetLocalSize(coordinatesNew, &gn)); in DMPlexCreateCoordinateSpace()
387 …gn == n * cdim, PetscObjectComm((PetscObject)dm), PETSC_ERR_ARG_WRONG, "Global coordinate size %" … in DMPlexCreateCoordinateSpace()
|