| /petsc/src/tao/leastsquares/tutorials/ |
| H A D | chwirut2.c | 36 PetscReal y[NOBSERVATIONS]; /* array of dependent variables */ member
166 PetscReal *t = user->t, *y = user->y; in InitializeData() local
170 y[i] = 92.9000; in InitializeData()
172 y[i] = 78.7000; in InitializeData()
174 y[i] = 64.2000; in InitializeData()
176 y[i] = 64.9000; in InitializeData()
178 y[i] = 57.1000; in InitializeData()
180 y[i] = 43.3000; in InitializeData()
182 y[i] = 31.1000; in InitializeData()
184 y[i] = 23.6000; in InitializeData()
[all …]
|
| H A D | chwirut1.c | 33 PetscReal y[NOBSERVATIONS]; /* array of dependent variables */ member
108 PetscReal *y = user->y, *f, *t = user->t; 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()
168 PetscReal *t = user->t, *y = user->y; in InitializeData() local
172 y[i] = 92.9000; in InitializeData()
174 y[i] = 78.7000; in InitializeData()
176 y[i] = 64.2000; in InitializeData()
178 y[i] = 64.9000; in InitializeData()
180 y[i] = 57.1000; in InitializeData()
182 y[i] = 43.3000; in InitializeData()
[all …]
|
| H A D | chwirut1f.F90 | 20 PetscReal y(0:213)
53 f_v(i + 1) = y(i) - exp(-x_v(1)*t(i))/(x_v(2) + x_v(3)*t(i))
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
82 y(i) = 43.3000; t(i) = 1.2500; i = i + 1
83 y(i) = 31.1000; t(i) = 1.7500; i = i + 1
84 y(i) = 23.6000; t(i) = 2.2500; i = i + 1
[all …]
|
| H A D | chwirut2f.F90 | 15 PetscReal y(0:213)
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
36 y(i) = 31.0500; t(i) = 1.7500; i = i + 1
[all …]
|
| /petsc/src/tao/leastsquares/tests/ |
| H A D | chwirut2.c | 34 PetscReal y[NOBSERVATIONS]; /* array of dependent variables */ member
130 PetscReal *y = user->y, *f, *t = user->t; 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()
191 PetscReal *t = user->t, *y = user->y; in InitializeData() local
195 y[i] = 92.9000; in InitializeData()
197 y[i] = 78.7000; in InitializeData()
199 y[i] = 64.2000; in InitializeData()
201 y[i] = 64.9000; in InitializeData()
203 y[i] = 57.1000; in InitializeData()
205 y[i] = 43.3000; in InitializeData()
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex31.c | 74 const PetscScalar *y; in RHSFunction_Hull1972A1() local
77 PetscCall(VecGetArrayRead(Y, &y)); in RHSFunction_Hull1972A1()
79 f[0] = -y[0]; in RHSFunction_Hull1972A1()
80 PetscCall(VecRestoreArrayRead(Y, &y)); in RHSFunction_Hull1972A1()
87 const PetscScalar *y; in RHSJacobian_Hull1972A1() local
92 PetscCall(VecGetArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
96 PetscCall(VecRestoreArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
102 const PetscScalar *y; in IFunction_Hull1972A1() local
106 PetscCall(VecGetArrayRead(Y, &y)); in IFunction_Hull1972A1()
108 f[0] = -y[0]; in IFunction_Hull1972A1()
[all …]
|
| H A D | ex36.c | 49 const PetscScalar *y, *ydot; in IFunctionImplicit() local
54 PetscCall(VecGetArrayRead(Y, &y)); in IFunctionImplicit()
58 f[0] = ydot[0] / 1.e6 - ydot[1] / 1.e6 - PetscSinReal(200 * PETSC_PI * t) / 2500. + y[0] / 1000.; in IFunctionImplicit()
59 …dot[1] / 1.e6 - 0.0006666766666666667 + PetscExpReal((500 * (y[1] - y[2])) / 13.) / 1.e8 + y[1] / … in IFunctionImplicit()
60 …f[2] = ydot[2] / 500000. + 1.e-6 - PetscExpReal((500 * (y[1] - y[2])) / 13.) / 1.e6 + y[2] / 9000.; in IFunctionImplicit()
61 …/ 1.e6 - 0.0006676566666666666 + (99 * PetscExpReal((500 * (y[1] - y[2])) / 13.)) / 1.e8 + y[3] / … in IFunctionImplicit()
62 f[4] = (3 * ydot[4]) / 1.e6 - (3 * ydot[3]) / 1.e6 + y[4] / 9000.; in IFunctionImplicit()
64 PetscCall(VecRestoreArrayRead(Y, &y)); in IFunctionImplicit()
76 const PetscScalar *y, *ydot; in IJacobianImplicit() local
80 PetscCall(VecGetArrayRead(Y, &y)); in IJacobianImplicit()
[all …]
|
| /petsc/src/dm/impls/plex/tests/output/ |
| H A D | ex37_sphere_iges.out | 15 (x, y, z) = (0.000000, 0.000000, 25.400051)
20 (x, y, z) = (0.000000, 25.400051, 0.000000)
22 (x, y, z) = (0.000000, 0.000000, 25.400051)
27 (x, y, z) = (0.000000, 25.400051, 0.000000)
29 (x, y, z) = (-25.400051, 0.000000, 0.000000)
34 (x, y, z) = (-25.400051, 0.000000, 0.000000)
36 (x, y, z) = (0.000000, 0.000000, 25.400051)
42 (x, y, z) = (0.000000, 0.000000, 25.400051)
47 (x, y, z) = (-25.400051, 0.000000, 0.000000)
49 (x, y, z) = (0.000000, 0.000000, 25.400051)
[all …]
|
| H A D | ex37_sphere_egads.out | 25 (x, y, z) = (0.000000, 0.000000, 0.025400)
31 (x, y, z) = (0.000000, 0.025400, 0.000000)
33 (x, y, z) = (0.000000, 0.000000, 0.025400)
39 (x, y, z) = (0.000000, 0.025400, 0.000000)
41 (x, y, z) = (-0.025400, 0.000000, 0.000000)
47 (x, y, z) = (-0.025400, 0.000000, 0.000000)
49 (x, y, z) = (0.000000, 0.000000, 0.025400)
62 (x, y, z) = (0.000000, 0.000000, 0.025400)
68 (x, y, z) = (-0.025400, 0.000000, 0.000000)
70 (x, y, z) = (0.000000, 0.000000, 0.025400)
[all …]
|
| H A D | ex37_sphere_step.out | 26 (x, y, z) = (0.000000, 0.000000, 0.025400)
32 (x, y, z) = (0.000000, 0.025400, 0.000000)
34 (x, y, z) = (0.000000, 0.000000, 0.025400)
40 (x, y, z) = (0.000000, 0.025400, 0.000000)
42 (x, y, z) = (-0.025400, 0.000000, 0.000000)
48 (x, y, z) = (-0.025400, 0.000000, 0.000000)
50 (x, y, z) = (0.000000, 0.000000, 0.025400)
63 (x, y, z) = (0.000000, 0.000000, 0.025400)
69 (x, y, z) = (-0.025400, 0.000000, 0.000000)
71 (x, y, z) = (0.000000, 0.000000, 0.025400)
[all …]
|
| /petsc/src/benchmarks/ |
| H A D | MPI_Wtime.c | 5 double x, y; in main() local
10 y = MPI_Wtime(); in main()
13 y = MPI_Wtime(); in main()
14 y = MPI_Wtime(); in main()
15 y = MPI_Wtime(); in main()
16 y = MPI_Wtime(); in main()
17 y = MPI_Wtime(); in main()
18 y = MPI_Wtime(); in main()
19 y = MPI_Wtime(); in main()
20 y = MPI_Wtime(); in main()
[all …]
|
| H A D | PetscTime.c | 6 PetscLogDouble x, y; in main() local
11 PetscTime(&y); in main()
15 PetscTime(&y); in main()
16 PetscTime(&y); in main()
17 PetscTime(&y); in main()
18 PetscTime(&y); in main()
19 PetscTime(&y); in main()
20 PetscTime(&y); in main()
21 PetscTime(&y); in main()
22 PetscTime(&y); in main()
[all …]
|
| H A D | PetscGetTime.c | 6 PetscLogDouble x, y; in main() local
11 PetscCall(PetscTime(&y)); in main()
15 PetscCall(PetscTime(&y)); in main()
16 PetscCall(PetscTime(&y)); in main()
17 PetscCall(PetscTime(&y)); in main()
18 PetscCall(PetscTime(&y)); in main()
19 PetscCall(PetscTime(&y)); in main()
20 PetscCall(PetscTime(&y)); in main()
21 PetscCall(PetscTime(&y)); in main()
22 PetscCall(PetscTime(&y)); in main()
[all …]
|
| /petsc/src/dm/impls/plex/tutorials/output/ |
| H A D | ex17_nozzle_stp.out | 28 (x, y, z) = (-3.732071, -65.122507, -1.523202)
30 (x, y, z) = (-61.516928, -264.236200, -34.928953)
36 (x, y, z) = (-3.732071, -65.122507, -1.523202)
38 (x, y, z) = (3.718856, -65.122507, 2.784221)
44 (x, y, z) = (3.718856, -65.122507, 2.784221)
46 (x, y, z) = (61.503713, -264.236200, 36.189972)
52 (x, y, z) = (-61.516928, -264.236200, -34.928953)
54 (x, y, z) = (61.503713, -264.236200, 36.189972)
68 (x, y, z) = (3.718856, -65.122507, 2.784221)
70 (x, y, z) = (3.544753, -64.352475, 2.683571)
[all …]
|
| H A D | ex17_nozzle_igs.out | 29 (x, y, z) = (-61.516928, -264.236200, -34.928953)
31 (x, y, z) = (61.503713, -264.236200, 36.189972)
37 (x, y, z) = (61.503713, -264.236200, 36.189972)
39 (x, y, z) = (-61.516928, -264.236200, -34.928953)
53 (x, y, z) = (61.503713, -264.236200, 36.189972)
55 (x, y, z) = (-61.516928, -264.236200, -34.928953)
61 (x, y, z) = (3.718856, -65.122507, 2.784221)
63 (x, y, z) = (61.503713, -264.236200, 36.189972)
69 (x, y, z) = (3.718856, -65.122507, 2.784221)
71 (x, y, z) = (-3.732071, -65.122507, -1.523202)
[all …]
|
| H A D | ex17_sphere_egads.out | 25 (x, y, z) = (0.000000, 0.000000, 0.025400)
31 (x, y, z) = (0.000000, 0.025400, 0.000000)
33 (x, y, z) = (0.000000, 0.000000, 0.025400)
39 (x, y, z) = (0.000000, 0.025400, 0.000000)
41 (x, y, z) = (-0.025400, 0.000000, 0.000000)
47 (x, y, z) = (-0.025400, 0.000000, 0.000000)
49 (x, y, z) = (0.000000, 0.000000, 0.025400)
62 (x, y, z) = (0.000000, 0.000000, 0.025400)
68 (x, y, z) = (-0.025400, 0.000000, 0.000000)
70 (x, y, z) = (0.000000, 0.000000, 0.025400)
[all …]
|
| H A D | ex17_sphere_egadslite.out | 24 (x, y, z) = (0.000000, 0.000000, 0.025400)
30 (x, y, z) = (0.000000, 0.025400, 0.000000)
32 (x, y, z) = (0.000000, 0.000000, 0.025400)
38 (x, y, z) = (0.000000, 0.025400, 0.000000)
40 (x, y, z) = (-0.025400, 0.000000, 0.000000)
46 (x, y, z) = (-0.025400, 0.000000, 0.000000)
48 (x, y, z) = (0.000000, 0.000000, 0.025400)
61 (x, y, z) = (0.000000, 0.000000, 0.025400)
67 (x, y, z) = (-0.025400, 0.000000, 0.000000)
69 (x, y, z) = (0.000000, 0.000000, 0.025400)
[all …]
|
| H A D | ex18_sphere_shapeMod.out | 27 (x, y, z) = (0.000000, 0.025400, 0.000000)
29 (x, y, z) = (-0.000000, 0.000000, 0.025400)
35 (x, y, z) = (0.000000, 0.025400, 0.000000)
37 (x, y, z) = (-0.025400, 0.000000, 0.000000)
43 (x, y, z) = (-0.025400, 0.000000, 0.000000)
45 (x, y, z) = (-0.000000, 0.000000, 0.025400)
50 (x, y, z) = (-0.000000, 0.000000, 0.025400)
64 (x, y, z) = (-0.025400, 0.000000, 0.000000)
66 (x, y, z) = (-0.000000, 0.000000, 0.025400)
72 (x, y, z) = (-0.025400, 0.000000, 0.000000)
[all …]
|
| /petsc/src/sys/classes/draw/impls/image/ |
| H A D | drawimage.h | 13 static inline void PetscImageSetClip(PetscImage img, int x, int y, int w, int h) in PetscImageSetClip() argument
16 img->clip[1] = PetscClipInterval(y, 0, img->h - 1); /* ymin */ in PetscImageSetClip()
18 img->clip[3] = PetscClipInterval(y + h, 0, img->h); /* ymax+1 */ in PetscImageSetClip()
24 int y, ys = img->clip[1], ye = img->clip[3]; in PetscImageClear() local
25 for (y = ys; y < ye; y++) in PetscImageClear()
26 for (x = xs; x < xe; x++) img->buffer[y * img->w + x] = 0; in PetscImageClear()
29 static inline void PetscImageDrawPixel(PetscImage img, int x, int y, int c) in PetscImageDrawPixel() argument
32 if (y < img->clip[1] || y >= img->clip[3]) return; in PetscImageDrawPixel()
33 img->buffer[y * img->w + x] = (unsigned char)c; in PetscImageDrawPixel()
75 static inline void PetscImageDrawRectangle(PetscImage img, int x, int y, int w, int h, int c) in PetscImageDrawRectangle() argument
[all …]
|
| /petsc/src/mat/impls/sell/seq/seqcuda/ |
| H A D | sellcuda.cu | 75 …cInt *acolidx, const MatScalar *aval, const PetscInt *sliidx, const PetscScalar *x, PetscScalar *y) in matmult_seqsell_basic_kernel() argument
86 y[row] = sum; in matmult_seqsell_basic_kernel()
90 …atScalar *aval, const PetscInt *sliidx, const PetscScalar *x, const PetscScalar *y, PetscScalar *z) in matmultadd_seqsell_basic_kernel() argument
101 z[row] = y[row] + sum; in matmultadd_seqsell_basic_kernel()
108 …cInt *acolidx, const MatScalar *aval, const PetscInt *sliidx, const PetscScalar *x, PetscScalar *y) in matmult_seqsell_tiled_kernel9() argument
112 int tid = threadIdx.x + threadIdx.y * 32; in matmult_seqsell_tiled_kernel9()
120 …for (i = sliidx[slice_id] + threadIdx.x + 32 * threadIdx.y; i < sliidx[slice_id + 1]; i += 32 * BL… in matmult_seqsell_tiled_kernel9()
125 if (threadIdx.x < sliceheight) shared[threadIdx.x][threadIdx.y] = t; in matmult_seqsell_tiled_kernel9()
132 if (row < nrows && threadIdx.y == 0 && threadIdx.x < sliceheight) y[row] = shared[0][threadIdx.x]; in matmult_seqsell_tiled_kernel9()
137 …atScalar *aval, const PetscInt *sliidx, const PetscScalar *x, const PetscScalar *y, PetscScalar *z) in matmultadd_seqsell_tiled_kernel9() argument
[all …]
|
| /petsc/src/vec/is/sf/tests/ |
| H A D | ex8.c | 8 Vec x, y; in main() local
37 PetscCall(VecScatterCreateToZero(x, &vscat, &y)); in main()
38 PetscCall(PetscObjectSetName((PetscObject)y, "y")); in main()
41 PetscCall(VecScatterBegin(vscat, x, y, INSERT_VALUES, SCATTER_FORWARD)); in main()
42 PetscCall(VecScatterEnd(vscat, x, y, INSERT_VALUES, SCATTER_FORWARD)); in main()
43 if (rank == 0) PetscCall(VecView(y, PETSC_VIEWER_STDOUT_SELF)); in main()
46 PetscCall(VecScatterBegin(vscat, x, y, ADD_VALUES, SCATTER_FORWARD)); in main()
47 PetscCall(VecScatterEnd(vscat, x, y, ADD_VALUES, SCATTER_FORWARD)); in main()
48 if (rank == 0) PetscCall(VecView(y, PETSC_VIEWER_STDOUT_SELF)); in main()
51 PetscCall(VecScatterBegin(vscat, y, x, INSERT_VALUES, SCATTER_REVERSE)); in main()
[all …]
|
| /petsc/src/mat/impls/sell/seq/seqhip/ |
| H A D | sellhip.hip.cxx | 78 …cInt *acolidx, const MatScalar *aval, const PetscInt *sliidx, const PetscScalar *x, PetscScalar *y) in matmult_seqsell_basic_kernel() argument
89 y[row] = sum; in matmult_seqsell_basic_kernel()
93 …atScalar *aval, const PetscInt *sliidx, const PetscScalar *x, const PetscScalar *y, PetscScalar *z) in matmultadd_seqsell_basic_kernel() argument
104 z[row] = y[row] + sum; in matmultadd_seqsell_basic_kernel()
112 …cInt *acolidx, const MatScalar *aval, const PetscInt *sliidx, const PetscScalar *x, PetscScalar *y) in matmult_seqsell_tiled_kernel9() argument
116 int tid = threadIdx.x + threadIdx.y * WARP_SIZE; in matmult_seqsell_tiled_kernel9()
124 …for (i = sliidx[slice_id] + threadIdx.x + WARP_SIZE * threadIdx.y; i < sliidx[slice_id + 1]; i += … in matmult_seqsell_tiled_kernel9()
129 if (threadIdx.x < sliceheight) shared[threadIdx.x][threadIdx.y] = t; in matmult_seqsell_tiled_kernel9()
136 if (row < nrows && threadIdx.y == 0 && threadIdx.x < sliceheight) y[row] = shared[0][threadIdx.x]; in matmult_seqsell_tiled_kernel9()
141 …atScalar *aval, const PetscInt *sliidx, const PetscScalar *x, const PetscScalar *y, PetscScalar *z) in matmultadd_seqsell_tiled_kernel9() argument
[all …]
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_pc_py.py | 17 def apply(self, pc, x, y): argument
20 def applyT(self, pc, x, y): argument
21 self.apply(pc, x, y)
23 def applyS(self, pc, x, y): argument
24 self.apply(pc, x, y)
26 def applySL(self, pc, x, y): argument
27 self.applyS(pc, x, y)
29 def applySR(self, pc, x, y): argument
30 self.applyS(pc, x, y)
32 def applyRich(self, pc, x, y, w, tols): argument
[all …]
|
| /petsc/src/sys/classes/draw/utils/ |
| H A D | lg.c | 23 PetscErrorCode PetscDrawLGAddCommonPoint(PetscDrawLG lg, const PetscReal x, const PetscReal *y) in PetscDrawLGAddCommonPoint() argument
34 PetscCall(PetscArraycpy(tmpy, lg->y, lg->len)); in PetscDrawLGAddCommonPoint()
35 PetscCall(PetscFree2(lg->x, lg->y)); in PetscDrawLGAddCommonPoint()
37 lg->y = tmpy; in PetscDrawLGAddCommonPoint()
43 if (y[i] > lg->ymax) lg->ymax = y[i]; in PetscDrawLGAddCommonPoint()
44 if (y[i] < lg->ymin) lg->ymin = y[i]; in PetscDrawLGAddCommonPoint()
47 lg->y[lg->loc++] = y[i]; in PetscDrawLGAddCommonPoint()
73 PetscErrorCode PetscDrawLGAddPoint(PetscDrawLG lg, const PetscReal *x, const PetscReal *y) in PetscDrawLGAddPoint() argument
85 PetscCall(PetscArraycpy(tmpy, lg->y, lg->len)); in PetscDrawLGAddPoint()
86 PetscCall(PetscFree2(lg->x, lg->y)); in PetscDrawLGAddPoint()
[all …]
|
| /petsc/src/ksp/ksp/utils/lmvm/blas_cyclic/ |
| H A D | blas_cyclic.c | 16 … m, PetscScalar alpha, const PetscScalar x[], PetscScalar beta, PetscScalar y[], PetscInt y_stride) in AXPBY_Private() argument
18 for (PetscInt i = 0; i < m; i++) y[i * y_stride] = alpha * x[i] + beta * y[i * y_stride]; in AXPBY_Private()
21 …xt, PetscScalar alpha, const PetscScalar x[], PetscScalar beta, PetscScalar y[], PetscInt y_stride) in AXPBYCylic_Private() argument
28 AXPBY_Private(m, alpha, x, beta, y, y_stride); in AXPBYCylic_Private()
30 AXPBY_Private(i_next - i_oldest, alpha, &x[i_oldest], beta, &y[i_oldest * y_stride], y_stride); in AXPBYCylic_Private()
32 AXPBY_Private(i_next, alpha, x, beta, y, y_stride); in AXPBYCylic_Private()
33 AXPBY_Private(m - i_oldest, alpha, &x[i_oldest], beta, &y[i_oldest * y_stride], y_stride); in AXPBYCylic_Private()
38 …e VecAXPBYCyclic(PetscInt oldest, PetscInt next, PetscScalar alpha, Vec x, PetscScalar beta, Vec y) in VecAXPBYCyclic() argument
48 PetscValidHeaderSpecific(y, VEC_CLASSID, 6); in VecAXPBYCyclic()
49 PetscCheckSameComm(x, 4, y, 6); in VecAXPBYCyclic()
[all …]
|