Home
last modified time | relevance | path

Searched refs:Lx (Results 1 – 12 of 12) sorted by relevance

/petsc/src/dm/tests/
H A Dex47.c14 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_3d() local
26 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
33 PetscScalar x = (Lx * i) / M; in test_3d()
36 …va[k][j][i] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2) + PetscPow… in test_3d()
57 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_2d() local
68 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_2d()
74 PetscScalar x = (Lx * i) / M; in test_2d()
76 va[j][i] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2); in test_2d()
96 const PetscScalar Lx = 1.0, Ly = 1.0; in test_2d_nocoord() local
112 PetscScalar x = (Lx * i) / M; in test_2d_nocoord()
[all …]
H A Dex42.c16 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_3d() local
28 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
35 PetscScalar x = (Lx * i) / M; in test_3d()
38 …va[k][j][i] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2) + PetscPow… in test_3d()
59 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_2d() local
70 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_2d()
76 PetscScalar x = (Lx * i) / M; in test_2d()
78 va[j][i] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y - 0.5 * Ly, 2); in test_2d()
98 const PetscScalar Lx = 1.0, Ly = 1.0; in test_2d_nocoord() local
114 PetscScalar x = (Lx * i) / M; in test_2d_nocoord()
[all …]
H A Dex48.c28 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_3d() local
41 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d()
48 const PetscScalar x = (Lx * i) / M; in test_3d()
51 …for (c = 0; c < dof; ++c) va[k][j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(… in test_3d()
71 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_2d() local
83 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_2d()
89 const PetscScalar x = (Lx * i) / M; in test_2d()
91 …for (c = 0; c < dof; ++c) va[j][i][c] = PetscPowScalarInt(x - 0.5 * Lx, 2) + PetscPowScalarInt(y -… in test_2d()
110 const PetscScalar Lx = 1.0, Ly = 1.0, Lz = 1.0; in test_3d_compat() local
123 PetscCall(DMDASetUniformCoordinates(da, 0.0, Lx, 0.0, Ly, 0.0, Lz)); in test_3d_compat()
[all …]
/petsc/src/snes/tutorials/
H A Dex48.c255 PetscReal Lx, Ly, Lz; /* Model domain */ member
304 …p->b = s - 1000 * units->meter + 500 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx)… in THIInitialize_HOM_A()
317 …p->beta2 = 1000 * (1 + PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / … in THIInitialize_HOM_C()
326 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_X()
337 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_Y()
350 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_Z()
460 thi->Lx = 10.e3; in THICreate()
473 L = thi->Lx; in THICreate()
475 if (flg) thi->Lx = thi->Ly = L; in THICreate()
476 PetscCall(PetscOptionsReal("-thi_Lx", "X Domain size (m)", "", thi->Lx, &thi->Lx, NULL)); in THICreate()
[all …]
H A Dex5.c227 const PetscReal Lx = 1., Ly = 1.; in MMSSolution4() local
231 …u[0] = (PetscPowReal(x, 4) - PetscSqr(Lx) * PetscSqr(x)) * (PetscPowReal(y, 4) - PetscSqr(Ly) * Pe… in MMSSolution4()
237 const PetscReal Lx = 1., Ly = 1.; in MMSForcing4() local
241 …Lx)) * (PetscSqr(Ly) - 6 * PetscSqr(y)) + 2 * PetscSqr(y) * (PetscSqr(Lx) - 6 * PetscSqr(x)) * (Pe… in MMSForcing4()
H A Dex55.c174 const PetscReal Lx = 1., Ly = 1.; in MMSSolution4() local
176 …u[0] = (PetscPowReal(x, 4) - PetscSqr(Lx) * PetscSqr(x)) * (PetscPowReal(y, 4) - PetscSqr(Ly) * Pe… in MMSSolution4()
182 const PetscReal Lx = 1., Ly = 1.; in MMSForcing4() local
184 …Lx)) * (PetscSqr(Ly) - 6 * PetscSqr(y)) + 2 * PetscSqr(y) * (PetscSqr(Lx) - 6 * PetscSqr(x)) * (Pe… in MMSForcing4()
H A Dex56.c205 PetscReal Lx, mdisp[10], err[10]; in main() local
214 Lx = 1.; /* or ne for rod */ in main()
218 PetscCall(PetscOptionsReal("-lx", "Length of domain", "", Lx, &Lx, NULL)); in main()
/petsc/src/ts/tutorials/
H A Dex14.c241 PetscReal Lx, Ly, Lz; /* Model domain */ member
343 …* units->meter + 500 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * … in THIInitialize_HOM_A()
356 …p->beta2 = 1000 * (1 + PetscSinReal(x * 2 * PETSC_PI / thi->Lx) * PetscSinReal(y * 2 * PETSC_PI / … in THIInitialize_HOM_C()
366 …* units->meter + 100 * units->meter * PetscSinReal(x * 2 * PETSC_PI / thi->Lx); /* * sin(y*2*PETSC… in THIInitialize_HOM_F()
368 …p->h = (1 - (atan((x - thi->Lx / 2) / 1.) + PETSC_PI / 2.) / PETSC_PI) * 500 * units->m… in THIInitialize_HOM_F()
378 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_X()
389 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_Y()
401 …PetscReal x = xx * 2 * PETSC_PI / thi->Lx - PETSC_PI, y = yy * 2 * PETSC_PI / thi->Ly - PETSC_PI; … in THIInitialize_HOM_Z()
524 thi->Lx = 10.e3; in THICreate()
543 L = thi->Lx; in THICreate()
[all …]
H A Dex11.c472 … initLinearWave(EulerNode *ux, const PetscReal gamma, const PetscReal coord[], const PetscReal Lx);
1623 …t initLinearWave(EulerNode *ux, const PetscReal gamma, const PetscReal coord[], const PetscReal Lx) in initLinearWave() argument
1645 wcp[2] = wc[2] + eps * PetscCosReal(coord[0] * 2. * twopi / Lx); in initLinearWave()
/petsc/src/dm/dt/tests/
H A Dex7.c9 PetscReal *Lstarw, *Lx, *Lstar, *Lstarwcheck, wLx, Lstarwx; in CheckPullback() local
26 PetscCall(PetscMalloc2(Nk, &Lstarw, M * k, &Lx)); in CheckPullback()
45 Lx[l * M + i] = sum; in CheckPullback()
75 PetscCall(PetscDTAltVApply(M, k, ww, Lx, &wLx)); in CheckPullback()
80 PetscCall(PetscFree2(Lstarw, Lx)); in CheckPullback()
/petsc/src/ksp/pc/impls/sor/
H A Dsor.tex79 For the AIJ format $t$ is updated a row at a time to contain $ (b - Lx^{1/2}).$
/petsc/src/binding/petsc4py/src/petsc4py/PETSc/
H A DMat.pyx5447 """Solve Lx = b, given a factored matrix A = LU.