Home
last modified time | relevance | path

Searched refs:ns (Results 1 – 25 of 33) sorted by relevance

12

/petsc/src/vec/is/tests/
H A Dex7.c11 PetscInt *p, *ns, **ids, *leaves, *mleaves; in main() local
94 PetscCall(ISLocalToGlobalMappingGetInfo(ltog, &np, &p, &ns, &ids)); in main()
97 … procs[%" PetscInt_FMT "] = %" PetscInt_FMT ", shared %" PetscInt_FMT "\n", rank, i, p[i], ns[i])); in main()
98 …for (j = 0; j < ns[i]; j++) PetscCall(PetscViewerASCIISynchronizedPrintf(viewer, "[%d] ids[%" … in main()
101 PetscCall(ISLocalToGlobalMappingRestoreInfo(ltog, &np, &p, &ns, &ids)); in main()
104 PetscCall(ISLocalToGlobalMappingGetNodeInfo(ltog, &n, &ns, &ids)); in main()
107 …ronizedPrintf(viewer, "[%d] sharedby[%" PetscInt_FMT "] = %" PetscInt_FMT "\n", rank, i, ns[i])); in main()
108 …for (j = 0; j < ns[i]; j++) PetscCall(PetscViewerASCIISynchronizedPrintf(viewer, "[%d] ids[%" … in main()
111 PetscCall(ISLocalToGlobalMappingRestoreNodeInfo(ltog, &n, &ns, &ids)); in main()
116 PetscCall(ISLocalToGlobalMappingGetBlockInfo(ltog, &np, &p, &ns, &ids)); in main()
[all …]
/petsc/src/ts/trajectory/impls/basic/
H A Dtrajbasic.c15 PetscInt ns, i; in TSTrajectorySet_Basic() local
26 PetscCall(TSGetStages(ts, &ns, &Y)); in TSTrajectorySet_Basic()
27 for (i = 0; i < ns; i++) { in TSTrajectorySet_Basic()
29 if (ts->stifflyaccurate && i == ns - 1) continue; in TSTrajectorySet_Basic()
42 PetscCall(TSForwardGetStages(ts, &ns, &S)); in TSTrajectorySet_Basic()
43 for (i = 0; i < ns; i++) PetscCall(MatView(S[i], tjbasic->viewer)); in TSTrajectorySet_Basic()
62 PetscInt ns, i; in TSTrajectoryGet_Basic() local
74 PetscCall(TSGetStages(ts, &ns, &Y)); in TSTrajectoryGet_Basic()
75 for (i = 0; i < ns; i++) { in TSTrajectoryGet_Basic()
77 if (ts->stifflyaccurate && i == ns - 1) continue; in TSTrajectoryGet_Basic()
[all …]
/petsc/src/mat/tests/
H A Dex169.c17 PetscInt ns = 2; in main() local
39 PetscCall(PetscOptionsGetInt(NULL, NULL, "-nsub", &ns, NULL)); in main()
40 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "Splitting in %" PetscInt_FMT " subcommunicators\n", ns)); in main()
42 PetscCall(PetscSubcommSetNumber(subc, ns)); in main()
H A Dex154.c10 PetscInt ns = 0, m, n; in main() local
54 PetscCall(PetscOptionsGetInt(NULL, NULL, "-ns", &ns, NULL)); in main()
61 PetscCall(ISCreateStride(PETSC_COMM_SELF, ns, m - ns, 1, &is_schur)); in main()
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Dperf_profile.m28 [nf,np,ns] = size(H); % Grab the dimensions
31 for j = 1:ns
42 T = zeros(np,ns);
45 for s = 1:ns
62 r = T./repmat(min(T,[],2),1,ns);
70 hl = zeros(ns,1);
71 for s = 1:ns
/petsc/src/ksp/pc/impls/mg/
H A Dgdsw.c5 …WSetUp(PC pc, PetscInt l, DM dm, KSP smooth, PetscInt Nc, Mat A, PetscInt *ns, Mat **sA_IG_n, KSP … in PCMGGDSWSetUp() argument
50 *ns = 1; in PCMGGDSWSetUp()
51 PetscCall(PetscMalloc1(*ns, &sA_IG)); in PCMGGDSWSetUp()
52 PetscCall(PetscMalloc1(*ns, &sksp)); in PCMGGDSWSetUp()
53 PetscCall(PetscMalloc1(*ns, &sI)); in PCMGGDSWSetUp()
54 PetscCall(PetscMalloc1(*ns, &sG)); in PCMGGDSWSetUp()
55 PetscCall(PetscMalloc1(*ns, &sGf)); in PCMGGDSWSetUp()
56 PetscCall(PetscMalloc1(*ns, &sGi)); in PCMGGDSWSetUp()
57 PetscCall(PetscMalloc1(*ns, &sGiM)); in PCMGGDSWSetUp()
223 PetscInt NG, ns, n, i, c, rbs, cbs[2]; in PCMGGDSWCreateCoarseSpace_Private() local
[all …]
/petsc/src/dm/impls/da/
H A Ddadd.c35 PetscInt ms = 0, ns = 0, ps = 0; in DMDACreatePatchIS() local
104 ns = 0; in DMDACreatePatchIS()
153 while (jj > ne - 1 || jj < ns) { in DMDACreatePatchIS()
155 ns = 0; in DMDACreatePatchIS()
160 ns = ne; in DMDACreatePatchIS()
178 ym = ne - ns; in DMDACreatePatchIS()
180 base = ms * ym * zm + ns * M * zm + ps * M * N; in DMDACreatePatchIS()
183 sj = jj - ns; in DMDACreatePatchIS()
198 PetscCall(DMDAGetCorners(da, &ms, &ns, &ps, &mw, &nw, &pw)); in DMDACreatePatchIS()
200 if (N > 1) ne = ns + nw; in DMDACreatePatchIS()
[all …]
H A Ddagetelem.c42 PetscInt cnt = 0, cell[4], ns = 2; in DMDAGetElements_2D() local
63 if (da->elementtype == DMDA_ELEMENT_P1) ns = 2; in DMDAGetElements_2D()
64 if (da->elementtype == DMDA_ELEMENT_Q1) ns = 1; in DMDAGetElements_2D()
73 da->ne = ns * (xe - xs - 1) * (ye - ys - 1); in DMDAGetElements_2D()
82 for (c = 0; c < ns * nn; c++) da->e[cnt++] = cell[split[c]]; in DMDAGetElements_2D()
85 for (c = 0; c < ns * nn; c++) da->e[cnt++] = cell[c]; in DMDAGetElements_2D()
108 PetscInt cnt = 0, cell[8], ns = 6; in DMDAGetElements_3D() local
129 if (da->elementtype == DMDA_ELEMENT_P1) ns = 6; in DMDAGetElements_3D()
130 if (da->elementtype == DMDA_ELEMENT_Q1) ns = 1; in DMDAGetElements_3D()
142 da->ne = ns * (xe - xs - 1) * (ye - ys - 1) * (ze - zs - 1); in DMDAGetElements_3D()
[all …]
/petsc/src/vec/vec/tests/
H A Dex48.c78 #define ns 4 macro
80 static const char *datasets[ns] = {"", "x", "nonExistingVec", "y"};
340 Vec vecs[nap][ns]; in testGroupsDatasets()
355 PetscCall(PetscMemzero(vecs, nap * ns * sizeof(Vec))); in testGroupsDatasets()
371 for (s = 0; s < ns; s++) { in testGroupsDatasets()
433 for (s = 0; s < ns; s++) { in testGroupsDatasets()
462 for (s = 0; s < ns; s++) PetscCall(VecDestroy(&vecs[p][s])); in testGroupsDatasets()
490 Capsule capsules[nap][ns], c = NULL, old = NULL; in testAttributesAbsolutePath()
505 PetscCall(PetscMemzero(capsules, nap * ns * sizeof(Capsule))); in testAttributesAbsolutePath()
510 for (s = 0; s < ns; s++) { in testAttributesAbsolutePath()
[all …]
/petsc/share/petsc/matlab/
H A DgeneratePetscTestFiles.m57 PetscBinaryWrite('ns-real-int32-float64',A,b,'indices','int32','precision','float64');
58 [At,bt] = PetscBinaryRead('ns-real-int32-float64','indices','int32','precision','float64');
61 '%s','Error in PetscBinaryWrite or/and PetscBinaryRead for ns-real-int32-float64');
63 PetscBinaryWrite('ns-real-int64-float64',A,b,'indices','int64','precision','float64');
64 [At,bt] = PetscBinaryRead('ns-real-int64-float64','indices','int64','precision','float64');
67 '%s','Error in PetscBinaryWrite or/and PetscBinaryRead for ns-real-int64-float64');
69 PetscBinaryWrite('ns-real-int32-float32',A,b,'indices','int32','precision','float32');
70 [At,bt] = PetscBinaryRead('ns-real-int32-float32','indices','int32','precision','float32');
73 '%s','Error in PetscBinaryWrite or/and PetscBinaryRead for ns-real-int32-float32');
75 PetscBinaryWrite('ns-real-int64-float32',A,b,'indices','int64','precision','float32');
[all …]
/petsc/src/ts/impls/eimex/
H A Deimex.c35 const PetscInt ns = ext->nstages; in TSEvaluateStep_EIMEX() local
38 PetscCall(VecCopy(ext->T[Map(ext->row_ind, ext->col_ind, ns)], X)); in TSEvaluateStep_EIMEX()
76 const PetscInt ns = ext->nstages; in TSStep_EIMEX() local
91 for (j = 0; j < ns; j++) { in TSStep_EIMEX()
96 for (i = 1; i < ns; i++) { in TSStep_EIMEX()
97 for (j = i; j < ns; j++) { in TSStep_EIMEX()
99 …PetscCall(VecAXPBYPCZ(T[Map(j, i, ns)], alpha, 1.0, 0, T[Map(j, i - 1, ns)], T[Map(j - 1, i - 1, n… in TSStep_EIMEX()
101 PetscCall(VecScale(T[Map(j, i, ns)], alpha)); in TSStep_EIMEX()
105 PetscCall(TSEvaluateStep(ts, ns, ts->vec_sol, NULL)); /*update ts solution */ in TSStep_EIMEX()
177 PetscInt ns; in TSReset_EIMEX() local
[all …]
/petsc/src/tao/pde_constrained/tutorials/
H A Delliptic.c9 PetscInt ns; /* Number of data samples (1<=ns<=8) member
112 user.ns = 6; in main()
113 …PetscCall(PetscOptionsInt("-ns", "Number of data samples (1<=ns<=8)", "", user.ns, &user.ns, NULL)… in main()
130 user.m = user.ns * user.mx * user.mx * user.mx; /* number of constraints */ in main()
303 if (user->ns == 1) { in StateMatMult()
306 for (i = 0; i < user->ns; i++) { in StateMatMult()
330 if (user->ns == 1) { in StateInvMatMult()
335 for (i = 0; i < user->ns; i++) { in StateInvMatMult()
353 if (user->ns == 1) { in QMatMult()
356 for (i = 0; i < user->ns; i++) { in QMatMult()
[all …]
H A Dparabolic.c9 PetscInt ns; /* Number of samples */ member
116 user.ns = 8; in main()
117 PetscCall(PetscOptionsInt("-ns", "Number of samples", "", user.ns, &user.ns, NULL)); in main()
234 for (j = 0; j < user->ns; j++) { in FormFunction()
238 PetscCall(Gather_i(user->dwork, user->di, user->di_scatter, user->ns)); in FormFunction()
263 for (j = 0; j < user->ns; j++) { in FormGradient()
267 PetscCall(Gather_i(user->dwork, user->di, user->di_scatter, user->ns)); in FormGradient()
269 PetscCall(Scatter_i(user->dwork, user->di, user->di_scatter, user->ns)); in FormGradient()
272 for (j = 0; j < user->ns; j++) { in FormGradient()
295 for (j = 0; j < user->ns; j++) { in FormFunctionGradient()
[all …]
/petsc/src/mat/impls/aij/seq/
H A Dinode.c10 static PetscErrorCode MatCreateColInode_Private(Mat A, PetscInt *size, PetscInt **ns) in MatCreateColInode_Private() argument
22 if (!ns) { in MatCreateColInode_Private()
46 *ns = ns_col; in MatCreateColInode_Private()
362 PetscInt *ns; in MatMult_SeqAIJ_Inode() local
367 ns = a->inode.size_csr; /* Node Size array */ in MatMult_SeqAIJ_Inode()
381 row = ns[i]; in MatMult_SeqAIJ_Inode()
382 nsz = ns[i + 1] - ns[i]; in MatMult_SeqAIJ_Inode()
571 const PetscInt *idx, *ns, *ii; in MatMultAdd_SeqAIJ_Inode() local
576 ns = a->inode.size_csr; /* Node Size array */ in MatMultAdd_SeqAIJ_Inode()
587 row = ns[i]; in MatMultAdd_SeqAIJ_Inode()
[all …]
/petsc/src/vec/vec/utils/
H A Dvinv.c992 PetscInt i, n, bs, ns; in VecStrideGather_Default() local
998 PetscCall(VecGetLocalSize(s, &ns)); in VecStrideGather_Default()
1003 …ns * bs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Subvector length * blocksize %" PetscInt_FMT … in VecStrideGather_Default()
1024 PetscInt i, n, bs, ns; in VecStrideScatter_Default() local
1030 PetscCall(VecGetLocalSize(s, &ns)); in VecStrideScatter_Default()
1035 …ns * bs, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Subvector length * blocksize %" PetscInt_FMT … in VecStrideScatter_Default()
1056 PetscInt i, j, n, bs, bss, ns; in VecStrideSubSetGather_Default() local
1062 PetscCall(VecGetLocalSize(s, &ns)); in VecStrideSubSetGather_Default()
1071 PetscCheck(n == ns / bss, PETSC_COMM_SELF, PETSC_ERR_ARG_SIZ, "Incompatible layout of vectors"); in VecStrideSubSetGather_Default()
1120 PetscInt j, i, n, bs, ns, bss; in VecStrideSubSetScatter_Default() local
[all …]
/petsc/src/ksp/pc/impls/bddc/
H A Dbddcgraph.c301 PetscInt maxc, ns; in PCBDDCGraphComputeConnectedComponents() local
305 for (ns = 1, maxc = 0; ns < n_neigh; ns++) maxc = PetscMax(maxc, n_shared[ns]); in PCBDDCGraphComputeConnectedComponents()
309 for (ns = 1; ns < n_neigh; ns++) { /* first proc is self */ in PCBDDCGraphComputeConnectedComponents()
311 PetscInt fst, j, k, d, cdim = graph->cdim, n = n_shared[ns]; in PCBDDCGraphComputeConnectedComponents()
319 PetscReal c = graph->coords[shared[ns][j] * cdim + d]; in PCBDDCGraphComputeConnectedComponents()
405 PetscCall(PetscBTSet(cornerp, shared[ns][point1])); in PCBDDCGraphComputeConnectedComponents()
406 PetscCall(PetscBTSet(cornerp, shared[ns][point2])); in PCBDDCGraphComputeConnectedComponents()
407 PetscCall(PetscBTSet(cornerp, shared[ns][point3])); in PCBDDCGraphComputeConnectedComponents()
408 PetscCall(PetscBTSet(cornerp, shared[ns][point4])); in PCBDDCGraphComputeConnectedComponents()
421 if (same[0] || same[1] || same[2] || same[3]) PetscCall(PetscBTSet(cornerp, shared[ns][j])); in PCBDDCGraphComputeConnectedComponents()
/petsc/lib/petsc/bin/maint/abi-compliance-checker/modules/Internals/Styles/
H A DSymbolsList.css37 span.ns {
H A DReport.css67 span.ns {
/petsc/src/ksp/pc/impls/tfs/
H A Dxxt.c44 PetscInt ns; member
98 xxt_handle->ns = 0; in XXT_factor()
204 m = j - xxt_handle->ns; in xxt_generate()
205 …null space exists %" PetscInt_FMT " %" PetscInt_FMT " %" PetscInt_FMT "\n", m, j, xxt_handle->ns)); in xxt_generate()
H A Dxyt.c46 PetscInt ns; member
100 xyt_handle->ns = 0; in XYT_factor()
260 m = j - xyt_handle->ns; in xyt_generate()
261 …null space exists %" PetscInt_FMT " %" PetscInt_FMT " %" PetscInt_FMT "\n", m, j, xyt_handle->ns)); in xyt_generate()
/petsc/src/ts/impls/implicit/theta/
H A Dtheta.c867 static PetscErrorCode TSForwardGetStages_Theta(TS ts, PetscInt *ns, Mat *stagesensip[]) in TSForwardGetStages_Theta() argument
872 if (ns) { in TSForwardGetStages_Theta()
873 if (!th->endpoint && th->Theta != 1.0) *ns = 1; /* midpoint form */ in TSForwardGetStages_Theta()
874 else *ns = 2; /* endpoint form */ in TSForwardGetStages_Theta()
1164 static PetscErrorCode TSGetStages_Theta(TS ts, PetscInt *ns, Vec *Y[]) in TSGetStages_Theta() argument
1169 if (ns) { in TSGetStages_Theta()
1170 if (!th->endpoint && th->Theta != 1.0) *ns = 1; /* midpoint form */ in TSGetStages_Theta()
1171 else *ns = 2; /* endpoint form */ in TSGetStages_Theta()
/petsc/src/ts/impls/implicit/discgrad/
H A Dtsdiscgrad.c269 static PetscErrorCode TSGetStages_DiscGrad(TS ts, PetscInt *ns, Vec **Y) in TSGetStages_DiscGrad() argument
274 if (ns) *ns = 1; in TSGetStages_DiscGrad()
/petsc/src/ts/impls/explicit/rk/
H A Drk.c724 static PetscErrorCode TSForwardGetStages_RK(TS ts, PetscInt *ns, Mat **stagesensip) in TSForwardGetStages_RK() argument
730 if (ns) *ns = tab->s; in TSForwardGetStages_RK()
1345 static PetscErrorCode TSGetStages_RK(TS ts, PetscInt *ns, Vec **Y) in TSGetStages_RK() argument
1350 if (ns) *ns = rk->tableau->s; in TSGetStages_RK()
/petsc/src/ts/impls/glee/
H A Dglee.c914 static PetscErrorCode TSGetStages_GLEE(TS ts, PetscInt *ns, Vec **Y) in TSGetStages_GLEE() argument
919 if (ns) *ns = glee->tableau->s; in TSGetStages_GLEE()
/petsc/src/ts/utils/dmplexlandau/tutorials/
H A Dex1.c93 static PetscErrorCode SetMaxwellians(DM dm, Vec X, PetscReal time, PetscReal temps[], PetscReal ns[… in SetMaxwellians() argument
104 data[i0].n = ns[ii]; in SetMaxwellians()

12