| /petsc/src/snes/tutorials/output/ |
| H A D | ex12_p4est_fas_q2_nonconformal_serial.out | 6 …0 KSP preconditioned resid norm 2.499198301910e+02 true resid norm 3.216249248292e+04 ||r(i)||/||b…
7 … resid norm 8.614011391108e+01 true resid norm 1.952597363694e+04 ||r(i)||/||b|| 6.071038694314e-01
8 … resid norm 2.745949105361e+01 true resid norm 1.930071384658e+04 ||r(i)||/||b|| 6.001000655293e-01
9 … resid norm 2.559253117247e+01 true resid norm 1.639100195786e+04 ||r(i)||/||b|| 5.096309611751e-01
10 … resid norm 2.282154921037e+01 true resid norm 9.122321919069e+03 ||r(i)||/||b|| 2.836323062932e-01
11 … resid norm 9.015546483100e+00 true resid norm 2.894802732061e+03 ||r(i)||/||b|| 9.000554710113e-02
12 … resid norm 6.311361207223e+00 true resid norm 1.903443566068e+03 ||r(i)||/||b|| 5.918209128471e-02
13 … resid norm 3.316691807429e+00 true resid norm 7.510656749157e+02 ||r(i)||/||b|| 2.335222232277e-02
14 … resid norm 2.729920508905e+00 true resid norm 6.941012731568e+02 ||r(i)||/||b|| 2.158107844177e-02
15 … resid norm 1.834338478600e+00 true resid norm 2.992015148139e+02 ||r(i)||/||b|| 9.302808697827e-03
[all …]
|
| H A D | ex12_mg_newton_coarse_2.out | 4 Number of 0-cells per rank: 4225
5 Number of 1-cells per rank: 12416
6 Number of 2-cells per rank: 8192
12 Number of 0-cells per rank: 1224
13 Number of 1-cells per rank: 3543
14 Number of 2-cells per rank: 2320
19 Number of 0-cells per rank: 408
20 Number of 1-cells per rank: 1154
21 Number of 2-cells per rank: 747
26 Number of 0-cells per rank: 136
[all …]
|
| H A D | ex12_mg_newton_coarse_1.out | 4 Number of 0-cells per rank: 4225
5 Number of 1-cells per rank: 12416
6 Number of 2-cells per rank: 8192
12 Number of 0-cells per rank: 1230
13 Number of 1-cells per rank: 3555
14 Number of 2-cells per rank: 2326
19 Number of 0-cells per rank: 414
20 Number of 1-cells per rank: 1171
21 Number of 2-cells per rank: 758
26 Number of 0-cells per rank: 144
[all …]
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexcheckinterface.c | 8 PetscMPIInt r; in ExchangeArrayByRank_Private() local
22 …for (r = 0; r < nrranks; r++) PetscCallMPI(MPIU_Irecv(&rsize[r], 1, MPIU_INT, rranks[r], tag, comm… in ExchangeArrayByRank_Private()
23 …for (r = 0; r < nsranks; r++) PetscCallMPI(MPIU_Isend(&ssize[r], 1, MPIU_INT, sranks[r], tag, comm… in ExchangeArrayByRank_Private()
28 for (r = 0; r < nrranks; r++) { in ExchangeArrayByRank_Private()
29 PetscCall(PetscMalloc(rsize[r] * unitsize, &rarr[r])); in ExchangeArrayByRank_Private()
30 PetscCallMPI(MPIU_Irecv(rarr[r], rsize[r], dt, rranks[r], tag, comm, &rreq[r])); in ExchangeArrayByRank_Private()
32 …for (r = 0; r < nsranks; r++) PetscCallMPI(MPIU_Isend(sarr[r], ssize[r], dt, sranks[r], tag, comm,… in ExchangeArrayByRank_Private()
45 PetscMPIInt r; in ExchangeVecByRank_Private() local
54 for (r = 0; r < nsranks; r++) { in ExchangeVecByRank_Private()
55 PetscCall(VecGetLocalSize(svecs[r], &ssize[r])); in ExchangeVecByRank_Private()
[all …]
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex18_2a.out | 4 [0] --------
5 r=0 ranks[r]=1 sntCoordinatesPerRank[r]:
6 [1] --------
7 r=0 ranks[r]=2 sntCoordinatesPerRank[r]:
8 [2] --------
27 r=1 ranks[r]=2 sntCoordinatesPerRank[r]:
38 ----------
39 ----------
40 r=0 iranks[r]=0 refCoordinatesPerRank[r]:
41 ----------
[all …]
|
| /petsc/src/sys/classes/random/interface/ |
| H A D | random.c | 17 PetscRandomGetValue - Generates a random number. Call this after first calling
23 . r - the random number generator context
26 . val - the value
42 PetscRandomCreate(PETSC_COMM_WORLD,&r);
43 PetscRandomGetValue(r,&value1);
44 PetscRandomGetValue(r,&value2);
45 PetscRandomGetValue(r,&value3);
46 PetscRandomDestroy(&r);
51 PetscErrorCode PetscRandomGetValue(PetscRandom r, PetscScalar *val) in PetscRandomGetValue() argument
54 PetscValidHeaderSpecific(r, PETSC_RANDOM_CLASSID, 1); in PetscRandomGetValue()
[all …]
|
| H A D | randomc.c | 21 PetscRandomDestroy - Destroys a `PetscRandom` object that was created by `PetscRandomCreate()`.
26 . r - the random number generator object
32 PetscErrorCode PetscRandomDestroy(PetscRandom *r) in PetscRandomDestroy() argument
35 if (!*r) PetscFunctionReturn(PETSC_SUCCESS); in PetscRandomDestroy()
36 PetscValidHeaderSpecific(*r, PETSC_RANDOM_CLASSID, 1); in PetscRandomDestroy()
37 if (--((PetscObject)*r)->refct > 0) { in PetscRandomDestroy()
38 *r = NULL; in PetscRandomDestroy()
41 PetscTryTypeMethod(*r, destroy); in PetscRandomDestroy()
42 PetscCall(PetscHeaderDestroy(r)); in PetscRandomDestroy()
47 PetscRandomGetSeed - Gets the random seed.
[all …]
|
| /petsc/src/ts/utils/ |
| H A D | tsconvest.c | 12 PetscCall(PetscObjectGetClassId(ce->solver, &id)); in PetscConvEstSetTS_Private()
14 PetscCall(TSGetDM((TS)ce->solver, &ce->idm)); in PetscConvEstSetTS_Private()
18 static PetscErrorCode PetscConvEstInitGuessTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u) in PetscConvEstInitGuessTS_Private() argument
21 PetscCall(TSComputeInitialCondition((TS)ce->solver, u)); in PetscConvEstInitGuessTS_Private()
25 static PetscErrorCode PetscConvEstComputeErrorTS_Private(PetscConvEst ce, PetscInt r, DM dm, Vec u,… in PetscConvEstComputeErrorTS_Private() argument
27 TS ts = (TS)ce->solver; in PetscConvEstComputeErrorTS_Private()
39 for (f = 1; f < ce->Nf; ++f) errors[f] = errors[0]; in PetscConvEstComputeErrorTS_Private()
45 PetscCall(DMComputeL2FieldDiff(dm, t, ce->exactSol, ce->ctxs, u, errors)); in PetscConvEstComputeErrorTS_Private()
52 TS ts = (TS)ce->solver; in PetscConvEstGetConvRateTS_Temporal_Private()
55 PetscInt Ns, oNs, Nf = ce->Nf, f, Nr = ce->Nr, r; in PetscConvEstGetConvRateTS_Temporal_Private() local
[all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | matrart.c | 3 C = R * A * R^T
15 PetscCall(MatTransposeColoringDestroy(&rart->matcoloring)); in MatProductCtxDestroy_SeqAIJ_RARt()
16 PetscCall(MatDestroy(&rart->Rt)); in MatProductCtxDestroy_SeqAIJ_RARt()
17 PetscCall(MatDestroy(&rart->RARt)); in MatProductCtxDestroy_SeqAIJ_RARt()
18 PetscCall(MatDestroy(&rart->ARt)); in MatProductCtxDestroy_SeqAIJ_RARt()
19 PetscCall(PetscFree(rart->work)); in MatProductCtxDestroy_SeqAIJ_RARt()
20 if (rart->destroy) PetscCall((*rart->destroy)(&rart->data)); in MatProductCtxDestroy_SeqAIJ_RARt()
25 PetscErrorCode MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart(Mat A, Mat R, PetscReal fill, Mat C) in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart() argument
36 …PetscCheck(!C->product->data, PetscObjectComm((PetscObject)C), PETSC_ERR_PLIB, "Product data not e… in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
38 PetscCall(MatTransposeSymbolic(R, &P)); in MatRARtSymbolic_SeqAIJ_SeqAIJ_colorrart()
[all …]
|
| /petsc/src/sys/classes/random/impls/sprng/ |
| H A D | sprng.c | 9 static PetscErrorCode PetscRandomSeed_Sprng(PetscRandom r) in PetscRandomSeed_Sprng() argument
12 init_sprng(r->seed, SPRNG_DEFAULT); in PetscRandomSeed_Sprng()
16 static PetscErrorCode PetscRandomGetValue_Sprng(PetscRandom r, PetscScalar *val) in PetscRandomGetValue_Sprng() argument
20 if (r->iset) { in PetscRandomGetValue_Sprng()
21 …*val = PetscRealPart(r->width) * sprng() + PetscRealPart(r->low) + (PetscImaginaryPart(r->width) *… in PetscRandomGetValue_Sprng()
26 if (r->iset) *val = r->width * sprng() + r->low; in PetscRandomGetValue_Sprng()
32 static PetscErrorCode PetscRandomGetValueReal_Sprng(PetscRandom r, PetscReal *val) in PetscRandomGetValueReal_Sprng() argument
36 if (r->iset) *val = PetscRealPart(r->width) * sprng() + PetscRealPart(r->low); in PetscRandomGetValueReal_Sprng()
39 if (r->iset) *val = r->width * sprng() + r->low; in PetscRandomGetValueReal_Sprng()
52 PETSCSPRNG - access to the publicly available random number generator sprng
[all …]
|
| /petsc/src/sys/classes/random/impls/rand/ |
| H A D | rand.c | 3 static PetscErrorCode PetscRandomSeed_Rand(PetscRandom r) in PetscRandomSeed_Rand() argument
6 srand((unsigned int)r->seed); in PetscRandomSeed_Rand()
11 static PetscErrorCode PetscRandomGetValue_Rand(PetscRandom r, PetscScalar *val) in PetscRandomGetValue_Rand() argument
15 …if (r->iset) *val = PetscRealPart(r->width) * RAND_WRAP + PetscRealPart(r->low) + (PetscImaginaryP… in PetscRandomGetValue_Rand()
18 if (r->iset) *val = r->width * RAND_WRAP + r->low; in PetscRandomGetValue_Rand()
24 static PetscErrorCode PetscRandomGetValueReal_Rand(PetscRandom r, PetscReal *val) in PetscRandomGetValueReal_Rand() argument
28 if (r->iset) *val = PetscRealPart(r->width) * RAND_WRAP + PetscRealPart(r->low); in PetscRandomGetValueReal_Rand()
31 if (r->iset) *val = r->width * RAND_WRAP + r->low; in PetscRandomGetValueReal_Rand()
48 PETSCRAND - access to the basic Unix random number generator
51 . -random_type <rand,rand48,sprng> - set the random number generator from the options database
[all …]
|
| /petsc/src/sys/classes/random/impls/rand48/ |
| H A D | rand48.c | 4 static PetscErrorCode PetscRandomSeed_Rand48(PetscRandom r) in PetscRandomSeed_Rand48() argument
7 srand48(r->seed); in PetscRandomSeed_Rand48()
11 static PetscErrorCode PetscRandomGetValue_Rand48(PetscRandom r, PetscScalar *val) in PetscRandomGetValue_Rand48() argument
15 if (r->iset) { in PetscRandomGetValue_Rand48()
16 …RealPart(r->width) * (PetscReal)drand48() + PetscRealPart(r->low) + (PetscImaginaryPart(r->width) … in PetscRandomGetValue_Rand48()
21 if (r->iset) *val = r->width * drand48() + r->low; in PetscRandomGetValue_Rand48()
27 static PetscErrorCode PetscRandomGetValueReal_Rand48(PetscRandom r, PetscReal *val) in PetscRandomGetValueReal_Rand48() argument
31 if (r->iset) *val = PetscRealPart(r->width) * drand48() + PetscRealPart(r->low); in PetscRandomGetValueReal_Rand48()
34 if (r->iset) *val = r->width * (PetscReal)drand48() + r->low; in PetscRandomGetValueReal_Rand48()
51 PETSCRAND48 - access to the basic Unix `drand48()` random number generator
[all …]
|
| /petsc/src/sys/classes/random/impls/random123/ |
| H A D | random123.c | 19 static PetscErrorCode PetscRandomSeed_Random123(PetscRandom r) in PetscRandomSeed_Random123() argument
22 PetscRandom123 *r123 = (PetscRandom123 *)r->data; in PetscRandomSeed_Random123()
25 ukey.v[0] = (R123_ULONG_LONG)r->seed; in PetscRandomSeed_Random123()
31 r123->key = threefry4x64keyinit(ukey); in PetscRandomSeed_Random123()
32 r123->counter.v[0] = 0; in PetscRandomSeed_Random123()
33 r123->counter.v[1] = 1; in PetscRandomSeed_Random123()
34 r123->counter.v[2] = 2; in PetscRandomSeed_Random123()
35 r123->counter.v[3] = 3; in PetscRandomSeed_Random123()
36 r123->result = threefry4x64(r123->counter, r123->key); in PetscRandomSeed_Random123()
37 r123->count = 0; in PetscRandomSeed_Random123()
[all …]
|
| /petsc/src/sys/classes/random/impls/rander48/ |
| H A D | rander48.c | 22 accu = (unsigned long)r48->mult[0] * (unsigned long)r48->seed[0] + (unsigned long)r48->add; in _dorander48()
25 …accu += (unsigned long)r48->mult[0] * (unsigned long)r48->seed[1] + (unsigned long)r48->mult[1] * … in _dorander48()
28 …u += (unsigned long)r48->mult[0] * r48->seed[2] + (unsigned long)r48->mult[1] * r48->seed[1] + (un… in _dorander48()
29 r48->seed[0] = temp[0]; in _dorander48()
30 r48->seed[1] = temp[1]; in _dorander48()
31 r48->seed[2] = (unsigned short)accu; in _dorander48()
32 …return ldexp((double)r48->seed[0], -48) + ldexp((double)r48->seed[1], -32) + ldexp((double)r48->se… in _dorander48()
35 static PetscErrorCode PetscRandomSeed_Rander48(PetscRandom r) in PetscRandomSeed_Rander48() argument
37 PetscRandom_Rander48 *r48 = (PetscRandom_Rander48 *)r->data; in PetscRandomSeed_Rander48()
40 r48->seed[0] = RANDER48_SEED_0; in PetscRandomSeed_Rander48()
[all …]
|
| /petsc/src/sys/tests/linter/output/ |
| H A D | testTypes.patch | 1 @@ -4,5 +4,5 @@
4 -PetscErrorCode testTypes(PetscRandom r, PetscViewer v, PetscObject o, PetscFunctionList f)
5 +static PetscErrorCode testTypes(PetscRandom r, PetscViewer v, PetscObject o, PetscFunctionList f)
8 @@ -7,5 +7,5 @@
11 - PetscValidType(r, -1);
12 + PetscValidType(r, 1);
13 PetscCheckSameType(r, -1, v, -1);
14 PetscCheckSameComm(o, -2, f, -2);
15 @@ -8,5 +8,5 @@
17 PetscValidType(r, -1);
[all …]
|
| /petsc/src/dm/impls/da/tests/output/ |
| H A D | ex1_0.out | 1 [0] p 0 (+0.4619,+0.1913) r=0.50 th=0.125*pi cellid=0
2 [0] p 0 (+0.1913,+0.4619) r=0.50 th=0.375*pi cellid=0
3 [0] p 0 (-0.1913,+0.4619) r=0.50 th=0.625*pi cellid=0
4 [0] p 0 (-0.4619,+0.1913) r=0.50 th=0.875*pi cellid=0
5 [0] p 0 (-0.4619,-0.1913) r=0.50 th=-0.875*pi cellid=0
6 [0] p 0 (-0.1913,-0.4619) r=0.50 th=-0.625*pi cellid=0
7 [0] p 0 (+0.1913,-0.4619) r=0.50 th=-0.375*pi cellid=0
8 [0] p 0 (+0.4619,-0.1913) r=0.50 th=-0.125*pi cellid=0
12 [0] p 0 (+0.4619,+0.1913) r=0.50 th=0.125*pi cellid=5
13 [0] p 0 (+0.1913,+0.4619) r=0.50 th=0.375*pi cellid=7
[all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | mesh-3d-box-innersphere_bin.msh | 11 …���i2E>ſ^s�̠��� m�ȿ+��������~��� x� ᧿�u����?,���x�ҝ(��1���������7���ο-����� =��_?����S2ϿGC�…
12 �*��?�Q���k�?U L�\ο8���_[f���?�R�}��˿��a��ͤ�9����O�$A��?(Jeܱ��?gz�j-E�?:�������ƿ �6|zƿ��
13 ��;�����*��ɿ��~8���?[�Y�h�?<����̒Ր�? ��M��ƿf��ږ�?=���L�Ɛ��ÿ%|��v��?�-Q�`��?>���+�
17 …-HD��?9�hI�?�?J���n��ދ�����:�?�?�ZS� οK�����_M'ڶ�q� ��[�?��t;��?L���<Mv �d�?�%�����?��.��?M��…
|
| H A D | gmsh-3d-binary-32.msh | 24 ���-���.���/���0���1���2���3���4���5���6���7���8���9�����������a������?a������?��������@������? ���…
32 ���G���H���I���J���K���L���M���N���O���P���Q���R���S��� ������?��������@������?@������?�������� ���…
44 ���n���o���p���q���r���s���t���u���v���w���x���y���z��� ������?@������?�������?������?������?����…
50 …���P��� ������!���P��������� ���Q������������Q���'���������R������������R���#�������…
54 …R������������O���������K���O������������S���K������"���L���!������!���L���I������I���P�…
56 …U���_���W���O���V���_���U���P���U���Z���V������2������(���Q������ ���i���R���������i���S��� �…
60 …r���b���a���f���s���b���d���e���t���d���b���h���u���c���b���f���v���b���c���g���w���k���b���g���x�…
|
| /petsc/src/sys/classes/draw/utils/ |
| H A D | cmap.c | 25 while (h > 1.0) h -= 1.0; in PetscHlsHelper()
27 if (h < 1 / 6.0) return m1 + (m2 - m1) * h * 6; in PetscHlsHelper()
29 if (h < 2 / 3.0) return m1 + (m2 - m1) * (2 / 3.0 - h) * 6; in PetscHlsHelper()
33 static inline void PetscHlsToRgb(double h, double l, double s, double *r, double *g, double *b) in PetscHlsToRgb() argument
36 double m2 = l <= 0.5 ? l * (1.0 + s) : l + s - (l * s); in PetscHlsToRgb()
37 double m1 = 2 * l - m2; in PetscHlsToRgb()
38 *r = PetscHlsHelper(m1, m2, h + 1 / 3.); in PetscHlsToRgb()
40 *b = PetscHlsHelper(m1, m2, h - 1 / 3.); in PetscHlsToRgb()
43 *r = *g = *b = l; in PetscHlsToRgb()
47 static inline void PetscGammaCorrect(double *r, double *g, double *b) in PetscGammaCorrect() argument
[all …]
|
| /petsc/src/ksp/ksp/impls/cg/groppcg/ |
| H A D | groppcg.c | 4 KSPSetUp_GROPPCG - Sets up the workspace needed by the GROPPCG method.
20 . ksp - the Krylov space object that was set to use conjugate gradient, by, for
28 Vec x, b, r, p, s, S, z, Z; in KSPSolve_GROPPCG() local
33 PetscCall(PCGetDiagonalScale(ksp->pc, &diagonalscale)); in KSPSolve_GROPPCG()
34 …ETSC_ERR_SUP, "Krylov method %s does not support diagonal scaling", ((PetscObject)ksp)->type_name); in KSPSolve_GROPPCG()
36 x = ksp->vec_sol; in KSPSolve_GROPPCG()
37 b = ksp->vec_rhs; in KSPSolve_GROPPCG()
38 r = ksp->work[0]; in KSPSolve_GROPPCG()
39 p = ksp->work[1]; in KSPSolve_GROPPCG()
40 s = ksp->work[2]; in KSPSolve_GROPPCG()
[all …]
|
| /petsc/src/sys/classes/random/impls/curand/ |
| H A D | curand.c | 9 static PetscErrorCode PetscRandomSeed_CURAND(PetscRandom r) in PetscRandomSeed_CURAND() argument
11 PetscRandom_CURAND *curand = (PetscRandom_CURAND *)r->data; in PetscRandomSeed_CURAND()
14 PetscCallCURAND(curandSetPseudoRandomGeneratorSeed(curand->gen, r->seed)); in PetscRandomSeed_CURAND()
20 static PetscErrorCode PetscRandomGetValuesReal_CURAND(PetscRandom r, PetscInt n, PetscReal *val) in PetscRandomGetValuesReal_CURAND() argument
22 PetscRandom_CURAND *curand = (PetscRandom_CURAND *)r->data; in PetscRandomGetValuesReal_CURAND()
23 size_t nn = n < 0 ? (size_t)(-2 * n) : (size_t)n; /* handle complex case */ in PetscRandomGetValuesReal_CURAND()
27 PetscCallCURAND(curandGenerateUniform(curand->gen, val, nn)); in PetscRandomGetValuesReal_CURAND()
29 PetscCallCURAND(curandGenerateUniformDouble(curand->gen, val, nn)); in PetscRandomGetValuesReal_CURAND()
31 if (r->iset) PetscCall(PetscRandomCurandScale_Private(r, nn, val, (PetscBool)(n < 0))); in PetscRandomGetValuesReal_CURAND()
35 static PetscErrorCode PetscRandomGetValues_CURAND(PetscRandom r, PetscInt n, PetscScalar *val) in PetscRandomGetValues_CURAND() argument
[all …]
|
| /petsc/src/mat/utils/ |
| H A D | zerorows.c | 7 PetscInt *owners = A->rmap->range; in MatZeroRowsMapLocal_Private()
8 PetscInt n = A->rmap->n; in MatZeroRowsMapLocal_Private()
13 PetscInt r, len = 0; in MatZeroRowsMapLocal_Private() local
19 for (r = 0; r < n; ++r) lrows[r] = -1; in MatZeroRowsMapLocal_Private()
20 if (!A->nooffproczerorows) PetscCall(PetscMalloc1(N, &rrows)); in MatZeroRowsMapLocal_Private()
21 for (r = 0; r < N; ++r) { in MatZeroRowsMapLocal_Private()
22 const PetscInt idx = rows[r]; in MatZeroRowsMapLocal_Private()
23 … >= 0 && A->rmap->N > idx, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Row %" PetscInt_FMT " out o… in MatZeroRowsMapLocal_Private()
24 …if (idx < owners[p] || owners[p + 1] <= idx) { /* short-circuit the search if the last p owns this… in MatZeroRowsMapLocal_Private()
25 PetscCall(PetscLayoutFindOwner(A->rmap, idx, &p)); in MatZeroRowsMapLocal_Private()
[all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_ksp_py.py | 1 # --------------------------------------------------------------------
7 # --------------------------------------------------------------------
29 def loop(self, ksp, r): argument
31 rnorm = r.norm()
47 r, z = self.work
49 A.mult(x, r)
50 r.aypx(-1, b)
51 P.apply(r, z)
54 A.mult(x, r)
55 r.aypx(-1, b)
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex87.c | 1 static char help[] = "Block-structured Nest matrix involving a HermitianTranspose block.\n\n"
3 " -n <n>, where <n> = dimension of the blocks.\n\n";
10 H = [ R C
11 -C^H -R^T ],
13 where R is Hermitian and C is complex symmetric. In particular, R and C have the
16 R = pentadiag{a,b,c,conj(b),conj(a)}
24 Mat H, R, C, block[4]; in main() local
25 Vec rhs, x, r; in main() local
37 PetscCall(PetscOptionsGetInt(NULL, NULL, "-n", &n, NULL)); in main()
38 PetscCall(PetscOptionsGetInt(NULL, NULL, "-M", &M, NULL)); in main()
[all …]
|
| /petsc/src/snes/utils/ |
| H A D | convest.c | 16 PetscConvEstDestroy - Destroys a PETSc convergence estimator `PetscConvEst` object
21 . ce - The `PetscConvEst` object
32 if (--((PetscObject)*ce)->refct > 0) { in PetscConvEstDestroy()
36 PetscCall(PetscFree3((*ce)->initGuess, (*ce)->exactSol, (*ce)->ctxs)); in PetscConvEstDestroy()
37 PetscCall(PetscFree2((*ce)->dofs, (*ce)->errors)); in PetscConvEstDestroy()
43 …PetscConvEstSetFromOptions - Sets a convergence estimator `PetscConvEst` object based on values in…
48 . ce - The `PetscConvEst` object
58 …scCall(PetscOptionsInt("-convest_num_refine", "The number of refinements for the convergence check… in PetscConvEstSetFromOptions()
59 …scCall(PetscOptionsReal("-convest_refine_factor", "The increase in resolution in each dimension", … in PetscConvEstSetFromOptions()
60 …tscCall(PetscOptionsBool("-convest_monitor", "Monitor the error for each convergence check", "Pets… in PetscConvEstSetFromOptions()
[all …]
|