Home
last modified time | relevance | path

Searched refs:pt (Results 1 – 25 of 31) sorted by relevance

12

/petsc/src/ts/tutorials/advection-diffusion-reaction/ !
H A Dshashi.F90104 PetscScalar temp, pt
108 pt = 0.1
197 PetscScalar temp, pt
206 pt = 0.1
300 part_p(ip) = (an_r(ip)/an_t)*pt
310 - (an_r(4)**i_cc)*(an_r(5)**i_h2o)*((pt/an_t)**idiff)
393 PetscScalar temp, pt
415 pt = 0.1
552 pt_cube = pt*pt*pt
553 pt_four = pt_cube*pt
[all …]
/petsc/src/dm/dt/space/impls/point/ !
H A Dspacepoint.c6 PetscSpace_Point *pt = (PetscSpace_Point *)sp->data; in PetscSpacePointView_Ascii() local
14 PetscCall(PetscQuadratureView(pt->quad, viewer)); in PetscSpacePointView_Ascii()
16 …space in dimension %" PetscInt_FMT " on %" PetscInt_FMT " points\n", sp->Nv, pt->quad->numPoints)); in PetscSpacePointView_Ascii()
34 PetscSpace_Point *pt = (PetscSpace_Point *)sp->data; in PetscSpaceSetUp_Point() local
37 if (!pt->quad->points && sp->degree >= 0) { in PetscSpaceSetUp_Point()
38 PetscCall(PetscQuadratureDestroy(&pt->quad)); in PetscSpaceSetUp_Point()
39 …PetscDTStroudConicalQuadrature(sp->Nv, sp->Nc, PetscMax(sp->degree + 1, 1), -1.0, 1.0, &pt->quad)); in PetscSpaceSetUp_Point()
46 PetscSpace_Point *pt = (PetscSpace_Point *)sp->data; in PetscSpaceDestroy_Point() local
49 PetscCall(PetscQuadratureDestroy(&pt->quad)); in PetscSpaceDestroy_Point()
50 PetscCall(PetscFree(pt)); in PetscSpaceDestroy_Point()
[all …]
/petsc/src/dm/dt/space/impls/ptrimmed/ !
H A Dspaceptrimmed.c5 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data; in PetscSpaceSetFromOptions_Ptrimmed() local
9 …"form degree of trimmed space", "PetscSpacePTrimmedSetFormDegree", pt->formDegree, &pt->formDegree… in PetscSpaceSetFromOptions_Ptrimmed()
16 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data; in PetscSpacePTrimmedView_Ascii() local
20 f = pt->formDegree; in PetscSpacePTrimmedView_Ascii()
40 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data; in PetscSpaceDestroy_Ptrimmed() local
45 if (pt->subspaces) { in PetscSpaceDestroy_Ptrimmed()
48 for (d = 0; d < sp->Nv; ++d) PetscCall(PetscSpaceDestroy(&pt->subspaces[d])); in PetscSpaceDestroy_Ptrimmed()
50 PetscCall(PetscFree(pt->subspaces)); in PetscSpaceDestroy_Ptrimmed()
51 PetscCall(PetscFree(pt)); in PetscSpaceDestroy_Ptrimmed()
57 PetscSpace_Ptrimmed *pt = (PetscSpace_Ptrimmed *)sp->data; in PetscSpaceSetUp_Ptrimmed() local
[all …]
/petsc/src/dm/dt/tests/ !
H A Dex13.c32 for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt]; in constructTabulationAndMass() local
111 for (PetscInt pt = 0; pt < npoints; pt++) v += p_i[pt] * p_j[pt] * weights[pt]; in test() local
190 for (PetscInt pt = 0; pt < npoints; pt++) p_i[pt] += p_j[pt] * points[pt * dim + k] * sign; in test() local
204 for (PetscInt pt = 0; pt < npoints; pt++) val += p_i[pt] * p_j[pt] * weights[pt]; in test() local
220 for (PetscInt pt = 0; pt < npoints; pt++) val += p_i[pt] * p_j[pt] * weights[pt]; in test() local
/petsc/doc/manual/anl_tech_report/ !
H A Dfirst.inc7 \vspace*{8pt}
9 \vspace*{8pt}
16 \vspace*{10pt}
19 \vspace*{20pt}
/petsc/src/dm/dt/space/impls/poly/ !
H A Dspacepoly.c123 … for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = 1.; in CoordinateBasis() local
125 …for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = points[… in CoordinateBasis() local
128 … for (PetscInt pt = 0; pt < npoints; pt++) pScalar[b * Njet * npoints + j * npoints + pt] = 1.; in CoordinateBasis() local
/petsc/src/dm/tests/ !
H A Dex33.c11 PetscInt stencil_width = 1, pt = 0, st = 0; in main() local
35 PetscCall(PetscOptionsGetInt(NULL, NULL, "-periodic", &pt, NULL)); in main()
36 if (pt == 1) bx = DM_BOUNDARY_PERIODIC; in main()
37 if (pt == 2) by = DM_BOUNDARY_PERIODIC; in main()
38 if (pt == 4) { in main()
115 …idth %" PetscInt_FMT " stencil_type %d periodic %d\n", stencil_width, (int)stencil_type, (int)pt)); in main()
H A Dex23.c9 …CIDE, n = PETSC_DECIDE, p = PETSC_DECIDE, M = 8, dof = 1, stencil_width = 1, P = 5, pt = 0, st = 0; in main() local
27 PetscCall(PetscOptionsGetInt(NULL, NULL, "-periodic", &pt, NULL)); in main()
29 if (pt == 1) bx = DM_BOUNDARY_PERIODIC; in main()
30 if (pt == 2) by = DM_BOUNDARY_PERIODIC; in main()
31 if (pt == 3) { in main()
35 if (pt == 4) bz = DM_BOUNDARY_PERIODIC; in main()
H A Dex15.c13 PetscBool pt; in main() local
23 PetscCall(PetscOptionsGetBool(NULL, NULL, "-periodic", (PetscBool *)&pt, NULL)); in main()
25 if (pt) { in main()
H A Dex7.c8 … i, start, end, P = 5, N = 6, m = PETSC_DECIDE, n = PETSC_DECIDE, p = PETSC_DECIDE, pt = 0, st = 0; in main() local
26 PetscCall(PetscOptionsGetInt(NULL, NULL, "-periodic", &pt, NULL)); in main()
28 periodic = (DMBoundaryType)pt; in main()
/petsc/src/mat/graphops/coarsen/impls/hem/ !
H A Dhem.c829 PetscInt *sbuff, *pt, scount; in MatCoarsenApply_HEM_private() local
840 sbuff = pt = sbuff + request_size; in MatCoarsenApply_HEM_private()
842 *pt++ = ndel; // number of deleted to send in MatCoarsenApply_HEM_private()
843 *pt++ = rank; // proc (not used) in MatCoarsenApply_HEM_private()
853 *pt++ = gid1; in MatCoarsenApply_HEM_private()
854 *pt++ = lid0 + Istart; // gid0 in MatCoarsenApply_HEM_private()
856 …PetscCheck(pt - sbuff == scount, PETSC_COMM_SELF, PETSC_ERR_SUP, "sbuff-pt != scount: %zu", pt - s… in MatCoarsenApply_HEM_private()
865 PetscInt *pt, *pt2, *pt3, *sbuff, tmp; in MatCoarsenApply_HEM_private() local
879 pt = rbuff; in MatCoarsenApply_HEM_private()
881 ndel = *pt++; // number of deleted to recv in MatCoarsenApply_HEM_private()
[all …]
/petsc/src/snes/impls/patch/ !
H A Dsnespatch.c18 PetscInt pt, size, i; in SNESPatchComputeResidual_Private() local
25 pt = pcpatch->currentPatch; in SNESPatchComputeResidual_Private()
26 PetscCall(ISGetSize(pcpatch->dofMappingWithoutToWithAll[pt], &size)); in SNESPatchComputeResidual_Private()
28 PetscCall(ISGetIndices(pcpatch->dofMappingWithoutToWithAll[pt], &indices)); in SNESPatchComputeResidual_Private()
36 PetscCall(ISRestoreIndices(pcpatch->dofMappingWithoutToWithAll[pt], &indices)); in SNESPatchComputeResidual_Private()
38 PetscCall(PCPatchComputeFunction_Internal(pc, pcpatch->patchStateWithAll, F, pt)); in SNESPatchComputeResidual_Private()
46 PetscInt pt, size, i; in SNESPatchComputeJacobian_Private() local
53 pt = pcpatch->currentPatch; in SNESPatchComputeJacobian_Private()
54 PetscCall(ISGetSize(pcpatch->dofMappingWithoutToWithAll[pt], &size)); in SNESPatchComputeJacobian_Private()
56 PetscCall(ISGetIndices(pcpatch->dofMappingWithoutToWithAll[pt], &indices)); in SNESPatchComputeJacobian_Private()
[all …]
/petsc/src/mat/impls/aij/seq/mkl_pardiso/ !
H A Dmkl_pardiso.c41 void pardiso_64init(void *pt, INT_TYPE *mtype, INT_TYPE iparm[]) in pardiso_64init() argument
46 pardisoinit(pt, &mtype_copy, iparm_copy); in pardiso_64init()
69 void *pt[IPARM_SIZE]; member
334 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatDestroy_MKL_PARDISO()
398 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatSolve_MKL_PARDISO()
403 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatSolve_MKL_PARDISO()
430 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatSolve_MKL_PARDISO()
456 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatForwardSolve_MKL_PARDISO()
481 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatBackwardSolve_MKL_PARDISO()
529 …PetscCallPardiso(MKL_PARDISO(mat_mkl_pardiso->pt, &mat_mkl_pardiso->maxfct, &mat_mkl_pardiso->mnum… in MatMatSolve_MKL_PARDISO()
[all …]
/petsc/src/dm/dt/interface/ !
H A Ddt.c497 PetscInt dim, Nc, imageNc, formSize, Ncopies, imageFormSize, Npoints, pt, i, j, c; in PetscQuadraturePushForward() local
518 for (pt = 0; pt < Npoints; pt++) { in PetscQuadraturePushForward()
519 const PetscReal *point = PetscSafePointerPlusOffset(points, pt * dim); in PetscQuadraturePushForward()
520 PetscReal *imagePoint = &imagePoints[pt * imageDim]; in PetscQuadraturePushForward()
529 const PetscReal *form = &weights[pt * Nc + c * formSize]; in PetscQuadraturePushForward()
530 PetscReal *imageForm = &imageWeights[pt * imageNc + c * imageFormSize]; in PetscQuadraturePushForward()
1099 PetscInt degidx, kidx, d, pt; in PetscDTPKDEvalJet() local
1163 for (pt = 0; pt < npoints; pt++) { in PetscDTPKDEvalJet()
1167 thetanm1x = dim - (d + 1) + 2. * points[pt * dim + d]; in PetscDTPKDEvalJet()
1168 for (e = d + 1; e < dim; e++) thetanm1x += points[pt * dim + e]; in PetscDTPKDEvalJet()
[all …]
/petsc/src/ksp/ksp/impls/gmres/agmres/ !
H A Dagmresleja.c6 static PetscErrorCode KSPAGMRESLejafmaxarray(PetscScalar *re, PetscInt pt, PetscInt n, PetscInt *po… in KSPAGMRESLejafmaxarray() argument
14 for (i = pt; i < n; i++) { in KSPAGMRESLejafmaxarray()
/petsc/src/ksp/ksp/tutorials/output/ !
H A Dbench_kspsolve_matmult.out3 Input matrix: 27-pt finite difference stencil
H A Dbench_kspsolve_ksp.out3 Input matrix: 27-pt finite difference stencil
/petsc/src/snes/tutorials/network/power/ !
H A DPFReadData.c106 double pg, qg, qt, qb, vs, mbase, pt, pb; in PFReadMatPowerData() local
108 …d %lf %lf %lf %lf %lf %lf %d %lf %lf", &bus_i, &pg, &qg, &qt, &qb, &vs, &mbase, &status, &pt, &pb); in PFReadMatPowerData()
117 Gen[geni].pt = pt; in PFReadMatPowerData()
H A Dpower.h92 PetscScalar pt; /* Gen max active power output: MW */ member
/petsc/src/mat/impls/aij/mpi/mkl_cpardiso/ !
H A Dmkl_cpardiso.c86 void *pt[IPARM_SIZE]; member
359 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatDestroy_MKL_CPARDISO()
389 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatSolve_MKL_CPARDISO()
412 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatForwardSolve_MKL_CPARDISO()
435 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatBackwardSolve_MKL_CPARDISO()
473 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatMatSolve_MKL_CPARDISO()
495 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatFactorNumeric_MKL_CPARDISO()
668 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatLUFactorSymbolic_AIJMKL_CPARDISO()
702 …PetscCallCluster(cluster_sparse_solver(mat_mkl_cpardiso->pt, &mat_mkl_cpardiso->maxfct, &mat_mkl_c… in MatCholeskyFactorSymbolic_AIJMKL_CPARDISO()
/petsc/src/ts/tutorials/output/ !
H A Dfig.tex.out9 \setlength{\topsep}{0pt}%
/petsc/src/benchmarks/results/ !
H A Dlap2d.m2 % Performance results for PETSc on the 2d Laplacian 5 pt stencil using GMRES(30) and no
/petsc/src/snes/tests/ !
H A Dex2.c47 PetscInt pt, fn; in ProcessOptions() local
55 pt = options->pointType; in ProcessOptions()
56 …point_type", "The point type", "ex2.c", pointTypes, 3, pointTypes[options->pointType], &pt, NULL)); in ProcessOptions()
57 options->pointType = (PointType)pt; in ProcessOptions()
/petsc/src/sys/classes/draw/impls/x/ !
H A Dxops.c173 XPoint pt[3]; in PetscDrawTriangle_X() local
175 pt[0].x = (short)XTRANS(draw, XiWin, X1); in PetscDrawTriangle_X()
176 pt[0].y = (short)YTRANS(draw, XiWin, Y_1); in PetscDrawTriangle_X()
177 pt[1].x = (short)XTRANS(draw, XiWin, X2); in PetscDrawTriangle_X()
178 pt[1].y = (short)YTRANS(draw, XiWin, Y2); in PetscDrawTriangle_X()
179 pt[2].x = (short)XTRANS(draw, XiWin, X3); in PetscDrawTriangle_X()
180 pt[2].y = (short)YTRANS(draw, XiWin, Y3); in PetscDrawTriangle_X()
181 …XFillPolygon(XiWin->disp, PetscDrawXiDrawable(XiWin), XiWin->gc.set, pt, 3, Convex, CoordModeOrigi… in PetscDrawTriangle_X()
/petsc/src/ksp/pc/impls/sor/ !
H A Dsor.tex1 \documentclass[11pt,english,pdftex]{article}

12