| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex68_1.out | 1 L^2 norm of the error 4 65.3605
2 L^2 norm of the error 6 4.07018
3 L^2 norm of the error 8 101.108
4 L^2 norm of the error 10 41.5333
5 L^2 norm of the error 12 4.98713
6 L^2 norm of the error 14 0.783661
7 L^2 norm of the error 16 0.193284
8 L^2 norm of the error 18 0.0419655
9 L^2 norm of the error 20 0.00744736
10 L^2 norm of the error 22 0.00107939
[all …]
|
| H A D | ex69_1.out | 1 L^2 norm of the error 4 0.0217259
2 L^2 norm of the error 6 0.000141342
3 L^2 norm of the error 8 1.51625e-06
4 L^2 norm of the error 10 1.23208e-08
5 L^2 norm of the error 12 7.45863e-11
6 L^2 norm of the error 14 3.4819e-13
7 L^2 norm of the error 16 1.36721e-13
8 L^2 norm of the error 18 1.02545e-13
9 L^2 norm of the error 20 8.73312e-15
10 L^2 norm of the error 22 2.93241e-13
[all …]
|
| H A D | ex69_2.out | 1 L^2 norm of the error 4 0.0217259
2 L^2 norm of the error 6 0.000141342
3 L^2 norm of the error 8 1.51625e-06
4 L^2 norm of the error 10 1.23208e-08
5 L^2 norm of the error 12 7.45863e-11
6 L^2 norm of the error 14 3.47475e-13
7 L^2 norm of the error 16 1.33098e-13
8 L^2 norm of the error 18 1.02825e-13
9 L^2 norm of the error 20 6.59372e-15
10 L^2 norm of the error 22 2.9704e-13
[all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex48.c | 96 #define L (0.5 * (1. - G)) macro
101 {H, 0, 0, 0, L, 0, 0, 0},
102 {0, H, 0, 0, 0, L, 0, 0},
103 {0, 0, H, 0, 0, 0, L, 0},
104 {0, 0, 0, H, 0, 0, 0, L},
105 {L, 0, 0, 0, H, 0, 0, 0},
106 {0, L, 0, 0, 0, H, 0, 0},
107 {0, 0, L, 0, 0, 0, H, 0},
108 {0, 0, 0, L, 0, 0, 0, H}
111 …0}, {0, 0, 0}, {0, P *H, 0}, {M * L, M *L, P}, {P * L, 0, 0}, {0, 0, 0}, {…
[all …]
|
| /petsc/src/mat/tests/ |
| H A D | ex241f.F90 | 11 Mat :: L
23 PetscCallA(MatCreate(PETSC_COMM_WORLD, L, ierr))
24 PetscCallA(MatSetType(L, MATAIJ, ierr))
25 PetscCallA(MatSetSizes(L, PETSC_DECIDE, PETSC_DECIDE, n, n, ierr))
27 PetscCallA(MatSeqAIJSetPreallocation(L, i1, PETSC_NULL_INTEGER_ARRAY, ierr))
28 …PetscCallA(MatMPIAIJSetPreallocation(L, i1, PETSC_NULL_INTEGER_ARRAY, i0, PETSC_NULL_INTEGER_ARRAY…
29 PetscCallA(MatSetOption(L, MAT_IGNORE_ZERO_ENTRIES, PETSC_TRUE, ierr))
30 PetscCallA(MatSetOption(L, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_TRUE, ierr))
31 PetscCallA(MatSetOption(L, MAT_NO_OFF_PROC_ENTRIES, PETSC_TRUE, ierr))
33 PetscCallA(MatGetOwnershipRange(L, istart, iend, ierr))
[all …]
|
| /petsc/src/ksp/ksp/tests/benchmarkscatters/ |
| H A D | Streams-Intel | 193 Package L#0
194 NUMANode L#0 (P#0 32GB) + L3 L#0 (23MB)
195 L2 L#0 (256KB) + L1d L#0 (32KB) + Core L#0
196 L2 L#1 (256KB) + L1d L#1 (32KB) + Core L#1
197 L2 L#2 (256KB) + L1d L#2 (32KB) + Core L#2
198 L2 L#3 (256KB) + L1d L#3 (32KB) + Core L#3
199 L2 L#4 (256KB) + L1d L#4 (32KB) + Core L#4
200 L2 L#5 (256KB) + L1d L#5 (32KB) + Core L#5
201 L2 L#6 (256KB) + L1d L#6 (32KB) + Core L#6
202 L2 L#7 (256KB) + L1d L#7 (32KB) + Core L#7
[all …]
|
| H A D | Streams-MPICH | 193 Socket L#0
194 NUMANode L#0 (P#0 32GB) + L3 L#0 (23MB)
195 L2 L#0 (256KB) + L1d L#0 (32KB) + Core L#0
196 L2 L#1 (256KB) + L1d L#1 (32KB) + Core L#1
197 L2 L#2 (256KB) + L1d L#2 (32KB) + Core L#2
198 L2 L#3 (256KB) + L1d L#3 (32KB) + Core L#3
199 L2 L#4 (256KB) + L1d L#4 (32KB) + Core L#4
200 L2 L#5 (256KB) + L1d L#5 (32KB) + Core L#5
201 L2 L#6 (256KB) + L1d L#6 (32KB) + Core L#6
202 L2 L#7 (256KB) + L1d L#7 (32KB) + Core L#7
[all …]
|
| H A D | Baseline-Intel-8 | 21 …L/home/bsmith/petsc/arch-opt/lib -lpetsc -Wl,-rpath,/soft/compilers/intel/compilers_and_libraries/…
23 …09/linux/mpi/intel64/include -L/soft/compilers/intel/compilers_and_libraries_2016.0.109/linux/mpi/…
|
| H A D | Baseline-MPICH-8 | 21 …L/home/bsmith/petsc/arch-opt-mpich/lib -lpetsc -Wl,-rpath,/home/bsmith/petsc/arch-opt-mpich/lib -L…
23 gcc -fPIC -g -O -I/home/bsmith/petsc/arch-opt-mpich/include -L/home/bsmith/petsc/arch-opt-mpich/lib…
|
| H A D | Baseline-Intel-16 | 21 …L/home/bsmith/petsc/arch-opt/lib -lpetsc -Wl,-rpath,/soft/compilers/intel/compilers_and_libraries/…
23 …09/linux/mpi/intel64/include -L/soft/compilers/intel/compilers_and_libraries_2016.0.109/linux/mpi/…
|
| H A D | Baseline-MPICH-16 | 21 …L/home/bsmith/petsc/arch-opt-mpich/lib -lpetsc -Wl,-rpath,/home/bsmith/petsc/arch-opt-mpich/lib -L…
23 gcc -fPIC -g -O -I/home/bsmith/petsc/arch-opt-mpich/include -L/home/bsmith/petsc/arch-opt-mpich/lib…
|
| /petsc/src/dm/dt/tests/output/ |
| H A D | ex7_1234.out | 44 L:
46 L*:
48 L*w:
69 L:
71 L*:
73 L*w:
106 L:
108 L*:
110 L*w:
113 L:
[all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex14.c | 79 #define L (0.5 * (1. - G)) macro
84 {H, 0, 0, 0, L, 0, 0, 0},
85 {0, H, 0, 0, 0, L, 0, 0},
86 {0, 0, H, 0, 0, 0, L, 0},
87 {0, 0, 0, H, 0, 0, 0, L},
88 {L, 0, 0, 0, H, 0, 0, 0},
89 {0, L, 0, 0, 0, H, 0, 0},
90 {0, 0, L, 0, 0, 0, H, 0},
91 {0, 0, 0, L, 0, 0, 0, H}
94 …0}, {0, 0, 0}, {0, P *H, 0}, {M * L, M *L, P}, {P * L, 0, 0}, {0, 0, 0}, {…
[all …]
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmresdeflation.c | 22 PetscInt i, L, R, j; in KSPAGMRESQuickSort() local
34 L = deb[i]; in KSPAGMRESQuickSort()
36 if (L < R) { in KSPAGMRESQuickSort()
37 pivot_r = val_r[L]; in KSPAGMRESQuickSort()
38 pivot_i = val_i[L]; in KSPAGMRESQuickSort()
40 ipivot = perm[L]; in KSPAGMRESQuickSort()
42 while (L < R) { in KSPAGMRESQuickSort()
44 while (abs_val >= abs_pivot && L < R) { in KSPAGMRESQuickSort()
48 if (L < R) { in KSPAGMRESQuickSort()
49 val_r[L] = val_r[R]; in KSPAGMRESQuickSort()
[all …]
|
| /petsc/interfaces/lua/src/ |
| H A D | luasrc.c | 4 int lua_VecCreate(lua_State *L) in lua_VecCreate() argument
10 lua_pushlightuserdata(L, vec); in lua_VecCreate()
14 int lua_VecSetSize(lua_State *L) in lua_VecSetSize() argument
21 vec = (Vec)lua_touserdata(L, 1); in lua_VecSetSize()
22 n = (PetscInt)lua_tointegerx(L, 2, &isnum); in lua_VecSetSize()
27 int luaopen_libpetsc(lua_State *L) in luaopen_libpetsc() argument
30 lua_register(L, "VecCreate", (lua_CFunction)lua_VecCreate); in luaopen_libpetsc()
31 lua_register(L, "VecSetSize", (lua_CFunction)lua_VecSetSize); in luaopen_libpetsc()
|
| /petsc/src/dm/interface/ |
| H A D | dmperiodicity.c | 22 …MGetPeriodicity(DM dm, const PetscReal *maxCell[], const PetscReal *Lstart[], const PetscReal *L[]) in DMGetPeriodicity() argument
28 if (L) *L = dm->L; in DMGetPeriodicity()
47 …e DMSetPeriodicity(DM dm, const PetscReal maxCell[], const PetscReal Lstart[], const PetscReal L[]) in DMSetPeriodicity() argument
55 if (L) PetscAssertPointer(L, 4); in DMSetPeriodicity()
71 if (L) { in DMSetPeriodicity()
72 if (!dm->L) PetscCall(PetscMalloc1(dim, &dm->L)); in DMSetPeriodicity()
73 for (d = 0; d < dim; ++d) dm->L[d] = L[d]; in DMSetPeriodicity()
75 PetscCall(PetscFree(dm->L)); in DMSetPeriodicity()
76 dm->L = NULL; in DMSetPeriodicity()
78 …PetscCheck((dm->maxCell && dm->L) || (!dm->maxCell && !dm->L), PetscObjectComm((PetscObject)dm), P… in DMSetPeriodicity()
[all …]
|
| /petsc/src/dm/impls/swarm/ |
| H A D | data_bucket.c | 38 …gistration_function[], const char name[], const size_t size, const PetscInt L, DMSwarmDataField *D… in DMSwarmDataFieldCreate() argument
47 df->L = L; in DMSwarmDataFieldCreate()
50 PetscCall(PetscMalloc(size * L, &df->data)); in DMSwarmDataFieldCreate()
51 PetscCall(PetscMemzero(df->data, size * L)); in DMSwarmDataFieldCreate()
79 db->L = -1; in DMSwarmDataBucketCreate()
196 *sum = df->L; in DMSwarmDataFieldGetNumEntries()
211 if (new_L == df->L) PetscFunctionReturn(PETSC_SUCCESS); in DMSwarmDataFieldSetSize()
212 if (new_L > df->L) { in DMSwarmDataFieldSetSize()
215 …PetscCall(PetscMemzero(((char *)df->data) + df->L * df->atomic_size, (new_L - df->L) * df->atomic_… in DMSwarmDataFieldSetSize()
220 df->L = new_L; in DMSwarmDataFieldSetSize()
[all …]
|
| /petsc/src/ksp/pc/impls/lsc/ |
| H A D | lsc.c | 7 Mat L; member
45 Mat L, Lp, Qscale; in PCSetUp_LSC() local
51 PetscCall(PetscObjectQuery((PetscObject)pc->mat, "LSC_L", (PetscObject *)&L)); in PCSetUp_LSC()
52 if (!L) PetscCall(PetscObjectQuery((PetscObject)pc->pmat, "LSC_L", (PetscObject *)&L)); in PCSetUp_LSC()
61 …PetscCheck(L || Lp, PetscObjectComm((PetscObject)pc), PETSC_ERR_USER, "The user must provide an L … in PCSetUp_LSC()
62 if (!L && Lp) L = Lp; in PCSetUp_LSC()
63 else if (L && !Lp) Lp = L; in PCSetUp_LSC()
72 if (!L) { in PCSetUp_LSC()
79 … if (!lsc->L) PetscCall(MatMatMult(CAdiaginv, B, MAT_INITIAL_MATRIX, PETSC_CURRENT, &lsc->L)); in PCSetUp_LSC()
80 else PetscCall(MatMatMult(CAdiaginv, B, MAT_REUSE_MATRIX, PETSC_CURRENT, &lsc->L)); in PCSetUp_LSC()
[all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dtaltv.c | 323 PetscErrorCode PetscDTAltVPullback(PetscInt N, PetscInt M, const PetscReal *L, PetscInt k, const Pe… in PetscDTAltVPullback() argument
339 for (j = 0; j < Mk; j++) sum += L[j * Nk + i] * w[j]; in PetscDTAltVPullback()
348 for (j = 0; j < Mk; j++) sum += L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * w[j] * mult[j]; in PetscDTAltVPullback()
360 …; j < Mk; j++) sum += (L[pairs[j][0] * N + pairs[i][0]] * L[pairs[j][1] * N + pairs[i][1]] - L[pai… in PetscDTAltVPullback()
375 … (L[pairs[offj + j][0] * N + pairs[offi + i][0]] * L[pairs[offj + j][1] * N + pairs[offi + i][1]] … in PetscDTAltVPullback()
379 …PetscReal detL = L[0] * (L[4] * L[8] - L[5] * L[7]) + L[1] * (L[5] * L[6] - L[3] * L[8]) + L[2] * … in PetscDTAltVPullback()
416 for (l = 0; l < k; l++) prod *= L[subsetw[perm[l]] * N + subsetv[l]]; in PetscDTAltVPullback()
452 PetscErrorCode PetscDTAltVPullbackMatrix(PetscInt N, PetscInt M, const PetscReal *L, PetscInt k, Pe… in PetscDTAltVPullbackMatrix() argument
472 for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[j * Nk + i]; in PetscDTAltVPullbackMatrix()
476 …for (j = 0; j < Mk; j++) Lstar[i * Mk + j] = L[(Mk - 1 - j) * Nk + (Nk - 1 - i)] * mult[i] * mult[… in PetscDTAltVPullbackMatrix()
[all …]
|
| /petsc/src/mat/graphops/order/ |
| H A D | spectral.c | 19 PetscErrorCode MatCreateLaplacian(Mat A, PetscReal tol, PetscBool weighted, Mat *L) in MatCreateLaplacian() argument
29 PetscCall(MatCreate(PetscObjectComm((PetscObject)A), L)); in MatCreateLaplacian()
32 PetscCall(MatSetSizes(*L, m, n, M, N)); in MatCreateLaplacian()
61 PetscCall(MatSetFromOptions(*L)); in MatCreateLaplacian()
62 PetscCall(MatXAIJSetPreallocation(*L, 1, dnnz, onnz, NULL, NULL)); in MatCreateLaplacian()
63 PetscCall(MatSetUp(*L)); in MatCreateLaplacian()
91 PetscCall(MatSetValues(*L, 1, &r, newcols, newCols, newVals, INSERT_VALUES)); in MatCreateLaplacian()
94 PetscCall(MatAssemblyBegin(*L, MAT_FINAL_ASSEMBLY)); in MatCreateLaplacian()
95 PetscCall(MatAssemblyEnd(*L, MAT_FINAL_ASSEMBLY)); in MatCreateLaplacian()
105 Mat L; in MatGetOrdering_Spectral() local
[all …]
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex84.c | 6 PetscErrorCode LSCLoadOperators(Mat *A, Mat *Q, Mat *L, Vec *rhs, IS *velocity, IS *pressure) in LSCLoadOperators() argument
15 if (L) PetscCall(MatCreate(PETSC_COMM_WORLD, L)); in LSCLoadOperators()
28 if (L) PetscCall(MatLoad(*L, viewer)); in LSCLoadOperators()
35 Mat A, Q, L = NULL; in port_lsd_bfbt() local
50 if (commute_lsc) PetscCall(LSCLoadOperators(&A, &Q, &L, &b, &isu, &isp)); in port_lsd_bfbt()
72 if (commute_lsc) PetscCall(PetscObjectCompose((PetscObject)lsc_pc_pmat, "LSC_L", (PetscObject)L)); in port_lsd_bfbt()
79 if (L) PetscCall(MatDestroy(&L)); in port_lsd_bfbt()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | 1d-2elems.e | 49 …���� A������ ���_Netcdf4Coordinates�������� ������������������ L��������DIMENSIO…
51 …�������dimension��������� ����������������<���������������ǑXFSSE�L������ ���������+…
61 …L��������DIMENSION_LIST������������������������������������������wOCHK�����…
66 ��� L��������DIMENSION_LIST����������������������������������������B��oOCHK |���…
68 ��� L��������DIMENSION_LIST����������������������������������������~9;�OCHK����…
72 …���� A������ ���_Netcdf4Coordinates�������� ������������������ L��������DIMENSIO…
73 …L��������DIMENSION_LIST����������������������������������������ʅ�OCHK������…
|
| H A D | gmsh-3d-binary-32.msh | 32 ���G���H���I���J���K���L���M���N���O���P���Q���R���S��� ������?��������@������?@������?�������� ���…
54 …L���!������!���L���I������I���P���!������#���L���"������G���L���#������#���R���G������'���N�…
56 …,���_���V���H���V���`���,���I���T���Y���U���J���U���^���T���K���U���Y���X���L���X���Z���U���M���W�…
60 …o���}�����������@���>���|�����������^���|���U���~�������e���b���~�����������L�������J�����������b�…
62 …L���/���G�������#��H���c���M�������$��.���3���d���~���%��:���V���_���{���&��B���r���:���{���'…
|
| H A D | gmsh-3d-binary-64.msh | 30 �������G�������H�������I�������J�������K�������L�������M�������N�������O�������P�������Q�������R���…
44 …���#���������������������S���������������������)�������S�������������� ������� �������L�������
46 �������H�������L������� ������� ������� �������P�������H������� �������
47 �������L������� �������
48 …L�������I��������������I�������Q������� ���������������������O����������������������������O�…
50 …����V�������J�������V�������Z�������T�������K�������U�������Y�������T�������L�������T�������[�����…
54 …������������-���������������������M�������~�������G���������������������L�������G�������������…
56 …L�������I���������������.������M�������0�������J�������~�������/������A�������>�������X���������…
|
| /petsc/src/ksp/pc/impls/sor/ |
| H A D | sor.tex | 11 Consider the matrix problem $ A x = b$, where $A = L + U + D$.
38 x^{1/2} = \omega (L + D)^{-1} b
41 x = (1 - \omega) x^{1/2} + \omega (U + D)^{-1}(b - L x^{1/2}) = x^{1/2} + \omega (U+D)^{-1}(b - A x…
58 …b $ and is updated a column of the matrix at a time to contain the value of $ (b - L x^{1/2})$ that
86 x = \omega (L + D)^{-1}b
97 x = (L + D/\omega)^{-1}b
100 Rather than applying the left preconditioner obtained by apply the two step process $ (L + D/\omega…
103 (L + D/\omega)^{-1} A (U + D/\omega)^{-1} y = (L + D/\omega)^{-1} b.
108 …(L + D/\omega)^{-1} A (U + D/\omega)^{-1} & = & (L + D/\omega)^{-1} (L + D + U) (U + D/\omega)^{-…
109 & = & (L + D/\omega)^{-1} (L + D/\omega + U + D/\omega + D - 2D/\omega) (U + D/\omega)^{-1} \\
[all …]
|