Home
last modified time | relevance | path

Searched refs:u (Results 1 – 25 of 536) sorted by relevance

12345678910>>...22

/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact8.c13 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering() local
56 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
58 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
59 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
60 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
61 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
62 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
64 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
65 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
66 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5_NaturalOrdering()
[all …]
H A Dsbaijfact7.c12 MatScalar *u, *d, *rtmp, *rtmp_ptr, work[25]; in MatCholeskyFactorNumeric_SeqSBAIJ_5() local
94 u = ba + ili * 25; in MatCholeskyFactorNumeric_SeqSBAIJ_5()
96 uik[0] = -(d[0] * u[0] + d[5] * u[1] + d[10] * u[2] + d[15] * u[3] + d[20] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
97 uik[1] = -(d[1] * u[0] + d[6] * u[1] + d[11] * u[2] + d[16] * u[3] + d[21] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
98 uik[2] = -(d[2] * u[0] + d[7] * u[1] + d[12] * u[2] + d[17] * u[3] + d[22] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
99 uik[3] = -(d[3] * u[0] + d[8] * u[1] + d[13] * u[2] + d[18] * u[3] + d[23] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
100 uik[4] = -(d[4] * u[0] + d[9] * u[1] + d[14] * u[2] + d[19] * u[3] + d[24] * u[4]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
102 uik[5] = -(d[0] * u[5] + d[5] * u[6] + d[10] * u[7] + d[15] * u[8] + d[20] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
103 uik[6] = -(d[1] * u[5] + d[6] * u[6] + d[11] * u[7] + d[16] * u[8] + d[21] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
104 uik[7] = -(d[2] * u[5] + d[7] * u[6] + d[12] * u[7] + d[17] * u[8] + d[22] * u[9]); in MatCholeskyFactorNumeric_SeqSBAIJ_5()
[all …]
H A Dsbaijfact5.c13 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering() local
57 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
59 uik[0] = -(diag[0] * u[0] + diag[4] * u[1] + diag[8] * u[2] + diag[12] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
60 uik[1] = -(diag[1] * u[0] + diag[5] * u[1] + diag[9] * u[2] + diag[13] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
61 uik[2] = -(diag[2] * u[0] + diag[6] * u[1] + diag[10] * u[2] + diag[14] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
62 uik[3] = -(diag[3] * u[0] + diag[7] * u[1] + diag[11] * u[2] + diag[15] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
64 uik[4] = -(diag[0] * u[4] + diag[4] * u[5] + diag[8] * u[6] + diag[12] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
65 uik[5] = -(diag[1] * u[4] + diag[5] * u[5] + diag[9] * u[6] + diag[13] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
66 uik[6] = -(diag[2] * u[4] + diag[6] * u[5] + diag[10] * u[6] + diag[14] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
67 uik[7] = -(diag[3] * u[4] + diag[7] * u[5] + diag[11] * u[6] + diag[15] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4_NaturalOrdering()
[all …]
H A Dsbaijfact6.c12 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_4() local
95 u = ba + ili * 16; in MatCholeskyFactorNumeric_SeqSBAIJ_4()
97 uik[0] = -(diag[0] * u[0] + diag[4] * u[1] + diag[8] * u[2] + diag[12] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
98 uik[1] = -(diag[1] * u[0] + diag[5] * u[1] + diag[9] * u[2] + diag[13] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
99 uik[2] = -(diag[2] * u[0] + diag[6] * u[1] + diag[10] * u[2] + diag[14] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
100 uik[3] = -(diag[3] * u[0] + diag[7] * u[1] + diag[11] * u[2] + diag[15] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
102 uik[4] = -(diag[0] * u[4] + diag[4] * u[5] + diag[8] * u[6] + diag[12] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
103 uik[5] = -(diag[1] * u[4] + diag[5] * u[5] + diag[9] * u[6] + diag[13] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
104 uik[6] = -(diag[2] * u[4] + diag[6] * u[5] + diag[10] * u[6] + diag[14] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
105 uik[7] = -(diag[3] * u[4] + diag[7] * u[5] + diag[11] * u[6] + diag[15] * u[7]); in MatCholeskyFactorNumeric_SeqSBAIJ_4()
[all …]
H A Dsbaijfact4.c13 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering() local
56 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
58 uik[0] = -(diag[0] * u[0] + diag[3] * u[1] + diag[6] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
59 uik[1] = -(diag[1] * u[0] + diag[4] * u[1] + diag[7] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
60 uik[2] = -(diag[2] * u[0] + diag[5] * u[1] + diag[8] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
62 uik[3] = -(diag[0] * u[3] + diag[3] * u[4] + diag[6] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
63 uik[4] = -(diag[1] * u[3] + diag[4] * u[4] + diag[7] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
64 uik[5] = -(diag[2] * u[3] + diag[5] * u[4] + diag[8] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
66 uik[6] = -(diag[0] * u[6] + diag[3] * u[7] + diag[6] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
67 uik[7] = -(diag[1] * u[6] + diag[4] * u[7] + diag[7] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3_NaturalOrdering()
[all …]
H A Dsbaijfact3.c12 MatScalar *u, *diag, *rtmp, *rtmp_ptr; in MatCholeskyFactorNumeric_SeqSBAIJ_3() local
94 u = ba + ili * 9; in MatCholeskyFactorNumeric_SeqSBAIJ_3()
96 uik[0] = -(diag[0] * u[0] + diag[3] * u[1] + diag[6] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
97 uik[1] = -(diag[1] * u[0] + diag[4] * u[1] + diag[7] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
98 uik[2] = -(diag[2] * u[0] + diag[5] * u[1] + diag[8] * u[2]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
100 uik[3] = -(diag[0] * u[3] + diag[3] * u[4] + diag[6] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
101 uik[4] = -(diag[1] * u[3] + diag[4] * u[4] + diag[7] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
102 uik[5] = -(diag[2] * u[3] + diag[5] * u[4] + diag[8] * u[5]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
104 uik[6] = -(diag[0] * u[6] + diag[3] * u[7] + diag[6] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
105 uik[7] = -(diag[1] * u[6] + diag[4] * u[7] + diag[7] * u[8]); in MatCholeskyFactorNumeric_SeqSBAIJ_3()
[all …]
/petsc/src/ts/tutorials/phasefield/
H A Dbiharmonic.c397 PetscScalar *u; in FormInitialSolution() local
411 PetscCall(DMDAVecGetArray(da, U, &u)); in FormInitialSolution()
435 if (r < .125) u[i] = 1.0; in FormInitialSolution()
436 else u[i] = -.5; in FormInitialSolution()
441 u[i] = f[scale * i]; in FormInitialSolution()
449 PetscCall(DMDAVecRestoreArray(da, U, &u)); in FormInitialSolution()
460 PetscScalar *u, l, r, c; in MyMonitor() local
486 PetscCall(DMDAVecGetArrayRead(da, localU, &u)); in MyMonitor()
512 …generate) yy[0] = PetscRealPart(.25 * (1. - u[i] * u[i]) * ctx->kappa * (u[i - 1] - u[i + 1]) * (u… in MyMonitor()
513 … else yy[0] = PetscRealPart(.25 * ctx->kappa * (u[i - 1] - u[i + 1]) * (u[i - 1] - u[i + 1]) * sx); in MyMonitor()
[all …]
/petsc/src/binding/petsc4py/demo/legacy/ode/
H A Dfastslowsplit.py25 def initialCondition(self, u): argument
26 u[0] = np.sqrt(2.0)
27 u[1] = np.sqrt(3.0)
28 u.assemble()
30 def evalRHSFunction(self, ts, t, u, f): argument
32 -2.0 * (-1.0 + u[0] * u[0] - np.cos(t)) / (2.0 * u[0])
33 + 0.05 * (-2.0 + u[1] * u[1] - np.cos(5.0 * t)) / (2.0 * u[1])
34 - np.sin(t) / (2.0 * u[0])
37 0.05 * (-1.0 + u[0] * u[0] - np.cos(t)) / (2.0 * u[0])
38 - (-2.0 + u[1] * u[1] - np.cos(5.0 * t)) / (2.0 * u[1])
[all …]
/petsc/src/ts/tutorials/hybrid/
H A Dex1fd.c31 const PetscScalar *u; in EventFunction() local
34 PetscCall(VecGetArrayRead(U, &u)); in EventFunction()
36 fvalue[0] = PetscRealPart(u[1] - actx->lambda1 * u[0]); in EventFunction()
38 fvalue[0] = PetscRealPart(u[1] - actx->lambda2 * u[0]); in EventFunction()
40 PetscCall(VecRestoreArrayRead(U, &u)); in EventFunction()
48 const PetscScalar *u; in ShiftGradients() local
54 PetscCall(VecGetArrayRead(U, &u)); in ShiftGradients()
57 denorm1 = -actx->lambda1 * (u[0] - 100. * u[1]) + 1. * (10. * u[0] + u[1]); in ShiftGradients()
58 denorm2 = -actx->lambda1 * (u[0] + 10. * u[1]) + 1. * (-100. * u[0] + u[1]); in ShiftGradients()
59 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm1 + 1.; in ShiftGradients()
[all …]
H A Dex1fwd.c33 PetscScalar *u; in MyMonitor() local
41 PetscCall(MatDenseGetColumn(sp, 2, &u)); in MyMonitor()
43 …5lf %20.15lf %20.15lf\n", (double)time, (double)PetscRealPart(u[0]), (double)PetscRealPart(u[1]))); in MyMonitor()
44 PetscCall(MatDenseRestoreColumn(sp, &u)); in MyMonitor()
53 const PetscScalar *u; in EventFunction() local
56 PetscCall(VecGetArrayRead(U, &u)); in EventFunction()
58 fvalue[0] = PetscRealPart(u[1] - actx->lambda1 * u[0]); in EventFunction()
60 fvalue[0] = PetscRealPart(u[1] - actx->lambda2 * u[0]); in EventFunction()
62 PetscCall(VecRestoreArrayRead(U, &u)); in EventFunction()
70 PetscScalar *u; in ShiftGradients() local
[all …]
H A Dex1adj.c31 const PetscScalar *u; in EventFunction() local
34 PetscCall(VecGetArrayRead(U, &u)); in EventFunction()
36 fvalue[0] = PetscRealPart(u[1] - actx->lambda1 * u[0]); in EventFunction()
38 fvalue[0] = PetscRealPart(u[1] - actx->lambda2 * u[0]); in EventFunction()
40 PetscCall(VecRestoreArrayRead(U, &u)); in EventFunction()
48 const PetscScalar *u; in ShiftGradients() local
54 PetscCall(VecGetArrayRead(U, &u)); in ShiftGradients()
57 denorm = -actx->lambda1 * (u[0] - 100. * u[1]) + 1. * (10. * u[0] + u[1]); in ShiftGradients()
58 A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm + 1.; in ShiftGradients()
59 A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm; in ShiftGradients()
[all …]
/petsc/src/ts/tutorials/multirate/
H A Dex3.c20 const PetscScalar *u; in RHSFunction() local
24 PetscCall(VecGetArrayRead(U, &u)); in RHSFunction()
26 …+ u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) + 0.05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.… in RHSFunction()
27 …u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) - (-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) /… in RHSFunction()
28 PetscCall(VecRestoreArrayRead(U, &u)); in RHSFunction()
35 const PetscScalar *u; in RHSFunctionslow() local
39 PetscCall(VecGetArrayRead(U, &u)); in RHSFunctionslow()
41 …+ u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) + 0.05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.… in RHSFunctionslow()
42 PetscCall(VecRestoreArrayRead(U, &u)); in RHSFunctionslow()
49 const PetscScalar *u; in RHSFunctionfast() local
[all …]
H A Dex3fastslowsplit.c25 const PetscScalar *u; in RHSFunctionslow() local
29 PetscCall(VecGetArrayRead(U, &u)); in RHSFunctionslow()
31 …+ u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) + 0.05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.… in RHSFunctionslow()
32 PetscCall(VecRestoreArrayRead(U, &u)); in RHSFunctionslow()
39 const PetscScalar *u; in RHSFunctionfast() local
43 PetscCall(VecGetArrayRead(U, &u)); in RHSFunctionfast()
45 …f[0] = 0.05 * (-1.0 + u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) - 5.0 * PetscSinScalar(5.0 *… in RHSFunctionfast()
46 PetscCall(VecRestoreArrayRead(U, &u)); in RHSFunctionfast()
53 const PetscScalar *u, *udot; in IFunctionfast() local
57 PetscCall(VecGetArrayRead(U, &u)); in IFunctionfast()
[all …]
/petsc/share/petsc/matlab/
H A Dlaplacian.m106 u = varargin{1}; variable
107 dim2 = size(u);
118 uint = round(u);
119 if max(uint~=u)
122 u = uint; clear uint variable
124 if max(u<=0 )
167 maxeigs = prod(u);
219 e1 = ones(u(1),1); %e1 = ones(u(1),1,'int8');
221 e2 = ones(u(2),1);
224 e3 = ones(u(3),1);
[all …]
/petsc/src/ts/tests/output/
H A Dex2_be_adapt.out1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00
5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00
7 At t = 4.000000e-03 u = 1.012111e+00 1.016155e+00 1.012111e+00
9 At t = 8.000000e-03 u = 1.024437e+00 1.032612e+00 1.024437e+00
11 At t = 1.256846e-02 u = 1.038732e+00 1.051712e+00 1.038732e+00
13 At t = 1.765770e-02 u = 1.054933e+00 1.073375e+00 1.054933e+00
15 At t = 2.275554e-02 u = 1.071441e+00 1.095468e+00 1.071441e+00
17 At t = 2.796519e-02 u = 1.088611e+00 1.118464e+00 1.088611e+00
19 At t = 3.312730e-02 u = 1.105922e+00 1.141669e+00 1.105922e+00
[all …]
H A Dex2_euler.out1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
3 At t = 1.000000e-03 u = 1.003000e+00 1.004000e+00 1.003000e+00
5 At t = 2.000000e-03 u = 1.006010e+00 1.008014e+00 1.006010e+00
7 At t = 3.000000e-03 u = 1.009030e+00 1.012042e+00 1.009030e+00
9 At t = 4.000000e-03 u = 1.012060e+00 1.016084e+00 1.012060e+00
11 At t = 5.000000e-03 u = 1.015100e+00 1.020140e+00 1.015100e+00
13 At t = 6.000000e-03 u = 1.018151e+00 1.024211e+00 1.018151e+00
15 At t = 7.000000e-03 u = 1.021211e+00 1.028296e+00 1.021211e+00
17 At t = 8.000000e-03 u = 1.024282e+00 1.032395e+00 1.024282e+00
19 At t = 9.000000e-03 u = 1.027363e+00 1.036508e+00 1.027363e+00
[all …]
H A Dex2_beuler.out1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00
5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00
7 At t = 3.000000e-03 u = 1.009060e+00 1.012084e+00 1.009060e+00
9 At t = 4.000000e-03 u = 1.012101e+00 1.016141e+00 1.012101e+00
11 At t = 5.000000e-03 u = 1.015151e+00 1.020212e+00 1.015151e+00
13 At t = 6.000000e-03 u = 1.018212e+00 1.024297e+00 1.018212e+00
15 At t = 7.000000e-03 u = 1.021283e+00 1.028396e+00 1.021283e+00
17 At t = 8.000000e-03 u = 1.024364e+00 1.032510e+00 1.024364e+00
19 At t = 9.000000e-03 u = 1.027456e+00 1.036638e+00 1.027456e+00
[all …]
/petsc/src/ts/tutorials/output/
H A Dex76_2d_tri_p2_p1_p1_tconv.out1 Timestep: 0000 time = 0 L_2 Error: [0.0156, 0.131, 0.0437] ||div u||: 1.27
2 Timestep: 0001 time = 0.1 L_2 Error: [0.0321, 31.2, 0.0346] ||div u||: 2.08
3 Timestep: 0002 time = 0.2 L_2 Error: [0.0332, 30.9, 0.0391] ||div u||: 2.2
4 Timestep: 0003 time = 0.3 L_2 Error: [0.0341, 30.1, 0.0399] ||div u||: 2.31
5 Timestep: 0004 time = 0.4 L_2 Error: [0.0351, 29.1, 0.0393] ||div u||: 2.44
6 Timestep: 0000 time = 0 L_2 Error: [0.0156, 0.131, 0.0437] ||div u||: 1.27
7 Timestep: 0001 time = 0.05 L_2 Error: [0.0317, 31.3, 0.0346] ||div u||: 2.03
8 Timestep: 0002 time = 0.1 L_2 Error: [0.0328, 31.2, 0.0328] ||div u||: 2.14
9 Timestep: 0003 time = 0.15 L_2 Error: [0.0332, 31, 0.0325] ||div u||: 2.19
10 Timestep: 0004 time = 0.2 L_2 Error: [0.0338, 30.8, 0.0324] ||div u||: 2.25
[all …]
H A Dex76_2d_tri_p2_p1_p1_tg_tconv.out1 Timestep: 0000 time = 0 L_2 Error: [0.0394, 0.174, 7.15e-17] ||div u||: 5.74
2 Timestep: 0001 time = 0.1 L_2 Error: [0.0113, 0.517, 2.89e-07] ||div u||: 0.681
3 Timestep: 0002 time = 0.2 L_2 Error: [0.00255, 0.0784, 8.85e-08] ||div u||: 0.0803
4 Timestep: 0003 time = 0.3 L_2 Error: [0.000358, 0.0155, 2.19e-08] ||div u||: 0.0112
5 Timestep: 0004 time = 0.4 L_2 Error: [3.67e-05, 0.00238, 5.21e-09] ||div u||: 0.00193
6 Timestep: 0000 time = 0 L_2 Error: [0.0394, 0.174, 7.15e-17] ||div u||: 5.74
7 Timestep: 0001 time = 0.05 L_2 Error: [0.0197, 0.427, 9.66e-07] ||div u||: 1.96
8 Timestep: 0002 time = 0.1 L_2 Error: [0.0097, 0.137, 7.05e-07] ||div u||: 0.679
9 Timestep: 0003 time = 0.15 L_2 Error: [0.00424, 0.0558, 3.45e-07] ||div u||: 0.243
10 Timestep: 0004 time = 0.2 L_2 Error: [0.00165, 0.0267, 1.46e-07] ||div u||: 0.089
[all …]
/petsc/src/vec/vec/tests/
H A Dex100cu.cu16 Vec u; in main() local
56 PetscCall(VecCreate(PETSC_COMM_WORLD, &u)); in main()
57 PetscCall(VecSetType(u, VECCUDA)); in main()
58 PetscCall(PetscObjectSetName((PetscObject)u, "Test_Vec")); in main()
59 PetscCall(VecSetSizes(u, PETSC_DECIDE, m)); in main()
60 PetscCall(VecSetFromOptions(u)); in main()
61 PetscCall(VecGetOwnershipRange(u, &low, &high)); in main()
62 PetscCall(VecGetLocalSize(u, &ldim)); in main()
66 PetscCall(VecSetValues(u, 1, &iglobal, &v, INSERT_VALUES)); in main()
68 PetscCall(VecAssemblyBegin(u)); in main()
[all …]
/petsc/src/ts/tutorials/autodiff/
H A Dadr_ex1.cxx43 const PetscScalar *u, *udot; in IFunctionPassive() local
47 PetscCall(VecGetArrayRead(U, &u)); in IFunctionPassive()
50 f[0] = udot[0] + ctx->k * u[0] * u[1]; in IFunctionPassive()
51 f[1] = udot[1] + ctx->k * u[0] * u[1]; in IFunctionPassive()
52 f[2] = udot[2] - ctx->k * u[0] * u[1]; in IFunctionPassive()
55 PetscCall(VecRestoreArrayRead(U, &u)); in IFunctionPassive()
65 const PetscScalar *u, *udot; in IFunctionActive1() local
72 PetscCall(VecGetArrayRead(U, &u)); in IFunctionActive1()
78 u_a[0] <<= u[0]; in IFunctionActive1()
79 u_a[1] <<= u[1]; in IFunctionActive1()
[all …]
/petsc/src/vec/vec/tutorials/
H A Dex15.c12 Vec u; in main() local
29 PetscCall(VecCreate(PETSC_COMM_WORLD, &u)); in main()
30 PetscCall(VecSetSizes(u, PETSC_DECIDE, m)); in main()
31 PetscCall(VecSetFromOptions(u)); in main()
32 PetscCall(VecGetOwnershipRange(u, &low, &high)); in main()
33 PetscCall(VecGetLocalSize(u, &ldim)); in main()
37 PetscCall(VecSetValues(u, 1, &iglobal, &v, INSERT_VALUES)); in main()
39 PetscCall(VecAssemblyBegin(u)); in main()
40 PetscCall(VecAssemblyEnd(u)); in main()
41 PetscCall(VecView(u, PETSC_VIEWER_STDOUT_WORLD)); in main()
[all …]
H A Dex10.c15 Vec u; in main() local
51 PetscCall(VecCreate(PETSC_COMM_WORLD, &u)); in main()
52 PetscCall(PetscObjectSetName((PetscObject)u, "Test_Vec")); in main()
53 PetscCall(VecSetSizes(u, PETSC_DECIDE, m)); in main()
54 PetscCall(VecSetFromOptions(u)); in main()
55 PetscCall(VecGetOwnershipRange(u, &low, &high)); in main()
56 PetscCall(VecGetLocalSize(u, &ldim)); in main()
60 PetscCall(VecSetValues(u, 1, &iglobal, &v, INSERT_VALUES)); in main()
62 PetscCall(VecAssemblyBegin(u)); in main()
63 PetscCall(VecAssemblyEnd(u)); in main()
[all …]
/petsc/src/ts/tutorials/advection-diffusion-reaction/
H A Dex2.c45 const PetscScalar *u, *udot; in IFunction() local
48 PetscCall(VecGetArrayRead(U, &u)); in IFunction()
51 f[0] = udot[0] - k1(ctx, t) * u[2] + ctx->k2 * u[0]; in IFunction()
52 f[1] = udot[1] - k1(ctx, t) * u[2] + ctx->k3 * u[1] * u[3] - ctx->sigma2; in IFunction()
53 f[2] = udot[2] - ctx->k3 * u[1] * u[3] + k1(ctx, t) * u[2]; in IFunction()
54 f[3] = udot[3] - ctx->k2 * u[0] + ctx->k3 * u[1] * u[3]; in IFunction()
55 PetscCall(VecRestoreArrayRead(U, &u)); in IFunction()
65 const PetscScalar *u, *udot; in IJacobian() local
68 PetscCall(VecGetArrayRead(U, &u)); in IJacobian()
75 J[1][1] = a + ctx->k3 * u[3]; in IJacobian()
[all …]
H A Dex1.c98 const PetscScalar *u, *udot; in IFunction() local
102 PetscCall(VecGetArrayRead(U, &u)); in IFunction()
105 f[0] = udot[0] + ctx->k * u[0] * u[1]; in IFunction()
106 f[1] = udot[1] + ctx->k * u[0] * u[1]; in IFunction()
107 f[2] = udot[2] - ctx->k * u[0] * u[1]; in IFunction()
108 PetscCall(VecRestoreArrayRead(U, &u)); in IFunction()
121 const PetscScalar *u, *udot; in IJacobian() local
124 PetscCall(VecGetArrayRead(U, &u)); in IJacobian()
126 J[0][0] = a + ctx->k * u[1]; in IJacobian()
127 J[0][1] = ctx->k * u[0]; in IJacobian()
[all …]

12345678910>>...22