| /petsc/src/tao/leastsquares/tutorials/ |
| H A D | chwirut1f.F90 | 42 PetscInt i
52 do i = 0, m - 1
53 f_v(i + 1) = y(i) - exp(-x_v(1)*t(i))/(x_v(2) + x_v(3)*t(i))
75 integer i local
76 i = 0
77 y(i) = 92.9000; t(i) = 0.5000; i = i + 1
78 y(i) = 78.7000; t(i) = 0.6250; i = i + 1
79 y(i) = 64.2000; t(i) = 0.7500; i = i + 1
80 y(i) = 64.9000; t(i) = 0.8750; i = i + 1
81 y(i) = 57.1000; t(i) = 1.0000; i = i + 1
[all …]
|
| H A D | chwirut2f.F90 | 26 PetscInt i
27 i = 0
28 y(i) = 92.9000; t(i) = 0.5000; i = i + 1
29 y(i) = 78.7000; t(i) = 0.6250; i = i + 1
30 y(i) = 64.2000; t(i) = 0.7500; i = i + 1
31 y(i) = 64.9000; t(i) = 0.8750; i = i + 1
32 y(i) = 57.1000; t(i) = 1.0000; i = i + 1
33 y(i) = 43.3000; t(i) = 1.2500; i = i + 1
34 y(i) = 31.1000; t(i) = 1.7500; i = i + 1
35 y(i) = 23.6000; t(i) = 2.2500; i = i + 1
[all …]
|
| H A D | chwirut1.c | 51 PetscInt i; /* iteration information */ in main() local
65 for (i = 0; i < NOBSERVATIONS; i++) user.idm[i] = i; in main()
67 for (i = 0; i < NPARAMETERS; i++) user.idn[i] = i; in main()
106 PetscInt i; in EvaluateFunction() local
114 …for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[… in EvaluateFunction()
126 PetscInt i; in EvaluateJacobian() local
133 for (i = 0; i < NOBSERVATIONS; i++) { in EvaluateJacobian()
134 base = PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]); in EvaluateJacobian()
136 user->j[i][0] = t[i] * base; in EvaluateJacobian()
137 user->j[i][1] = base / (x[1] + x[2] * t[i]); in EvaluateJacobian()
[all …]
|
| H A D | chwirut2.c | 46 PetscErrorCode RunSimulation(PetscReal *x, PetscInt i, PetscReal *f, AppCtx *user);
103 PetscInt i; in EvaluateFunction() local
111 for (i = 0; i < NOBSERVATIONS; i++) PetscCall(RunSimulation(x, i, &f[i], user)); in EvaluateFunction()
167 PetscInt i = 0; in InitializeData() local
170 y[i] = 92.9000; in InitializeData()
171 t[i++] = 0.5000; in InitializeData()
172 y[i] = 78.7000; in InitializeData()
173 t[i++] = 0.6250; in InitializeData()
174 y[i] = 64.2000; in InitializeData()
175 t[i++] = 0.7500; in InitializeData()
[all …]
|
| /petsc/src/dm/tutorials/output/ |
| H A D | ex5_1.out | 1 i 0 e[3*i] 0 1 10
2 i 1 e[3*i] 11 10 1
3 i 2 e[3*i] 1 2 11
4 i 3 e[3*i] 12 11 2
5 i 4 e[3*i] 2 3 12
6 i 5 e[3*i] 13 12 3
7 i 6 e[3*i] 3 4 13
8 i 7 e[3*i] 14 13 4
9 i 8 e[3*i] 4 5 14
10 i 9 e[3*i] 15 14 5
[all …]
|
| H A D | ex5_2.out | 1 i 0 e[3*i] 0 1 6
2 i 1 e[3*i] 7 6 1
3 i 2 e[3*i] 1 2 7
4 i 3 e[3*i] 8 7 2
5 i 4 e[3*i] 2 3 8
6 i 5 e[3*i] 9 8 3
7 i 6 e[3*i] 3 4 9
8 i 7 e[3*i] 10 9 4
9 i 8 e[3*i] 6 7 12
10 i 9 e[3*i] 13 12 7
[all …]
|
| /petsc/src/tao/leastsquares/tests/ |
| H A D | chwirut2.c | 54 PetscInt i; /* iteration information */ in main() local
75 for (i = 0; i < NOBSERVATIONS; i++) { in main()
76 w_row[i] = i; in main()
77 w_col[i] = i; in main()
78 w_vals[i] = 1.0; in main()
84 for (i = 0; i < NOBSERVATIONS; i++) user.idm[i] = i; in main()
86 for (i = 0; i < NPARAMETERS; i++) user.idn[i] = i; in main()
129 PetscInt i; in EvaluateFunction() local
137 …for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[… in EvaluateFunction()
149 PetscInt i; in EvaluateJacobian() local
[all …]
|
| /petsc/src/ksp/pc/tests/output/ |
| H A D | ex6_1.out | 2 i = 0,flg = 1
3 i = 1,flg = 1
4 i = 2,flg = 1
5 i = 3,flg = 1
6 i = 4,flg = 1
7 i = 5,flg = 1
8 i = 6,flg = 1
9 i = 7,flg = 1
10 i = 8,flg = 1
11 i = 9,flg = 1
[all …]
|
| /petsc/src/dm/impls/stag/tests/ |
| H A D | ex2.c | 11 PetscInt dim, dof[4], i, elx, ely, elz; in main() local
16 for (i = 0; i < 4; ++i) dof[i] = 1; in main()
36 PetscInt i, low, high, n; in Test_3d_4x4x4_3x3x3() local
46 …for (i = 0; i < 3; ++i) PetscCheck(nRanks[i] == 3, PetscObjectComm((PetscObject)dmstag), PETSC_ERR… in Test_3d_4x4x4_3x3x3()
48 …for (i = 0; i < 3; ++i) PetscCheck(n[i] == 4, PetscObjectComm((PetscObject)dmstag), PETSC_ERR_SUP,… in Test_3d_4x4x4_3x3x3()
56 for (i = 0; i < n; ++i) arr[i] = (PetscScalar)(i + low); in Test_3d_4x4x4_3x3x3()
86 for (i = 0; i < entriesGhost; ++i) arrLocalExpected[i] = arrLocalExpectedHere[i]; in Test_3d_4x4x4_3x3x3()
95 for (i = 0; i < entriesGhost; ++i) arrLocalExpected[i] = arrLocalExpectedHere[i]; in Test_3d_4x4x4_3x3x3()
104 for (i = 0; i < entriesGhost; ++i) arrLocalExpected[i] = arrLocalExpectedHere[i]; in Test_3d_4x4x4_3x3x3()
113 for (i = 0; i < entriesGhost; ++i) arrLocalExpected[i] = arrLocalExpectedHere[i]; in Test_3d_4x4x4_3x3x3()
[all …]
|
| /petsc/src/ts/tutorials/multirate/ |
| H A D | finitevolume1d.c | 42 PetscInt i; in Limit_Upwind() local
43 for (i = 0; i < info->m; i++) lmt[i] = 0; in Limit_Upwind()
47 PetscInt i; in Limit_LaxWendroff() local
48 for (i = 0; i < info->m; i++) lmt[i] = jR[i]; in Limit_LaxWendroff()
52 PetscInt i; in Limit_BeamWarming() local
53 for (i = 0; i < info->m; i++) lmt[i] = jL[i]; in Limit_BeamWarming()
57 PetscInt i; in Limit_Fromm() local
58 for (i = 0; i < info->m; i++) lmt[i] = 0.5 * (jL[i] + jR[i]); in Limit_Fromm()
62 PetscInt i; in Limit_Minmod() local
63 for (i = 0; i < info->m; i++) lmt[i] = MinMod2(jL[i], jR[i]); in Limit_Minmod()
[all …]
|
| /petsc/src/ts/tutorials/phasefield/ |
| H A D | biharmonic.c | 200 PetscInt i, Mx, xs, xm; in FormFunction() local
238 for (i = xs; i < xs + xm; i++) { in FormFunction()
240 c = (1. - x[i] * x[i]) * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
241 r = (1. - x[i + 1] * x[i + 1]) * (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
242 l = (1. - x[i - 1] * x[i - 1]) * (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
244 c = (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
245 r = (x[i] + x[i + 2] - 2.0 * x[i + 1]) * sx; in FormFunction()
246 l = (x[i - 2] + x[i] - 2.0 * x[i - 1]) * sx; in FormFunction()
248 f[i] = -ctx->kappa * (l + r - 2.0 * c) * sx; in FormFunction()
252 …f[i] += 6. * .25 * x[i] * (x[i + 1] - x[i - 1]) * (x[i + 1] - x[i - 1]) * sx + (3. * x[i] * x[i] -… in FormFunction()
[all …]
|
| /petsc/src/benchmarks/ |
| H A D | Index.c | 21 int i, ierr, *z, *zi, intval; in test1() local
37 for (i = 0; i < 2000; i++) { in test1()
40 z[i] = intval; in test1()
43 for (i = 0; i < 2000; i++) { in test1()
46 zi[i] = intval; in test1()
53 for (i = 0; i < 2000; i++) x[i] = y[i]; in test1()
60 for (i = 0; i < 500; i += 4) { in test1()
61 x[i] = y[z[i]]; in test1()
62 x[1 + i] = y[z[1 + i]]; in test1()
63 x[2 + i] = y[z[2 + i]]; in test1()
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex64.c | 8 PetscInt i, j; in main() local
27 for (i = 0; i < 16; i++) values[i] = i; in main()
28 for (i = 0; i < 4; i++) values[4 * i + i] += 5; in main()
29 i = 0; in main()
31 PetscCall(MatSetValuesBlocked(A, 1, &i, 1, &j, values, INSERT_VALUES)); in main()
33 for (i = 0; i < 16; i++) values[i] = i; in main()
34 i = 0; in main()
36 PetscCall(MatSetValuesBlocked(A, 1, &i, 1, &j, values, INSERT_VALUES)); in main()
38 for (i = 0; i < 16; i++) values[i] = i; in main()
39 i = 1; in main()
[all …]
|
| /petsc/src/ksp/pc/impls/mg/ |
| H A D | smg.c | 8 PetscInt i, l = mglevels[0]->levels; in PCMGACycle_Private() local
12 for (i = l - 1; i > 0; i--) { in PCMGACycle_Private()
13 …if (mglevels[i]->eventinterprestrict) PetscCall(PetscLogEventBegin(mglevels[i]->eventinterprestric… in PCMGACycle_Private()
15 … if (matapp) PetscCall(MatMatRestrict(mglevels[i]->restrct, mglevels[i]->B, &mglevels[i - 1]->B)); in PCMGACycle_Private()
16 else PetscCall(MatRestrict(mglevels[i]->restrct, mglevels[i]->b, mglevels[i - 1]->b)); in PCMGACycle_Private()
18 …if (matapp) PetscCall(MatMatRestrict(mglevels[i]->interpolate, mglevels[i]->B, &mglevels[i - 1]->B… in PCMGACycle_Private()
19 else PetscCall(MatRestrict(mglevels[i]->interpolate, mglevels[i]->b, mglevels[i - 1]->b)); in PCMGACycle_Private()
21 …if (mglevels[i]->eventinterprestrict) PetscCall(PetscLogEventEnd(mglevels[i]->eventinterprestrict,… in PCMGACycle_Private()
24 for (i = 0; i < l; i++) { in PCMGACycle_Private()
26 if (!mglevels[i]->X) { in PCMGACycle_Private()
[all …]
|
| H A D | fmg.c | 8 PetscInt i, l = mglevels[0]->levels; in PCMGFCycle_Private() local
13 for (i = l - 1; i > 0; i--) { in PCMGFCycle_Private()
14 …if (mglevels[i]->eventinterprestrict) PetscCall(PetscLogEventBegin(mglevels[i]->eventinterprestric… in PCMGFCycle_Private()
15 … if (matapp) PetscCall(MatMatRestrict(mglevels[i]->restrct, mglevels[i]->B, &mglevels[i - 1]->B)); in PCMGFCycle_Private()
16 else PetscCall(MatRestrict(mglevels[i]->restrct, mglevels[i]->b, mglevels[i - 1]->b)); in PCMGFCycle_Private()
17 …if (mglevels[i]->eventinterprestrict) PetscCall(PetscLogEventEnd(mglevels[i]->eventinterprestrict,… in PCMGFCycle_Private()
30 for (i = 0; i < l - 1; i++) { in PCMGFCycle_Private()
31 PetscCall(PCMGMCycle_Private(pc, &mglevels[i], transpose, matapp, NULL)); in PCMGFCycle_Private()
32 …if (mglevels[i + 1]->eventinterprestrict) PetscCall(PetscLogEventBegin(mglevels[i + 1]->eventinter… in PCMGFCycle_Private()
33 …if (matapp) PetscCall(MatMatInterpolate(mglevels[i + 1]->interpolate, mglevels[i]->X, &mglevels[i … in PCMGFCycle_Private()
[all …]
|
| /petsc/src/vec/is/sf/tests/ |
| H A D | ex24.c | 17 for (PetscInt i = 0; i < nroots; i++) roota[i] = 0; in testOverlappingCommunication() local
18 for (PetscInt i = 0; i < nroots; i++) rootb[i] = 0; in testOverlappingCommunication() local
19 for (PetscInt i = 0; i < maxleaf + 1; i++) leafa[i] = (i + 1); in testOverlappingCommunication() local
20 for (PetscInt i = 0; i < maxleaf + 1; i++) leafb[i] = -(i + 1); in testOverlappingCommunication() local
26 …for (PetscInt i = 0; i < nroots; i++) PetscCheck(roota[i] >= 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "… in testOverlappingCommunication() local
27 …for (PetscInt i = 0; i < nroots; i++) PetscCheck(rootb[i] <= 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "… in testOverlappingCommunication() local
29 for (PetscInt i = 0; i < nroots; i++) roota[i] = 0; in testOverlappingCommunication() local
30 for (PetscInt i = 0; i < nroots; i++) rootb[i] = 0; in testOverlappingCommunication() local
37 …for (PetscInt i = 0; i < nroots; i++) PetscCheck(roota[i] >= 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "… in testOverlappingCommunication() local
38 …for (PetscInt i = 0; i < nroots; i++) PetscCheck(rootb[i] <= 0, PETSC_COMM_SELF, PETSC_ERR_PLIB, "… in testOverlappingCommunication() local
[all …]
|
| /petsc/src/sys/utils/ |
| H A D | mpimesg.c | 35 PetscMPIInt size, rank, *recv_buf, i, *iflags_local = NULL, *iflags_localm; in PetscGatherNumberOfMessages() local
47 for (i = 0; i < size; i++) { in PetscGatherNumberOfMessages()
48 if (ilengths[i]) iflags_local[i] = 1; in PetscGatherNumberOfMessages()
49 else iflags_local[i] = 0; in PetscGatherNumberOfMessages()
92 PetscMPIInt size, rank, tag, i, j; in PetscGatherMessageLengths() local
107 …for (i = 0; i < nrecvs; i++) PetscCallMPI(MPIU_Irecv((*olengths) + i, 1, MPI_INT, MPI_ANY_SOURCE, … in PetscGatherMessageLengths()
110 for (i = 0, j = 0; i < size; ++i) { in PetscGatherMessageLengths()
111 if (ilengths[i]) { in PetscGatherMessageLengths()
112 PetscCallMPI(MPIU_Isend((void *)(ilengths + i), 1, MPI_INT, i, tag, comm, s_waits + j)); in PetscGatherMessageLengths()
122 for (i = 0; i < nrecvs; ++i) { in PetscGatherMessageLengths()
[all …]
|
| /petsc/src/mat/impls/aij/seq/bas/ |
| H A D | spbas.c | 96 PetscInt i; in spbas_allocate_data() local
111 for (i = 1; i < nrows; i++) result->icols[i] = result->icols[i - 1] + result->row_nnz[i - 1]; in spbas_allocate_data()
120 … for (i = 1; i < nrows; i++) result->values[i] = result->values[i - 1] + result->row_nnz[i - 1]; in spbas_allocate_data()
123 for (i = 0; i < nrows; i++) { in spbas_allocate_data()
124 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
125 PetscCall(PetscMalloc1(r_nnz, &result->icols[i])); in spbas_allocate_data()
128 for (i = 0; i < nrows; i++) { in spbas_allocate_data()
129 r_nnz = result->row_nnz[i]; in spbas_allocate_data()
130 PetscCall(PetscMalloc1(r_nnz, &result->values[i])); in spbas_allocate_data()
182 PetscInt i, i1, i2; /* Loop counters for (partly) sorted arrays */ in spbas_mergesort_icols() local
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex10.c | 23 PetscInt i, its, m = 3, rdim, cdim, rstart, rend; in main() local
47 for (i = rstart; i < rend; i++) { in main()
48 v = (PetscScalar)(i - rstart + 100 * rank); in main()
49 PetscCall(VecSetValues(u, 1, &i, &v, INSERT_VALUES)); in main()
85 PetscInt i, j, k, i1, i2, j_1, j2, k1, k2, h1, h2, shiftx, shifty, shiftz; in GetElasticityMatrix() local
99 for (i = 0; i < 81; i++) PetscCall(PetscMalloc1(81, &K[i])); in GetElasticityMatrix()
108 for (i = 0; i < m; i++) { in GetElasticityMatrix()
110 base = 2 * k * shiftz + 2 * j * shifty + 2 * i * shiftx; in GetElasticityMatrix()
133 for (i = 0; i < 81; i++) PetscCall(PetscFree(K[i])); in GetElasticityMatrix()
143 for (i = nstart; i < N; i++) { in GetElasticityMatrix()
[all …]
|
| /petsc/src/vec/is/sf/tutorials/ |
| H A D | ex1.c | 21 PetscInt i, nroots, nleaves; in PetscSFViewCustomLocals_Private() local
31 …i = 0; i < nleaves; i++) PetscCall(PetscViewerASCIISynchronizedPrintf(viewer, "[%d] %" PetscInt_FM… in PetscSFViewCustomLocals_Private()
40 PetscInt i, bs = 2, nroots, nrootsalloc, nleaves, nleavesalloc, *mine, stride; in main() local
127 PetscInt i; in main() local
130 for (i = 0; i < nleaves; i++) mine[i] = stride * i; in main()
177 for (i = 0; i < nrootsalloc; i++) rootdata[i] = -1; in main()
178 for (i = 0; i < nroots; i++) rootdata[i * stride] = 100 * (rank + 1) + i; in main()
180 for (i = 0; i < nleavesalloc; i++) leafdata[i] = -1; in main()
184 for (PetscInt i = 0; i < 10; i++) { in main() local
202 for (i = 0; i < nrootsalloc; i++) rootdata[i] = '*'; in main()
[all …]
|
| /petsc/src/mat/tests/output/ |
| H A D | ex88_2.out | 5 -6.2902938000000000e+07 + 4.5741953400000000e+08i -8.6359366000000000e+07 + 8.2368680200000000e+08i
6 2.1562260800000000e+08 + 8.0739766600000000e+08i 1.0828994620000000e+09 + 1.2955916360000000e+09i
10 -5.68471e+19 + 3.02618e+19 i
11 -7.30611e+19 + 9.22255e+19 i
15 -5.68471e+19 + 3.02618e+19 i
16 -7.30611e+19 + 9.22255e+19 i
20 1.49262e+08 - 1.28111e+09 i
21 -1.29852e+09 - 2.10299e+09 i
24 1.49262e+08 - 1.28111e+09 i
25 -1.29852e+09 - 2.10299e+09 i
[all …]
|
| /petsc/src/ksp/pc/impls/pbjacobi/ |
| H A D | pbjacobi.c | 7 PetscInt i, ib, jb; in PCApply_PBJacobi() local
19 for (i = 0; i < m; i++) yy[i] = diag[i] * xx[i]; in PCApply_PBJacobi()
22 for (i = 0; i < m; i++) { in PCApply_PBJacobi()
23 x0 = xx[2 * i]; in PCApply_PBJacobi()
24 x1 = xx[2 * i + 1]; in PCApply_PBJacobi()
25 yy[2 * i] = diag[0] * x0 + diag[2] * x1; in PCApply_PBJacobi()
26 yy[2 * i + 1] = diag[1] * x0 + diag[3] * x1; in PCApply_PBJacobi()
31 for (i = 0; i < m; i++) { in PCApply_PBJacobi()
32 x0 = xx[3 * i]; in PCApply_PBJacobi()
33 x1 = xx[3 * i + 1]; in PCApply_PBJacobi()
[all …]
|
| /petsc/src/vec/vec/utils/tagger/tutorials/output/ |
| H A D | ex1_12.out | 3 0.720032 + 0.061794 i
4 0.0100223 + 0.144639 i
5 0.397778 + 0.0730366 i
6 0.103866 + 0.25078 i
7 0.984052 + 0.738516 i
8 0.401147 + 0.527649 i
9 0.144054 + 0.996504 i
10 0.106773 + 0.989053 i
11 0.731313 + 0.30846 i
12 0.180556 + 0.986045 i
[all …]
|
| H A D | ex1_13.out | 3 0.720032 + 0.061794 i
4 0.0100223 + 0.144639 i
5 0.397778 + 0.0730366 i
6 0.103866 + 0.25078 i
7 0.984052 + 0.738516 i
8 0.401147 + 0.527649 i
9 0.144054 + 0.996504 i
10 0.106773 + 0.989053 i
11 0.731313 + 0.30846 i
12 0.180556 + 0.986045 i
[all …]
|
| /petsc/src/mat/interface/ |
| H A D | matnull.c | 93 PetscInt n, N, dim, nmodes, i, j; in MatNullSpaceCreateRigidBody() local
114 for (i = 1; i < nmodes; i++) PetscCall(VecDuplicate(vec[0], &vec[i])); in MatNullSpaceCreateRigidBody()
115 for (i = 0; i < nmodes; i++) PetscCall(VecGetArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody()
117 for (i = 0; i < n; i++) { in MatNullSpaceCreateRigidBody()
119 v[0][i * 2 + 0] = sN; in MatNullSpaceCreateRigidBody()
120 v[0][i * 2 + 1] = 0.; in MatNullSpaceCreateRigidBody()
121 v[1][i * 2 + 0] = 0.; in MatNullSpaceCreateRigidBody()
122 v[1][i * 2 + 1] = sN; in MatNullSpaceCreateRigidBody()
124 v[2][i * 2 + 0] = -x[i * 2 + 1]; in MatNullSpaceCreateRigidBody()
125 v[2][i * 2 + 1] = x[i * 2 + 0]; in MatNullSpaceCreateRigidBody()
[all …]
|