Home
last modified time | relevance | path

Searched full:y (Results 1 – 25 of 1114) sorted by relevance

12345678910>>...45

/petsc/share/petsc/datafiles/meshes/
H A Dhybrid_hexwedge.msh7 …Y���������@�@��������������Y�������(\��@�������?������Y���������ʥ@��������������Y���������@�@�…
8 ���������@�(\��?�@������Y� ��������ǥ@�z�GA�@������Y� ��������@�@�����@�@������Y�
/petsc/src/ts/tutorials/
H A Dex31.c71 PetscErrorCode RHSFunction_Hull1972A1(TS ts, PetscReal t, Vec Y, Vec F, void *s) in RHSFunction_Hull1972A1() argument
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()
85 PetscErrorCode RHSJacobian_Hull1972A1(TS ts, PetscReal t, Vec Y, Mat A, Mat B, void *s) in RHSJacobian_Hull1972A1() argument
87 const PetscScalar *y; in RHSJacobian_Hull1972A1() local
92 PetscCall(VecGetArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
96 PetscCall(VecRestoreArrayRead(Y, &y)); in RHSJacobian_Hull1972A1()
100 PetscErrorCode IFunction_Hull1972A1(TS ts, PetscReal t, Vec Y, Vec Ydot, Vec F, void *s) in IFunction_Hull1972A1() argument
[all …]
H A Dex36.c4 …` This example illustrates the implementation of an implicit DAE index-1 of form M y'=f(t,y) with …
12 [ -(U(t) - y[0])/1000 ]
13 [ -6/R + y[1]/4500 + 0.01 * h(y[1]-y[2]) ]
14 f(t,y)= [ y[2]/R - h(y[1]-y[2]) ]
15 [ (y[3]-6)/9000 + 0.99 * h([y1]-y[2]) ]
16 [ y[4]/9000 ]
47 static PetscErrorCode IFunctionImplicit(TS ts, PetscReal t, Vec Y, Vec Ydot, Vec F, PetscCtx ctx) in IFunctionImplicit() argument
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()
[all …]
/petsc/src/tao/leastsquares/tutorials/
H A Dchwirut2.c24 y = exp[-b1*x]/(b2+b3*x) + e \n";
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()
[all …]
H A Dchwirut1.c24 y = exp[-b1*x]/(b2+b3*x) + e \n";
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()
[all …]
H A Dchwirut1f.F9020 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 Dchwirut2f.F9015 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 Dchwirut2.c25 y = exp[-b1*x]/(b2+b3*x) + e \n";
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()
[all …]
/petsc/src/mat/utils/
H A Daxpy.c3 static PetscErrorCode MatTransposeAXPY_Private(Mat Y, PetscScalar a, Mat X, MatStructure str, Mat T) in MatTransposeAXPY_Private() argument
21 A = Y; in MatTransposeAXPY_Private()
23 …PetscCall(PetscInfo(NULL, "Transposing X because Y of type MATTRANSPOSEVIRTUAL to perform MatAXPY(… in MatTransposeAXPY_Private()
31 A = Y; in MatTransposeAXPY_Private()
33 …PetscCall(PetscInfo(NULL, "Hermitian transposing X because Y of type MATHERMITIANTRANSPOSEVIRTUAL … in MatTransposeAXPY_Private()
43 static PetscErrorCode MatAXPY_BasicWithTypeCompare(Mat Y, PetscScalar a, Mat X, MatStructure str) in MatAXPY_BasicWithTypeCompare() argument
48 PetscCall(MatIsShell(Y, &flg)); in MatAXPY_BasicWithTypeCompare()
52 PetscCall(MatGetOperation(Y, MATOP_AXPY, (PetscErrorCodeFn **)&f)); in MatAXPY_BasicWithTypeCompare()
54 PetscCall((*f)(Y, a, X, str)); in MatAXPY_BasicWithTypeCompare()
58 /* no need to preallocate if Y is dense */ in MatAXPY_BasicWithTypeCompare()
[all …]
/petsc/src/dm/impls/plex/tutorials/output/
H A Dex17_nozzle_stp.out28 (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 Dex17_nozzle_igs.out29 (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 Dex17_sphere_egads.out25 (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 Dex17_sphere_egadslite.out24 (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 Dex18_sphere_shapeMod.out27 (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/binding/petsc4py/test/
H A Dtest_pc_py.py17 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/ksp/ksp/utils/lmvm/blas_cyclic/
H A Dblas_cyclic.c16 … 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 …]
/petsc/src/mat/impls/sell/seq/seqcuda/
H A Dsellcuda.cu75 …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/mat/impls/sell/seq/seqhip/
H A Dsellhip.hip.cxx78 …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/benchmarks/
H A DMPI_Wtime.c5 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 DPetscTime.c6 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 …]
/petsc/src/snes/tutorials/output/
H A Dex58_3.out5 0 Line search: fty/||y|| = -0.033171, lambda = 0.5
6 1 Line search: fty/||y|| = 0.0180532, lambda = 0.25
7 2 Line search: fty/||y|| = 0.0126471, lambda = 0.375
8 3 Line search: fty/||y|| = 0.000436031, lambda = 0.4375
9 4 Line search: fty/||y|| = -0.0166589, lambda = 0.46875
10 5 Line search: fty/||y|| = -0.00698054, lambda = 0.453125
11 6 Line search: fty/||y|| = -0.00295682, lambda = 0.445312
12 7 Line search: fty/||y|| = -0.00118288, lambda = 0.441406
13 8 Line search: fty/||y|| = -0.000354409, lambda = 0.439453
14 9 Line search: fty/||y|| = 4.55115e-05, lambda = 0.438477
[all …]
/petsc/src/sys/classes/draw/utils/
H A Dlg.c12 - y - the new y coordinate point for each curve.
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()
62 - y - array containing the y coordinate for the point on each curve
73 PetscErrorCode PetscDrawLGAddPoint(PetscDrawLG lg, const PetscReal *x, const PetscReal *y) in PetscDrawLGAddPoint() argument
[all …]
/petsc/src/vec/is/sf/tests/
H A Dex8.c8 Vec x, y; in main() local
37 PetscCall(VecScatterCreateToZero(x, &vscat, &y)); in main()
38 PetscCall(PetscObjectSetName((PetscObject)y, "y")); in main()
40 /* Test PetscSFBcastAndOp with op = MPI_REPLACE, which does y = x on rank 0 */ 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()
45 /* Test PetscSFBcastAndOp with op = MPI_SUM, which does y += x */ 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()
[all …]
H A Dex2.c15 Vec x, y; in main() local
26 /* Create two CUDA vectors x, y. Though we only care y's memory on host, we make y a CUDA vector, in main()
27 … since we want to have y's memory on host pinned (i.e.,non-pagable), to really trigger asynchronous in main()
32 PetscCall(VecCreateSeq(PETSC_COMM_WORLD, n, &y)); in main()
33 PetscCall(VecSetFromOptions(y)); in main()
35 /* Init x, y, and push them to GPU (their offloadmask = PETSC_OFFLOAD_GPU) */ in main()
40 PetscCall(VecSet(y, 314)); in main()
42 /* Pull y to CPU (make its offloadmask = PETSC_OFFLOAD_CPU) */ in main()
43 PetscCall(VecGetArray(y, &val)); in main()
44 PetscCall(VecRestoreArray(y, &val)); in main()
[all …]
/petsc/src/mat/tests/
H A Dex201f.F9011 subroutine mymatmult(A, x, y, ierr) argument
13 Vec x, y
19 subroutine mymatmultadd(A, x, y, z, ierr) argument
21 Vec x, y, z
27 subroutine mymatmulttranspose(A, x, y, ierr) argument
29 Vec x, y
35 subroutine mymatmulthermitiantranspose(A, x, y, ierr) argument
37 Vec x, y
43 subroutine mymatmulttransposeadd(A, x, y, z, ierr) argument
45 Vec x, y, z
[all …]

12345678910>>...45